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_s.f90

    r2001 r2037  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! Anelastic approximation implemented
    2323!
    2424! Former revisions:
     
    113113
    114114       USE arrays_3d,                                                          &
    115            ONLY:  ddzu, ddzw, kh, tend
     115           ONLY:  ddzu, ddzw, kh, tend, drho_air, rho_air_zw
    116116       
    117117       USE control_parameters,                                                 &
     
    190190                                       + 0.5_wp * (                            &
    191191            ( 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)    &
    192193          - ( 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)
    194196             ENDDO
    195197
     
    205207                                                  * ( s(k+1,j,i)-s(k,j,i) )    &
    206208                                                  * ddzu(k+1)                  &
     209                                                  * rho_air_zw(k)              &
    207210                                           + s_flux_b(j,i)                     &
    208                                          ) * ddzw(k)
     211                                         ) * ddzw(k) * drho_air(k)
    209212
    210213             ENDIF
     
    222225                                                    * ( s(k,j,i)-s(k-1,j,i) )  &
    223226                                                    * ddzu(k)                  &
    224                                          ) * ddzw(k)
     227                                                    * rho_air_zw(k-1)          &
     228                                         ) * ddzw(k) * drho_air(k)
    225229
    226230             ENDIF
     
    240244
    241245       USE arrays_3d,                                                          &
    242            ONLY:  ddzu, ddzw, kh, tend
     246           ONLY:  ddzu, ddzw, kh, tend, drho_air, rho_air_zw
    243247           
    244248       USE control_parameters,                                                 &
     
    324328                                       + 0.5_wp * (                            &
    325329            ( 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)    &
    326331          - ( 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)
    328334                ENDIF
    329335             ENDDO
     
    338344                                                     * ( s(k+1,j,i)-s(k,j,i) ) &
    339345                                                     * ddzu(k+1)               &
     346                                                     * rho_air_zw(k)           &
    340347                                              + s_flux_b(j,i)                  &
    341                                             ) * ddzw(k)
     348                                            ) * ddzw(k) * drho_air(k)
    342349                ENDIF
    343350
     
    351358                                                       * ( s(k,j,i)-s(k-1,j,i) )  &
    352359                                                       * ddzu(k)                  &
    353                                             ) * ddzw(k)
     360                                                       * rho_air_zw(k-1)          &
     361                                            ) * ddzw(k) * drho_air(k)
    354362                ENDIF
    355363             ENDDO
     
    370378
    371379       USE arrays_3d,                                                          &
    372            ONLY:  ddzu, ddzw, kh, tend
     380           ONLY:  ddzu, ddzw, kh, tend, drho_air, rho_air_zw
    373381           
    374382       USE control_parameters,                                                 &
     
    446454                                       + 0.5_wp * (                            &
    447455            ( 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)    &
    448457          - ( 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)
    450460       ENDDO
    451461
     
    459469                                               * ( s(k+1,j,i)-s(k,j,i) )       &
    460470                                               * ddzu(k+1)                     &
     471                                               * rho_air_zw(k)                 &
    461472                                        + s_flux_b(j,i)                        &
    462                                       ) * ddzw(k)
     473                                      ) * ddzw(k) * drho_air(k)
    463474
    464475       ENDIF
     
    474485                                               * ( s(k,j,i)-s(k-1,j,i) )       &
    475486                                               * ddzu(k)                       &
    476                                       ) * ddzw(k)
     487                                               * rho_air_zw(k-1)               &
     488                                      ) * ddzw(k) * drho_air(k)
    477489
    478490       ENDIF
Note: See TracChangeset for help on using the changeset viewer.