def main(): ofname = "fitting_results_%s.csv" % ("Mercado") results_dir = "results" data_dir = "data" plot_dir = "plots" # Normal way model = FarquharC3(peaked_Jmax=True, peaked_Vcmax=True, model_Q10=True) F = FitJmaxVcmaxRd(model, ofname, results_dir, data_dir, plot_dir) F.main(print_to_screen=True) # Using Rdark measurement, assuming Rd = Rdark * 0.6 ofname = "fitting_results_%s.csv" % ("Mercado_Rdark") model = FarquharC3(peaked_Jmax=True, peaked_Vcmax=True) F = FitJmaxVcmaxKnownRdark(model, ofname, results_dir, data_dir, plot_dir) F.main(print_to_screen=True)
Example fit to some synthetic data... That's all folks. """ __author__ = "Martin De Kauwe" __version__ = "1.0 (13.08.2012)" __email__ = "*****@*****.**" import os import sys import glob import numpy as np from fit_farquhar_model.farquhar_model import FarquharC3 from fit_farquhar_model.fit_model import FitJmaxVcmaxRd ############################## # Fit Jmax, Vcmax + Rd ############################## ofname = "fitting_results.csv" results_dir = "results" data_dir = "data2" plot_dir = "plots" Egamma = 37830.0 model = FarquharC3() ############################## F = FitJmaxVcmaxRd(model, ofname, results_dir, data_dir, plot_dir) F.main(print_to_screen=False)
delimiter=delimiter, names=True, case_sensitive=True) return data ############################## # Fit Jmax, Vcmax + Rd ############################## ofname = "fitting_results.csv" results_dir = "results" data_dir = "data" plot_dir = "plots" model = FarquharC3(peaked_Jmax=True, peaked_Vcmax=True) ############################## F = FitJmaxVcmaxRd(model, ofname, results_dir, data_dir, plot_dir) F.main(print_to_screen=False) # OK what are the real values?? fit = read_data("results/fitting_results.csv") deg2kelvin = 273.15 index = 0 for Tleaf in np.arange(15.0, 40.0, 5.0): Tleaf += deg2kelvin Jmax25 = 150.0 Vcmax25 = Jmax25 / 1.6 Eaj = 30000.0 Eav = 60000.0 deltaSj = 650.0 deltaSv = 650.0 Hdv = 200000.0