238 |
enddo |
enddo |
239 |
enddo |
enddo |
240 |
enddo |
enddo |
241 |
ncldras = ncldras + 1 |
if ( (bi.eq.1) .and. (bj.eq.1) ) ncldras = ncldras + 1 |
242 |
endif |
endif |
243 |
|
|
244 |
if( icldtot.gt.0 ) then |
if( icldtot.gt.0 ) then |
250 |
enddo |
enddo |
251 |
enddo |
enddo |
252 |
enddo |
enddo |
253 |
ncldtot = ncldtot + 1 |
if ( (bi.eq.1) .and. (bj.eq.1) ) ncldtot = ncldtot + 1 |
254 |
endif |
endif |
255 |
|
|
256 |
if( icldlow.gt.0 ) then |
if( icldlow.gt.0 ) then |
259 |
qdiag(i,j,icldlow,bi,bj) = qdiag(i,j,icldlow,bi,bj) + cldlow(i,j) |
qdiag(i,j,icldlow,bi,bj) = qdiag(i,j,icldlow,bi,bj) + cldlow(i,j) |
260 |
enddo |
enddo |
261 |
enddo |
enddo |
262 |
ncldlow = ncldlow + 1 |
if ( (bi.eq.1) .and. (bj.eq.1) ) ncldlow = ncldlow + 1 |
263 |
endif |
endif |
264 |
|
|
265 |
if( icldmid.gt.0 ) then |
if( icldmid.gt.0 ) then |
268 |
qdiag(i,j,icldmid,bi,bj) = qdiag(i,j,icldmid,bi,bj) + cldmid(i,j) |
qdiag(i,j,icldmid,bi,bj) = qdiag(i,j,icldmid,bi,bj) + cldmid(i,j) |
269 |
enddo |
enddo |
270 |
enddo |
enddo |
271 |
ncldmid = ncldmid + 1 |
if ( (bi.eq.1) .and. (bj.eq.1) ) ncldmid = ncldmid + 1 |
272 |
endif |
endif |
273 |
|
|
274 |
if( icldhi.gt.0 ) then |
if( icldhi.gt.0 ) then |
277 |
qdiag(i,j,icldhi,bi,bj) = qdiag(i,j,icldhi,bi,bj) + cldhi(i,j) |
qdiag(i,j,icldhi,bi,bj) = qdiag(i,j,icldhi,bi,bj) + cldhi(i,j) |
278 |
enddo |
enddo |
279 |
enddo |
enddo |
280 |
ncldhi = ncldhi + 1 |
if ( (bi.eq.1) .and. (bj.eq.1) ) ncldhi = ncldhi + 1 |
281 |
endif |
endif |
282 |
|
|
283 |
if( ilzrad.gt.0 ) then |
if( ilzrad.gt.0 ) then |
289 |
enddo |
enddo |
290 |
enddo |
enddo |
291 |
enddo |
enddo |
292 |
nlzrad = nlzrad + 1 |
if ( (bi.eq.1) .and. (bj.eq.1) ) nlzrad = nlzrad + 1 |
293 |
endif |
endif |
294 |
|
|
295 |
c Albedo Diagnostics |
c Albedo Diagnostics |
301 |
. albvisdr(i,j) |
. albvisdr(i,j) |
302 |
enddo |
enddo |
303 |
enddo |
enddo |
304 |
nalbvisdr = nalbvisdr + 1 |
if ( (bi.eq.1) .and. (bj.eq.1) ) nalbvisdr = nalbvisdr + 1 |
305 |
endif |
endif |
306 |
|
|
307 |
if( ialbvisdf.gt.0 ) then |
if( ialbvisdf.gt.0 ) then |
311 |
. albvisdf(i,j) |
. albvisdf(i,j) |
312 |
enddo |
enddo |
313 |
enddo |
enddo |
314 |
nalbvisdf = nalbvisdf + 1 |
if ( (bi.eq.1) .and. (bj.eq.1) ) nalbvisdf = nalbvisdf + 1 |
315 |
endif |
endif |
316 |
|
|
317 |
if( ialbnirdr.gt.0 ) then |
if( ialbnirdr.gt.0 ) then |
321 |
. albnirdr(i,j) |
. albnirdr(i,j) |
322 |
enddo |
enddo |
323 |
enddo |
enddo |
324 |
nalbnirdr = nalbnirdr + 1 |
if ( (bi.eq.1) .and. (bj.eq.1) ) nalbnirdr = nalbnirdr + 1 |
325 |
endif |
endif |
326 |
|
|
327 |
if( ialbnirdf.gt.0 ) then |
if( ialbnirdf.gt.0 ) then |
331 |
. albnirdf(i,j) |
. albnirdf(i,j) |
332 |
enddo |
enddo |
333 |
enddo |
enddo |
334 |
nalbnirdf = nalbnirdf + 1 |
if ( (bi.eq.1) .and. (bj.eq.1) ) nalbnirdf = nalbnirdf + 1 |
335 |
endif |
endif |
336 |
|
|
337 |
C Compute Optical Thicknesses and Diagnostics |
C Compute Optical Thicknesses and Diagnostics |
356 |
enddo |
enddo |
357 |
enddo |
enddo |
358 |
enddo |
enddo |
359 |
ntauave = ntauave + 1 |
if ( (bi.eq.1) .and. (bj.eq.1) ) ntauave = ntauave + 1 |
360 |
endif |
endif |
361 |
|
|
362 |
if( itaucld.gt.0 ) then |
if( itaucld.gt.0 ) then |