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