Ignore:
Timestamp:
Oct 26, 2016 11:15:40 AM (5 years ago)
Author:
knoop
Message:

Anelastic approximation implemented

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/diffusion_w.f90

    r2001 r2037  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Anelastic approximation implemented
    2323!
    2424! Former revisions:
     
    119119
    120120       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
    122122           
    123123       USE control_parameters,                                                 &
     
    184184                      & + 2.0_wp * (                                           &
    185185                      &   km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1)  &
     186                      &               * rho_air(k+1)                           &
    186187                      & - 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)
    188190             ENDDO
    189191
     
    227229                                 + 2.0_wp * (                                  &
    228230                           km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) &
     231                                       * rho_air(k+1)                          &
    229232                         - 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)
    231235                ENDDO
    232236             ENDIF
     
    246250
    247251       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
    249253           
    250254       USE control_parameters,                                                 &
     
    316320                         & + 2.0_wp * (                                          &
    317321                         &   km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) &
     322                         &               * rho_air(k+1)                          &
    318323                         & - 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)
    320326                ENDIF
    321327             ENDDO
     
    361367                                 + 2.0_wp * (                                  &
    362368                           km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) &
     369                                       * rho_air(k+1)                          &
    363370                         - 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)
    365373                ENDIF
    366374             ENDDO
     
    381389
    382390       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
    384392           
    385393       USE control_parameters,                                                 &
     
    430438                      & + 2.0_wp * (                                           &
    431439                      &   km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1)  &
     440                      &               * rho_air(k+1)                           &
    432441                      & - 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)
    434444       ENDDO
    435445
     
    485495                                 + 2.0_wp * (                                  &
    486496                           km(k+1,j,i) * ( w(k+1,j,i) - w(k,j,i) ) * ddzw(k+1) &
     497                                       * rho_air(k+1)                          &
    487498                         - 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)
    489501          ENDDO
    490502       ENDIF
Note: See TracChangeset for help on using the changeset viewer.