source: trunk/SOURCES/Fichiers-parametres/mism3_5_param_list.dat @ 10

!  Parametres du run : mismip_3D en 5 km


!___________________________________________________________
&runpar                       ! nom du bloc parametres du run
 runname      =  "MSM5_015"   ! 8 caracteres 
 icompteur    =  0     
 iout         =  2            
 reprcptr     = ""
 itracebug    = 1
 num_tracebug = 167     
 comment_run  = "mismip_3D :  avec le nouveau depart de Vincent : test opposx  "
! 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
 /
! igrdline :  1 ligne d echouage fixée, sinon 0
! paleo grounding line -> 2
!___________________________________________________________
&timesteps                ! bloc timestep

 tend      =   5.
 tbegin    =   0.    !1.e10      ! si tbegin > 1.e9 on prend le temps du fichier cptr
 dtmin     =   2.e-3
 dtmax     =   1.        ! 10.   ! 0.1
 dtt       =   1   ! 5
 testdiag  =   0.016      ! 0.025      ! 0.016
 /
! 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-2, dtmax=1., dtt=5., tesdiag=0.02
 
!___________________________________________________________
&mass_conserv            ! conservation de la masse
adv_frac        =  2.    ! 2-> advection seule
V_limit         = 30000. ! 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 =  F
 hwatermax      =    0.000    
 infiltr        =   1.0000
 /
! 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.000000E-06
/ 
! hmax_till (m) : epaisseur max du sediment 
! poro_till : porosite du sediment 
! conductivite du sediment :  kond0 (m/s)

!____________________________________________________________
! loi de deformation 1      module deformation_mod_2lois_isoth 
&loidef_1_iso
exposant_1       = 3.
enhanc_fact_1    = 1.         
coef_iso_1       = 1.E-25    ! en Pa m-3 s-1 
/
! pour MISMIP  : coef_iso_1 = 2* A * secyear**(1/3)
!  exposant (glen), coefficient et enhancement factor (sf) 
!________________________________________________________ 
! loi de deformation 2       module deformation_mod_2lois_isoth 
&loidef_2_iso
exposant_2       = 1.
enhanc_fact_2    = 1.         
coef_iso_2       = 0.
/
!  exposant (glen), coefficient et enhancement factor (sf) 
!____________________________________________________________
! loi de deformation 1          module deformation_mod_2lois
&loidef_1

exposant_1       = 3.
temp_trans_1     = -6.5
enhanc_fact_1    = 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    = 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           =        1      !     0.125  ! 1/8
sf03           =        1      !     0.125
pvimin         =        1.e4   ! 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
!___________________________________________________________ 

! glissement                   module  sliding_Bindschadler 
&slid_bindsh                    ! nom du bloc
 
kweert       =  0.             ! 5.e-11
loigliss     =  2
coefbmax     =  10. 

! kweert : coefficent, loigliss le type de loi
! coefbmax : facteur de normalisation pour influence eau
!___________________________________________________________ 
&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'
/
!----------------------------------------------------------
& 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 

!----------------------------------------------------------
! PDD base Tann et Tjuly                module ablation_ann    

&ablation_ann
Cice          = 0.008
Csnow         = 0.005
Sigma         = 5.
csi           = 0.6
/
! Cice and Csnow, melting factors for ice and snow 
! sigma variabilite Tday
! csi proportion of melted water that can refreeze 

!----------------------------------------------------------
!  module  lect_topo_mismip3
& topo_mismip3
topo_file = 'topo_mis3_stnd_LXly_05_y23_new.nc'   !'topo_mis3_stnd_LXly_05_y23.nc'
/
!----------------------------------------------------------
!  module  lect_clim_acc_T_ant_gen  
&climat_acc_T_gen

precip_file    = 'acca_15km.grd'                                   ! precipitation
coef_dens      = 1.                                                ! direct. en m glace 
temp_annual_file ='temp_15km.grd'                      ! annual surface temperature
/
!----------------------------------------------------------
&bmelt_seuil                     ! module  bmelt_seuil_prof
bm_grz           = 0.
bmshelf_plateau  = 0.
bmshelf_abysses  = 0..
depth_talus      = -3000.
/                      
! 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
!___________________________________________________________ 
&vitbil_upwind             ! nom du bloc vitbil calcule sur le mailles staggered
balance_Ux_file = 'Ux-balvel-01.grd'
balance_Uy_file = 'Uy-balvel-01.grd'
/
!balance_vel_file = 'Umoy-for-calc-beta_LBq15-06_ZBL_15km.dat'
! balance_vel_file = 'Uslid-for-calc-beta_LBq15-06_ZBL_15km.dat'
! balance_vel_file = 'Umoy-for-calc-beta_LBq15-06_ZBL_15km.dat'
! balance_vel_file = 'Vitbil-stag-pour-spinup_merge_runs_Lbq_ZBL_15km.dat'
! balance velocities on staggered grid Ux, Uy 
!___________________________________________________________

! module dragging_vit_bil_LBq_gen
 &drag_vit_bil_LBq_gen
 HWATSTREAM              =   50.00000    ,
 CF                      =  9.9999997E-06,
 TOBMAX                  =   5000.000    ,
 TOBLIM                  =   70000.00    ,
 SEUIL_VEL               =   200.0000    ,
 BALANCE_VEL_FILE        = 'gwavel-LeBrocq-grounded_ZBL_15km.dat' 
 /
! hwatstream (m) :  critere de passage en stream vitesses de bilan
!  si hwater > hwatstream 
! cf coefficient de la loi de frottement fonction Neff
! seulement pour les points cotiers
! tobmax : (Pa) frottement maxi sous les streams 
! toblim : (Pa)  pour les iles 
! seuil_vel (m/an) : seuil sur les vitesses pour definir le masque stream 
! balance_vel_file : nom du fichier qui contient les vitesse de bilan
!___________________________________________________________
&beta_prescr                                   ! dragging_prescr_beta'
beta_c_file = 'beta_first_optim.grd'           ! 'beta-driving-over-sliding.grd' 
beta_limgz = .5e5                              ! -1.e12 pour plastic ?    ! .5e6 en visqueux   
beta_min = 10.                                 ! for grounded ice (Pa)
beta_mult = 1                                  ! coefficient multiplicateur    
/
! read beta on centered grid
! beta_file : nom des fichiers qui contiennent les betamx et betamy
! beta-estime-run07.dat calcul direct beta sans modif 
! beta-staggered-15km-opt1.dat premier test en utilisant beta calcule  avec iterations
! above beta_limgz, gzmx is false if negative : no sliding
! betamx_file = not read anymore
! betamy_file = not read anymore
!___________________________________________________________
&beta_stag                           !read beta on staggered grid'
betamx_file = ''
betamy_file = ''
beta_limgz = 2.e3        
/

! above beta_limgz, gzmx is false
!___________________________________________________________
&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'
!
!___________________________________________________________