Ignore:
Timestamp:
Oct 7, 2019 1:29:08 PM (2 years ago)
Author:
suehring
Message:

Land-surface model: Revise limitation for soil moisture in case it exceeds its saturation value; Revise initialization of soil moisture and temperature in a nested run in case dynamic input information is available. This case, the soil within the child domains can be initialized separately; As part of this revision, migrate the netcdf input of soil temperature / moisture to this module, as well as the routine to inter/extrapolate soil profiles between different grids.; Plant-canopy: Check if any LAD is prescribed when plant-canopy model is applied, in order to avoid crashes in the radiation transfer model; Surface-layer fluxes: Initialization of Obukhov length also at vertical surfaces (if allocated); Urban-surface model: Add checks to ensure that relative fractions of walls, windowns and green surfaces sum-u to one; Revise message calls dealing with local checks

File:
1 edited

Legend:

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

    r4226 r4258  
    2727! -----------------
    2828! $Id$
     29! Check if any LAD is prescribed when plant-canopy model is applied.
     30!
     31! 4226 2019-09-10 17:03:24Z suehring
    2932! Bugfix, missing initialization of heating rate
    3033!
     
    941944       INTEGER(iwp) ::  m   !< running index
    942945
    943        REAL(wp) ::  int_bpdf        !< vertical integral for lad-profile construction
     946       REAL(wp) ::  canopy_height   !< canopy height for lad-profile construction
    944947       REAL(wp) ::  gradient        !< gradient for lad-profile construction
    945        REAL(wp) ::  canopy_height   !< canopy height for lad-profile construction
     948       REAL(wp) ::  int_bpdf        !< vertical integral for lad-profile construction     
     949       REAL(wp) ::  lad_max         !< maximum LAD value in the model domain, used to perform a check
    946950
    947951       IF ( debug_output )  CALL debug_message( 'pcm_init', 'start' )
     
    11151119 
    11161120       END SELECT
     1121!
     1122!--    Check that at least one grid point has an LAD /= 0, else this may
     1123!--    cause errors in the radiation model.
     1124       lad_max = MAXVAL( lad_s )
     1125#if defined( __parallel )
     1126       CALL MPI_ALLREDUCE( MPI_IN_PLACE, lad_max, 1, MPI_REAL, MPI_MAX,        &
     1127                           comm2d, ierr)
     1128#endif
     1129       IF ( lad_max <= 0.0_wp )  THEN
     1130          message_string = 'Plant-canopy model is switched-on but no ' //      &
     1131                           'plant canopy is present in the model domain.'
     1132          CALL message( 'pcm_init', 'PA0685', 1, 2, 0, 6, 0 )
     1133       ENDIF
     1134   
    11171135!
    11181136!--    Initialize 2D index array indicating canopy top index.
     
    13011319
    13021320!
    1303 !--    Try to find radiation model package
     1321!--    Try to find plant-canopy model package
    13041322       REWIND ( 11 )
    13051323       line = ' '
     
    13141332
    13151333!
    1316 !--    Set flag that indicates that the radiation model is switched on
     1334!--    Set flag that indicates that the plant-canopy model is switched on
    13171335       plant_canopy = .TRUE.
    13181336
     
    13411359
    13421360!
    1343 !--    Set flag that indicates that the radiation model is switched on
     1361!--    Set flag that indicates that the plant-canopy model is switched on
    13441362       plant_canopy = .TRUE.
    13451363
Note: See TracChangeset for help on using the changeset viewer.