3.3.2. Setting the radiation field

The user must supply the radiation field he/she is interested in by editing the body of function URAD in module POPRATIO.

For example, in order to include the FGAL and FUVB radiation fields defined in module FIELDINTENSITIES the user should change the body of this function to:

urad = Fgal(lambda) 

if (z/=-1.0_WP) then 
    urad = urad + Fuvb(lambda,z,3) 

end if

If the user does not intend to use any of the radiation fields defined in module FIELDINTENSITIES then it does not need to be included in the project, and the corresponding use FieldIntensities statement may be deleted from module POPRATIO.

The default is a null radiation field:

urad = 0.0_WP ;