Changeset 2037 for palm/trunk/SOURCE/diffusion_s.f90
- Timestamp:
- Oct 26, 2016 11:15:40 AM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/diffusion_s.f90
r2001 r2037 20 20 ! Current revisions: 21 21 ! ------------------ 22 ! 22 ! Anelastic approximation implemented 23 23 ! 24 24 ! Former revisions: … … 113 113 114 114 USE arrays_3d, & 115 ONLY: ddzu, ddzw, kh, tend 115 ONLY: ddzu, ddzw, kh, tend, drho_air, rho_air_zw 116 116 117 117 USE control_parameters, & … … 190 190 + 0.5_wp * ( & 191 191 ( kh(k,j,i) + kh(k+1,j,i) ) * ( s(k+1,j,i)-s(k,j,i) ) * ddzu(k+1) & 192 * rho_air_zw(k) & 192 193 - ( kh(k,j,i) + kh(k-1,j,i) ) * ( s(k,j,i)-s(k-1,j,i) ) * ddzu(k) & 193 ) * ddzw(k) 194 * rho_air_zw(k-1) & 195 ) * ddzw(k) * drho_air(k) 194 196 ENDDO 195 197 … … 205 207 * ( s(k+1,j,i)-s(k,j,i) ) & 206 208 * ddzu(k+1) & 209 * rho_air_zw(k) & 207 210 + s_flux_b(j,i) & 208 ) * ddzw(k) 211 ) * ddzw(k) * drho_air(k) 209 212 210 213 ENDIF … … 222 225 * ( s(k,j,i)-s(k-1,j,i) ) & 223 226 * ddzu(k) & 224 ) * ddzw(k) 227 * rho_air_zw(k-1) & 228 ) * ddzw(k) * drho_air(k) 225 229 226 230 ENDIF … … 240 244 241 245 USE arrays_3d, & 242 ONLY: ddzu, ddzw, kh, tend 246 ONLY: ddzu, ddzw, kh, tend, drho_air, rho_air_zw 243 247 244 248 USE control_parameters, & … … 324 328 + 0.5_wp * ( & 325 329 ( kh(k,j,i) + kh(k+1,j,i) ) * ( s(k+1,j,i)-s(k,j,i) ) * ddzu(k+1) & 330 * rho_air_zw(k) & 326 331 - ( kh(k,j,i) + kh(k-1,j,i) ) * ( s(k,j,i)-s(k-1,j,i) ) * ddzu(k) & 327 ) * ddzw(k) 332 * rho_air_zw(k-1) & 333 ) * ddzw(k) * drho_air(k) 328 334 ENDIF 329 335 ENDDO … … 338 344 * ( s(k+1,j,i)-s(k,j,i) ) & 339 345 * ddzu(k+1) & 346 * rho_air_zw(k) & 340 347 + s_flux_b(j,i) & 341 ) * ddzw(k) 348 ) * ddzw(k) * drho_air(k) 342 349 ENDIF 343 350 … … 351 358 * ( s(k,j,i)-s(k-1,j,i) ) & 352 359 * ddzu(k) & 353 ) * ddzw(k) 360 * rho_air_zw(k-1) & 361 ) * ddzw(k) * drho_air(k) 354 362 ENDIF 355 363 ENDDO … … 370 378 371 379 USE arrays_3d, & 372 ONLY: ddzu, ddzw, kh, tend 380 ONLY: ddzu, ddzw, kh, tend, drho_air, rho_air_zw 373 381 374 382 USE control_parameters, & … … 446 454 + 0.5_wp * ( & 447 455 ( kh(k,j,i) + kh(k+1,j,i) ) * ( s(k+1,j,i)-s(k,j,i) ) * ddzu(k+1) & 456 * rho_air_zw(k) & 448 457 - ( kh(k,j,i) + kh(k-1,j,i) ) * ( s(k,j,i)-s(k-1,j,i) ) * ddzu(k) & 449 ) * ddzw(k) 458 * rho_air_zw(k-1) & 459 ) * ddzw(k) * drho_air(k) 450 460 ENDDO 451 461 … … 459 469 * ( s(k+1,j,i)-s(k,j,i) ) & 460 470 * ddzu(k+1) & 471 * rho_air_zw(k) & 461 472 + s_flux_b(j,i) & 462 ) * ddzw(k) 473 ) * ddzw(k) * drho_air(k) 463 474 464 475 ENDIF … … 474 485 * ( s(k,j,i)-s(k-1,j,i) ) & 475 486 * ddzu(k) & 476 ) * ddzw(k) 487 * rho_air_zw(k-1) & 488 ) * ddzw(k) * drho_air(k) 477 489 478 490 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.