The file 'cntrlcomparison.csv', which is in comma-separated-value format, contains outputs from PISCAL that 
facilitates examination of how well the fitting is. First model prediction is calculated at each sampling 
point. Then model prediction is calculated at many selected levels of intercelluar CO2 partial pressure 
levels under three limitation states to enable plotting curves. The same structure is then repeated for 
each fitting of the same curve with different parameters to be estimated and with different curves. Note 
that some common information about a curve is repeated for each data point to make the structure of the 
file as regular as possible.

!The first section compares model calculation at the optimized parameters with measurements
CurveID - (NA), the curve identifier, repeated for each point
ChlFlUse? -(0/1), whether or not chlorophyll fluorescence data are used for identifying the limitation states
        0 = not used, 1 = used. (Currently this choice is still under evaluation and therefore ChlFlUse? = 0).
FitGi? _(0/1), whether or not the internal conductance (gi) is fitted for. 
         0 = not fitted and gi is either infinite or fixed at the value provided by the user. 
         1 = gi is estiamted, repeated for each point
FitGamma*? - (0/1), whether or not the chloroplastic CO2 partial pressure photocompensation point is
         fitted for. 0= not fitted, a prescribed value is used, repeated for each point
FitKco? - (0/1), whether or not the apparent Michaelis - Menten constant Kco = Kc(1+O/Ko) is fitted for.
         0= not fitted, a prescribed value is calculated from Kc, Ko and the oxygen partial pressure.
         1= Fitted for, repeated for each point
FitRd? - (0/1), whether the dark respiration is fitted for. 0= not fitted for, a prescribed value is used.
         1= fitted for. repeated for each point
FitAlpha? - (0/1), whether alpha (the non-returned fraction of the glycolate carbon
         recycled in the photorespiratory cycle) is fitted for.
         0= not fitted for and alpha = 0
         1= fitted for, repeated for each point
LimitCombina - (NA), the combination of limitation states present in the A/Ci dataset. Parameters 
         are estimated for those limitation states that occur in the measured curve, repeated for each point
         Rubisco: Rubisco limitation (Vcmax, Kco)
         RuBP: RuBP regeneration limitation (J)
         Tpu: export limitation (TPU, alpha)
PCO2i - (Pa), intercellular CO2 partial pressure
PCO2c - (Pa), chloroplastic CO2 partial pressure corresponding to the PCO2i of a point. 
         PCO2c=PCO2i-AnetCal/internal conductance
AnetMeas - (umol/m2/s), measured net assimilation rate
AnetCal  - (umol/m2/s), calculated net assimilation rate 
weitedrms - (umol/m2/s), root mean square error of the fitting (weiting coefficient =1)
PointLimitType - (1,2,3). 1 = point limited by rubisco
                          2 = point limited by RuBP regeneration
                          3 = point limited by TPU

!The second section gives photosynthesis for each of the three limitation states at selected values of intercellular
!partial pressure. Note that the corresponding values of chloraplastic CO2 partial pressure depend on the limitation state.
CO2i - (Pa), intercellular CO2 partial pressure
CO2cc- (Pa), Chloraplastic CO2 partial pressure for Rubisco limited photosynthesis
Ac- (umol/m2/s), Rubisco-limited net assimilation rate calculated with the optimized parameters
CO2cj- (Pa), Chloraplastic CO2 partial pressure for RuBP regeneration limited photosynthesis
Aj- (umol/m2/s), RuBP regeneration-limited net assimilation rate calculated with estimated parameters
CO2ct- (Pa), Chloraplastic CO2 partial pressure for export limited photosynthesis
At- (umol/m2/s), TPU-limited net assimilation rate calculated with estimated parameters
