wrudisill_2357
New member
Hello,
This question does not refer specifically to WRF-Solar but regular WRF, but since it's related to radiation I might find a good audience here.
I'm trying to understand the equation that WRF uses to compute the diffuse radiation component of shortwave. Specifically his line in module_ra_sw.
https://github.com/wrf-model/WRF/blob/eed56d74b865af4ce9f98ea029226acfe52b0569/phys/module_ra_sw.F#L495
The raw fortran code is:
diffuse_frac = min(1.,1/(max(0.1,2.1-2.8*log(log(SDOWN(kts)/max(SDOWN(kte+1),1.e-3))))))
The code comment says "! Parameterize diffuse fraction of global solar radiation as a function of the ratio between TOA radiation and surface global radiation".
Two questions:
1) Does anyone know where this parameterization comes from? I haven't been able to find a source for this.
2) Can anyone explain the max(SDOWN(kts) term? Is this the "surface global radiation"? I don't fully understand the timesteps either (kte/kts). I think kte is the current model timestep. Maybe kts is the starting timestep(?)
The reason that I'm asking this is that I want to further downscale WRF radiation outputs (SWDOMW) to a finer resolution grid to do offline-land modelling. So I want to re-compute terrain normal shortwave radiation (SWNORM) given the fine resolution topography, and I'd like to use the same method that WRF uses (meaning I need to compute the diffuse frac).
Thanks and I hope that made sense,
-Wil
This question does not refer specifically to WRF-Solar but regular WRF, but since it's related to radiation I might find a good audience here.
I'm trying to understand the equation that WRF uses to compute the diffuse radiation component of shortwave. Specifically his line in module_ra_sw.
https://github.com/wrf-model/WRF/blob/eed56d74b865af4ce9f98ea029226acfe52b0569/phys/module_ra_sw.F#L495
The raw fortran code is:
diffuse_frac = min(1.,1/(max(0.1,2.1-2.8*log(log(SDOWN(kts)/max(SDOWN(kte+1),1.e-3))))))
The code comment says "! Parameterize diffuse fraction of global solar radiation as a function of the ratio between TOA radiation and surface global radiation".
Two questions:
1) Does anyone know where this parameterization comes from? I haven't been able to find a source for this.
2) Can anyone explain the max(SDOWN(kts) term? Is this the "surface global radiation"? I don't fully understand the timesteps either (kte/kts). I think kte is the current model timestep. Maybe kts is the starting timestep(?)
The reason that I'm asking this is that I want to further downscale WRF radiation outputs (SWDOMW) to a finer resolution grid to do offline-land modelling. So I want to re-compute terrain normal shortwave radiation (SWNORM) given the fine resolution topography, and I'd like to use the same method that WRF uses (meaning I need to compute the diffuse frac).
Thanks and I hope that made sense,
-Wil