Home


About us
PALM
PALM-4U
Installation
Release notes
Documentation
Getting Help
Gallery
Developers area

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

Anelastic approximation implemented

File:
1 edited

Legend:

Unmodified
Added
Removed

  • palm/trunk/SOURCE/diffusion_e.f90

    r2001 r2037  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Anelastic approximation implemented
    2323!
    2424! Former revisions:
    … …  
    128128
    129129       USE arrays_3d,                                                          &
    130            ONLY:  dd2zu, ddzu, ddzw, diss, e, km, l_grid, tend, zu, zw
     130           ONLY:  dd2zu, ddzu, ddzw, diss, e, km, l_grid, tend, zu, zw,        &
     131                  drho_air, rho_air_zw
    131132           
    132133       USE control_parameters,                                                 &
    … …  
    222223                                        + (                                    &
    223224               ( km(k,j,i)+km(k+1,j,i) ) * ( e(k+1,j,i)-e(k,j,i) ) * ddzu(k+1) &
     225                                                             * rho_air_zw(k)   &
    224226             - ( km(k,j,i)+km(k-1,j,i) ) * ( e(k,j,i)-e(k-1,j,i) ) * ddzu(k)   &
    225                                           ) * ddzw(k)                          &
     227                                                             * rho_air_zw(k-1) &
     228                                          ) * ddzw(k) * drho_air(k)            &
    226229                             - dissipation(k,j)
    227230
    … …  
    294297                                        + (                                    &
    295298               ( km(k,j,i)+km(k+1,j,i) ) * ( e(k+1,j,i)-e(k,j,i) ) * ddzu(k+1) &
     299                                                             * rho_air_zw(k)   &
    296300             - ( km(k,j,i)+km(k-1,j,i) ) * ( e(k,j,i)-e(k-1,j,i) ) * ddzu(k)   &
    297                                           ) * ddzw(k)                          &
     301                                                             * rho_air_zw(k-1) &
     302                                          ) * ddzw(k) * drho_air(k)            &
    298303                             - dissipation(k,j)
    299304
    … …  
    339344
    340345       USE arrays_3d,                                                          &
    341            ONLY:  dd2zu, ddzu, ddzw, diss, e, km, l_grid, tend, zu, zw
     346           ONLY:  dd2zu, ddzu, ddzw, diss, e, km, l_grid, tend, zu, zw,        &
     347                  drho_air, rho_air_zw
    342348         
    343349       USE control_parameters,                                                 &
    … …  
    430436                                        + (                                    &
    431437               ( km(k,j,i)+km(k+1,j,i) ) * ( e(k+1,j,i)-e(k,j,i) ) * ddzu(k+1) &
     438                                                             * rho_air_zw(k)   &
    432439             - ( km(k,j,i)+km(k-1,j,i) ) * ( e(k,j,i)-e(k-1,j,i) ) * ddzu(k)   &
    433                                           ) * ddzw(k)                          &
     440                                                             * rho_air_zw(k-1) &
     441                                          ) * ddzw(k) * drho_air(k)            &
    434442                                  - dissipation
    435443
    … …  
    497505                                        + (                                    &
    498506               ( km(k,j,i)+km(k+1,j,i) ) * ( e(k+1,j,i)-e(k,j,i) ) * ddzu(k+1) &
     507                                                             * rho_air_zw(k)   &
    499508             - ( km(k,j,i)+km(k-1,j,i) ) * ( e(k,j,i)-e(k-1,j,i) ) * ddzu(k)   &
    500                                           ) * ddzw(k)                          &
     509                                                             * rho_air_zw(k-1) &
     510                                          ) * ddzw(k) * drho_air(k)            &
    501511                                  - dissipation
    502512
    … …  
    542552
    543553       USE arrays_3d,                                                          &
    544            ONLY:  dd2zu, ddzu, ddzw, diss, e, km, l_grid, tend, zu, zw
     554           ONLY:  dd2zu, ddzu, ddzw, diss, e, km, l_grid, tend, zu, zw,        &
     555                  drho_air, rho_air_zw
    545556         
    546557       USE control_parameters,                                                 &
    … …  
    625636                                       + (                                    &
    626637              ( km(k,j,i)+km(k+1,j,i) ) * ( e(k+1,j,i)-e(k,j,i) ) * ddzu(k+1) &
     638                                                            * rho_air_zw(k)   &
    627639            - ( km(k,j,i)+km(k-1,j,i) ) * ( e(k,j,i)-e(k-1,j,i) ) * ddzu(k)   &
    628                                          ) * ddzw(k)                          &
     640                                                            * rho_air_zw(k-1) &
     641                                         ) * ddzw(k) * drho_air(k)            &
    629642                                       - dissipation(k)
    630643
Note: See TracChangeset for help on using the changeset viewer.
                                                                                                                                                                                                                                                                                                                                                                               
  | Impressum | ©Leibniz Universität Hannover |