41 subroutine pdd(temp_mm, s_stat, ET)
47 real(dp),
dimension(12),
intent(in) :: temp_mm
48 real(dp),
intent(in) :: s_stat
50 real(dp),
intent(out) :: et
53 real(dp) :: inv_sqrt2pi, inv_s_stat, inv_sqrt2
56 real(dp),
parameter :: pi=3.141592653589793_dp
57 real(dp),
parameter :: time_year = 1.0_dp, &
58 time_year_inv = 1.0_dp/time_year, &
59 d_time = 1.0_dp/12.0_dp
61 inv_sqrt2pi = 1.0_dp/sqrt(2.0_dp*pi)
62 inv_s_stat = 1.0_dp/s_stat
63 inv_sqrt2 = 1.0_dp/sqrt(2.0_dp)
69 + ( s_stat*inv_sqrt2pi*exp(-0.5_dp*(temp_mm(n)*inv_s_stat)**2) &
70 + 0.5_dp*temp_mm(n)*
erfcc(-temp_mm(n)*inv_s_stat*inv_sqrt2) ) &
74 et = pdd_sum*time_year_inv
Declarations of kind types for SICOPOLIS.
subroutine pdd(temp_mm, s_stat, ET)
Computation of the positive degree days (PDD) with statistical temperature fluctuations; based on sem...
real(dp) function erfcc(x)
Computation of the complementary error function erfc(x) = 1-erf(x) with a fractional error everywhere...