Changeset 4302 for palm


Ignore:
Timestamp:
Nov 22, 2019 1:15:56 PM (4 years ago)
Author:
suehring
Message:

Omit tall canopy mapped on top of buildings, which may happen due to topography filtering

File:
1 edited

Legend:

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

    r4279 r4302  
    2727! -----------------
    2828! $Id$
     29! Omit tall canopy mapped on top of buildings
     30!
     31! 4279 2019-10-29 08:48:17Z scharf
    2932! unused variables removed
    3033!
     
    932935           ONLY: message_string, ocean_mode
    933936
     937       USE indices,                                                            &
     938           ONLY:  wall_flags_0
     939
    934940       USE netcdf_data_input_mod,                                              &
    935941           ONLY:  leaf_area_density_f
     942
     943       USE pegrid
    936944
    937945       USE surface_mod,                                                        &
     
    10461054          DO k = 1, pch_index
    10471055             lad(k) = 0.5 * ( pre_lad(k-1) + pre_lad(k) )
    1048           ENDDO         
     1056          ENDDO
    10491057
    10501058       ENDIF
     
    10531061!--    Allocate 3D-array for the leaf area density (lad_s).
    10541062       ALLOCATE( lad_s(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    1055 
    10561063!
    10571064!--    Initialize canopy parameters cdc (canopy drag coefficient),
     
    10961103                         lad_s(k,j,i) = leaf_area_density_f%var(k,j,i)
    10971104                      ENDDO
    1098                    ENDDO
    1099                 ENDDO
    1100                
     1105!
     1106!--                   Check if resolved vegetation is mapped onto buildings.
     1107!--                   Even if the user defines LAD only for non-building covered
     1108!--                   grid points, the filter process of topography may cause
     1109!--                   that LAD overlaps with buildings. This case, an
     1110!--                   informative message is given and the LAD field is set to
     1111!--                   zero at these
     1112                      IF ( ANY( lad_s(:,j,i) /= 0.0_wp )  .AND.                &
     1113                           ANY( BTEST( wall_flags_0(:,j,i), 6 ) ) )  THEN
     1114                         lad_s(:,j,i) = 0.0_wp
     1115                         WRITE( message_string, * )                            &
     1116                                          'Resolved plant-canopy is ' //       &
     1117                                          'is defined on top of a building ' //&
     1118                                          '- lad is omitted at this grid '   //&
     1119                                          'point: (i,j) = ', i, j
     1120                         CALL message( 'pcm_init', 'PA0313', 0, 0, myid, 6, 0 )
     1121                      ENDIF
     1122                   ENDDO
     1123                ENDDO
    11011124                CALL exchange_horiz( lad_s, nbgp )
    11021125!
Note: See TracChangeset for help on using the changeset viewer.