Ignore:
Timestamp:
Mar 19, 2008 9:41:30 AM (14 years ago)
Author:
steinfeld
Message:

Update for the plant canopy model

File:
1 edited

Legend:

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

    r142 r153  
    44! Actual revisions:
    55! -----------------
     6! Consideration of heat sources within the forest canopy which represent the
     7! rate of heat input into the air from the forest leaves
     8!
     9! Evaluation of sinks and sources for scalar concentration due to canopy
     10! elements
     11!
    612! Bugfix: remove IF statement in plant_canopy_model_ij
    713!
     
    1521! Description:
    1622! ------------
    17 ! Evaluation of the drag due to vegetation
     23! Evaluation of sinks and sources of momentum, heat and scalar concentration
     24! due to canopy elements
    1825!------------------------------------------------------------------------------!
    1926
     
    5966                                              v(k,j,i)        +  &
    6067                                              v(k,j+1,i)      +  &
    61                                               v(k,j+1,i+1) )     &
     68                                              v(k,j+1,i-1) )     &
    6269                                            / 4.0 )**2        +  &
    6370                                          ( ( w(k-1,j,i-1)    +  &
     
    120127
    121128!
     129!--       potential temperature
     130          CASE ( 4 )
     131             DO  i = nxl, nxr
     132                DO  j = nys, nyn
     133                   DO  k = nzb_s_inner(j,i)+1, pch_index
     134                      tend(k,j,i) = tend(k,j,i) +                     &
     135                                    ( canopy_heat_flux(k,j,i) -     &
     136                                      canopy_heat_flux(k-1,j,i) ) / &
     137                                      dzw(k)
     138                   ENDDO
     139                ENDDO
     140             ENDDO
     141
     142!
     143!--       scalar concentration
     144          CASE ( 5 )
     145             DO  i = nxl, nxr
     146                DO  j = nys, nyn
     147                   DO  k = nzb_s_inner(j,i)+1, pch_index
     148                      tend(k,j,i) = tend(k,j,i) -                     &
     149                                    sec(k,j,i) * lad_s(k,j,i) *       &
     150                                    SQRT( ( ( u(k,j,i)        +       &
     151                                              u(k,j,i+1) )            &
     152                                            / 2.0 )**2        +       &
     153                                          ( ( v(k,j,i)        +       &
     154                                              v(k,j+1,i) )            &
     155                                            / 2.0 )**2        +       &
     156                                          ( ( w(k-1,j,i)      +       &
     157                                              w(k,j,i) )              &
     158                                            / 2.0 )**2 )      *       &
     159                                    ( q(k,j,i) - sls(k,j,i) )
     160                   ENDDO
     161                ENDDO
     162             ENDDO
     163
     164!
    122165!--       sgs-tke
    123           CASE ( 4 )
     166          CASE ( 6 )
    124167             DO  i = nxl, nxr
    125168                DO  j = nys, nyn
     
    180223                                        v(k,j,i)    +        &
    181224                                        v(k,j+1,i)  +        &
    182                                         v(k,j+1,i+1) )       &
     225                                        v(k,j+1,i-1) )       &
    183226                                      / 4.0 )**2    +        &
    184227                                    ( ( w(k-1,j,i-1) +       &
     
    232275
    233276!
     277!--       potential temperature
     278          CASE ( 4 )
     279             DO  k = nzb_s_inner(j,i)+1, pch_index
     280                tend(k,j,i) = tend(k,j,i) +                     &
     281                              ( canopy_heat_flux(k,j,i) -     &
     282                                canopy_heat_flux(k-1,j,i) ) / &
     283                                dzw(k)
     284             ENDDO
     285
     286
     287!
     288!--       scalar concentration
     289          CASE ( 5 )
     290             DO  k = nzb_s_inner(j,i)+1, pch_index
     291                tend(k,j,i) = tend(k,j,i) -                     &
     292                              sec(k,j,i) * lad_s(k,j,i) *       &
     293                              SQRT( ( ( u(k,j,i)        +       &
     294                                        u(k,j,i+1) )            &
     295                                      / 2.0 )**2        +       &
     296                                    ( ( v(k,j,i)        +       &
     297                                        v(k,j+1,i) )            &
     298                                      / 2.0 )**2        +       &
     299                                    ( ( w(k-1,j,i)      +       &
     300                                        w(k,j,i) )              &
     301                                      / 2.0 )**2 )      *       &
     302                              ( q(k,j,i) - sls(k,j,i) )
     303             ENDDO   
     304
     305!
    234306!--       sgs-tke
    235        CASE ( 4 )
     307       CASE ( 6 )
    236308          DO  k = nzb_s_inner(j,i)+1, pch_index   
    237309             tend(k,j,i) = tend(k,j,i) -                     &
Note: See TracChangeset for help on using the changeset viewer.