Default WRF shortwave code

Improvement of WRF shortwave code

REAL :: diffuse_frac, corr_fac, csza_slp

GSW=0.0

bext340=5.E-6

bexth2o=5.E-6

SOLTOP=SOLCON

if (present(coszen)) then

csza=coszen

else

xt24=mod(xtime+radfr

*0.5,1440.)+eot

tloctm = gmt + xt24/60. + xlong/15.

hrang = 15. * (tloctm-

12.) * degrad

xxlat = xlat * degrad

csza = sin(xxlat) *

sin(declin) +

cos(xxlat) * cos(declin)

* cos(hrang)

end if

XMU=csza

SDOWN(1)=SOLTOP*XMU

GSW=(1-ALBEDO)*SDOWN(kte+1)

REAL:: diffuse_frac, corr_fac, csza_slp

GSW=0.0

bext340=5.E-6

bexth2o=5.E-6

SOLTOP=SOLCON

if (present(coszen)) then

csza=coszen

else

xt24=mod(xtime+radfrq

*0.5,1440.)+eot

tloctm = gmt + xt24/60. + xlong/15.

hrang = 15. * (tloctm-

12.) * degrad

xxlat = xlat * degrad

csza=sin(xxlat) * sin(declin)

+ cos(xxlat) *cos(declin)

* cos(hrang)

cszm=(cos(declin)*cos(xxlat)*tan(declin)+sin(xxlat)*cos(hrang))/cos(csza)

end if

XMU=csza

XMM=cszm

SDOWN(1)=SOLTOP*(1+0.03344*cos(0.9856*julian-2.72))*sin(XMU) - cos(XMM)

GSW=(1-ALBEDO)*SDOWN(kte+1)