Ignore:
Timestamp:
Jul 12, 2016 4:34:24 PM (5 years ago)
Author:
suehring
Message:

Separate balance equations for humidity and passive_scalar

File:
1 edited

Legend:

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

    r1954 r1960  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Separate humidity and passive scalar
    2222!
    2323! Former revisions:
     
    8686 
    8787    USE arrays_3d,                                                             &
    88         ONLY:  dzu, dzw, e, q, shf, tend, u, v, w, zu, zw
     88        ONLY:  dzu, dzw, e, q, s, shf, tend, u, v, w, zu, zw
    8989
    9090    USE indices,                                                               &
     
    818818
    819819!
    820 !--       scalar concentration
     820!--       humidity
    821821          CASE ( 5 )
    822822             DO  i = nxl, nxr
     
    866866                ENDDO
    867867             ENDDO
     868!
     869!--       scalar concentration
     870          CASE ( 7 )
     871             DO  i = nxl, nxr
     872                DO  j = nys, nyn
     873                   DO  k = nzb_s_inner(j,i)+1, nzb_s_inner(j,i)+pch_index
     874                      kk = k - nzb_s_inner(j,i)  !- lad arrays are defined flat
     875                      tend(k,j,i) = tend(k,j,i) -                              &
     876                                       lsec *                                  &
     877                                       lad_s(kk,j,i) *                         &
     878                                       SQRT( ( 0.5_wp * ( u(k,j,i) +           &
     879                                                          u(k,j,i+1) )         &
     880                                             )**2 +                            &
     881                                             ( 0.5_wp * ( v(k,j,i) +           &
     882                                                          v(k,j+1,i) )         &
     883                                             )**2 +                            &
     884                                             ( 0.5_wp * ( w(k-1,j,i) +         &
     885                                                          w(k,j,i) )           &
     886                                             )**2                              &
     887                                           ) *                                 &
     888                                       ( s(k,j,i) - lsc )
     889                   ENDDO
     890                ENDDO
     891             ENDDO   
     892
    868893
    869894
     
    11461171
    11471172!
    1148 !--       scalar concentration
     1173!--       humidity
    11491174          CASE ( 5 )
    11501175             DO  k = nzb_s_inner(j,i)+1, nzb_s_inner(j,i)+pch_index
     
    11861211                                 e(k,j,i)
    11871212             ENDDO
     1213             
     1214!
     1215!--       scalar concentration
     1216          CASE ( 7 )
     1217             DO  k = nzb_s_inner(j,i)+1, nzb_s_inner(j,i)+pch_index
     1218                kk = k - nzb_s_inner(j,i)  !- lad arrays are defined flat
     1219                tend(k,j,i) = tend(k,j,i) -                                    &
     1220                                 lsec *                                        &
     1221                                 lad_s(kk,j,i) *                               &
     1222                                 SQRT( ( 0.5_wp * ( u(k,j,i) +                 &
     1223                                                    u(k,j,i+1) )               &
     1224                                       )**2  +                                 &
     1225                                       ( 0.5_wp * ( v(k,j,i) +                 &
     1226                                                    v(k,j+1,i) )               &
     1227                                       )**2 +                                  &
     1228                                       ( 0.5_wp * ( w(k-1,j,i) +               &
     1229                                                    w(k,j,i) )                 &
     1230                                       )**2                                    &
     1231                                     ) *                                       &
     1232                                 ( s(k,j,i) - lsc )
     1233             ENDDO               
    11881234
    11891235       CASE DEFAULT
Note: See TracChangeset for help on using the changeset viewer.