Changeset 4361 for palm/trunk/SOURCE/plant_canopy_model_mod.f90
- Timestamp:
- Jan 7, 2020 12:22:38 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/plant_canopy_model_mod.f90
r4360 r4361 27 27 ! ----------------- 28 28 ! $Id$ 29 ! - Remove unused arrays in pmc_rrd_local 30 ! - Remove one exchange of ghost points 31 ! 32 ! 4360 2020-01-07 11:25:50Z suehring 29 33 ! - Bugfix, read restart data for time-averaged pcm output quantities 30 34 ! - Output of plant-canopy quantities will fill values … … 1152 1156 pch_index_ji = 0 1153 1157 1154 DO i = nxl , nxr1155 DO j = nys , nyn1158 DO i = nxlg, nxrg 1159 DO j = nysg, nyng 1156 1160 DO k = 0, pch_index 1157 1161 IF ( lad_s(k,j,i) /= 0 ) pch_index_ji(j,i) = k … … 1168 1172 ENDDO 1169 1173 ENDDO 1170 1171 CALL exchange_horiz_2d_int( pch_index_ji, nys, nyn, nxl, nxr, nbgp )1172 1174 ! 1173 1175 !-- Calculate global pch_index value (index of top of plant canopy from ground) 1174 1176 pch_index = MAXVAL( pch_index_ji ) 1175 1176 1177 1177 ! 1178 1178 !-- Exchange pch_index from all processors … … 1475 1475 SUBROUTINE pcm_rrd_local( k, nxlf, nxlc, nxl_on_file, nxrf, nxrc, & 1476 1476 nxr_on_file, nynf, nync, nyn_on_file, nysf, & 1477 nysc, nys_on_file, tmp_2d, tmp_3d,found )1477 nysc, nys_on_file, found ) 1478 1478 1479 1479 INTEGER(iwp) :: k !< … … 1492 1492 1493 1493 LOGICAL, INTENT(OUT) :: found 1494 1495 REAL(wp), DIMENSION(nys_on_file-nbgp:nyn_on_file+nbgp, &1496 nxl_on_file-nbgp:nxr_on_file+nbgp) :: tmp_2d !< temporary 2D array1497 1498 REAL(wp), DIMENSION(nzb:nzt+1, &1499 nys_on_file-nbgp:nyn_on_file+nbgp, &1500 nxl_on_file-nbgp:nxr_on_file+nbgp) :: tmp_3d !< temporary 3D array with entire vertical dimension1501 1494 1502 1495 REAL(wp), DIMENSION(0:pch_index, &
Note: See TracChangeset
for help on using the changeset viewer.