jhan_temp_cbl_init_rad

** cbl_init_radiation.F90 **

The subroutine init_radiation is called from cable_rad_driver() and also cbl_model_driver() (Formerly cbm) on the explicit call to CABLE.

It does potentially two things for CABLE:

In offline/standalone simulations only.

1. It computes the direct beam fraction of downward shortwave given the total downward shortwave.

In all simulations.

2. It calculates the extinction coefficients for both VIS and NIR bands, separately for both direct beam and diffuse components.

The subroutine init_radiation is shown below as it appears in the code - less the argument lists.

The title at the top of the page is a link to the file in the repository containing the init_radiation and module.

The subroutines called in this subroutine are also linked here to wiki versions of those subroutines.

subroutine init_radiation(...)

...

...

! Compute common scaling co-efficients used throughout init_radiation. Uses only parameters and constants
! with only model dependence being on Effctive LAI considering snow coverage

call Common_InitRad_Scalings(...)

! Define Raw extinction co-efficients for direct beam/diffuse radiation
! Largely parametrized per PFT. Does depend on zenith angle and effective LAI 
! [Formerly rad%extkb, rad%extkd]  

call ExtinctionCoeff(...)

! Define effective Extinction co-efficient for direct beam/diffuse radiation
! Extincion Co-eff defined by parametrized leaf reflect(transmit)ance - used in
! canopy transmitance calculations (cbl_albeo)
! [Formerly rad%extkbm, rad%extkdm ]

call EffectiveExtinctCoeffs(...)

! Offline/standalone forcing gives us total downward Shortwave. We have
! previosuly, arbitratily split this into NIR/VIS (50/50). We use 
!Spitter function to split these bands into direct beam and diffuse components

IF( jls_standalone .AND. .NOT. jls_radiation ) & CALL BeamFraction(...)

End subroutine