@@ -283,60 +283,63 @@ SUBROUTINE MISR_simulator(
283283
284284 !
285285 ! Modify MISR CTH for satellite spatial / pattern matcher effects
286+ !
287+ ! Code in this region added by roj 5 / 2006 to account
288+ ! for spatial effect of the MISR pattern matcher.
289+ ! Basically, if a column is found between two neighbors
290+ ! at the same CTH, and that column has no hieght or
291+ ! a lower CTH, THEN misr will tend to but place the
292+ ! odd column at the same height as it neighbors.
293+ !
294+ ! This setup assumes the columns represent a about a 1 to 4 km scale
295+ ! it will need to be modified significantly, otherwise
296+
297+ ! Commented out . Should only be used with subcolumns are not randomly ordered. Roj 4 / 2018
286298 !
287- ! Code in this region added by roj 5 / 2006 to account
288- ! for spatial effect of the MISR pattern matcher.
289- ! Basically, if a column is found between two neighbors
290- ! at the same CTH, and that column has no hieght or
291- ! a lower CTH, THEN misr will tend to but place the
292- ! odd column at the same height as it neighbors.
299+ ! if (ncol.eq. 1 ) then
300+ !
301+ ! adjust based on neightboring points ... i.e. only 2D grid was input
302+ ! do j= 2 ,npoints-1
303+ !
304+ ! if (box_MISR_ztop(j-1 ,1 ).gt. 0 .and.
305+ ! & box_MISR_ztop(j+1 ,1 ).gt. 0 ) then
293306 !
294- ! This setup assumes the columns represent a about a 1 to 4 km scale
295- ! it will need to be modified significantly, otherwise
296- if (ncol.eq. 1 ) then
297-
298- ! adjust based on neightboring points ... i.e. only 2D grid was input
299- do j= 2 ,npoints-1
300-
301- if (box_MISR_ztop(j-1 ,1 ).gt. 0 .and.
302- & box_MISR_ztop(j+1 ,1 ).gt. 0 ) then
303-
304- if ( abs ( box_MISR_ztop(j-1 ,1 ) -
305- & box_MISR_ztop(j+1 ,1 ) ) .lt. 500
306- & .and.
307- & box_MISR_ztop(j,1 ) .lt.
308- & box_MISR_ztop(j+1 ,1 ) ) then
309-
310- box_MISR_ztop(j,1 ) =
311- & box_MISR_ztop(j+1 ,1 )
312- endif
313-
314- endif
315- enddo
316- else
317-
307+ ! if ( abs ( box_MISR_ztop(j-1 ,1 ) -
308+ ! & box_MISR_ztop(j+1 ,1 ) ) .lt. 500
309+ ! & .and.
310+ ! & box_MISR_ztop(j,1 ) .lt.
311+ ! & box_MISR_ztop(j+1 ,1 ) ) then
312+ !
313+ ! box_MISR_ztop(j,1 ) =
314+ ! & box_MISR_ztop(j+1 ,1 )
315+ ! endif
316+ !
317+ ! endif
318+ ! enddo
319+ ! else
320+ !
318321 ! adjust based on neighboring subcolumns ....
319- do ibox= 2 ,ncol-1
320-
321- if (box_MISR_ztop(1 ,ibox-1 ).gt. 0 .and.
322- & box_MISR_ztop(1 ,ibox+1 ).gt. 0 ) then
323-
324- if ( abs ( box_MISR_ztop(1 ,ibox-1 ) -
325- & box_MISR_ztop(1 ,ibox+1 ) ) .lt. 500
326- & .and.
327- & box_MISR_ztop(1 ,ibox) .lt.
328- & box_MISR_ztop(1 ,ibox+1 ) ) then
329-
330- box_MISR_ztop(1 ,ibox) =
331- & box_MISR_ztop(1 ,ibox+1 )
332- endif
333-
334- endif
335- enddo
336-
337- endif
322+ ! do ibox= 2 ,ncol-1
323+ !
324+ ! if (box_MISR_ztop(1 ,ibox-1 ).gt. 0 .and.
325+ ! & box_MISR_ztop(1 ,ibox+1 ).gt. 0 ) then
326+ !
327+ ! if ( abs ( box_MISR_ztop(1 ,ibox-1 ) -
328+ ! & box_MISR_ztop(1 ,ibox+1 ) ) .lt. 500
329+ ! & .and.
330+ ! & box_MISR_ztop(1 ,ibox) .lt.
331+ ! & box_MISR_ztop(1 ,ibox+1 ) ) then
332+ !
333+ ! box_MISR_ztop(1 ,ibox) =
334+ ! & box_MISR_ztop(1 ,ibox+1 )
335+ ! endif
336+ !
337+ ! endif
338+ ! enddo
339+ !
340+ ! endif
338341
339- !
342+ !
340343 ! DETERMINE CLOUD TYPE FREQUENCIES
341344 !
342345 ! Now that ztop and tau have been determined,
0 commit comments