MCRCDefines.h

#ifndef MCRC_DEFINES_H
#define MCRC_DEFINES_H
 
// VERSION and DATE of KERNEL LIBRARY
// See Parameters.cpp
 
// Facies //
// Simulation regular discretization default values for unavailable sample (below or above simulation)
#define UNAVAILABLE_FACIES     255   
#define UNAVAILABLE_AGE          0   
#define UNAVAILABLE_GRAIN_SIZE  -1.  
 
// Following family values correspond to the position of each facies in the FaciesRegister (See Facies.cpp)
// These defines are intended to be used by plugins.
//                          Family             Abbr.     Type of deposit
#define FACIES_UDF_ID       0             
#define FACIES_CL_ID        1             
#define FACIES_PB_ID        2             
#define FACIES_SP_ID        3             
#define FACIES_CSI_ID       4             
#define FACIES_CCH_ID       5             
#define FACIES_CSII_ID      6             
#define FACIES_LV_ID        7             
#define FACIES_OB_ID        8             
#define FACIES_MP_ID        9             
#define FACIES_CF_ID        10            
#define FACIES_WL_ID        11            
#define FACIES_DR_ID        12            
#define FACIES_PL_ID        13            
 
#define FACIES_MIN_ID       0             
#define FACIES_MAX_ID       13            
 
#define FACIES_UDF_LABEL  "Undefined"
#define FACIES_CL_LABEL   "Channel Lag"
#define FACIES_PB_LABEL   "Point Bar / LAPs"
#define FACIES_SP_LABEL   "Sand Plug"
#define FACIES_CSI_LABEL  "Crevasse Splay I"
#define FACIES_CCH_LABEL  "Splay II Channels"
#define FACIES_CSII_LABEL "Crevasse Splay II"
#define FACIES_LV_LABEL   "Levee"
#define FACIES_OB_LABEL   "Overbank"
#define FACIES_MP_LABEL   "Mud Plug"
#define FACIES_CF_LABEL   "Channel Fill"
#define FACIES_WL_LABEL   "Wetland"
#define FACIES_DR_LABEL   "Draping"
#define FACIES_PL_LABEL   "Pelagic"
 
#define FACIES_UDF_BRIEF  "UDF"
#define FACIES_CL_BRIEF   "CL"
#define FACIES_PB_BRIEF   "PB"
#define FACIES_SP_BRIEF   "SP"
#define FACIES_CSI_BRIEF  "CSI"
#define FACIES_CCH_BRIEF  "CCh"
#define FACIES_CSII_BRIEF "CSII"
#define FACIES_LV_BRIEF   "LV"
#define FACIES_OB_BRIEF   "OB"
#define FACIES_MP_BRIEF   "MP"
#define FACIES_CF_BRIEF   "CF"
#define FACIES_WL_BRIEF   "WL"
#define FACIES_DR_BRIEF   "DR"
#define FACIES_PL_BRIEF   "PL"
 
#ifndef SWIG
#define FACIES_UDF_COLOR  0.6,0.6,0.9
#define FACIES_CL_COLOR   1.,0.5,0.
#define FACIES_PB_COLOR   1.,1.,0.
#define FACIES_SP_COLOR   0.75,0.75,0.55
#define FACIES_CSI_COLOR  0.8,0.5,0.2
#define FACIES_CCH_COLOR  0.8,1.,0.2
#define FACIES_CSII_COLOR 0.8,0.8,0.2
#define FACIES_LV_COLOR   0.4,0.8,0.2
#define FACIES_OB_COLOR   0.,1.,0.
#define FACIES_MP_COLOR   0.,0.8,0.5
#define FACIES_CF_COLOR   0.5,0.5,0.5
#define FACIES_WL_COLOR   0.85,0.45,0.85
#define FACIES_DR_COLOR   0.6,0.8,0.9
#define FACIES_PL_COLOR   1.,0.8,1.
#endif
 
// Grain size //
// Grain size class {0,15}
#define FACIES_UDF_GRAIN  15  // Arbitrary
#define FACIES_CL_GRAIN   13
#define FACIES_PB_GRAIN   10
#define FACIES_SP_GRAIN   9
#define FACIES_CSI_GRAIN  9
#define FACIES_CCH_GRAIN  8
#define FACIES_CSII_GRAIN 7
#define FACIES_LV_GRAIN   6
#define FACIES_OB_GRAIN   2
#define FACIES_MP_GRAIN   1
#define FACIES_CF_GRAIN   1
#define FACIES_WL_GRAIN   0
#define FACIES_DR_GRAIN   0
#define FACIES_PL_GRAIN   0
 
#define MD_MAX_GRANULO    1.0000 // 15
//                        0.9375 // 14
#define MD_CL_GRANULO     0.8750 // 13
//                        0.8125 // 12
//                        0.7500 // 11
#define MD_PB_GRANULO     0.6875 // 10
#define MD_SP_GRANULO     0.6250 //  9
#define MD_CSI_GRANULO    0.6250 //  9
#define MD_CCH_GRANULO    0.5625 //  8
#define MD_CSII_GRANULO   0.5000 //  7
#define MD_LV_GRANULO     0.4375 //  6
//                        0.3750 //  5
//                        0.3125 //  4
//                        0.2500 //  3
#define MD_OB_GRANULO     0.1875 //  2
#define MD_MP_GRANULO     0.1250 //  1
#define MD_CF_GRANULO     0.1250 //  1
#define MD_WL_GRANULO     0.0625 //  0
#define MD_DR_GRANULO     0.0625 //  0
#define MD_PL_GRANULO     0.0625 //  0
#define MD_MIN_GRANULO    0.0625 //  0
// Grain size 0. value is forbidden
 
#define MD_CL_GRANULO_MAX    1.0000
#define MD_CL_GRANULO_MIN    0.8125
#define MD_PB_GRANULO_MAX    0.8125
#define MD_PB_GRANULO_MIN    0.5625
#define MD_SP_GRANULO_MAX    0.7500
#define MD_SP_GRANULO_MIN    0.5625
#define MD_CSI_GRANULO_MAX   0.7500
#define MD_CSI_GRANULO_MIN   0.5625
#define MD_CCH_GRANULO_MAX   0.6875
#define MD_CCH_GRANULO_MIN   0.5000
#define MD_CSII_GRANULO_MAX  0.5625
#define MD_CSII_GRANULO_MIN  0.3750
#define MD_LV_GRANULO_MAX    0.5625
#define MD_LV_GRANULO_MIN    0.3750 
#define MD_OB_GRANULO_MAX    0.3125
#define MD_OB_GRANULO_MIN    0.0625
#define MD_MP_GRANULO_MAX    0.1250
#define MD_MP_GRANULO_MIN    0.1250
#define MD_CF_GRANULO_MAX    0.1250
#define MD_CF_GRANULO_MIN    0.1250
#define MD_WL_GRANULO_MAX    0.0625
#define MD_WL_GRANULO_MIN    0.0625
#define MD_DR_GRANULO_MAX    0.0625
#define MD_DR_GRANULO_MIN    0.0625
#define MD_PL_GRANULO_MAX    0.0625
#define MD_PL_GRANULO_MIN    0.0625
 
#define MD_NB_CLASS_GR        16          
#define MD_NB_TEMPLATES        4          
 
#define MD_GR_CLASSES {  " 0 - Clay (13-14)",         \
                         " 1 - Clay (11-12)",         \
                         " 2 - Clay (9-10)",          \
                         " 3 - Clay (8)",             \
                         " 4 - Silt (7)",             \
                         " 5 - Silt (6)",             \
                         " 6 - Silt (5)",             \
                         " 7 - Silt (4)",             \
                         " 8 - Very Fine Sand (3)",   \
                         " 9 - Fine Sand (2)",        \
                         "10 - Medium Sand (1)",      \
                         "11 - Coarse Sand (0)",      \
                         "12 - Very Coarse Sand (-1)",\
                         "13 - Gravel (-2)",          \
                         "14 - Pebble (-3,-5)",       \
                         "15 - Cobble (-6,-8)"        \
                      }
 
//                                Clay    =>     Silt    =>       Sand       => Gravel/Pebble
//                     Classes:   0  1  2  3  4  5    6   7   8   9  10  11  12 13 14 15
 
#define MD_FLUV_GR_TEMPLATES  { { 0, 0, 0, 0, 0, 2,   3,  6, 13, 27, 26, 10,  8, 5, 0, 0 }, /*Standard*/     \
                                { 0, 0, 0, 0, 0, 2,   2,  4, 13, 20, 18, 16, 12, 8, 3, 2 }, /*Sandy*/        \
                                { 0, 0, 0, 0, 0, 0,  10, 13, 33, 18, 10,  8,  5, 3, 0, 0 }, /*Fine grained*/ \
                                { 0, 0, 2, 3, 5, 10, 24, 29, 22,  2,  2,  1,  0, 0, 0, 0 }, /*Muddy*/        \
                               }
 
#define MD_TURBI_GR_TEMPLATES { { 0, 0, 1, 2, 3, 4,   8, 17, 25, 22,  8,  5,  3, 2, 0, 0 }, /*Standard*/     \
                                { 0, 0, 0, 0, 0, 2,   3,  6, 13, 27, 26, 10,  8, 5, 0, 0 }, /*Sandy*/        \
                                { 0, 0, 0, 0, 0, 0,  10, 13, 50, 22,  5,  0,  0, 0, 0, 0 }, /*Fine grained*/ \
                                { 1, 2, 2, 3, 5, 10, 24, 29, 22,  2,  0,  0,  0, 0, 0, 0 }, /*Muddy*/        \
                               }
 
 
#define MD_GR_TEMPLATE_NAMES {  "Standard",     \
                                "Sandy",        \
                                "Fine grained"  \
                                "Muddy",        \
                             }
 
#define MD_CLASS_SAND_MEDIUM_IDX 10   
#define MD_CLASS_SAND_FINE_IDX    9   
 
#define MD_STD_IDX          0   
#define MD_SANDY_IDX        1   
#define MD_FINE_IDX         2   
#define MD_MUDDY_IDX        3   
 
 
// Erodibility Factor //
#define MD_E_FAC   1.e-8   
 
 
// User classes //
#define MD_MAX_USER_CLASS    26 
 
 
// Menu selection strings //
 
#define MD_MENU_FLUVIAL    "Fluvial"
#define MD_MENU_TURBIDITE  "Turbidite"
#define MD_MENU_TURBIDITE2 "Turbidite (Research)" // Kept for retro compatibility
 
#define MD_MENU_NONE       "None"
 
#define MD_MENU_UDF        FACIES_UDF_LABEL
#define MD_MENU_DR         FACIES_DR_LABEL
#define MD_MENU_PL         FACIES_PL_LABEL
#define MD_MENU_PB         FACIES_PB_LABEL
#define MD_MENU_LV         FACIES_LV_LABEL
#define MD_MENU_OB         FACIES_OB_LABEL
#define MD_MENU_CF         FACIES_CF_LABEL
#define MD_MENU_MP         FACIES_MP_LABEL
#define MD_MENU_WL         FACIES_WL_LABEL
 
#define MD_MENU_FL         "Overbank Flow"
#define MD_MENU_ER         "Erosive"
#define MD_MENU_GH         "Ghost Migration"
 
#define MD_MENU_NEVER      "Never"
#define MD_MENU_PER        "Periodic"
#define MD_MENU_POIS       "Poisson"
 
#define MD_MENU_CONST      "Constant"
#define MD_MENU_UNIF       "Uniform"
#define MD_MENU_NORM       "Normal"
#define MD_MENU_LOGN       "Lognormal"
 
#define MD_MENU_FILE       "File"
 
typedef enum { FLUVIAL_SYSTEM   = 0, 
               TURBIDITE_SYSTEM = 1, 
} SysType;
 
typedef enum { AGG_NONE,  
               AGG_DR,    
               AGG_OB,    
               AGG_GH,    
               AGG_ER     
} AggType;
 
typedef enum { FREQ_NEVER,    
               FREQ_PERIODIC, 
               FREQ_POISSON   
} FreqType;
 
typedef enum { DIST_CONSTANT, 
               DIST_UNIFORM,  
               DIST_NORMAL,   
               DIST_LOGNORMAL 
} DistType;
 
typedef enum { FILL_UNDEFINED               = FACIES_UDF_ID,  
               FILL_CHANNEL_LAG             = FACIES_CL_ID,   
               FILL_POINT_BAR               = FACIES_PB_ID,   
               FILL_SAND_PLUG               = FACIES_SP_ID,   
               FILL_CREVASSE_SPLAY_I        = FACIES_CSI_ID,  
               FILL_CREVASSE_SPLAY_CHANNELS = FACIES_CCH_ID,  
               FILL_CREVASSE_SPLAY_II       = FACIES_CSII_ID, 
               FILL_LEVEE                   = FACIES_LV_ID,   
               FILL_OVERBANK                = FACIES_OB_ID,   
               FILL_MUD_PLUG                = FACIES_MP_ID,   
               FILL_CHANNEL_FILL            = FACIES_CF_ID,   
               FILL_WETLAND                 = FACIES_WL_ID,   
               FILL_DRAPING                 = FACIES_DR_ID,   
               FILL_PELAGIC                 = FACIES_PL_ID    
} FillingFacies;
 
typedef enum { DRAP_UNDEFINED               = FACIES_UDF_ID,  
               DRAP_DRAPING                 = FACIES_DR_ID,   
               DRAP_PELAGIC                 = FACIES_PL_ID    
} DrapingFacies;
 
typedef enum { ZUL_NONE,    // No upper limit surface
               ZUL_NEUTRAL, // When upper limit has been reach, the simulation stops and do nothing more
               ZUL_EROSIVE, // When upper limit has been reach, the simulation stops and erodes all exceeding deposits
               ZUL_STRATI   // The simulation topography will be conditioned to fit the upper limit surface at the end
} ZulType;
 
typedef enum { SURF_NONE,  
               SURF_CONST, 
               SURF_FILE   
} SurfType;
 
typedef enum { EMAP_CONST,     
               EMAP_LOAD,      
               EMAP_LOAD_CONV, 
               EMAP_BUILD      
} EmapType;
 
typedef enum { SURF_EXCD      = 0, 
               SURF_REPL      = 1, 
               SURF_AGGR      = 2, 
               SURF_EROD_DOWN = 3, 
               SURF_EROD_UP   = 4  
} SurfMode;
 
typedef enum { STE_NO_ERROR,            
               STE_UNKNOWN_ERROR,       
               STE_INIT_ERROR,          
               STE_NOT_READY,           
               STE_CHANNEL_ERROR,       
               STE_CONSISTENCY_WARN,    
               STE_CONSISTENCY_ERROR,   
               STE_ZUL_REACHED,         
               STE_STOP_HONORED,        
               STE_STOPPED              
} StatusError;
 
typedef enum { TL_NONE    = 0,
               TL_FATAL   = 1,
               TL_ERROR   = 2,
               TL_WARNING = 3,
               TL_INFO    = 4,
               TL_DEBUG   = 5,
               TL_ALL     = 6
} TraceLevel;
 
typedef enum { STEEP_UNKNOWN    = -2, 
               STEEP_CUTOFF     = -1, 
               STEEP_OK         = 0,  
               STEEP_SELF_CUT   = 1,  
               STEEP_UP_OR_DOWN = 2,  
               STEEP_SHORT      = 3,  
               STEEP_TRAP       = 4,  
               STEEP_OUT_EXT    = 5,  
               STEEP_MARGINS    = 6,  
               STEEP_COMPLETE   = 7   
} SteepestError;
 
typedef enum { NO_DISTRIB,
               FREE_DISTRIB,
               PREMIUM_DISTRIB,
               RESEARCH_DISTRIB
} MCRCDistrib;
 
#define NO_DISTRIB_STR       "No valid license"
#define FREE_DISTRIB_STR     "Free"
#define PREMIUM_DISTRIB_STR  "Premium"
#define RESEARCH_DISTRIB_STR "Research"
 
#ifndef SWIG
typedef enum { COL_FACIES_CODING,                 
               COL_GRANULO_CODING,                
               COL_FACIES_AND_GRANULO_CODING,     
               COL_FACIES_GRANULO_AND_AGE_CODING, 
               COL_AGE_CODING,                    
               COL_ERODIBILITY_CODING             
} ColorsCoding;
 
typedef enum {
  ONLY_SAND,
  ONLY_CHANNELIZED,
  EACH_FACIES
} ComputedFacies;
#endif
 
// Trace function pointer declaration //
typedef void (*TraceFunc)(const char*, TraceLevel); 
typedef void (*DumpFunc)(const char*);              
 
// Progress event function pointer declaration //
typedef bool (*ProgressFunc)(int); 
 
// Parameters keys //
 
// Others parameters
#define KEY_SIM_SEED                 "SIM_SEED"                 
#define KEY_SIM_TYPE                 "SIM_TYPE"                 
#define KEY_DOMAIN_OX                "DOMAIN_OX"                
#define KEY_DOMAIN_OY                "DOMAIN_OY"                
#define KEY_DOMAIN_NX                "DOMAIN_NX"                
#define KEY_DOMAIN_NY                "DOMAIN_NY"                
#define KEY_DOMAIN_DX                "DOMAIN_DX"                
#define KEY_DOMAIN_DY                "DOMAIN_DY"                
#define KEY_DOMAIN_GRID_DIR          "DOMAIN_GRID_DIR"          
#define KEY_DOMAIN_ZREF              "DOMAIN_ZREF"              
#define KEY_DOMAIN_SLOPE             "DOMAIN_SLOPE"             
#define KEY_CHNL_WIDTH               "CHNL_WIDTH"               
#define KEY_CHNL_MAX_DEPTH           "CHNL_MAX_DEPTH"           
#define KEY_CHNL_WAVELENGTH          "CHNL_WAVELENGTH"          
#define KEY_CHNL_MARGIN              "CHNL_MARGIN"              
#define KEY_CHNL_FLW_DIR             "CHNL_FLW_DIR"             
#define KEY_CHNL_SCALE_DIST          "CHNL_SCALE_DIST"          
#define KEY_CHNL_SCALE_CONST         "CHNL_SCALE_CONST"         
#define KEY_CHNL_SCALE_MIN           "CHNL_SCALE_MIN"           
#define KEY_CHNL_SCALE_MAX           "CHNL_SCALE_MAX"           
#define KEY_CHNL_SCALE_NORM_MEAN     "CHNL_SCALE_NORM_MEAN"     
#define KEY_CHNL_SCALE_NORM_STDEV    "CHNL_SCALE_NORM_STDEV"    
#define KEY_CHNL_SCALE_LOGNORM_MEAN  "CHNL_SCALE_LOGNORM_MEAN"  
#define KEY_CHNL_SCALE_LOGNORM_STDEV "CHNL_SCALE_LOGNORM_STDEV" 
#define KEY_DYN_EMAP_ACTIVE          "DYN_EMAP_ACTIVE"          
#define KEY_EROD_COEF                "EROD_COEF"                
#define KEY_EROD_OFFSET              "EROD_OFFSET"              
#define KEY_ZUL_TYPE                 "ZUL_TYPE"                 
#define KEY_AV_LV_OB                 "AV_LV_OB"                 
#define KEY_AV_NO_SPLAY              "AV_NO_SLAY"               
#define KEY_AV_LOC_FREQ              "AV_LOC_FREQ"              
#define KEY_AV_LOC_PERIOD            "AV_LOC_PERIOD"            
#define KEY_AV_LOC_POISSON           "AV_LOC_POISSON"           
#define KEY_AV_LOC_PROB1             "AV_LOC_PROB1"             
#define KEY_AV_LOC_PROB2             "AV_LOC_PROB2"             
#define KEY_AV_REG_FREQ              "AV_REG_FREQ"              
#define KEY_AV_REG_PERIOD            "AV_REG_PERIOD"            
#define KEY_AV_REG_POISSON           "AV_REG_POISSON"           
#define KEY_AG_TYPE                  "AG_TYPE"                  
#define KEY_AG_DRAP                  "AG_DRAP"                  
#define KEY_AG_EXP_DEC_THICK         "AG_EXP_DEC_THICK"         
#define KEY_AG_LV_WIDTH              "AG_LV_WIDTH"              
#define KEY_AG_RAND_LAMBDA           "AG_RAND_LAMBDA"           
#define KEY_AG_OB_FREQ               "AG_OB_FREQ"               
#define KEY_AG_OB_PERIOD             "AG_OB_PERIOD"             
#define KEY_AG_OB_POISSON            "AG_OB_POISSON"            
#define KEY_AG_OB_DIST               "AG_OB_DIST"               
#define KEY_AG_OB_CONST              "AG_OB_CONST"              
#define KEY_AG_OB_MIN                "AG_OB_MIN"                
#define KEY_AG_OB_MAX                "AG_OB_MAX"                
#define KEY_AG_OB_NORM_MEAN          "AG_OB_NORM_MEAN"          
#define KEY_AG_OB_NORM_STDEV         "AG_OB_NORM_STDEV"         
#define KEY_AG_OB_LOGNORM_MEAN       "AG_OB_LOGNORM_MEAN"       
#define KEY_AG_OB_LOGNORM_STDEV      "AG_OB_LOGNORM_STDEV"      
#define KEY_AG_OB_PEAT               "AG_OB_PEAT"               
#define KEY_AG_OB_PEL_DRAP_INT       "AG_OB_PEL_DRAP_INT"       
#define KEY_AG_EP_INIT_ELEV          "AG_EP_INIT_ELEV"          
#define KEY_AG_EP_FREQ               "AG_EP_FREQ"               
#define KEY_AG_EP_PERIOD             "AG_EP_PERIOD"             
#define KEY_AG_EP_POISSON            "AG_EP_POISSON"            
#define KEY_AG_EP_DIST               "AG_EP_DIST"               
#define KEY_AG_EP_CONST              "AG_EP_CONST"              
#define KEY_AG_EP_MIN                "AG_EP_MIN"                
#define KEY_AG_EP_MAX                "AG_EP_MAX"                
#define KEY_AG_EP_NORM_MEAN          "AG_EP_NORM_MEAN"          
#define KEY_AG_EP_NORM_STDEV         "AG_EP_NORM_STDEV"         
#define KEY_AG_EP_LOGNORM_MEAN       "AG_EP_LOGNORM_MEAN"       
#define KEY_AG_EP_LOGNORM_STDEV      "AG_EP_LOGNORM_STDEV"      
#define KEY_COND_WELL_RELAX          "COND_WELL_RELAX"          
#define KEY_COND_WELL_CLOSURE_ON     "COND_WELL_CLOSURE_ON"     
#define KEY_COND_WELL_STOP_HONORED   "COND_WELL_STOP_HONORED"   
#define KEY_COND_WELL_CLOSURE        "COND_WELL_CLOSURE"        
#define KEY_SED_LOAD_ACTIVE          "SED_LOAD_ACTIVE"          
#define KEY_SED_LOAD_MEAN            "SED_LOAD_MEAN"            
#define KEY_SED_LOAD_VOL             "SED_LOAD_VOL"             
 
 
// Parameters default values //
// WARNING: in accordance with nexus default parameters !
 
// Parameters default values
#define DEF_SIM_SEED                 165426111
#define DEF_SIM_TYPE                 MD_MENU_FLUVIAL
#define DEF_DOMAIN_OX                0.    // 0 m
#define DEF_DOMAIN_OY                0.    // 0 m
#define DEF_DOMAIN_NX                251   // 251 # of nodes along Ox axis
#define DEF_DOMAIN_NY                201   // 201 # of nodes along Oy axis
#define DEF_DOMAIN_DX                10.   // 10 m (w / 3)
#define DEF_DOMAIN_DY                10.   // 10 m (w / 3)
#define DEF_DOMAIN_GRID_DIR          0.    // 0 degrees
#define DEF_DOMAIN_ZREF              0.    // 0 m
#define DEF_DOMAIN_SLOPE             0.001 // 1 m / 1 km
#define DEF_CHNL_WIDTH               30.   // 30 m
#define DEF_CHNL_MAX_DEPTH           3.    // 3 m (Which gives mean depth = 2 m)
#define DEF_CHNL_WAVELENGTH          375.  // 375 m (Which gives a default Cf = 0.16)
#define DEF_CHNL_MARGIN              12.   // 12 x w
#define DEF_CHNL_FLW_DIR             90.   // 90 degrees
#define DEF_CHNL_SCALE_DIST          DIST_CONSTANT
#define DEF_CHNL_SCALE_CONST         100.  // 100 %
#define DEF_CHNL_SCALE_MIN           80.   //  80 %
#define DEF_CHNL_SCALE_MAX           120.  // 120 %
#define DEF_CHNL_SCALE_NORM_MEAN     100.  // 100 %
#define DEF_CHNL_SCALE_NORM_STDEV    20.   //  20 %
#define DEF_CHNL_SCALE_LOGNORM_MEAN  100.  // 100 %
#define DEF_CHNL_SCALE_LOGNORM_STDEV 20.   //  20 %
#define DEF_EROD_COEF                4*MD_E_FAC // 4.e-8 (TODO : [Cleaning] : Which gives MigRate = XXXm / it)
#define DEF_EROD_OFFSET              0     // 0 m
#define DEF_ZUL_TYPE                 0     // None
#define DEF_AV_LV_OB                 false
#define DEF_AV_NO_SPLAY              false
#define DEF_AV_LOC_FREQ              FREQ_POISSON
#define DEF_AV_LOC_PERIOD            280   // 280 it
#define DEF_AV_LOC_POISSON           280   // 280 it
#define DEF_AV_LOC_PROB1             0.5   // 50 % of chances
#define DEF_AV_LOC_PROB2             0.9   // 90 % of chances
#define DEF_AV_REG_FREQ              FREQ_POISSON
#define DEF_AV_REG_PERIOD            500   // 500 it
#define DEF_AV_REG_POISSON           500   // 500 it
#define DEF_AG_TYPE                  AGG_OB
#define DEF_AG_EXP_DEC_THICK         919.  // 919 m (~ w x 30)
#define DEF_AG_LV_WIDTH              6.    // 6w m
#define DEF_AG_RAND_LAMBDA           true
#define DEF_COND_WELL_RELAX          20.   // 20 %
#define DEF_COND_WELL_CLOSURE_ON     false
#define DEF_COND_WELL_STOP_HONORED   false
#define DEF_COND_WELL_CLOSURE        0.1   // 10 cm
#define DEF_AG_EP_INIT_ELEV          9999. // 9999 m
#define DEF_AG_EP_FREQ               FREQ_NEVER
#define DEF_AG_EP_PERIOD             70   // 70 it
#define DEF_AG_EP_POISSON            70   // 70 it
#define DEF_AG_EP_DIST               DIST_CONSTANT
#define DEF_AG_EP_CONST              0.1  // 10 cm
#define DEF_AG_EP_MIN                0.07 //  7 cm
#define DEF_AG_EP_MAX                0.13 // 13 cm
#define DEF_AG_EP_NORM_MEAN          0.1  // 10 cm
#define DEF_AG_EP_NORM_STDEV         0.03 //  3 cm
#define DEF_AG_EP_LOGNORM_MEAN       0.1  // 10 cm
#define DEF_AG_EP_LOGNORM_STDEV      0.03 //  3 cm
#define DEF_AG_OB_PEAT               0.   // 0 m
#define DEF_AG_OB_PEL_DRAP_INT       0.   // 0 m / 10000 it
#define DEF_AG_OB_FREQ               FREQ_POISSON
#define DEF_AG_OB_PERIOD             70   // 70 it
#define DEF_AG_OB_POISSON            70   // 70 it
#define DEF_AG_OB_DIST               DIST_NORMAL
#define DEF_AG_OB_CONST              0.1  // 10 cm
#define DEF_AG_OB_MIN                0.07 //  7 cm
#define DEF_AG_OB_MAX                0.13 // 13 cm
#define DEF_AG_OB_NORM_MEAN          0.1  // 10 cm
#define DEF_AG_OB_NORM_STDEV         0.03 //  3 cm
#define DEF_AG_OB_LOGNORM_MEAN       0.1  // 10 cm
#define DEF_AG_OB_LOGNORM_STDEV      0.03 //  3 cm
#define DEF_AG_DRAP                  DRAP_DRAPING
#define DEF_SED_LOAD_ACTIVE          false
#define DEF_SED_LOAD_MEAN            0.01 // TODO : Really ?
#define DEF_NEXUS_ISBX_TEMPLATE      { 80,  50,  80, 110 }
#define DEF_NEXUS_NG_TEMPLATE        { 50., 30., 50., 70.}
 
// Default values for new turbidites scenario (keep Lemay 2018 (PhD) values)
//
// Lemay 2020 proposes another scenario but with too big turbidites
// See : Geomorphic variability of submarine channelized systems along continental margins:
//       Comparison with fluvial meandering channels
//       Fig 7 and 8 - Sediment Conduit : Leveed Channels
// So we choose to keep the Martin's PhD scenario.
//
// The default Erodibility Coefficient 5.10-9 can be increased to reduce simulation time.
// In that case, one iteration in turbidites doesn't correspond to 1 year anymore.
// Default avulsion and overbank periods must be updated accordingly.
// Keep in mind that for a same migration rate between fluvial and turbidites,
// the aggradation rate should be 10 times bigger in turbidites.
//
#define DEF_TURBI_DOMAIN_NX            251   // 251 # of nodes along Ox axis
#define DEF_TURBI_DOMAIN_NY            201   // 201 # of nodes along Oy axis
#define DEF_TURBI_DOMAIN_DX            200.  // 200 m (w/4)
#define DEF_TURBI_DOMAIN_DY            200.  // 200 m (w/4)
#define DEF_TURBI_DOMAIN_SLOPE         0.005 // 5 m / 1 km
#define DEF_TURBI_CHNL_WIDTH           800.  // 800 m (1200m in Lemay 2020)
#define DEF_TURBI_CHNL_WAVELENGTH      4000. // 4000 m (4200m in Lemay 2020)
#define DEF_TURBI_CHNL_MAX_DEPTH       40.   // 40 m (mean depth = 27m, 34m in Lemay 2020)
#define DEF_TURBI_CHNL_MARGIN          5.    // 5w m
#define DEF_TURBI_EROD_COEF            4*MD_E_FAC // 4.e-8
#define DEF_TURBI_COND_WELL_CLOSURE    1.4   // 1.4 m
#define DEF_TURBI_AG_LV_WIDTH          2.    // 2w m
#define DEF_TURBI_AG_EXP_DEC_THICK     8107. // 8107 m
#define DEF_TURBI_AG_EP_PERIOD         24    // 24 it
#define DEF_TURBI_AG_EP_POISSON        24    // 24 it
#define DEF_TURBI_AG_EP_CONST          2.0   // 2.0 m
#define DEF_TURBI_AG_EP_MIN            1.4   // 1.4 m
#define DEF_TURBI_AG_EP_MAX            2.6   // 2.6 m
#define DEF_TURBI_AG_EP_NORM_MEAN      2.0   // 2.0 m
#define DEF_TURBI_AG_EP_NORM_STDEV     0.6   // 60 cm
#define DEF_TURBI_AG_EP_LOGNORM_MEAN   2.0   // 2.0 m
#define DEF_TURBI_AG_EP_LOGNORM_STDEV  0.6   // 60 cm
#define DEF_TURBI_AG_OB_CONST          2.0   // 2.0 m
#define DEF_TURBI_AG_OB_MIN            1.4   // 1.4 m
#define DEF_TURBI_AG_OB_MAX            2.6   // 2.6 m
#define DEF_TURBI_AG_OB_NORM_MEAN      2.0   // 2.0 m
#define DEF_TURBI_AG_OB_NORM_STDEV     0.6   // 60 cm
#define DEF_TURBI_AG_OB_LOGNORM_MEAN   2.0   // 2.0 m
#define DEF_TURBI_AG_OB_LOGNORM_STDEV  0.6   // 60 cm
#define DEF_TURBI_AG_OB_PERIOD         24    // 24 it
#define DEF_TURBI_AG_OB_POISSON        24    // 24 it
#define DEF_TURBI_AG_DRAP              DRAP_PELAGIC
#define DEF_TURBI_AV_LOC_PERIOD        288   // 288 it
#define DEF_TURBI_AV_LOC_POISSON       288   // 288 it
#define DEF_TURBI_AV_REG_PERIOD        470   // 470 it
#define DEF_TURBI_AV_REG_POISSON       470   // 470 it
#define DEF_TURBI_NEXUS_ISBX_TEMPLATE  { 60,  30, 60, 90 }
#define DEF_TURBI_SED_LOAD_MEAN        0.1
#define DEF_TURBI_NEXUS_NG_TEMPLATE    {  6., 3., 6., 10.} // Lemay 2016
 
 
// Parameters range values //
 
// Parameters range values
#define MIN_SIM_SEED                 1
#define MIN_DOMAIN_NX                20            // 20 nodes
#define MIN_DOMAIN_NY                20            // 20 nodes
#define MAX_DOMAIN_NX                5000          // 5000 nodes
#define MAX_DOMAIN_NY                5000          // 5000 nodes
#define MIN_DOMAIN_DX                1             // 1 m
#define MIN_DOMAIN_DY                1             // 1 m
#define MIN_DOMAIN_GRID_DIR          -360          // degrees
#define MAX_DOMAIN_GRID_DIR          +360          // degrees
#define MIN_DOMAIN_SLOPE             0.0001        // 10 cm / 1 km
#define MAX_DOMAIN_SLOPE             0.1           // 100 m / 1 km
#define MIN_CHNL_WIDTH               5.            // 5 m
#define MAX_CHNL_WIDTH               500.          // 500 m
#define MIN_CHNL_WAVELENGTH          10.           // 10 m
#define MAX_CHNL_WAVELENGTH          10000.        // 10 km
#define MIN_CHNL_MAX_DEPTH           0.5           // 0.5 m (Mean depth of 0.33 m)
#define MAX_CHNL_MAX_DEPTH           15.           // 15 m (Mean depth of 10 m)
#define MIN_CHNL_MARGIN              1             // 1 time channel width
#define MAX_CHNL_MARGIN              100           // 100 times channel width
#define MIN_CHNL_FLW_DIR             -360.         // degrees
#define MAX_CHNL_FLW_DIR             +360.         // degrees
#define MIN_CHNL_SCALE_CONST         80.           // 80 %
#define MAX_CHNL_SCALE_CONST         120.          // 120 %
#define MIN_CHNL_SCALE_MIN           50.           // 50 %
#define MAX_CHNL_SCALE_MIN           100.          // 100 %
#define MIN_CHNL_SCALE_MAX           100.          // 100 %
#define MAX_CHNL_SCALE_MAX           150.          // 150 %
#define MIN_CHNL_SCALE_NORM_MEAN     90.           // 90 %
#define MAX_CHNL_SCALE_NORM_MEAN     110.          // 110 %
#define MIN_CHNL_SCALE_NORM_STDEV    0.            // No variation
#define MAX_CHNL_SCALE_NORM_STDEV    40.           // +-40%
#define MIN_CHNL_SCALE_LOGNORM_MEAN  90.           // 90 %
#define MAX_CHNL_SCALE_LOGNORM_MEAN  110.          // 110 %
#define MIN_CHNL_SCALE_LOGNORM_STDEV 0.            // No variation
#define MAX_CHNL_SCALE_LOGNORM_STDEV 40.           // +-40%
#define MIN_EROD_COEF                0.            // No erodibility
#define MAX_EROD_COEF                20*MD_E_FAC   // 2.e-7
#define MIN_AV_LOC_PERIOD            1             // 1 iteration
#define MIN_AV_LOC_POISSON           1             // 1 iteration
#define MIN_AV_LOC_PROB1             0             // 0 % of chance
#define MAX_AV_LOC_PROB1             1             // 100 % of chances
#define MIN_AV_LOC_PROB2             0             // 0 % of chances
#define MAX_AV_LOC_PROB2             1             // 100 % of chances
#define MIN_AV_REG_PERIOD            1             // 1 iteration
#define MIN_AV_REG_POISSON           1             // 1 iteration
#define MIN_AG_EXP_DEC_THICK         5.            // 5 m
#define MIN_AG_LV_WIDTH              0.05          // 0.05 times channel width
#define MIN_COND_WELL_RELAX          0             // 0 %
#define MAX_COND_WELL_RELAX          100           // 100 %
#define MIN_COND_WELL_CLOSURE        0.01          // 1 cm
#define MAX_COND_WELL_CLOSURE        1.            // 1 m
#define MIN_AG_EP_PERIOD             1             // 1 iteration
#define MIN_AG_EP_POISSON            1             // 1 iteration
#define MIN_AG_EP_CONST              0.02          // 2 cm
#define MIN_AG_EP_MIN                0.02          // 2 cm
#define MIN_AG_EP_MAX                0.03          // 3 cm
#define MIN_AG_EP_NORM_MEAN          0.02          // 2 cm
#define MIN_AG_EP_NORM_STDEV         0             // No variation
#define MIN_AG_EP_LOGNORM_MEAN       0.02          // 2 cm
#define MIN_AG_EP_LOGNORM_STDEV      0             // No variation
#define MIN_AG_OB_PEAT               0             // 0 %
#define MAX_AG_OB_PEAT               99            // 99 %
#define MIN_AG_OB_PEL_DRAP_INT       0.            // 0. m / 10000 it (no pelagic / draping)
#define MAX_AG_OB_PEL_DRAP_INT       100           // 100 m / 10000 it
#define MIN_AG_OB_PERIOD             1             // 1 iteration
#define MIN_AG_OB_POISSON            1             // 1 iteration
#define MIN_AG_OB_CONST              0.            // 0 cm (Authorized for steady migration)
#define MIN_AG_OB_MIN                0.            // 0 cm (Authorized for steady migration)
#define MIN_AG_OB_MAX                0.03          // 3 cm
#define MIN_AG_OB_NORM_MEAN          0.            // 0 cm (Authorized for steady migration)
#define MIN_AG_OB_NORM_STDEV         0.            // No variation
#define MIN_AG_OB_LOGNORM_MEAN       0.            // 0 cm (Authorized for steady migration)
#define MIN_AG_OB_LOGNORM_STDEV      0.            // No variation
#define MIN_SED_LOAD_MEAN            0.001         // 0.1% (default 1%)
#define MIN_SED_LOAD_VOL             0.            // 0% (minimal volume in a class)
 
// Ranges for turbidites parameters // range has been modified for sensitivity analysis
#define MAX_TURBI_DOMAIN_SLOPE        0.5   // 0.5 km/1 km
#define MIN_TURBI_CHNL_WIDTH          20.   // 20 m
#define MAX_TURBI_CHNL_WIDTH          4000. // 4 km
#define MIN_TURBI_CHNL_WAVELENGTH     100.  // 100 m
#define MAX_TURBI_CHNL_WAVELENGTH     20000.// 20 km
#define MIN_TURBI_CHNL_MAX_DEPTH      10.    // 10 m
#define MAX_TURBI_CHNL_MAX_DEPTH      500.  // 500 m
#define MIN_TURBI_COND_WELL_CLOSURE   0.1   // 10 cm
#define MAX_TURBI_COND_WELL_CLOSURE   10.   // 10 m
#define MIN_TURBI_AG_EXP_DEC_THICK    100.  // 100 m
 
// Non expert user parameters range values
#define MIN_NEXUS_ISBX    20         // Minimum index: ribbon type
#define MAX_NEXUS_ISBX    160        // Maximum index: sheet type
#define MIN_NEXUS_NG      1.         // 1%
#define MAX_NEXUS_NG      99.        // 99%
 
#endif // MCRC_DEFINES_H