sf_surf_irr_scheme (max_dom) = 0: no irrigation scheme
= 1: channel scheme (in Noah only)
= 2: drip scheme
= 3: sprinkler scheme
irr_daily_amount (max_dom) = daily equivalent irrigation water amount in mm (float, e.g. 5.7)
irr_start_hour (max_dom) = UTC start hour for irrigation (int, e.g. 5)
irr_num_hours (max_dom) = number of consecutive hours for irrigation (int, e.g. 3)
irr_start_julianday (max_dom) = Julian start day of irrigation (int, e.g. 135)
irr_end_julianday (max_dom) = Julian end day of irrigation (int, e.g. 255)
irr_ph (max_dom) = 0: all the static IRRIGATION field is activated
= 1: pseudo random activation field: ensures repeatability across compilers
= 2: random activation field: results might depend on the compiler fortran RANDOM_SEED
irr_freq (max_dom) = irrigation frequency in days (int, e.g. 3)
185 constants_irrigation=irr_freq*irr_daily_amount/(irr_num_hours*3600*rho(a,kms,b)*dz8w(a,kms,b)*100)
186 IF (irrigation(a,b).GE.0.1 .AND. irr_on.eq.1 ) THEN
187 CALL irr_calc_phase(irr_ph,phase,irr_rand_field_val(a,b),a,b,irrigation(a,b),irr_freq)
188 IF(irr_ph.EQ.0) THEN
189 qr_curr(a,kms,b)=qr_curr(a,kms,b)+irrigation(a,b)*constants_irrigation*dt
190 ELSE
191 IF(timing.EQ.int(phase)) THEN
192 qr_curr(a,kms,b)=qr_curr(a,kms,b)+irrigation(a,b)*constants_irrigation*dt
193 END IF
194 END IF
195 END IF