Changeset 3022 for palm/trunk/SOURCE/plant_canopy_model_mod.f90
- Timestamp:
- May 18, 2018 11:12:35 AM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/plant_canopy_model_mod.f90
r3014 r3022 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Bugfix in allocation of transpiration rate 28 ! 29 ! 3014 2018-05-09 08:42:38Z maronga 27 30 ! Bugfix: nzb_do and nzt_do were not used for 3d data output 28 31 ! Added pc_transpiration_rate … … 894 897 895 898 ALLOCATE( cum_lai_hf(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 896 pc_heating_rate(nzb:nzt+1,nysg:nyng,nxlg:nxrg), & 897 pc_transpiration_rate(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 899 pc_heating_rate(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 898 900 ! 899 901 !-- Piecewise calculation of the cumulative leaf area index by vertical … … 962 964 DO k = 1, pch_index_ji(j,i) 963 965 IF ( cum_lai_hf(0,j,i) /= 0.0_wp ) THEN 964 pc_heating_rate(k,j,i) = cthf * 965 ( exp(-ext_coef*cum_lai_hf(k,j,i)) - 966 pc_heating_rate(k,j,i) = cthf * & 967 ( exp(-ext_coef*cum_lai_hf(k,j,i)) - & 966 968 exp(-ext_coef*cum_lai_hf(k-1,j,i) ) ) / dzw(k) 967 969 ENDIF … … 971 973 972 974 ENDIF 975 ! 976 !-- Allocate transpiration rate 977 IF ( humidity ) & 978 ALLOCATE( pc_transpiration_rate(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ) 973 979 974 980
Note: See TracChangeset
for help on using the changeset viewer.