source: branches/GRISLIv3/Param/anteis1_param_list_Cyc-Scho.dat @ 496

 !  Parametres du run :CY-S-ANT    anteis1
 !  sources dans  ../../SOURCES
 
!___________________________________________________________
&runpar                  ! nom du bloc  parametres du run
 
 runname      =  "RUNNAME"
 icompteur    =            1
 iout         =            2
 reprcptr     = "../../../Fichier-CPTR/TB40S252-k126-cptr.nc"
 itracebug    =            0
 num_tracebug  =          163
 comment_run   = "Cycle Antarctique Schoof"
 /
! runname    : nom de l experience (8 caracteres)  
! icompteur  : reprise dans un fichier  0 -> non, 1 -> oui, 2 -> T et Hw
! icompteur  : 3-> T seulement
! iout   : 1-> sortie cptr pour reprise,2 -> sortie nc pour reprise 
! reprcptr   : nom du fichier restart
! itracebug: 1-> ecriture de traces au debut des routines
! num_tracebug   numero d unite ecriture traces
! comment_run: commentaire court sur le run
 
!___________________________________________________________
&grdline                 ! bloc grounding line
 
 igrdline     =            0
 Schoof       =            1
 ibmelt_inv   =            0
 /
! igrdline   : 1 ligne d echouage fixée, sinon 0
! Schoof     : 0 pas de Schoof, 1 flux de Schoof
! ibmelt_inv : 0 cas std, 1 inversion du bmelt (avec igrdline=1)
 
!___________________________________________________________
&timesteps                 ! bloc time steps
 
 tend      =   RUNTEND    
 tbegin    =   RUNTBEGIN    
 dtmin     =   2.0000001E-03
 dtmax     =    1.000000    
 dtt       =    5.000000    
 testdiag  =   1.6000001E-02
 /
! tous les temps en annees. tbegin et tend : debut et fin du run 
! pour equation masse, pas de temps mini -> dtmin, maxi -> dtmax
! dtt : pas de temps long
! testdiag, pour gérer le pas de temps dynamique dt
! ordres de grandeur (a moduler selon dx) : 
! 40 km dtmin=2.e-2, dtmax=1., dtt=5., tesdiag=0.02
 
!___________________________________________________________
&topo_file                                  !  input_topo  
topo_ref = SHB_bedmap2_Ant40_relax20ans.nc                                                                     
topo_dep = SHB_bedmap2_Ant40_relax20ans.nc                                                                     
grid_topo =coord-Ant-40km.dat                                                                                  
ghf_fich = ghf_Ant_Lebrocq_40km.nc                                                                             
 /
! topo_ref= topo ref isostasie
! topo_dep= topo de depart
! grid_topo : fichier i,j,x,y,lon,lat
! ghf_fich  : fichier flux geothermique
 
!  module climat_Grice2sea_years_perturb_mod                     
 &LAPSE_RATES
 T_LAPSE_RATE    = -8.0000004E-03
 /
!T_lapse_rate       = lapse rate temp annuelle                   
!________________________________________________________________
! module climat_Grice2sea_years_perturb_mod                      
 &CLIM_PERT_MASSB
 COEFT   =   1.000000    ,
 FILFORC = forcage-429k-2007-Jouzel-Waelbroeck_AMOC.dat                                                                              
                             ,
 PERTSMB =           1,
 RAPSMB  =  7.0000000E-02,
 PERTBMB =           1,
 COEFBMB =   10.00000    
 /
! coefT : coef amplitude perturbation T
! filforc : fichier de forcage temporel
! pertsmb : 0 SMB fixe, 1 SMB exp temperature
! rapsmb : coef variation SMB
! pertbmb : 0 bmelt fixe, 1 bmelt variable
! coefbmb : coef variation bmelt
!________________________________________________________________
!  module climat_Grice2sea_years_perturb_mod                     
 &CLIM_SMB_T_GEN
 SMB_FILE        = RACMO2.3p2_ANT27_smb_ltm_1979_2016_ant40.grd                                                        ,
 COEF_SMB_UNIT   =   1.000000    ,
 TEMP_ANNUAL_FILE        = RACMO2.3p2_ANT27_t2m_ltm_1979_2016_ant40.grd                                                        
 /
! smb_file          = fichier SMB (kg/m2/an)                     
! coef_smb_unit     = coef passage m glace/an  (1/910 ou 1/918)  
! temp_annual_file  = Temp moy annuelle  (°C)                    
!________________________________________________________________
!  module climat_Grice2sea_years_perturb_mod                            
 &CLIM_SNAP
 NB_SNAP =         200,
 TIME_DEPART_SNAPS       =   2000000.    ,
 ECART_SNAP      =   1.000000    ,
 FILE_SMB_SNAP   = smb_ECHAM5_A1B_MAR_highres_all_time.nc                                                              ,
 MASSB_TIME      =           0
 /
! nb_snap           = nombre de snapshots                               
! time_depart_snaps = debut du forçage                                  
! ecart_snap        = ecart entre les snapshots                         
! file_smb_snap     = fichier serie temp anomalie SMB de GCM            
! massb_time        = 0:fixe, 1:interp snapshots, 2:snapsh+interp vert  
!_______________________________________________________________________
!___________________________________________________________
&eaubasale1              ! nom du premier bloc eau basale 
 
 ecoulement_eau =  T
 hwatermax      =    5000.000    
 infiltr        =   1.0000000E-03
 /
! ecoulement eau : .false. -> modele bucket, sinon equ. diffusion
! hwatermax :  hauteur d eau basale maximum dans le sediment (m)
! infiltr est la quantite d eau qui peut s infiltrer dans le sol (m/an)
 
!___________________________________________________________
&param_hydr             ! nom du  bloc parametres hydrauliques 
 
 hmax_till      =    20.00000    
 poro_till      =   0.5000000    
 kond0          =   1.7748050E-04
 /
! hmax_till (m) : epaisseur max du sediment 
! poro_till : porosite du sediment 
! conductivite du sediment :  kond0 (m/s)
 
!  module bmelt_ant_regions_initmip 
 &BMELT_ANT_REG_INITMIP
 BMELT_REGIONS   = 2*0.5000000      ,  0.1000000    ,   1.000000    ,  0.3000000    , 2*1.000000       ,  0.5000000    ,
    1.000000    ,   2.000000    ,   1.000000    ,  0.8000000    ,  0.9000000    ,  0.2200000    ,  0.1200000    ,  0.2000000    ,
  2*0.5000000      ,
 BMGRZ_REGIONS   =   6.000000    , 2*1.500000       ,   6.000000    ,   3.000000    ,   6.500000    ,   7.000000    ,   4.000000    ,
    6.000000    ,   5.000000    ,   6.000000    ,   8.000000    ,   3.500000    ,  0.5400000    ,  0.4500000    ,  0.2500000    ,
    4.700000    ,  0.1000000    ,
 BMELT_TALUS     =   20.00000    ,
 BMGRZ_TALUS     =   20.00000    ,
 BMELT_COEF      =  0.9688509    ,
 FILE_NUMBER_SHELVES     = shelves-from-initMIP-40km.nc                                                                        ,
 FLAG_DIST       =           1,
 FILE_DIST_TALUS = distance_talus_1500m.grd                                                                            ,
 BMELT_DIST0     =   20.00000    
 /
! bmelt_regions    :  fonte basale sous shelves pour les 18 regions initMIP
! bmgrz_regions    :  fonte basale grounding zone pour les 18 regions initMIP
! bmelt_talus & bmgrz_talus  :  fonte basale apres talus cont
! bmelt_coef                 :  coef fonte (1 pour conserver val)
! file_numer_ice-shelves     : fichier zones ice shelves
! flag_dist : flag pour bmelt fnct distance talus
! file_dist_talus : fichier de distance talus
! bmelt_dist0 : coef de bmelt au talus
!_______________________________________________________________________
!  module bmelt_ant_regions_initmip                                     
 &BMELT_ANOM_INITMIP
 FILE_BMELT_ANOM = basal_melt_anomaly_16km.nc                                                                          ,
 BMELT_TIME      =           0
 /
! file_bmelt_anom   = fichier anomalie bmelt                            
! bmelt_time        = 0:fixe, 1:anomalies                               
!_______________________________________________________________________
 &LOIDEF_1
 EXPOSANT_1      =   3.000000    ,
 TEMP_TRANS_1    =  -6.500000    ,
 ENHANC_FACT_1   =   2.732374    ,
 COEF_COLD_1     =  1.6599999E-16,
 Q_COLD_1        =   78200.00    ,
 COEF_WARM_1     =  2.0000000E-16,
 Q_WARM_1        =   95450.00    
 /
 &LOIDEF_2
 EXPOSANT_2      =   1.000000    ,
 TEMP_TRANS_2    =  -10.00000    ,
 ENHANC_FACT_2   =   1.000000    ,
 COEF_COLD_2     =  8.3129997E-08,
 Q_COLD_2        =   40000.00    ,
 COEF_WARM_2     =  8.3129997E-08,
 Q_WARM_2        =   60000.00    
 /
 !___________________________________________________________
 ! loi de deformation            module deformation_mod_2lois
 ! exposant (glen), temperature de transition (ttrans)
 ! enhancement factor (sf)
 ! pour les temperatures inf. a Temp_trans :
 !                            coef_cold (Bat1) et Q_cold (Q1)
 ! pour les temperatures sup. a Temp_trans :
 !                            coef_warm (Bat2) et Q_warm (Q2)
 !___________________________________________________________
!___________________________________________________________
&meca_SIA_L1             ! bloc resol_meca 
 
 i_resolmeca =            2
 /
! i_resolmeca type d association entre SIA et L1
! i_resolmeca=0  chacun dans sa zone
! i_resolmeca=1  dans les zones stream, addition si uxdef > uxL1 (MIS11 Cairns)
! i_resolmeca=2  addition systematique dans les zones stream 
 
!___________________________________________________________
&furst_schoof                                ! nom du bloc 
 
 frot_coef  =   0.6000000    
 gr_select  =            2
 /
 ! frot_coef : solid friction law coef 0.6 in GMD 2018
 ! gr_select = 1 : Tsai , 2 : Schoof
!___________________________________________________________
&diagno_rheol              ! nom du bloc  diagno_rheol
 
 sf01           =   0.1250000    
 sf03           =   0.1250000    
 sffunc         =   0
 pvimin         =    1500.000    
 /
! coefficients par rapport a la loi glace posee 
! sf01 : coefficient viscosite loi lineaire 
! sf03 : coefficient viscosite loi n=3 
! sffunc : 1 for visco coef that depends on enh factor
! pvimin : valeur de pvi pour les noeuds fictifs ~ 1.e3
! tres petit par rapport aux valeurs standards ~ 1.e10
 
 !___________________________________________________________
 !        spinup          module no_spinup               
 &SPINUP
 ISPINUP =           0
 /
 
 
 ! ispinup = 0     run standard
 ! ispinup = 1     equilibre temperature avec vitesses grisli 
 ! ispinup = 2     conservation de la masse avec vitesses bilan 
 ! ispinup = 3     equilibre temperature avec vitesses bilan
 
!___________________________________________________________
&drag_param_beta          ! nom du bloc dragging param beta
 
 beta_slope            =   1.1707746E-03
 beta_expo             =    1.000000    
 betamax               =    500000.0    
 betamin               =    1.000000    
 /
! Slope & expo of beta = - slope x Neff ** expo
 &MASS_CONSERV
 ADV_FRAC        =   2.000000    ,
 V_LIMIT =   3000.000    
 /
 ! Conservation de la masse avec equation advection-diffusion 
 ! la repartition depend de adv_frac
 ! >1  -> advection seule
 ! 0   -> diffusion seule
 ! 0<*<1   -> fraction de l advection
 ! -1 -> zones diffusion + zones advection
 ! V_limit depend de la calotte : 
 ! typiquement 3000 m/an en Antarctique, 10000 m/an au Groenland
 !___________________________________________________________
!___________________________________________________________
&calving              ! nom du bloc calving méthode Vincent 
 
 Hcoup_plateau  =    250.0000    
 Hcoup_abysses  =    250.0000    
 prof_plateau  =    50.00000    
 prof_abysses  =    1000.000    
 ifrange    =            4
 meth_hcoup =            0
 /
! Hcoup epaisseurs de coupure pour les zones peu prodondes et profondes
! Hcoup_plateau<Hcoup_abysses && prof_plateau<prof_abysses
! prof profondeur delimitant les zones peu prodondes et profondes
! ifrange=0 -> pas de traitement particulier sur les bords
! ifrange=1 -> traitement de Vincent avec ice shelves frangeants partout
! ifrange=2 -> ice shelves frangeants seulement si bm-bmelt positif
 ! meth_hcoup pour faire eventuellement varier Hcoup avec le climat