!
! This is an initial function for running WEPS model
! read in IFC (soil) file
! read in Management file

      subroutine weps_init(isr)
      include 'p1werm.inc'
      include 'wpath.inc'
      include 'm1subr.inc'
      include 'm1sim.inc'

      include 'm1geo.inc'
      include 'm1flag.inc'
      include 'm1dbug.inc'
      include 's1layr.inc'
      include 's1surf.inc'
      include 's1phys.inc'
      include 's1agg.inc'
      include 's1dbh.inc'
      include 's1dbc.inc'
      include 's1sgeo.inc'
      include 'h1hydro.inc'
      include 'd1glob.inc'
      include 'b1glob.inc'
      include 'c1glob.inc'


      include 'h1scs.inc'
      include 'h1db1.inc'
      include 'file.fi'
      include 'command.inc'
      include 'precision.inc' !declaration for portable math range checking
! local variabl
      integer pd
      integer isr 
! read in the soil file and manage file
!       isr = 1
!      call inputFile()
! initialize variables defined from cmdline
                                !Note that all paths MUST end with a "/" for now.
      wc_type = 4         !default soil ifc file type (estimate Brooks and Corey from FC and WP)
      ifc_format = 0      !default soil ifc file format type (uses IFC specified slope value)
      report_debug = 0    !default report debug printing (0=off, 1=on)
      saeinp_daysim = 0   !0 value skips creation of a stand alone erosion input file.
      saeinp_jday = 0     !0 value skips creation of a stand alone erosion input file.
      init_cycle = 1      !default is to do one initialization cycle before reporting
      run_erosion = 1     !default is to run erosion submodel
      calibrate_crops = 0 !default is to NOT run in crop calibration mode
      calibrate_rotcycles = 0 ! default is to run specified number of rotation cycles during calibration
      cook_yield = 1      !default is to use partitioned Yield/Residue ratio   
!      cook_yield = 0      !default is to grow using full partitioning
      growth_stress = 3   !default is to turn on all stress values
      water_stress_max = 0.0 ! default is to allow full stress
      layer_scale = 2     !default is set to 2 mm
      layer_infla = 25    !default is set to 25 percent increase per layer
      layer_weighting = 0 !default is set to use the no weighting ".5" method
      puddle_warm = 1     !default is set to use warm puddling
      winter_ann_root = 1 !default is set to use no delay fall root growth

      wind_max_value = 0.0 !default is set to not use max wind speed capping
      wind_max_flag = 0   !default is set to not use max wind speed capping
      wepp_hydro = 0      !default is set to use Darcy method for hydrology
      soil_cond = 1       !default is set to output soil conditioning index file
      resurf_roots = 1    !default is set to resurface buried roots (process 26)
      calc_confidence = 0 !default is set to not calculate confidence intervals.
      
      am0csr = 1  ! set global current subregion variable?  
      aslrrc(isr) = 10.
      as0rrk(isr) = 0.9
    
      ! initialize variables
         abmf(isr) = 0.0
         abmst(isr) = 0.0 
         abffcv(isr)= 0.0
         abfscv(isr) = 0.0
         acfscv(isr) = 0.0
         adfscvtot(isr) = 0.0
         adrsaitot(isr) = 0.0
         abfcancov(isr) = 0.0
         admftot(isr) = 0.0
         admsttot(isr) = 0.0 
! initialize for the space
        amxsim(1,1) = 0.0
        amxsim(2,1) = 0.0
        amxsim(1,2) = 1000.0
        amxsim(2,2) = 1000.0 
        

! initialize
! manage file initialization will be added here!  
!    call init_report_vars(nperiods, maxper, n_rot_cycles)
!      call mandates(1)
    
      call asdini()
!     move the erodinit to the web_driver for daily initial     
!      call erodinit(isr)
      am0gdf = .true.
    
!     Initializations unique to particular submodels
   
       call decoinit(isr)
       call cropinit(isr)
       ! initialize all dependent variables
!       call soilinit(isr)
     
       call updres(isr)       
       call sumbio(isr)        
       call sci_init(isr)
         
!       Initialize the water holding capacity variable
       call hydrinit(isr)

!       initialize soil depth to bottom of layers (mm) from layer thickness (mm)
!       and initialize applied NO3
     
       call soilinit(isr)        
!read "yearly average info" from cligen header
       call cliginit          
      end