***********************************************************************************************
* Filename: SoyBeanD.crp
* Contents: SWAP 3.2 - Data for detailed crop model
* Compiled by Ana Paula Schwantes, 2015-16
***********************************************************************************************
* Soybean (Glycine max (L.) Merrill)
***********************************************************************************************

*** PLANT GROWTH SECTION ***

***********************************************************************************************
* Part 1: Crop factor or crop height

  SWCF = 2 ! choice between crop factor [=1] or crop height [=2]
* Choose crop factor if ETref is used, either from meteo input file (SWETR = 1) or with Penman-Monteith
* Choose crop height if Penman-Monteith should be used with actual crop height, albedo and resistance

* If SWCF = 1, list crop factor CF [0.5..1.5, R],   as function of dev. stage DVS [0..2 -,R]:
* If SWCF = 2, list crop height CH [0..1000 cm, R], as function of dev. stage DVS [0..2 -,R]:
* (maximum 36 records)

    DVS       CH     CF
    0.0      25.0    1.0
    1.0      60.0    1.2
    2.0      80.0    1.2
* End of Table

* If SWCF = 2, list crop specifi values for:
  ALBEDO =   0.23 ! crop reflection coefficient [0..1.0 -, R]                    
  RSC    =   70.0 ! Minimum canopy resistance [0..10^6 s/m, R]                    
  RSW    =    0.0 ! Canopy resistance of intercepted water [0..10^6 s/m, R]
***********************************************************************************************

***********************************************************************************************
* Part 2 : Crop development
*
  IDSL   =        0 ! Switch:
*                     0 = Crop development before anthesis depends on temperature only
*                     1 = Crop development before anthesis depends on daylength only
*                     2 = Crop development before anthesis depends on both
*
* If IDSL = 1 or 2, specify:
  DLO    =  14.0     ! Optimum daylength for crop development [0..24 h, R]
  DLC    =   8.0     ! Minimum daylength, [0..24 h, R]
*
* If IDSL = 0 or 2 specify:
  TSUMEA =   860.00 ! Temperature sum from emergence to anthesis, [0..10000 C, R]
  TSUMAM =   720.00 ! Temperature sum from anthesis to maturity  [0..10000 C, R]
*
* List increase in temperature sum [0..60 C, R] as function of daily average temp. [0..100 C, R]
*         TAV  DTSM    (maximum 15 records)
  DTSMTB =
            0.00   0.00
           10.00   0.00
           20.00  10.00 
           35.00  25.00
           60.00  25.00
* End of Table
*
  DVSEND =      2.00 ! development stage at harvest [-]
***********************************************************************************************

***********************************************************************************************
* Part 3: Initial values
*
  TDWI   = 120.0 ! Initial total crop dry weight [0..10000 kg/ha, R]
  LAIEM  =  0.0163 ! Leaf area index at emergence [0..10 m2/m2, R]
  RGRLAI = 0.00700 ! Maximum relative increase in LAI [0..1 m2/m2/d, R]
***********************************************************************************************

***********************************************************************************************
* Part 4: Green surface area
*
* List specific leaf area [0..1 ha/kg, R] as function of devel. stage [0..2, R]
*         DVS  SLA    (maximum 15 records)
  SLATB =
           0.00 0.0014
           0.45 0.0025
           0.90 0.0025
           2.00 0.0007
* End of Table 
*
  SPA    =  0.0000 ! Specific pod area  [0..1 ha/kg, R]
  SSA    =  0.0000 ! Specific stem area [0..1 ha/kg, R]
  SPAN   =   23.00 ! Life span under leaves under optimum conditions, [0..366 d, R]
  TBASE  =    15.00 ! Lower threshold temperature for ageing of leaves ,[-10..30 C, R]
***********************************************************************************************

***********************************************************************************************
* Part 5: Assimilation
*
  KDIF   =    0.50 ! Extinction coefficient for diffuse visible light, [0..2 -, R]
  KDIR   =    0.75 ! Extinction coefficient for direct visible light, [0..2 -, R]
  EFF    =    0.40 ! Light use efficiency for real leaf [0..10 kg/ha/hr/(Jm2s), R]
*
* List max CO2 assimilation rate [0..100 kg/ha/hr, R] as function of development stage [0..2 -, R]
*          DVS    AMAX   (maximum 15 records)
  AMAXTB =
           0.00 37.000
           1.60 37.000
           2.00  0.000
* End of table 
*
* List reduction factor of AMAX [-, R] as function of average day temp. [-10..50 C, R]
*          TAVD   TMPF  (maximum 15 records)
  TMPFTB =
           0.00  0.000
          10.00  0.300
          20.00  0.600
          25.00  0.800
          30.00  1.000
          35.00  1.000
* End of table 
*
* List reduction factor of AMAX [-, R] as function of minimum day temp. [-10..50 C, R]
*          TMNR    TMNF  (maximum 15 records)
  TMNFTB = 
           0.00  0.000
           3.00  1.000
* End of table 
***********************************************************************************************

***********************************************************************************************
* Part 6: Conversion of assimilates into biomass
*
  CVL    =  0.6800 ! Efficiency of conversion into leaves,         [0..1 kg/kg, R]
  CVO    =  0.7200 ! Efficiency of conversion into storage organs, [0..1 kg/kg, R]
  CVR    =  0.7200 ! Efficiency of conversion into roots,          [0..1 kg/kg, R]
  CVS    =  0.6900 ! Efficiency of conversion into stems,          [0..1 kg/kg, R]
***********************************************************************************************

***********************************************************************************************
* Part 7: Maintenance respiration
*
  Q10    =  2.0000 ! Rel. increase in respiration rate with temperature, [0..5 /10 C, R]
  RML    =  0.0300 ! Rel. maintenance respiration rate of leaves,  [0..1 kgCH2O/kg/d, R]
  RMO    =  0.0100 ! Rel. maintenance respiration rate of st. org.,[0..1 kgCH2O/kg/d, R]
  RMR    =  0.0100 ! Rel. maintenance respiration rate of roots,   [0..1 kgCH2O/kg/d, R]
  RMS    =  0.0150 ! Rel. maintenance respiration rate of stems,   [0..1 kgCH2O/kg/d, R]
*
* List reduction factor of senescence [-, R] as function of dev. stage [0..2 -, R]
*          DVS    RFSE  (maximum 15 records) - o que significa este reduction factor of senescence?
  RFSETB = 
           0.00   1.00
           2.00   1.00
* End of table 
***********************************************************************************************

***********************************************************************************************
* Part 8: Partitioning
*
* List fraction of total dry matter increase partitioned to the roots [kg/kg, R]
* as function of development stage [0..2 -, R]
*          DVS     FR    (maximum 15 records)
  FRTB = 
           0.00   0.50
           0.75   0.10
           1.50   0.00
           2.00   0.00
* End of table 
*
* List fraction of total above ground dry matter incr. part. to the leaves [kg/kg, R]
* as function of development stage [0..2 -, R]
*          DVS     FL   (maximum 15 records)
  FLTB = 
           0.00   0.75
           1.00   0.75
           1.15   0.60
           1.30   0.46
           1.50   0.27
           1.70   0.00
           2.00   0.00
* End of table 
*
* List fraction of total above ground dry matter incr. part. to the stems [kg/kg, R]
* as function of development stage [0..2 -, R]
*          DVS    FS   (maximum 15 records)
  FSTB = 
           0.00   0.25
           1.00   0.25
           1.15   0.27
           1.30   0.27
           1.50   0.28
           1.70   0.00
           2.00   0.00
* End of table 
*
* List fraction of total above ground dry matter incr. part. to the st. organs [kg/kg, R]
* as function of development stage [0..2 -, R]
*          DVS    FO    (maximum 15 records)
  FOTB = 
           0.00   0.00
           1.00   0.00
           1.15   0.13
           1.30   0.27
           1.50   0.45
           1.70   1.00
           2.00   1.00
* End of table
***********************************************************************************************

***********************************************************************************************
* Part 9: Death rates
*
  PERDL =   0.030 ! Maximum rel. death rate of leaves due to water stress [0..3 /d, R]
*
* List relative death rates of roots [kg/kg/d] as function of dev. stage [0..2 -, R]
*          DVS    RDRR    (maximum 15 records)
  RDRRTB = 
          0.0000 0.0000
          1.5000 0.0000
          1.5001 0.0200
          2.0000 0.0200
* End of table
*
* List relative death rates of stems [kg/kg/d] as function of dev. stage [0..2 -, R]
*          DVS     RDRS    (maximum 15 records)
  RDRSTB = 
          0.0000 0.0000
          1.5000 0.0000
          1.5001 0.0200
          2.0000 0.0200
* End of table
***********************************************************************************************

***********************************************************************************************
* Part 10: Crop water use                                            
  swroottyp = 1      ! Switch for type root water extraction  [1,2 -, I] 
*                    ! (1 = Feddes et al., 1978; 2 = De Jong van Lier et al., 2006)
* if swroottyp=1 then enter HLIM1 - ADCRL
* if swroottyp=2 then enter wiltpoint, rootradius, rootcoefa
*
*
  HLIM1  =     -10.0 ! No water extraction at higher pressure heads, [-100..100 cm, R]
  HLIM2U =     -25.0 ! h below which optimum water extr. starts for top layer, [-1000..100 cm, R]
  HLIM2L =     -25.0 ! h below which optimum water extr. starts for sub layer, [-1000..100 cm, R]
  HLIM3H =    -200.0 ! h below which water uptake red. starts at high Tpot, [-10000..100 cm, R]
  HLIM3L =    -350.0 ! h below which water uptake red. starts at low Tpot, [-10000..100 cm, R]
  HLIM4  =   -5000.0 ! No water extraction at lower pressure heads, [-16000..100 cm, R]
  ADCRH  =       0.5 ! Level of high atmospheric demand, [0..5 cm/d, R]     
  ADCRL  =       0.1 ! Level of low atmospheric demand,  [0..5 cm/d, R]     
***********************************************************************************************

***********************************************************************************************
* Part 11: salt stress                                            

* only when solutes are simulated (SWSOLU=1 in SWP-file)

* relation between ECsat and crop reduction
  ECMAX  =       5.0 ! ECsat level at which salt stress starts, [0..20 dS/m, R] 
  ECSLOP =      20.0 ! Decline of rootwater uptake above ECMAX [0..40 %/dS/m, R] 

* relation between concentration and ECsat
  C2ECa  =     4.21  ! coefficient a to convert concentration to EC [0.0..1000.0 -, R]
  C2ECb  =     0.763 ! exponent b to convert concentration to EC [0.0..10.0 -, R]
* Switch to enter factor f (SWC2ECF ) per profile or per soil layer/horizon [1,2 -, I]
*    if SWC2ECF = 1 then enter one C2ECf-value for whole model profile
*    if SWC2ECF = 2 then enter one C2ECf-value for each model/soil layer/horizon
  SWC2ECF = 1
* factor f to convert concentration to EC [0.0..10.0 -, R]; 
*    dependent on SWC2ECF one value for model profile or a value for each soil horizon
  C2ECf  = 1.7
***********************************************************************************************

***********************************************************************************************
* Part 12: Interception                                            
*
  COFAB  =      0.25 ! Interception coefficient Von Hoyningen-Hune and Braden, [0..1 cm, R]
***********************************************************************************************

***********************************************************************************************
* Part 13: Root density distribution and root growth               
*
* List relative root density [0..1 -, R], as function of rel. rooting depth [0..1 -, R]:
*          RD     RDC   (maximum 11 records)
  RDCTB = 
           0.00   1.00
           1.00   1.00
* End of table
*
  RDI    =   5.00 ! Initial rooting depth, [0..1000 cm, R]
  RRI    =    1.20 ! Maximum daily increase in rooting depth, [0..100 cm/d, R]
  RDC    =  100.00 ! Maximum rooting depth crop/cultivar, [0..1000 cm, R]
*
*********************************************************************************

*** IRRIGATION SCHEDULING SECTION ***

**********************************************************************************
* Part 1: General

  SCHEDULE = 0  ! Switch for application irrigation scheduling [Y=1, N=0] 

* If SCHEDULE = 0, no more information is required in this input file! 
* If SCHEDULE = 1, continue ....

  STARTIRR = 30 3 ! Specify day and month after which irrigation scheduling is allowed [dd mm]
  ENDIRR = 31 12  ! Specify day and month after which irrigation scheduling is NOT allowed [dd mm]
  CIRRS = 0.0     ! solute concentration of scheduled irrig. water, [0..100 mg/cm3, R]
  ISUAS = 1       ! Switch for type of irrigation method: 
                  ! 0 = sprinkling irrigation
                  ! 1 = surface irrigation

* Specify pressure head at field capacity
* required for timing options  TCS = 2, 3, or 4 and depth option DCS = 1, else dummy 
  phFieldCapacity = -100.0   ! soil hydraulic pressure head [-1000.0 .. 0.0,cm, R] 

**********************************************************************************


**********************************************************************************
* Part 2: Irrigation time criteria

*** Choose one of the following 5 timing options:
  TCS = 1  ! Switch, timing criterion  [1..6, I]]
!            TCS = 1   :  Daily Stress
!            TCS = 2   :  Depletion of Readily Available Water
!            TCS = 3   :  Depletion of Totally Available Water
!            TCS = 4   :  Depletion Water Amount
!            TCS = 5   :  Pressure head or moisture content
!            TCS = 6   :  Fixed weekly irrigation, rootzone to field capacity
 

*** Daily stress criterion (TCS = 1)
* If TCS = 1, specify mimimum of ratio actual/potential transpiration Trel [0..1, R],
* as function of development stage DVS_tc1 [0..2, R], maximum 7 records:
  DVS_tc1  Trel
      0.0  0.95
      2.0  0.95
* End of table


*** Depletion of Readily Available Water (TCS = 2) 
* If TCS = 2, specify minimal fraction of readily available water RAW [0..1, R],
* as function of development stage DVS_tc2 [0..2, R], maximum 7 records:
  DVS_tc2   RAW
      0.0  0.95
      2.0  0.95
* End of table


*** Depletion of Totally Available Water (TCS = 3)
* If TCS = 3, specify minimal fraction of totally available water TAW [0..1, R],
* as function of development stage DVS_tc3 [0..2, R], maximum 7 records:
  DVS_tc3   TAW
      0.0  0.50
      2.0  0.50
* End of table


*** Depletion Water Amount (TCS = 4)
* If TCS = 4, specify maximum amount of water depleted below field cap. DWA [0..500 mm, R],
* as function of development stage DVS_tc4 [0..2, R], maximum 7 records:
  DVS_tc4   DWA
      0.0  40.0
      2.0  40.0
* End of table


*** Pressure head or Moisture content (TCS = 5)
* If TCS = 5, specify:
  PHORMC = 0   ! Switch, use pressure head (PHORMC=0) or water content (PHORMC=1)
  DCRIT = -30.0! Depth of the sensor [-100..0 cm, R]
* Also specify critical pressure head [-1.d6..-100 cm, R] or moisture content 
* [0..1.0 cm3/cm3, R], as function of development stage DVS_tc5 [0..2, R]:
  DVS_tc5  Value_tc5
      0.0    -1000.0
      2.0    -1000.0
* End of table


*** Fixed weekly irrigation, rootzone to field capacity (TCS = 6)
* If TCS = 6, specify: 
* Threshold for weekly irrigation; only irrigate when deficit is higher than threshold
  irgthreshold = 1.0       ! threshold value  [0.0..20.0 mm, R]


*** Select (optional) fixed time interval:
  tcsfix = 0  ! Switch, fixed timing criterion  [0 or 1, I]]
* If tcsfix = 1, specify:
  irgdayfix = 7    !   length of interval (number of days) [1..365, I]

**********************************************************************************


**********************************************************************************
* Part 3: Irrigation depth criteria

*** Choose one of the following 2 options for irrigation depth:
  DCS = 1      ! Switch, depth criterion  [1..2, I]]
!                DCS = 1 :  Back to Field Capacity
!                DCS = 2 :  Fixed Irrigation Depth


*** Back to Field Capacity (DCS = 1)  
* If DCS = 1, specify amount of under (-) or over (+) irrigation dI [-100..100 mm, R],
* as function of development stage DVS_dc1 [0..2, R], maximum 7 records:
  DVS_dc1   dI
     0.0  10.0
     2.0  10.0
* End of table


*** Fixed Irrigation Depth (DCS = 2)
* If DCS = 2, specify fixed irrigation depth FID [0..400 mm, R],
* as function of development stage DVS_dc2 [0..2, R], maximum 7 records:
  DVS_dc2   FID
      0.0  60.0
      2.0  60.0
* End of table


*** Select (optional) limitations of irrigation depth:
  dcslim = 0  ! Switch, limited irrigation depth  [0=No, 1=Yes]   [0..1, I]
* If dcslim = 1, specify:
  irgdepmin = 0.0    !   minimum irrigation depth [0.0d0 .. 100.0d0, mm, I]
  irgdepmax = 0.0    !   maximum irrigation depth [irgdepmin .. 1.0d7, mm, I]


* End of .crp file !