source: trunk/Param/alps2_param_list.dat @ 436

!  Parametres du run : geoplace GI2S_5T Groenland Ice2sea 5 km cut Tamsin

!___________________________________________________________
&runpar                       ! nom du bloc parametres du run
 runname      =  "ALPS2-PI"   ! 8 caracteres 
 icompteur    =    0
 iout         =    2  
 reprcptr     =  "../../ALPS2-PI-part1/bin/ALPS2-PI_1.nc"
 itracebug    = 0
 num_tracebug = 166     
 comment_run  = "Alps 2 km"
/
! runname    : nom de l experience (8 caracteres)
! icompteur  : reprise dans un fichier  0 -> non, 1 -> oui, 2 -> T et Hwat
!              3-> T seulement
! iout       : 1-> sortie cptr pour reprise,2 -> sortie nc pour reprise 
! reprcptr   : nom du fichier
!___________________________________________________________
&grdline                 ! bloc grounding line

 igrdline     =            0
 Schoof       =            1
 ibmelt_inv   =            0
 /
! igrdline :  1 ligne d echouage fixée, 3 recul forcé, 0 sans traitement specifique, 2 paleo gr line
! Schoof   :  0 pas de Schoof, 1 flux de Schoof

!___________________________________________________________
&timesteps                ! bloc timestep

 tend      =   10000.
 tbegin    =   0.
 dtmin     =   1e-3
 dtmax     =   0.5
 dtt         =   1.
 testdiag  =   0.008
 /
! 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 gerer le pas de temps dynamique dt
! ordres de grandeur (a moduler selon dx) : 
! 40 km dtmin=2.e-3, dtmax=1., dtt=5., tesdiag=0.02

!___________________________________________________________
&topo_file
 topo_ref        = "SHB_SRTM_alps2.nc"
 topo_dep        = "SHB_SRTM_alps2.nc"
 grid_topo       = "coord-alps2.dat"
/
! topo_ref= topo ref isostasie
! topo_dep= topo de depart
! grid_topo : fichier i,j,x,y,lon,lat
!___________________________________________________________
&mass_conserv            ! conservation de la masse
adv_frac        =  2.    ! 2-> advection seule
V_limit         = 3000.  ! maximum velocity for mass conservation : depend de la calotte
/
! 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 advecttion'
! V_limit depend de la calotte : typiquement 3000 en Antarctique, 10000 au Groenland

!--------------------------------------------------------------------
&eaubasale1              ! nom du premier bloc eau basale 
 
 ecoulement_eau =  T
 hwatermax      =    5000.000    
 infiltr        =   1.e-3   !1.0000001E-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          =   0.000193342957286257   !10.000000E-05
/
! hmax_till (m) : epaisseur max du sediment 
! poro_till : porosite du sediment 
! conductivite du sediment :  kond0 (m/s)

!____________________________________________________________
&calving                        !   nom du bloc calving méthode Vincent

Hcoup_plateau   =      250. !50     !   tres petit quand shelves fixes sinon 250
Hcoup_abysses   =       250. !400     !   tres petit quand shelves fixes sinon 250
prof_plateau   =       50. !100     !   tres petit quand shelves fixes sinon 250
prof_abysses   =       1000. !800     !   tres petit quand shelves fixes sinon 250
ifrange         =       4
meth_Hcoup      =       0
/

! Hcoup epaisseurs de coupure pour les zones peu prodondes et profondes
! Hcoup_plateau<Hcoup_abysses && prof_plateau<prof_abysse
! 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
! ifrange=2 -> ice shelves frangeant seulement si bm-bmelt positif
! meth_hcoup pour faire eventuellement varier Hcoup avec le climat
!____________________________________________________________
! loi de deformation 1          module deformation_mod_2lois
&loidef_1

exposant_1       = 3.
temp_trans_1     = -6.5
enhanc_fact_1    = 5.0          ! Grice2sea 3.         
coef_cold_1      = 1.660E-16
Q_cold_1         = 7.820E+04
coef_warm_1      = 2.000E-16
Q_warm_1         = 9.545E+04
/
! 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)
!________________________________________________________ 
! loi de deformation 2          module deformation_mod_2lois
&loidef_2

exposant_2       = 1.
temp_trans_2     = -10.
enhanc_fact_2    = 1.            ! Grice2sea 3.         
coef_cold_2      = 8.313E-08
Q_cold_2         = 4.000E+04
coef_warm_2      = 8.313E-08
Q_warm_2         = 6.000E+04
/
! 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)
!___________________________________________________________ 
 
&diagno_rheol              ! nom du bloc  diagno_rheol

sf01           =        0.125   ! 0.125
sf03           =        0.125   ! 0.125
pvimin         =        1.5e3   ! 1.e3
/  
! coefficients par rapport a la loi glace posee                        
! sf01 : coefficient viscosite loi lineaire 
! sf03 : coefficient viscosite loi n=3 
! pvimin : valeur de pvi pour les noeuds fictifs ~ 1.e3
! tres petit par rapport aux valeurs standards ~ 1.e10
!___________________________________________________________ 

! furst_schoof
&furst_schoof

frot_coef = 0.6   !0.035
gr_select = 1
/
! frot_coef  = ',frot_coef
! gr_select  = ',gr_select
! gr_select = 1 : Tsai , 2 : Schoof'

!___________________________________________________________
&ablation                                     !  module ablation_mod'

pdd_type        = 0 
annual          = .false.
Cice            = 0.008
Csnow           = 0.003
csi             = 0.6
sigma_ice       = 5.
/
! pdd_type : 0 reeh, 1 Fausto, 2 Tarasov
! annual : T = annuel, F = mensuel
! Cice and Csnow, melting factors for ice and snow 
! sigma variabilite Tday
! csi proportion of melted water that can refreeze

!___________________________________________________________ 
&clim_pert                                    ! nom du bloc

coefT           =       1.
rappact         =       0.07                    ! 0.05->0.75  0.07 -> 0.5
retroac         =       1
rapbmshelf      =       5.  
mincoefbmelt    =       0.
maxcoefbmelt    =       2.
filforc         =     'forcage-900k-2007-Parrenin-Bassinot.dat' !'LGM_permanent.dat'
/
!___________________________________________________________ 
&clim_pert_massb                                   ! nom du bloc

coefT           =     0. ! 1.
filforc         =   'forcage-429k-2007-Jouzel-Waelbroeck_AMOC.dat' ! 'quiquet-et-al_CP2013_forcage.dat' !'forcage-900k-2007-Parrenin-Bassinot.dat' !'LGM_permanent.dat'
pertsmb         =     0  !1  !1
rapsmb          =     0. !0.07 !7
pertbmb         =     1
coefbmb   =     10.
/
! pertsmb = 1, smb pert via tpert / pertsmb = 0, smb = cst
!----------------------------------------------------------
&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 

!----------------------------------------------------------
! pour le module tracer (suivi d'isotopes)
&tracer
file_tr_dat     = "../INPUT/tracer.dat"
file_tr_out     = "../RESULTATS/tracer.out"
file_tr_dep     = "../RESULTATS/deparchive.out"
lambdab         = 0.0001
agemax          = 1000000
coefT_tra       = 1.
rappact_tra     = 0.07
filforc         = "quiquet-et-al_CP2013_forcage.dat" !"forcage-900k-2007-Parrenin-Bassinot.dat"
pert_type       = 0
/
!----------------------------------------------------------
! module climat_forcage
&clim_forcage
clim_ref_file = 'climato_monthly_t2m_precip_orog_ERA5-1979-2008-alps2.nc'
ntr = 2
ttr_temp = -21000, 0.
forcage_file = 'climato_monthly_t2m_precip_orog_ERA5-1979-2008-alps2_forclim.nc' , 'climato_monthly_t2m_precip_orog_ERA5-1979-2008-alps2_forclim.nc'
typerun = 1
lapserate = 0.007
rappact = 0.07 !0.11
filforc = 'forcage-429k-quiquet-et-al_CP2013_glacial-index_Waelbroeck_AMOC.dat'
pertbmb = 1
coefT   = 1.
coefbmb = 10. !10.
r_atmvar = 0.25
/
!----------------------------------------------------------
&imposed_regions_rsl               ! module  imposed_regions_rsl
nb_reg           = 1
filreg           = Hemin40_Hudson.grd
filforc          = rsl_hudson_forc.txt
/
!
!___________________________________________________________
&bmelt_seuil                     ! module  bmelt_seuil_prof
bm_grz           = 0.3      !3
bmshelf_plateau  = 0.2        !0.4
bmshelf_abysses  = 10.
depth_talus      = -1000.
/
! Pour l actuel : bm_grz a la grounding line
!                 bmshelf_plateau sur le plateau continental
!                 bmshelf_abysses pour les grandes profondeurs
!       depth_talus, negative, separation entre les 2 domaines
!___________________________________________________________ 
&drag_param_beta                              ! drag_beta_param_mod

beta_slope     = 0.001584297807517     !1.e-3
beta_expo      = 1.      !1.0
betamax    = 5.e5     ! maximum value of beta
betamin    = 100. !10.   ! minimum value of beta
coef_ile   = 1.
/
!
!___________________________________________________________
&drag_param_beta_sedim                       ! drag_beta_param_sedim_mod

beta_slope    = 0.001584297807517    !1.e-3
beta_expo     = 1.     !1.0
betamax    = 5.e5     ! maximum value of beta
betamin    = 100. !10.   ! minimum value of beta
coef_ile   = 1
file_sedim  = 'sedmap_laure16.nc'
seuil_sedim = 200.
coef_sedim  = 0.5e-1
/
!
!___________________________________________________________
&drag_coulomb_friction                        ! drag_coulomb_friction

cf            = 0.001584297807517    !1.e-3
m_nolin       = 1.
niter_nolin   = 1
betamax    = 5.e5     ! maximum value of beta
betamin    = 100.   ! minimum value of beta
bool_sedim = .F.
file_sedim  = 'sedmap_laure16.nc'
seuil_sedim = 200.
coef_sedim  = 0.5e-1
/
! cf: a friction coefficient (to be tuned)
! m_nolin: non-linear exponent, from 1 to infinity (put -1 for infinity)
!          m_nolin=1/q in Pattyn TC 2017, q in [0:1]
! niter_nolin: number of iterations to solve the non-linearity (expensive!)
!
!___________________________________________________________
&drag_coulomb_friction_simplhydro           ! drag_coulomb_friction_simplhydro

till_c0       = 0.
TFA_coeff     = 20.
Ob_coeff      = 0.001584297807517    !1.e-3
m_nolin       = 1.
niter_nolin   = 1
betamax    = 5.e5     ! maximum value of beta
betamin    = 100.   ! minimum value of beta
bool_sedim = .F.
file_sedim  = 'sedmap_laure16.nc'
seuil_sedim = 200.
coef_sedim  = 0.5e-1
/
! till_c0: sediment cohesion (usually = 0)
! TFA_coeff: till friction angle
! Ob_coeff: a friction coefficient
! m_nolin: non-linear exponent, from 1 to infinity (put -1 for infinity)
!          m_nolin=1/q in Pattyn TC 2017, q in [0:1]
! niter_nolin: number of iterations to solve the non-linearity (expensive!)
!
!___________________________________________________________
&spinup                                ! warning : 2 different modules
ispinup = 0
/
! with module no_spinup
! ispinup = 0     run standard ou calcul du beta
! ispinup = 1     temperature equilibrium with grisli velocities'
!
! with module spinup_vitbil
! ispinup = 2     conservation de la masse avec vitesses bilan '
! ispinup = 3     equilibre temperature avec vitesses bilan'
!
!___________________________________________________________