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/prognostic_equations.f90

    r139 r153  
    44! Actual revisions:
    55! -----------------
    6 !
     6! add call of plant_canopy_model in the prognostic equation for
     7! the potential temperature
     8! add call of plant_canopy_model in the prognostic equation for
     9! the concentration of a passive scalar
    710!
    811! Former revisions:
     
    422425             CALL impact_of_latent_heat( i, j )
    423426          ENDIF
     427
     428!
     429!--       Consideration of heat sources within the plant canopy
     430          IF ( plant_canopy .AND. ( cthf /= 0.0 ) ) THEN
     431             CALL plant_canopy_model( i, j, 4 )
     432          ENDIF
     433
    424434          CALL user_actions( i, j, 'pt-tendency' )
    425435
     
    612622                CALL calc_precipitation( i, j )
    613623             ENDIF
     624
     625!
     626!--          Sink or source of scalar concentration due to canopy elements
     627             IF ( plant_canopy ) CALL plant_canopy_model( i, j, 5 )
     628
    614629             CALL user_actions( i, j, 'q-tendency' )
    615630
     
    754769!
    755770!--          Additional sink term for flows through plant canopies
    756              IF ( plant_canopy )  CALL plant_canopy_model( i, j, 4 )         
     771             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 6 )         
    757772
    758773             CALL user_actions( i, j, 'e-tendency' )
     
    10341049             CALL impact_of_latent_heat( i, j )
    10351050          ENDIF
     1051
     1052!
     1053!--       Consideration of heat sources within the plant canopy
     1054          IF ( plant_canopy .AND. ( cthf /= 0.0 ) ) THEN
     1055             CALL plant_canopy_model( i, j, 4 )
     1056          ENDIF
     1057
    10361058          CALL user_actions( i, j, 'pt-tendency' )
    10371059
     
    11421164                CALL calc_precipitation( i, j )
    11431165             ENDIF
     1166
     1167!
     1168!--          Sink or source of scalar concentration due to canopy elements
     1169             IF ( plant_canopy ) CALL plant_canopy_model( i, j, 5 )
     1170
     1171
    11441172             CALL user_actions( i, j, 'q-tendency' )
    11451173
     
    12191247!
    12201248!--          Additional sink term for flows through plant canopies
    1221              IF ( plant_canopy )  CALL plant_canopy_model( i, j, 4 )
     1249             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 6 )
    12221250
    12231251             CALL user_actions( i, j, 'e-tendency' )
     
    15821610       CALL impact_of_latent_heat
    15831611    ENDIF
     1612
     1613!
     1614!-- Consideration of heat sources within the plant canopy
     1615    IF ( plant_canopy .AND. ( cthf /= 0.0 ) ) THEN
     1616       CALL plant_canopy_model( 4 )
     1617    ENDIF
     1618
     1619
    15841620    CALL user_actions( 'pt-tendency' )
    15851621
     
    17821818          CALL calc_precipitation
    17831819       ENDIF
     1820
     1821!
     1822!--    Sink or source of scalar concentration due to canopy elements
     1823       IF ( plant_canopy ) CALL plant_canopy_model( 5 )
     1824
    17841825       CALL user_actions( 'q-tendency' )
    17851826
     
    19191960!
    19201961!--    Additional sink term for flows through plant canopies
    1921        IF ( plant_canopy )  CALL plant_canopy_model( 4 )
     1962       IF ( plant_canopy )  CALL plant_canopy_model( 6 )
    19221963       CALL user_actions( 'e-tendency' )
    19231964
Note: See TracChangeset for help on using the changeset viewer.