!  Parametres du run : Antar  Anteis1


!___________________________________________________________
&runpar                  ! nom du bloc  parametres du run

 runname      =  "M11An129"
 icompteur    =    1
 reprcptr     =  "../Fichier-CPTR/M11An095-k430.CPTR"    ! "../Fichier-CPTR/M11An095-k400.CPTR"
 itracebug    =  0
 num_tracebug = 166
 /
! runname    : nom de l experience (8 caracteres)
! icompteur  : reprise dans un fichier  0 -> non, 1 -> oui, 2 -> T et Hwat
!              3-> T seulement
! reprcptr   : nom du fichier

!___________________________________________________________
&netcdf                  ! bloc netcdf

 dtncdf      =       500.
 /
! dtncdf :espace en temps entre les sorties
 
!___________________________________________________________
&grdline                 ! bloc grounding line

 igrdline     =            0
 /
! igrdline :  1 ligne d echouage fixée, sinon 0

!___________________________________________________________
&timesteps                ! bloc timestep

 tend      =   -380000.
 tbegin    =  -650000.	
 dtmin     =   2.e-3
 dtmax     =   1.   ! 1.
 dtt	   =   10.   ! 5
 testdiag  =   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 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
 
!___________________________________________________________
&eaubasale1              ! nom du premier bloc eau basale 
 
 ecoulement_eau =  T
 hwatermax      =    5000.000    
 infiltr        =   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          =   1.000000E-06
 
! hmax_till (m) : epaisseur max du sediment 
! poro_till : porosite du sediment 
! conductivite du sediment :  kond0 (m/s)

!____________________________________________________________
&drag_hwat_cont              ! nom du bloc dragging hwater contigu

hwatstream	=	50.  
cf		=	1.e-6             ! 1.e-5             ! 1.e-4  
tobmax		=	200.              ! 1000.  
toblim		=	0.7e5              ! 0.25e5 
/
! hwatstream (m) :  critere de passage en stream en partant de la cote'
! si hwater > hwatstream '
! cf coefficient de la loi de frottement fonction Neff'
! tobmax : (Pa) frottement maxi '
! toblim : (Pa) tes pour les iles '

!___________________________________________________________ 
&diagno_rheol              ! nom du bloc  diagno_rheol

sf01           = 	0.2
sf03           = 	0.0
/  
! coefficients par rapport a la loi glace posee                        
! sf01 : coefficient viscosite loi lineaire 
! sf03 : coefficient viscosite loi n=3 
!___________________________________________________________ 
&bmelt_ant_reg                                 ! nom du bloc 

bmelt_Ross	=	0.1	     ! 0.
bmgrz_Ross  	=	1.           ! 0.8      
bmelt_FRis	=	0.2          ! 0.3
bmgrz_FRis  	=	3.           ! 2.   
bmelt_Amery	=	0.5          ! 1.
bmgrz_Amery 	=	60.          ! 70.
bmelt_PIG	=	1.           ! 5.
bmgrz_PIG	=	7.           ! 10.           
bmelt_Pen	=	0.5	     !	1.	
bmgrz_Pen	=	3.	     ! 4           
bmelt_other	=	0.1	     ! 0.5
bmgrz_other     = 	0.1          ! 0.5 
bmelt_talus	=	20.
bmgrz_talus     =	20.  
ajout_bmelt 	= 	0.          ! ajoute 0.5 m/degre au dessus de Tafor_ajout
Tafor_ajout  	= 	-2.
/
!___________________________________________________________ 
&clim_pert                                    ! nom du bloc

rappact		=	0.05                    ! 0.05
retroac		=	1
rapbmshelf      =       5.  
mincoefbmelt 	= 	0.
maxcoefbmelt	 = 	2.
/
!----------------------------------------------------------