= MICT Reference Simulations =

== Protocol to launch MICT simulations with the spinup Carbon (LibIGCM config) ==

[https://forge.ipsl.jussieu.fr/orchidee/attachment/wiki/DevelopmentActivities/Branches/ORCHIDEE-MICT-IMBALANCE-P/ReferenceSimulations/spinup_MICT.pdf]

== Main features of the run.def ==

This options belong to revision [2916/branches/ORCHIDEE-MICT/ORCHIDEE]:

{{{
### General flags ###

### added by Ye ###
# Activate STOMATE?
# Def   = n
STOMATE_OK_STOMATE=y

# Activate CO2?
# set to TRUE if photosynthesis is to be activated
# Def   = n
STOMATE_OK_CO2=y
### end Ye ###

# Activate the DGVM?
STOMATE_OK_DGVM = n

# Activate grassland management?
ENABLE_GRAZING = n

# Activate Spitfire?
FIRE_DISABLE = y
ALLOW_DEFOREST_FIRE = n

# Activate the soil freezing scheme?
OK_FREEZE = y

# Activate permafrost carbon?
OK_PC = y

# Activate the explicit snow scheme?
OK_EXPLICITSNOW = y

# Activate computation of flux densities of CH4 from wetlands?
CH4_CALCUL = n

# Activate Topmodel?
TOPM_CALCUL = n

# Activate CO2 downregulation?
DOWNREGULATION_CO2 = n

### Options related to soil carbon ###

#Activating this option allows the soil carbon to go below 0.5 m everywhere.
new_carbinput_intdepzlit=y

#Model for the temperature dependence of the respiration related to microbial activity
frozen_respiration_func=1 (default value is 1) 

#No insolation effect of the organic top soil layer
USE_SOILC_TEMPDIFF=n


### Vertical discretization scheme ###
# recommended values to get ndeep=32 (default values)
DEPTH_TMAX=38.
RATIO_GEOM_BELOW=1.05

### Options related to the routing scheme ###
RIVER_ROUTING = TRUE
DO_FLOODPLAINS =  y
DO_FLOODINFILT =  TRUE
DO_SWAMPS =  TRUE

### Added by Dan ###

# These are all current default values #
cryoturbate = y
cryoturbation_diff_k = 0.001
bioturbation_diff_k = 0.0001
bioturbation_depth = 2.
adjust_method_new = y

### End Dan ###

# Use age class and activate gross land use change?
GLUC_USE_AGE_CLASS = n [default n]

}}}

 == Spitfire ==

Main options for Spitfire module when it is active

{{{
#Chao add
  
## IF FIRE_DISABLE = n
# Read population density? [always set yes, otherwise the model may break]    
READ_POPDENS = y [always yes]                                                 
POPDENS_FILE = /home/orchidee01/ychao/FIRE_INPUT_DATA/HYDE/05deg/popd_1700.nc 
# Read human ignition parameter? [always yes]                                 
READ_HUMIGN = y                                                               
HUMIGN_FILE = /home/orchidee01/ychao/FIRE_INPUT_DATA/humign.nc                
# Lightning flashes input file                                                
LIGHTNING_FILE = /home/orchidee01/ychao/FIRE_INPUT_DATA/LISOTD_HRMC_V2.3.2011.lightn.nc
                                                                                
# Read the ratio flag? [always yes]                                           
READ_RATIO_FLAG = y                                                           
# When the flag is greater than 1, BA and CF are forced, otherwise            
# they are prognositcally determined.                                         
RATIO_FLAG_FILE = /home/orchidee01/ychao/FIRE_INPUT_DATA/FORCE_BA/flag_minus_ones.nc
# Read the ratio? [always y]                                                  
READ_RATIO = y                                                                
# Default ratio file with all ratios as 1, i.e., the BA is prognostically     
# determined. Current GFED3 ratio file is derived by calibrating with GFED BA.
RATIO_FILE = /home/orchidee01/ychao/FIRE_INPUT_DATA/FORCE_BA/ratio_GFED3_BvmaxT9_Opt3.nc
                                                                                
# Force the combustion fraction for the coarse fuel [default: n]              
READ_CF_COARSE = n                                                            
CF_COARSE_FILE = /home/orchidee01/ychao/FIRE_DATA/GFED3/CF_coarse.nc          
# Force the combustion fraction for the fine fuel [default: n]                
READ_CF_FINE = n                                                              
CF_FINE_FILE = /home/orchidee01/ychao/FIRE_DATA/GFED3/CF_fine.nc   
}}}


 == Gross land use change with age classes ==

Input for gross land use change when it is activated

{{{
#Chao add

##IF GLUC_USE_AGE_CLASS = y

GLUC_NAGEC_TREE = 6     # number of age class for forest                                      
GLUC_NAGEC_HERB = 4     # number of age class for crop, pasture and grassland
 ​ ​
#Total number of ​metaclasses​
GLUC_​NVMAP = 15              

​GLUC_​AGEC_GROUP__​xx: look-up table to link MTC to PFTs​
GLUC_AGE_CLASS_BOUND__xx: boundaries used to distinguish age class
 
# allow forestry harvest ? [default n]
​GLUC_ALLOW_FORESTRY_HARVEST = n    
​
# Assume there is a single age class in each MTC?​ ​[default n]
GLUC_SINGLE_AGE_CLASS = n      

​# Use a spital file for input as the boundaries of age class? [default n]
GLUC_USE_BOUND_SPA​ = n
                                               
  ​# If yes, file name of the spatial boundary input file:
  GLUC_AGE_THRESHOLD_FILE​

# File containing matrix for Net land cover change
GLUC_NET_LCC_FILE​
# File containing matrix for shifting-cultivation land cover change (from primary and secondary land, and forestry harvest)
GLUC_SHIFT_PRI_FILE
​GLUC_SHIFT_SEC_FILE
​# File containing matrix for forestry wood harvest
GLUC_FORESTRY_HARVEST_FILE
 
}}}