Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer, | intent(in) | :: | mb | |||
integer, | intent(in) | :: | meap | |||
integer, | intent(in) | :: | subbasin_input_file_id |
subroutine evapotranspiration_initialise(mb, meap, subbasin_input_file_id)
use input, only: get_config_fid
use output, only: output_register_hydrotope_var
integer, intent(in) :: mb, meap, subbasin_input_file_id
read(get_config_fid(), EVAPOTRANSPIRATION_PARAMETERS)
if (iemeth.eq.0) call log_debug('evapotranspiration_initialise', &
'Evaporation - Priestley-Taylor, iemeth =', int=iemeth)
if (iemeth.eq.1) call log_debug('evapotranspiration_initialise', &
'Evaporation - Turc-Ivanov, iemeth =', int=iemeth)
if (idvwk.eq.1) call log_debug('evapotranspiration_initialise', &
'Evaporation variables after DVWK-M 238, idvwk =', int=idvwk)
if (idvwk.eq.0) call log_debug('evapotranspiration_initialise', &
'Evaporation variables from legacy method, idvwk =', int=idvwk)
eta_output_id = output_register_hydrotope_var("eta", .false.)
etp_output_id = output_register_hydrotope_var("etp", .false.)
soil_evaporation_output_id = &
output_register_hydrotope_var("soil_evaporation", .false.)
call evapotranspiration_allocate(mb, meap)
call evapotranspiration_read_input(subbasin_input_file_id)
yls = sin(lat / 57.296)
ylc = cos(lat / 57.296)
end subroutine evapotranspiration_initialise