Changeset 4302 for palm/trunk/SOURCE/plant_canopy_model_mod.f90
- Timestamp:
- Nov 22, 2019 1:15:56 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/plant_canopy_model_mod.f90
r4279 r4302 27 27 ! ----------------- 28 28 ! $Id$ 29 ! Omit tall canopy mapped on top of buildings 30 ! 31 ! 4279 2019-10-29 08:48:17Z scharf 29 32 ! unused variables removed 30 33 ! … … 932 935 ONLY: message_string, ocean_mode 933 936 937 USE indices, & 938 ONLY: wall_flags_0 939 934 940 USE netcdf_data_input_mod, & 935 941 ONLY: leaf_area_density_f 942 943 USE pegrid 936 944 937 945 USE surface_mod, & … … 1046 1054 DO k = 1, pch_index 1047 1055 lad(k) = 0.5 * ( pre_lad(k-1) + pre_lad(k) ) 1048 ENDDO 1056 ENDDO 1049 1057 1050 1058 ENDIF … … 1053 1061 !-- Allocate 3D-array for the leaf area density (lad_s). 1054 1062 ALLOCATE( lad_s(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 1055 1056 1063 ! 1057 1064 !-- Initialize canopy parameters cdc (canopy drag coefficient), … … 1096 1103 lad_s(k,j,i) = leaf_area_density_f%var(k,j,i) 1097 1104 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 1101 1124 CALL exchange_horiz( lad_s, nbgp ) 1102 1125 !
Note: See TracChangeset
for help on using the changeset viewer.