test_data/write_eco/WHAPS048_ECO.R

In palderman/DSSAT: A Comprehensive R Interface for the DSSAT Cropping Systems Model

structure(list(`ECO#` = c("IB0001", "IB0002", "IB0003", "IB0004", 
"IB0005", "IB0006", "IB0007", "IB0008", "DFAULT", "IB0099", "IB0098"
), ECONAME = c("GENERIC", "GENERIC2", "GENERIC3", "GENERIC4", 
"GENERIC5", "GENERIC6", "GENERIC7", "GENERIC8", "DEFAULT", "Cold Hardy (test)", 
"Cold Sens. (test)"), TBASE = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0), TOPT = c(26, 26, 26, 26, 26, 26, 26, 26, 34, 26, 26), ROPT = c(26, 
26, 26, 26, 26, 26, 26, 26, 34, 26, 26), TTOP = c(34, 34, 34, 
34, 34, 34, 34, 34, 34, 34, 34), P20 = c(12.5, 12.5, 12.5, 12.5, 
12.5, 12.5, 12.5, 12.5, 12.5, 12.5, 12.5), VREQ = c(50, 50, 50, 
50, 50, 50, 50, 50, 50, 50, 50), GDDE = c(10.2, 10.2, 10.2, 10.2, 
10.2, 10.2, 10.2, 10.2, 6, 6, 6), DSGFT = c(200, 200, 200, 200, 
200, 200, 200, 200, 170, 170, 170), RUE1 = c(3.8, 3.8, 3.8, 3.8, 
3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.8), RUE2 = c(3.8, 3.8, 3.8, 3.8, 
3.8, 3.8, 3.8, 3.8, 3.8, 3.8, 3.8), KVAL1 = c(0.63, 0.63, 0.63, 
0.63, 0.63, 0.63, 0.63, 0.63, 0.63, 0.63, 0.63), KVAL2 = c(0.63, 
0.63, 0.63, 0.63, 0.63, 0.63, 0.63, 0.63, 0.63, 0.63, 0.63), 
    SLAP2 = c(300, 270, 270, 225, 270, 270, 270, 270, 270, 300, 
    300), TC1P1 = c(2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 
    2.5, 2.5), TC1P2 = c(0.6, 0.6, 1, 1, 0.6, 0.6, 0.6, 0.6, 
    0.6, 0.6, 0.6), DTNP1 = c(0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005), PLGP1 = c(1400, 
    1400, 1400, 1400, 1400, 1400, 1400, 1400, 1400, 1400, 1400
    ), PLGP2 = c(0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 
    0.6, 0.6), P2AF = c(0.6, 0.6, 0.6, 0.6, 0.6, 0.1, 0.6, 0.1, 
    0.6, 0.6, 0.6), P3AF = c(50, 50, 50, 50, 50, 50, 50, 50, 
    50, 50, 50), P4AF = c(3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3), P5AF = c(1, 
    1, 1, 1, 3, 3, 1, 3, 1, 1, 1), P6AF = c(3, 3, 3, 3, 3, 3, 
    3, 3, 3, 3, 3), ADLAI = c(1, 1, 1, 1, 1, 1, 0.7, 0.5, 1, 
    1, 1), ADTIL = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), ADPHO = c(1, 
    1, 1, 1, 1, 1, 1, 1, 1, 1, 1), STEMN = c(0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0), MXNUP = c(0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 0.6, 
    0.6, 0.6, 0.6, 0.6), MXNCR = c(0.04, 0.04, 0.04, 0.04, 0.04, 
    0.04, 0.04, 0.04, 0.04, 0.04, 0.04), WFNU = c(2, 2, 2, 2, 
    2, 2, 2, 2, 2, 2, 2), PNUPR = c(0.45, 0.45, 0.45, 0.45, 0.45, 
    0.45, 0.45, 0.45, 0.45, 0.45, 0.45), EXNO3 = c(6.75, 6.75, 
    6.75, 6.75, 6.75, 6.75, 6.75, 6.75, 6.75, 6.75, 6.75), MNNO3 = c(0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0), EXNH4 = c(6.5, 6.5, 6.5, 6.5, 
    6.5, 6.5, 6.5, 6.5, 6.5, 6.5, 6.5), MNNH4 = c(0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0), INGWT = c(3.5, 3.5, 3.5, 3.5, 3.5, 
    3.5, 3.5, 3.5, 3.5, 3.5, 3.5), INGNC = c(0.03, 0.03, 0.03, 
    0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03, 0.03), FREAR = c(0.25, 
    0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25
    ), MNNCR = c(1.23, 1.23, 1.23, 1.23, 1.23, 1.23, 1.23, 1.23, 
    1.23, 1.23, 1.23), GPPSS = c(2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 
    2), GPPES = c(5, 5, 5, 5, 5, 5, 5, 5, 5, 5, 5), MXGWT = c(55, 
    55, 55, 55, 55, 55, 55, 55, 55, 55, 55), MNRTN = c(4.5, 4.5, 
    4.5, 4.5, 4.5, 4.5, 4.5, 4.5, 4.5, 4.5, 4.5), NOMOB = c(0.25, 
    0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25
    ), RTDP1 = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), 
    RTDP2 = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), FOZ1 = c(0.06, 
    0.06, 0.06, 0.06, 0.06, 0.06, 0.06, 0.06, 0.06, 0.06, 0.06
    ), SFOZ1 = c(0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 0.08, 
    0.08, 0.08, 0.08), TSEN = c(NA, NA, NA, NA, NA, NA, NA, NA, 
    NA, 0.5, 10), CDAY = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, 
    20L, 5L)), class = c("DSSAT_tbl", "tbl_df", "tbl", "data.frame"
), row.names = c(NA, -11L), v_fmt = c(`ECO#` = "%-7s", ECONAME = "%-18s", 
TBASE = "%5.0f", TOPT = "%6.0f", ROPT = "%6.0f", TTOP = "%6.0f", 
P20 = "%6.1f", VREQ = "%6.0f", GDDE = "%6.1f", DSGFT = "%6.0f", 
RUE1 = "%6.1f", RUE2 = "%6.1f", KVAL1 = "%6.2f", KVAL2 = "%6.2f", 
SLAP2 = "%6.0f", TC1P1 = "%6.1f", TC1P2 = "%6.1f", DTNP1 = "%6.3f", 
PLGP1 = "%6.0f", PLGP2 = "%6.1f", P2AF = "%6.1f", P3AF = "%6.0f", 
P4AF = "%6.0f", P5AF = "%6.0f", P6AF = "%6.0f", ADLAI = "%6.1f", 
ADTIL = "%6.0f", ADPHO = "%6.0f", STEMN = "%6.0f", MXNUP = "%6.1f", 
MXNCR = "%6.2f", WFNU = "%6.0f", PNUPR = "%6.2f", EXNO3 = "%6.2f", 
MNNO3 = "%6.0f", EXNH4 = "%6.1f", MNNH4 = "%6.0f", INGWT = "%6.1f", 
INGNC = "%6.2f", FREAR = "%6.2f", MNNCR = "%6.2f", GPPSS = "%6.0f", 
GPPES = "%6.0f", MXGWT = "%6.0f", MNRTN = "%6.1f", NOMOB = "%6.2f", 
RTDP1 = "%6.0f", RTDP2 = "%6.0f", FOZ1 = "%6.2f", SFOZ1 = "%6.2f", 
TSEN = "%6.1f", CDAY = "%6.0f"), first_line = "*WHEAT ECOTYPE COEFFICIENTS: WHAPS048 MODEL", comments = c("                            1     2     3     4     5     6     7     8     9    10    11    12    13    14    15    16    17    18    19    20    21    22    23    24    25    26    27    28    29    30    31    32    33    34    35    36    37    38    39    40    41    42    43    44    45    46    47    48    49    50", 
"GENERIC with different SLAP2", "GENERIC with different SLAP2 and TC1P2", 
"GENERIC with different SLAP2 and P5AF", "GENERIC with different SLAP2, P2AF, and P5AF", 
"GENERIC with different SLAP2 and ADLAI", "GENERIC with different SLAP2, P2AF, P5AF, and ADLAI", 
"chp added this for test of optional temperature input for cold damage", 
"TBASE - base temperature below which no development occurs, C", 
"TOPT  - temperature at which maximum development rate occurs during vegetative stages, C", 
"ROPT  - temperature at which maximum development rate occurs for reproductive stages, C (no effect, WHAPS)", 
"TTOP  - temperature above which no development occurs, C", "P2O   - Daylength below which daylength does not affect development rate, hours\t(no effect, WHAPS)", 
"VREQ  - Vernalization required for max.development rate (VDays)", 
"GDDE  - Growing degree days per cm seed depth required for emergence, GDD/cm", 
"DSGFT - GDD from End Ear Growth to Start Grain Filling period, C", 
"RUE1  - Pre-anthesis radiation use efficiency, g plant dry matter/MJ PAR", 
"RUE2  - Post-anthesis radiation use efficiency, g plant dry matter/MJ PAR", 
"KVAL1 - Pre-anthesis canopy light extinction coefficient for daily PAR", 
"KVAL2 - Post-anthesis canopy light extinction coefficient for daily PAR", 
"SLAP2 - ratio of leaf area to mass at end of leaf growth (cm2/g)", 
"TC1P1 - for calculating tc1: tiller number from emerg. to term. spik.(=stem elongation)", 
"TC1P2 - tc1 = tc1_p1 + tc1_p2 *cumph_nw(istage)", "DTNP1 - for calculating dtiln: tiller number after term. spik.(=stem elongation)", 
"PLGP1 - for calculating plag: potential leaf growth.  plag= plag_p1*cumph(istage)**plag_p2", 
"PLGP2 - for calculating plag: potential leaf growth.  plag= plag_p1*cumph(istage)**plag_p2", 
"P2AF  - threshold AD in a layer becoming effective on root growth", 
"P3AF  - length of downwards root not effected under aeration deficit", 
"P4AF  - days to be accumulated before aeration deficit effects root growth", 
"P5AF  - power term at af1", "P6AF  - days to be accumulated before aeration deficit effects crop growth", 
"ADLAI - threshold aeration deficit (AF2) affecting LAI        (set to 1.0 for no stress run)", 
"ADTIL - threshold aeration deficit (AF2) effecting tillering  (set to 1.0 for no stress run)", 
"ADPHO - threshold aeration deficit (AF2) effecting photosyn.  (set to 1.0 for no stress run)", 
"STEMN - 0=original C to grain translocation, >0 to 1.0 sets % of C of stem to be transloc. to grain", 
"MXNUP - max N uptake per day", "MXNCR - 0.035=20%,   .04=23% protein, max n:c ratio of grain growth", 
"WFNU  - power term for water effect on N supply", "PNUPR - potential uptake rate (mg/meter/day)", 
"EXNO3 - exponent for NO3 supply factor", "MNNO3 - minimum for NO3 supply factor", 
"EXNH4 - exponent for NH4 supply factor", "MNNH4 - minimum for NH4 supply factor", 
"INGWT - initial grain weight  (mg/100 grain)", "INGNC - % protein, initial grain N conc   (init_grain_nconc or p_init_grain_nconc from APSIM Nwheat)", 
"FREAR - fraction between end ear and begin grainfilling, setting min stem weight for remobilisation", 
"MNNCR - % protein, min n:c ratio of grain growth", "GPPSS - gpp_start_stage - Grain per plant: 2=stem elong., 3=end leaf stage at which to start accumo. stem for gpp calc.", 
"GPPES - gpp_end_stage - Start grainfilling stage at which to end accumulation stem for Grain per plant calc", 
"MXGWT - maximum kernal weight  [100 = no effect]", "MNRTN - min root n due to grain n initialisation [0 = off] (root_n_min or p_root_n_min from APSIM Nwheat)", 
"NOMOB - fraction of accum stem weight that is not mobile [0 = original]", 
"RTDP1 - sw effect, 0 = old version", "RTDP2 - crop stress effect, 0 = old version", 
"FOZ1  - ozone effect rate on photosynthesis: 0.06 = ozone tolerant, 0.10 = intermediate, 0.50 = ozone sensitive", 
"FOZ2  - ozone effect initialization on photosynthesis (dependent on FOZ1): 1.015 = ozone tolerant, 1.025 = intermediate, 1.125 = ozone sensitive", 
"SFOZ1 - ozone effect rate on leaf senescence: 0.08 = ozone tolerant, 0.25 = intermediate, 0.40 = ozone sensitive", 
"SFOZ2 - ozone effect initialization on leaf senescence (dependent on SFOZ1): 0.80 = ozone tolerant, 0.38 = intermediate, 0.00 = ozone sensitive", 
"TSEN  - Critical temperature below which leaf damage occurs (default 6°C)", 
"CDAY  - cold days parameter (default 15.0)"))