-
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() -> 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] (formerly known as rad%reffdf and rad%reffbm).
The direct beam and diffuse components describing the Effective Surface Reflectance are computed slightly differently [see cbl_albdo.F90], however
In general we require the Canopy Reflectance, and the Canopy Transmittance (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 therefore compute the Canopy Reflectance, and the Canopy Transmittance (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 ).