Python FitJmaxVcmaxRd Exemples

Exemple #1

Fichier : fit_jmax_vcmax_using_measured_Rdark.py Projet : zhujiedong/FitFarquharModel

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)

Exemple #2

Fichier : fit_jmax_vcmax_using_measured_Rdark.py Projet : mdekauwe/FitFarquharModel

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)

Exemple #3

Fichier : fit_jmax_vcmax_rd_PAR_example.py Projet : PaulPGauthier/FitFarquharModel

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)     

Exemple #4

                         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