Ignore:
Timestamp:
Nov 28, 2007 10:03:58 AM (14 years ago)
Author:
letzel
Message:

Plant canopy model of Watanabe (2004,BLM 112,307-341) added.

File:
1 edited

Legend:

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

    r129 r138  
    44! Actual revisions:
    55! -----------------
     6! add call of subroutines that evaluate the canopy drag terms
    67! add wall_*flux to parameter list of calls of diffusion_s
    78!
     
    8081    USE diffusion_w_mod
    8182    USE impact_of_latent_heat_mod
     83    USE plant_canopy_model_mod
    8284    USE production_e_mod
    8385    USE user_actions_mod
     
    160162          CALL coriolis( i, j, 1 )
    161163          IF ( sloping_surface )  CALL buoyancy( i, j, pt, pt_reference, 1, 4 )
     164
     165!
     166!--       Drag by plant canopy
     167          IF ( plant_canopy )  CALL plant_canopy_model( i, j, 1 )
    162168          CALL user_actions( i, j, 'u-tendency' )
    163169
     
    227233          ENDIF
    228234          CALL coriolis( i, j, 2 )
     235
     236!
     237!--       Drag by plant canopy
     238          IF ( plant_canopy )  CALL plant_canopy_model( i, j, 2 )     
     239
    229240          CALL user_actions( i, j, 'v-tendency' )
    230241
     
    303314             ENDIF
    304315          ENDIF
     316
     317!
     318!--       Drag by plant canopy
     319          IF ( plant_canopy )  CALL plant_canopy_model( i, j, 3 )
     320
    305321          CALL user_actions( i, j, 'w-tendency' )
    306322
     
    732748             ENDIF
    733749             CALL production_e( i, j )
     750
     751!
     752!--          Additional sink term for flows through plant canopies
     753             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 4 )         
     754
    734755             CALL user_actions( i, j, 'e-tendency' )
    735756
     
    833854             IF ( sloping_surface )  CALL buoyancy( i, j, pt, pt_reference, 1, &
    834855                                                    4 )
     856
     857!
     858!--          Drag by plant canopy
     859             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 1 )
     860
    835861             CALL user_actions( i, j, 'u-tendency' )
    836862
     
    882908             ENDIF
    883909             CALL coriolis( i, j, 2 )
     910
     911!
     912!--          Drag by plant canopy
     913             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 2 )       
     914
    884915             CALL user_actions( i, j, 'v-tendency' )
    885916
     
    938969             ENDIF
    939970          ENDIF
     971
     972!
     973!--       Drag by plant canopy
     974          IF ( plant_canopy )  CALL plant_canopy_model( i, j, 3 )
     975
    940976          CALL user_actions( i, j, 'w-tendency' )
    941977
     
    11771213             ENDIF
    11781214             CALL production_e( i, j )
     1215
     1216!
     1217!--          Additional sink term for flows through plant canopies
     1218             IF ( plant_canopy )  CALL plant_canopy_model( i, j, 4 )
     1219
    11791220             CALL user_actions( i, j, 'e-tendency' )
    11801221
     
    12691310    CALL coriolis( 1 )
    12701311    IF ( sloping_surface )  CALL buoyancy( pt, pt_reference, 1, 4 )
     1312
     1313!
     1314!-- Drag by plant canopy
     1315    IF ( plant_canopy )  CALL plant_canopy_model( 1 )
     1316
    12711317    CALL user_actions( 'u-tendency' )
    12721318
     
    13401386    ENDIF
    13411387    CALL coriolis( 2 )
     1388
     1389!
     1390!-- Drag by plant canopy
     1391    IF ( plant_canopy )  CALL plant_canopy_model( 2 )
    13421392    CALL user_actions( 'v-tendency' )
    13431393
     
    14201470       ENDIF
    14211471    ENDIF
     1472
     1473!
     1474!-- Drag by plant canopy
     1475    IF ( plant_canopy )  CALL plant_canopy_model( 3 )
     1476
    14221477    CALL user_actions( 'w-tendency' )
    14231478
     
    18581913       ENDIF
    18591914       CALL production_e
     1915
     1916!
     1917!--    Additional sink term for flows through plant canopies
     1918       IF ( plant_canopy )  CALL plant_canopy_model( 4 )
    18601919       CALL user_actions( 'e-tendency' )
    18611920
Note: See TracChangeset for help on using the changeset viewer.