!  Parametres du run : Heino201    heino_1


!___________________________________________________________
&runpar                  ! nom du bloc  parametres du run

 runname      =  "Hudson15"
 icompteur    =            0
 reprcptr     = "Hudson15.CPTR"
 itracebug    =  0
 num_tracebug = 166

 comment_run  = "exactmnt=hudson14 , mais parametrisation de la fusion basale type jalvarez "
! runname    : nom de l experience (8 caracteres)
! icompteur  : reprise dans un fichier  0 -> non, 1 -> oui, 2 -> T seulement
! reprcptr   : nom du fichier

!___________________________________________________________
&netcdf                  ! bloc netcdf

 dtncdf      =       100
 /
! 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      =   100000.
 tbegin    =   0	
 dtmin     =   2.e-3
 dtmax     =   1.
 dtt	   =   10.
 testdiag  =   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        =   2.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          =   3.0000001E-06
 /
! 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           =       250
ifrange         =       4
meth_Hcoup      =       0
/

! Hcoup epaisseur de coupure
! 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
! ifrange=4 -> jalv: conseille par Catherine
!____________________________________________________________

!____________________________________________________________
 &diagno_rheol              ! nom du bloc  diagno_rheol

sf01           = 	0.2
sf03           = 	0.0
pvimin         =        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
!___________________________________________________________
&sliding_Heino    ! nom du  bloc parametres glissement Heino

 seuil_gliss    =       0.0001
 coefmax        =       5000.   ! meme valeur que hwatermax
 Cr             =       1.e5
! Cr             =       5.e5
 Cs             =       500.
 /
! seuil_gliss (m) : glissement si eau basale des 6 voisins > seuil_gliss
! coefmax : hauteur d eau max dans le sediment (=hwatermax)
! Cr : coeficient glissemnt sur hard rock (loi en puissance 3)
! Cs : coefficient glissemnt sur sediments (loi lineaire)

!___________________________________________________________
&drag_hudson    ! nom du  bloc parametres frottement Heino

 seuil_sedim    =       0.0001  ! 0.0001
 coefmax        =       5000.   ! meme valeur que hwatermax
 Cfs            =       500.
 /
! seuil_sedim (m) : frottement sediment si eau basale des 6 voisins > seuil_sedim
! coefmax : hauteur d eau max dans le sediment (=hwatermax)
! Cfs : coefficient frottement sur sediments

!___________________________________________________________
&drag_hwat_stream     ! nom du bloc dragging hwatermax stream

 hwatstream     =       49.9
 cf             =       -10.e-5
 tobmax         =       200.
 toblim         =       0.25e5
 Cfs            =       1000.
 seuil_sedim    =       0.0001
 coefmax        =       5000.   ! meme valeur que hwatermax
 /
! 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 '
! Cfs : coef de frottement sediment
! seuil_sedim : seuil pour passer en stream si zone sediment
! coefmax : hauteur d eau max dans le sediment (=hwatermax)

!----------------------------------------------------------
&drag_hudson_jorge              ! nom du bloc dragging hwater contigu

hwatstream      =       50.
cf              =       1.e-5             ! 1.e-4
tobmax          =       1000.
toblim          =       0.7e5              ! 0.25e5
seuil_sedim    =       0.00001             ! 0.0001
/
! 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 '

!___________________________________________________________

& meca_SIA_L1                        ! bloc resol_meca 

i_resolmeca 	=	0
/                            
! i_resolmeca type d association entre SIA et L1'
! i_resolmeca=0  chacun dans sa zone'
! i_resolmeca=1  dans ls zones stream, addition si uxdef > uxL1 (MIS11 Cairns)
! i_resolmeca=2  addition 0systematique dans les zones stream