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) |