Changeset 2037 for palm/trunk/SOURCE/diffusion_w.f90
- Timestamp:
- Oct 26, 2016 11:15:40 AM (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/diffusion_w.f90
r2001 r2037 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! Anelastic approximation implemented 23 23 ! 24 24 ! Former revisions: … … 119 119 120 120 USE arrays_3d, & 121 ONLY : ddzu, ddzw, km, tend, u, v, w 121 ONLY : ddzu, ddzw, km, tend, u, v, w, drho_air_zw, rho_air 122 122 123 123 USE control_parameters, & … … 184 184 & + 2.0_wp * ( & 185 185 & km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) & 186 & * rho_air(k+1) & 186 187 & - km(k,j,i) * ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k) & 187 & ) * ddzu(k+1) 188 & * rho_air(k) & 189 & ) * ddzu(k+1) * drho_air_zw(k) 188 190 ENDDO 189 191 … … 227 229 + 2.0_wp * ( & 228 230 km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) & 231 * rho_air(k+1) & 229 232 - km(k,j,i) * ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k) & 230 ) * ddzu(k+1) 233 * rho_air(k) & 234 ) * ddzu(k+1) * drho_air_zw(k) 231 235 ENDDO 232 236 ENDIF … … 246 250 247 251 USE arrays_3d, & 248 ONLY : ddzu, ddzw, km, tend, u, v, w 252 ONLY : ddzu, ddzw, km, tend, u, v, w, drho_air_zw, rho_air 249 253 250 254 USE control_parameters, & … … 316 320 & + 2.0_wp * ( & 317 321 & km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) & 322 & * rho_air(k+1) & 318 323 & - km(k,j,i) * ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k) & 319 & ) * ddzu(k+1) 324 & * rho_air(k) & 325 & ) * ddzu(k+1) * drho_air_zw(k) 320 326 ENDIF 321 327 ENDDO … … 361 367 + 2.0_wp * ( & 362 368 km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) & 369 * rho_air(k+1) & 363 370 - km(k,j,i) * ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k) & 364 ) * ddzu(k+1) 371 * rho_air(k) & 372 ) * ddzu(k+1) * drho_air_zw(k) 365 373 ENDIF 366 374 ENDDO … … 381 389 382 390 USE arrays_3d, & 383 ONLY : ddzu, ddzw, km, tend, u, v, w 391 ONLY : ddzu, ddzw, km, tend, u, v, w, drho_air_zw, rho_air 384 392 385 393 USE control_parameters, & … … 430 438 & + 2.0_wp * ( & 431 439 & km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) & 440 & * rho_air(k+1) & 432 441 & - km(k,j,i) * ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k) & 433 & ) * ddzu(k+1) 442 & * rho_air(k) & 443 & ) * ddzu(k+1) * drho_air_zw(k) 434 444 ENDDO 435 445 … … 485 495 + 2.0_wp * ( & 486 496 km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) & 497 * rho_air(k+1) & 487 498 - km(k,j,i) * ( w(k,j,i) - w(k-1,j,i) ) * ddzw(k) & 488 ) * ddzu(k+1) 499 * rho_air(k) & 500 ) * ddzu(k+1) * drho_air_zw(k) 489 501 ENDDO 490 502 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.