-
Notifications
You must be signed in to change notification settings - Fork 0
jhan_HaC1.4
It was intended the HaC 1.3 only include certain changes. However it was found that in order to accomodate these changes - many more were required.
Generally the flow through JULES looks like:
JULES - IO -> control() -> surf_couple_*() -> CABLE
The first of these surf_couple_*() routines is surf_couple_radiation().
==surf_couple_radiation()==
surf_couple_radiation() -> cable_rad_main() - etc. The required variables to be returned are:
land_albedo, & ! GridBoxMean albedo per rad band (row_length,rows,4) [land_albedo_ij]
alb_surft, & ! albedo per rad band per tile (land_pts, ntiles, 4) [alb_surft]
actually we can aggregate alb_surft to create land_albedo anywhere so we only need one thing back from CABLE really, alb_surft. However we do this aggregate within CABLE.
So to get alb_surft, alb_tile it is called in CABLE, we UNPACK from the calculated variables describing the effective surface reflectance, EffSurfRefl_beam(dif) l (formerly known as rad%reffdf and rad%reffbm).
The direct beam and diffuse components describing the Effective Surface Reflectance are computed slightly differently than they were previously [see cbl_albdo.F90], however scientifically identical to the previous calculation.
In general we require the Canopy Reflectance [CanopyRefl_beam(dif)], and the Canopy Transmittance [CanopyTransmit_beam(dif)] (diffuse/beam),
The surface albedo (considering that their may be snow present. AlbSnow (formerly snow%albsoilsn) and knowledge of the sunlit/vegetated status of the site (sunlit_veg_mask ).
We compute the Canopy Reflectance (formerly rad%rhocbm [beam] rad%rhocdf [dif] in cbl_albedo.F90)
We compute the Canopy Transmittance (formerly rad%cexpkbm [beam] rad%cexpkdm [dif] in cbl_albedo.F90)
The surface albedo (considering that their may be snow present (AlbSnow - formerly snow%albsoilsn) and knowledge of the sunlit/vegetated status of the site (sunlit_veg_mask ).