Validate the externalized version

For the moment, we validate the externalization of parameters on a fluxnet site. We consider the following reference simulation (named REF-TEST) :

# - Fluxnet files path
FluxnetPath=/home/orchidee01/vuichard/AR5/FORCING

# - Number of PFTs :
NbPFTs= 13

#   * ORCHIDEE name for physical parameters on each sites
#   PFT (IMPOSE_VEG), \
#   initial LAI (IMPOSE_VEG)
NameSitesParam= ( SECHIBA_VEGMAX, SECHIBA_LAI )
# by Default : 
# 1) first line is for PFT
# 2) second line is for LAI default for SLOWPROC lai model with : 
#    llaimax = 0.,  8.,  8.,  4., 4.5, 4.5,  4., 4.5,  4.,  2.,  2.,  2.,  2.)

#   * Name of component for each physical parameter described in NameSitesParam
#     (in SECHIBA, STOMATE, DRIVER)
CompSitesParam= ( SECHIBA, SECHIBA )

# Sites descriptions
#       Abbrv,  Filename ,      Inital year (for gregorian calendar) , Length (Y),  \
#param 1,   2,   3,   4,   5,   6,   7,   8,   9,  10,  11,  12,  13
Sites= ( NL-Loo, NL-Loo.nc   ,          1996,                         11     , \
     0.1, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.9, 0.0, 0.0, 0.0, \
      0.,  8.,  8.,  4., 4.5, 4.5,  4., 4.5,  4.,  2.,  2.,  2.,  2. ) \ 

Then we tried to compare the outputs files of the reference job to the following tests :

1. REF-TEST on 2PFTS : we want to check if the matching array PFT_TO_MTC works correctly. So we have in fluxnet.card :

# - Number of PFTs :
NbPFTs= 2

# Sites descriptions
#       Abbrv,  Filename ,      Inital year (for gregorian calendar) , Length (Y),  \
#param 1,   2,   3,   4,   5,   6,   7,   8,   9,  10,  11,  12,  13
Sites= ( NL-Loo, NL-Loo.nc   ,          1996,                         11     , \
      0.1, 0.9, \
      0.0,  2. ) \ 

and in PARAM/driver.def :

NVM = 2

PFT_TO_MTC__01 = 01
PFT_TO_MTC__02 = 10

2. The second test consists to test if the re-parametrization works correctly.

As there are some few differences between C3 crop and C3 grass, we change the parameters of the C3 crop and we give the values of the C3 grass.

So we turn a C3 crop in a C3 grass. Now we have in fluxnet.card :

# - Number of PFTs :
NbPFTs= 13

# Sites descriptions
#       Abbrv,  Filename ,      Inital year (for gregorian calendar) , Length (Y),  \
#param 1,   2,   3,   4,   5,   6,   7,   8,   9,  10,  11,  12,  13
Sites= ( NL-Loo, NL-Loo.nc   ,          1996,                         11     , \
     0.1, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.9, 0.0, \
      0.,  8.,  8.,  4., 4.5, 4.5,  4., 4.5,  4.,  2.,  2.,  2.,  2. ) \ 

In driver.def, we only got :

VCMAX_FIX__12 = 60.

VJMAX_FIX__12 = 120.

NATURAL__12 = TRUE

VJMAX_OPT__12 = 160.

TPHOTO_MIN_A__12 = 0.0025  

TPHOTO_MIN_B__12 = 0.1

TPHOTO_MIN_C__12 = -3.25

TPHOTO_OPT_A__12 = 0.0025

TPHOTO_OPT_B__12 = 0.25

TPHOTO_OPT_C__12 = 27.25

TPHOTO_MAX_A__12 =  0.00375

TPHOTO_MAX_B__12 =  0.35

TPHOTO_MAX_C__12 = 41.125

FLAM__12 = .25

LAI_MAX__12 = 2.5

PHENO_TYPE__12 = 4

PHENO_GDD_CRIT_A__12 = 0.03125 

PHENO_GDD_CRIT_B__12 = 6.25

PHENO_GDD_CRIT_C__12 = 270.

HUM_MIN_TIME__12 = 35.

LEAFAGECRIT__12 = 120.

SENESCENCE_HUM__12 = .2

SENESCENCE_TEMP_C__12 = -1.375 

SENESCENCE_TEMP_B__12 = .1

SENESCENCE_TEMP_A__12 = .00375

RSTRUCT_CONST__12 = 2.5

To be sure that we have the real C3 crop, we add this test consisting to write in driver.def ( with the same fluxnet.card) :

NVM = 13

PFT_TO_MTC__01 = 1
PFT_TO_MTC__02 = 2
PFT_TO_MTC__03 = 3
PFT_TO_MTC__04 = 4
PFT_TO_MTC__05 = 5
PFT_TO_MTC__06 = 6
PFT_TO_MTC__07 = 7
PFT_TO_MTC__08 = 8
PFT_TO_MTC__09 = 9
PFT_TO_MTC__10 = 10
PFT_TO_MTC__11 = 11
PFT_TO_MTC__12 = 10
PFT_TO_MTC__13 = 13

in order to see that the system of Metaclasses work well.