82 subroutine rrtmgp_sw_main_run(doSWrad, doSWclrsky, top_at_1, doGP_sgs_cnv, doGP_sgs_pbl, &
83 nCol, nDay, nLay, nGases, rrtmgp_phys_blksz, idx, icseed_sw, iovr, iovr_convcld, &
84 iovr_max, iovr_maxrand, iovr_rand, iovr_dcorr, iovr_exp, iovr_exprand, isubc_sw, &
85 iSFC, sfc_alb_nir_dir, sfc_alb_nir_dif, sfc_alb_uvvis_dir, sfc_alb_uvvis_dif, coszen,&
86 p_lay, p_lev, t_lay, t_lev, vmr_o2, vmr_h2o, vmr_o3, vmr_ch4, vmr_n2o, vmr_co2, &
87 cld_frac, cld_lwp, cld_reliq, cld_iwp, cld_reice, cld_swp, cld_resnow, cld_rwp, &
88 cld_rerain, precip_frac, cld_cnv_lwp, cld_cnv_reliq, cld_cnv_iwp, cld_cnv_reice, &
89 cld_pbl_lwp, cld_pbl_reliq, cld_pbl_iwp, cld_pbl_reice, cloud_overlap_param, &
90 active_gases_array, aersw_tau, aersw_ssa, aersw_g, solcon, scmpsw, &
91 fluxswUP_allsky, fluxswDOWN_allsky, fluxswUP_clrsky, fluxswDOWN_clrsky, cldtausw, &
95 logical,
intent(in) :: &
101 integer,
intent(in) :: &
117 integer,
intent(in),
dimension(:) :: &
119 integer,
intent(in),
dimension(:),
optional :: &
121 real(kind_phys),
dimension(:),
intent(in) :: &
127 real(kind_phys),
dimension(:,:),
intent(in),
optional :: &
138 real(kind_phys),
dimension(:,:),
intent(in) :: &
148 real(kind_phys),
dimension(:,:),
intent(in),
optional :: &
159 real(kind_phys),
dimension(:,:,:),
intent(in) :: &
163 character(len=*),
dimension(:),
intent(in),
optional :: &
165 real(kind_phys),
intent(in) :: &
169 character(len=*),
intent(out) :: &
171 integer,
intent(out) :: &
173 real(kind_phys),
dimension(:,:),
intent(inout) :: &
175 real(kind_phys),
dimension(:,:),
intent(inout),
optional :: &
180 type(
cmpfsw_type),
dimension(:),
intent(inout) :: &
190 type(
cmpfsw_type),
dimension(rrtmgp_phys_blksz) :: scmpsw_clrsky, scmpsw_allsky
191 type(ty_fluxes_byband) :: flux_allsky, flux_clrsky
192 real(kind_phys) :: tau_rain, tau_snow, ssa_rain, ssa_snow, asy_rain, asy_snow, &
193 tau_prec, asy_prec, ssa_prec, asyw, ssaw, za1, za2, flux_dir, flux_dif
194 real(kind_phys),
dimension(rrtmgp_phys_blksz) :: zcf0, zcf1
195 real(kind_dbl_prec),
dimension(sw_gas_props%get_ngpt()) :: rng1d
196 real(kind_dbl_prec),
dimension(sw_gas_props%get_ngpt(),nLay,rrtmgp_phys_blksz) :: rng3d,rng3d2
197 real(kind_dbl_prec),
dimension(sw_gas_props%get_ngpt()*nLay) :: rng2d
198 logical,
dimension(rrtmgp_phys_blksz,nLay,sw_gas_props%get_ngpt()) :: maskmcica
199 logical :: cloudy_column, clear_column
200 real(kind_phys),
dimension(sw_gas_props%get_nband(),rrtmgp_phys_blksz) :: &
201 sfc_alb_dir, sfc_alb_dif
202 real(kind_phys),
dimension(rrtmgp_phys_blksz,nLay+1,sw_gas_props%get_nband()),
target :: &
203 fluxsw_up_allsky, fluxsw_up_clrsky, fluxsw_dn_dir_clrsky, fluxsw_dn_allsky, &
204 fluxsw_dn_clrsky, fluxsw_dn_dir_allsky
205 integer :: iband, ibd, ibd_uv, icol, igas, ilay, ix, ix2, iblck
206 integer,
dimension(rrtmgp_phys_blksz) :: ipseed_sw, icols
208 real(kind_phys),
dimension(2,sw_gas_props%get_nband()) :: bandlimits
209 real(kind_phys),
dimension(2),
parameter :: &
210 nir_uvvis_bnd = (/12850,16000/), &
211 uvb_bnd = (/29000,38000/)
212 real(kind_phys),
dimension(rrtmgp_phys_blksz,sw_gas_props%get_ngpt()) :: toa_src_sw
214 type(ty_gas_concs) :: gas_concs
215 type(ty_optical_props_2str) :: sw_optical_props_accum, sw_optical_props_aerosol_local, &
216 sw_optical_props_cloudsbyband, sw_optical_props_cnvcloudsbyband, &
217 sw_optical_props_pblcloudsbyband, sw_optical_props_precipbyband, &
218 sw_optical_props_clouds
224 if (.not. doswrad)
return
227 call check_error_msg(
'rrtmgp_sw_main_gas_concs_init',gas_concs%init(active_gases_array))
230 call check_error_msg(
'rrtmgp_sw_main_accumulated_optics_init',&
231 sw_optical_props_accum%alloc_2str(rrtmgp_phys_blksz, nlay, sw_gas_props))
232 call check_error_msg(
'rrtmgp_sw_main_cloud_optics_init',&
233 sw_optical_props_cloudsbyband%alloc_2str(rrtmgp_phys_blksz, nlay, sw_gas_props%get_band_lims_wavenumber()))
234 call check_error_msg(
'rrtmgp_sw_main_precip_optics_init',&
235 sw_optical_props_precipbyband%alloc_2str(rrtmgp_phys_blksz, nlay, sw_gas_props%get_band_lims_wavenumber()))
236 call check_error_msg(
'rrtmgp_sw_mian_cloud_sampling_init', &
237 sw_optical_props_clouds%alloc_2str(rrtmgp_phys_blksz, nlay, sw_gas_props))
238 call check_error_msg(
'rrtmgp_sw_main_aerosol_optics_init',&
239 sw_optical_props_aerosol_local%alloc_2str(rrtmgp_phys_blksz, nlay, sw_gas_props%get_band_lims_wavenumber()))
240 if (dogp_sgs_cnv)
then
241 call check_error_msg(
'rrtmgp_sw_main_cnv_cloud_optics_init',&
242 sw_optical_props_cnvcloudsbyband%alloc_2str(rrtmgp_phys_blksz, nlay, sw_gas_props%get_band_lims_wavenumber()))
244 if (dogp_sgs_pbl)
then
245 call check_error_msg(
'rrtmgp_sw_main_pbl_cloud_optics_init',&
246 sw_optical_props_pblcloudsbyband%alloc_2str(rrtmgp_phys_blksz, nlay, sw_gas_props%get_band_lims_wavenumber()))
249 if (nday .gt. 0)
then
251 bandlimits = sw_gas_props%get_band_lims_wavenumber()
257 do icol=1,nday,rrtmgp_phys_blksz
260 icols = idx(icol:icol + rrtmgp_phys_blksz - 1)
263 zcf0(:) = 1._kind_phys
264 zcf1(:) = 1._kind_phys
265 do iblck = 1, rrtmgp_phys_blksz
267 zcf0(iblck) = min(zcf0(iblck), 1._kind_phys - cld_frac(icols(iblck),ilay))
269 if (zcf0(iblck) <= ftiny) zcf0(iblck) = 0._kind_phys
270 if (zcf0(iblck) > oneminus) zcf0(iblck) = 1._kind_phys
271 zcf1(iblck) = 1._kind_phys - zcf0(iblck)
273 cloudy_column = any(zcf1 .gt. eps)
274 clear_column = .true.
275 if (cloudy_column) clear_column = .false.
282 sw_optical_props_clouds%tau = 0._kind_phys
283 sw_optical_props_clouds%ssa = 0._kind_phys
284 sw_optical_props_clouds%g = 0._kind_phys
285 sw_optical_props_accum%tau = 0._kind_phys
286 sw_optical_props_accum%ssa = 0._kind_phys
287 sw_optical_props_accum%g = 0._kind_phys
288 sw_optical_props_cloudsbyband%tau = 0._kind_phys
289 sw_optical_props_cloudsbyband%ssa = 0._kind_phys
290 sw_optical_props_cloudsbyband%g = 0._kind_phys
291 sw_optical_props_precipbyband%tau = 0._kind_phys
292 sw_optical_props_precipbyband%ssa = 0._kind_phys
293 sw_optical_props_precipbyband%g = 0._kind_phys
294 if (dogp_sgs_cnv)
then
295 sw_optical_props_cnvcloudsbyband%tau = 0._kind_phys
296 sw_optical_props_cnvcloudsbyband%ssa = 0._kind_phys
297 sw_optical_props_cnvcloudsbyband%g = 0._kind_phys
299 if (dogp_sgs_pbl)
then
300 sw_optical_props_pblcloudsbyband%tau = 0._kind_phys
301 sw_optical_props_pblcloudsbyband%ssa = 0._kind_phys
302 sw_optical_props_pblcloudsbyband%g = 0._kind_phys
304 scmpsw_clrsky=
cmpfsw_type( 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 )
305 scmpsw_allsky=
cmpfsw_type( 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 )
306 cldtausw = 0._kind_phys
309 fluxsw_up_allsky = 0._kind_phys
310 fluxsw_dn_allsky = 0._kind_phys
311 fluxsw_dn_dir_allsky = 0._kind_phys
312 fluxsw_up_clrsky = 0._kind_phys
313 fluxsw_dn_clrsky = 0._kind_phys
314 flux_allsky%bnd_flux_up => fluxsw_up_allsky
315 flux_allsky%bnd_flux_dn => fluxsw_dn_allsky
316 flux_allsky%bnd_flux_dn_dir => fluxsw_dn_dir_allsky
317 flux_clrsky%bnd_flux_up => fluxsw_up_clrsky
318 flux_clrsky%bnd_flux_dn => fluxsw_dn_clrsky
325 call check_error_msg(
'rrtmgp_sw_main_set_vmr_o2', &
326 gas_concs%set_vmr(trim(active_gases_array(istr_o2)), vmr_o2(icols,:)))
327 call check_error_msg(
'rrtmgp_sw_main_set_vmr_co2', &
328 gas_concs%set_vmr(trim(active_gases_array(istr_co2)),vmr_co2(icols,:)))
329 call check_error_msg(
'rrtmgp_sw_main_set_vmr_ch4', &
330 gas_concs%set_vmr(trim(active_gases_array(istr_ch4)),vmr_ch4(icols,:)))
331 call check_error_msg(
'rrtmgp_sw_main_set_vmr_n2o', &
332 gas_concs%set_vmr(trim(active_gases_array(istr_n2o)),vmr_n2o(icols,:)))
333 call check_error_msg(
'rrtmgp_sw_main_set_vmr_h2o', &
334 gas_concs%set_vmr(trim(active_gases_array(istr_h2o)),vmr_h2o(icols,:)))
335 call check_error_msg(
'rrtmgp_sw_main_set_vmr_o3', &
336 gas_concs%set_vmr(trim(active_gases_array(istr_o3)), vmr_o3(icols,:)))
344 call check_error_msg(
'rrtmgp_sw_main_gas_optics',sw_gas_props%gas_optics(&
349 sw_optical_props_accum, &
353 do iblck = 1, rrtmgp_phys_blksz
354 toa_src_sw(iblck,:) = toa_src_sw(iblck,:)*solcon / sum(toa_src_sw(iblck,:))
366 do iblck = 1, rrtmgp_phys_blksz
367 do iband=1,sw_gas_props%get_nband()
368 if (bandlimits(1,iband) .lt. nir_uvvis_bnd(1))
then
369 sfc_alb_dir(iband,iblck) = sfc_alb_nir_dir(icols(iblck))
370 sfc_alb_dif(iband,iblck) = sfc_alb_nir_dif(icols(iblck))
372 if (bandlimits(1,iband) .eq. nir_uvvis_bnd(1))
then
373 sfc_alb_dir(iband,iblck) = 0.5_kind_phys*(sfc_alb_nir_dir(icols(iblck)) + &
374 sfc_alb_uvvis_dir(icols(iblck)))
375 sfc_alb_dif(iband,iblck) = 0.5_kind_phys*(sfc_alb_nir_dif(icols(iblck)) + &
376 sfc_alb_uvvis_dif(icols(iblck)))
379 if (bandlimits(1,iband) .ge. nir_uvvis_bnd(2))
then
380 sfc_alb_dir(iband,iblck) = sfc_alb_uvvis_dir(icols(iblck))
381 sfc_alb_dif(iband,iblck) = sfc_alb_uvvis_dif(icols(iblck))
383 if (bandlimits(1,iband) .eq. uvb_bnd(1)) ibd_uv = iband
392 if (cloudy_column)
then
394 call check_error_msg(
'rrtmgp_sw_main_cloud_optics',sw_cloud_props%cloud_optics(&
397 cld_reliq(icols,:), &
398 cld_reice(icols,:), &
399 sw_optical_props_cloudsbyband))
401 cldtausw(icols,:) = sw_optical_props_cloudsbyband%tau(:,:,11)
404 if (dogp_sgs_cnv)
then
406 call check_error_msg(
'rrtmgp_sw_main_cnv_cloud_optics',sw_cloud_props%cloud_optics(&
407 cld_cnv_lwp(icols,:), &
408 cld_cnv_iwp(icols,:), &
409 cld_cnv_reliq(icols,:), &
410 cld_cnv_reice(icols,:), &
411 sw_optical_props_cnvcloudsbyband))
414 call check_error_msg(
'rrtmgp_sw_main_increment_cnvclouds_to_clouds',&
415 sw_optical_props_cnvcloudsbyband%increment(sw_optical_props_cloudsbyband))
419 if (dogp_sgs_pbl)
then
421 call check_error_msg(
'rrtmgp_sw_main_pbl_cloud_optics',sw_cloud_props%cloud_optics(&
422 cld_pbl_lwp(icols,:), &
423 cld_pbl_iwp(icols,:), &
424 cld_pbl_reliq(icols,:), &
425 cld_pbl_reice(icols,:), &
426 sw_optical_props_pblcloudsbyband))
429 call check_error_msg(
'rrtmgp_sw_main_increment_pblclouds_to_clouds',&
430 sw_optical_props_pblcloudsbyband%increment(sw_optical_props_cloudsbyband))
434 do iblck = 1, rrtmgp_phys_blksz
436 if (cld_frac(icols(iblck),ilay) .gt. ftiny)
then
438 tau_rain = cld_rwp(icols(iblck),ilay)*a0r
439 if (cld_swp(icols(iblck),ilay) .gt. 0. .and. cld_resnow(icols(iblck),ilay) .gt. 10._kind_phys)
then
440 tau_snow = cld_swp(icols(iblck),ilay)*1.09087*(a0s + a1s/(1.0315*cld_resnow(icols(iblck),ilay)))
442 tau_snow = 0._kind_phys
446 do iband=1,sw_gas_props%get_nband()
448 ssa_rain = tau_rain*(1.-b0r(iband))
449 asy_rain = ssa_rain*c0r(iband)
450 ssa_snow = tau_snow*(1.-(b0s(iband)+b1s(iband)*1.0315*cld_resnow(icols(iblck),ilay)))
451 asy_snow = ssa_snow*c0s(iband)
453 tau_prec = max(1.e-12_kind_phys, tau_rain + tau_snow)
454 ssa_prec = max(1.e-12_kind_phys, ssa_rain + ssa_snow)
455 asy_prec = max(1.e-12_kind_phys, asy_rain + asy_snow)
456 asyw = asy_prec/max(1.e-12_kind_phys, ssa_prec)
457 ssaw = min(1._kind_phys-0.000001, ssa_prec/tau_prec)
460 sw_optical_props_precipbyband%tau(iblck,ilay,iband) = (1._kind_phys - za2) * tau_prec
461 sw_optical_props_precipbyband%ssa(iblck,ilay,iband) = (ssaw - za2) / (1._kind_phys - za2)
462 sw_optical_props_precipbyband%g(iblck,ilay,iband) = asyw/(1+asyw)
468 call check_error_msg(
'rrtmgp_sw_main_increment_precip_to_clouds',&
469 sw_optical_props_precipbyband%increment(sw_optical_props_cloudsbyband))
477 if(isubc_sw == 1)
then
478 do iblck = 1, rrtmgp_phys_blksz
479 ipseed_sw(iblck) = sw_gas_props%get_ngpt() + icols(iblck)
481 elseif (isubc_sw == 2)
then
482 do iblck = 1, rrtmgp_phys_blksz
483 ipseed_sw(iblck) = icseed_sw(icols(iblck))
488 do iblck = 1, rrtmgp_phys_blksz
491 if (iovr == iovr_max)
then
494 rng3d(:,ilay,iblck) = rng1d
499 rng3d(:,ilay,iblck) = rng1d
506 if (iovr == iovr_maxrand .or. iovr == iovr_rand .or. iovr == iovr_max)
then
507 call sampled_mask(real(rng3d, kind=kind_phys), cld_frac(icols,:), maskmcica)
510 if (iovr == iovr_dcorr)
then
511 do iblck = 1, rrtmgp_phys_blksz
515 rng3d2(:,:,iblck) = reshape(source = rng2d,shape=[sw_gas_props%get_ngpt(),nlay])
518 call sampled_mask(real(rng3d, kind=kind_phys), cld_frac(icols,:), maskmcica, &
519 overlap_param = cloud_overlap_param(icols,1:nlay-1), randoms2 = real(rng3d2, kind=kind_phys))
522 if (iovr == iovr_exp .or. iovr == iovr_exprand)
then
523 call sampled_mask(real(rng3d, kind=kind_phys), cld_frac(icols,:), maskmcica, &
524 overlap_param = cloud_overlap_param(icols,1:nlay-1))
527 call check_error_msg(
'rrtmgp_sw_main_cloud_sampling',&
528 draw_samples(maskmcica, .true., &
529 sw_optical_props_cloudsbyband, sw_optical_props_clouds))
538 sw_optical_props_aerosol_local%tau = aersw_tau(icols,:,:)
539 sw_optical_props_aerosol_local%ssa = aersw_ssa(icols,:,:)
540 sw_optical_props_aerosol_local%g = aersw_g(icols,:,:)
541 call check_error_msg(
'rrtmgp_sw_main_increment_aerosol_to_clrsky', &
542 sw_optical_props_aerosol_local%increment(sw_optical_props_accum))
545 if (clear_column .or. doswclrsky)
then
546 call check_error_msg(
'rrtmgp_sw_main_rte_sw_clrsky',rte_sw( &
547 sw_optical_props_accum, &
556 fluxswup_clrsky(icols,:) = sum(flux_clrsky%bnd_flux_up, dim=3)
557 fluxswdown_clrsky(icols,:) = sum(flux_clrsky%bnd_flux_dn, dim=3)
560 do iblck = 1, rrtmgp_phys_blksz
561 do iband=1,sw_gas_props%get_nband()
562 flux_dir = flux_clrsky%bnd_flux_dn(iblck,isfc,iband)
563 flux_dif = 0._kind_phys
565 if (iband < ibd)
then
566 scmpsw_clrsky(iblck)%nirbm = scmpsw_clrsky(iblck)%nirbm + flux_dir
567 scmpsw_clrsky(iblck)%nirdf = scmpsw_clrsky(iblck)%nirdf + flux_dif
570 if (iband == ibd)
then
571 scmpsw_clrsky(iblck)%nirbm = scmpsw_clrsky(iblck)%nirbm + flux_dir*0.5_kind_phys
572 scmpsw_clrsky(iblck)%nirdf = scmpsw_clrsky(iblck)%nirdf + flux_dif*0.5_kind_phys
573 scmpsw_clrsky(iblck)%visbm = scmpsw_clrsky(iblck)%visbm + flux_dir*0.5_kind_phys
574 scmpsw_clrsky(iblck)%visdf = scmpsw_clrsky(iblck)%visdf + flux_dif*0.5_kind_phys
577 if (iband > ibd)
then
578 scmpsw_clrsky(iblck)%visbm = scmpsw_clrsky(iblck)%visbm + flux_dir
579 scmpsw_clrsky(iblck)%visdf = scmpsw_clrsky(iblck)%visdf + flux_dif
582 scmpsw_clrsky(iblck)%uvbfc = flux_clrsky%bnd_flux_dn(iblck,isfc,ibd_uv)
586 fluxswup_clrsky(icols,:) = 0._kind_phys
587 fluxswdown_clrsky(icols,:) = 0._kind_phys
588 scmpsw =
cmpfsw_type( 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 )
596 if (cloudy_column)
then
601 call check_error_msg(
'rrtmgp_sw_main_increment_clouds_to_clrsky', &
602 sw_optical_props_clouds%increment(sw_optical_props_accum))
605 call check_error_msg(
'rrtmgp_sw_main_rte_sw_allsky',rte_sw( &
606 sw_optical_props_accum, &
615 fluxswup_allsky(icols,:) = sum(flux_allsky%bnd_flux_up, dim=3)
616 fluxswdown_allsky(icols,:) = sum(flux_allsky%bnd_flux_dn, dim=3)
619 do iblck = 1, rrtmgp_phys_blksz
621 do iband=1,sw_gas_props%get_nband()
622 flux_dir = flux_allsky%bnd_flux_dn_dir(iblck,isfc,iband)
623 flux_dif = flux_allsky%bnd_flux_dn(iblck,isfc,iband) - flux_allsky%bnd_flux_dn_dir(iblck,isfc,iband)
625 if (iband < ibd)
then
626 scmpsw_allsky(iblck)%nirbm = scmpsw_allsky(iblck)%nirbm + flux_dir
627 scmpsw_allsky(iblck)%nirdf = scmpsw_allsky(iblck)%nirdf + flux_dif
630 if (iband == ibd)
then
631 scmpsw_allsky(iblck)%nirbm = scmpsw_allsky(iblck)%nirbm + flux_dir*0.5_kind_phys
632 scmpsw_allsky(iblck)%nirdf = scmpsw_allsky(iblck)%nirdf + flux_dif*0.5_kind_phys
633 scmpsw_allsky(iblck)%visbm = scmpsw_allsky(iblck)%visbm + flux_dir*0.5_kind_phys
634 scmpsw_allsky(iblck)%visdf = scmpsw_allsky(iblck)%visdf + flux_dif*0.5_kind_phys
637 if (iband > ibd)
then
638 scmpsw_allsky(iblck)%visbm = scmpsw_allsky(iblck)%visbm + flux_dir
639 scmpsw_allsky(iblck)%visdf = scmpsw_allsky(iblck)%visdf + flux_dif
642 scmpsw_allsky(iblck)%uvbfc = flux_allsky%bnd_flux_dn(iblck,isfc,ibd_uv)
645 if (zcf1(iblck) .gt. eps)
then
646 scmpsw(icols(iblck))%nirbm = scmpsw_allsky(iblck)%nirbm
647 scmpsw(icols(iblck))%nirdf = scmpsw_allsky(iblck)%nirdf
648 scmpsw(icols(iblck))%visbm = scmpsw_allsky(iblck)%visbm
649 scmpsw(icols(iblck))%visdf = scmpsw_allsky(iblck)%visdf
650 scmpsw(icols(iblck))%uvbfc = flux_allsky%bnd_flux_dn(iblck,isfc,ibd_uv)
652 scmpsw(icols(iblck))%nirbm = scmpsw_clrsky(iblck)%nirbm
653 scmpsw(icols(iblck))%nirdf = scmpsw_clrsky(iblck)%nirdf
654 scmpsw(icols(iblck))%visbm = scmpsw_clrsky(iblck)%visbm
655 scmpsw(icols(iblck))%visdf = scmpsw_clrsky(iblck)%visdf
656 scmpsw(icols(iblck))%uvbfc = flux_clrsky%bnd_flux_dn(iblck,isfc,ibd_uv)
658 scmpsw(icols(iblck))%uvbf0 = flux_clrsky%bnd_flux_dn(iblck,isfc,ibd_uv)
661 fluxswup_allsky(icols,:) = sum(flux_clrsky%bnd_flux_up, dim=3)
662 fluxswdown_allsky(icols,:) = sum(flux_clrsky%bnd_flux_dn, dim=3)
663 do iblck = 1, rrtmgp_phys_blksz
664 scmpsw(icols(iblck))%nirbm = scmpsw_clrsky(iblck)%nirbm
665 scmpsw(icols(iblck))%nirdf = scmpsw_clrsky(iblck)%nirdf
666 scmpsw(icols(iblck))%visbm = scmpsw_clrsky(iblck)%visbm
667 scmpsw(icols(iblck))%visdf = scmpsw_clrsky(iblck)%visdf
668 scmpsw(icols(iblck))%uvbfc = flux_clrsky%bnd_flux_dn(iblck,isfc,ibd_uv)
669 scmpsw(icols(iblck))%uvbf0 = flux_clrsky%bnd_flux_dn(iblck,isfc,ibd_uv)
675 fluxswup_allsky(:,:) = 0._kind_phys
676 fluxswdown_allsky(:,:) = 0._kind_phys
677 fluxswup_clrsky(:,:) = 0._kind_phys
678 fluxswdown_clrsky(:,:) = 0._kind_phys
679 scmpsw =
cmpfsw_type( 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 )