source: trunk/SOURCES/Fichiers-parametres/Templates-param_list/greeneem45_param_list.dat @ 10

  Parametres du run : greeneem45


!___________________________________________________________
&runpar                  ! nom du bloc  parametres du run

 runname      =  "nm45ts17"   ! 8 caracteres experiences greeneem
 icompteur    =    0        
 reprcptr     =  "../Fichier-CPTR/Gmint-15-steady.CPTR"   ! a voir comment je le change (aurel)
 itracebug    = 0
 num_tracebug = 166
 comment_run  = "neem v17 test execution"
 /
! 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      =      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      =   100.
 tbegin    =   0.
 dtmin     =   2.e-3
 dtmax     =   1.  
 dtt       =   5.   ! 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 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
!___________________________________________________________
&toponeem                            ! module lect_topo_greeneem
topo_ref  = MAR_alti_45km.dat
topo_dep  = groen_grisli_45km.dat
grid_topo = groen_grisli_grid_xylonlat_45km.dat
ghf_fich      = ghf_cism_45km.dat     !50.   
/                     
! topo_ref= topo actuelle  topo_dep = celle de depart du run
!___________________________________________________________
&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)
!___________________________________________________________

&bmelt_seuil                     ! module  bmelt_seuil_prof
bm_grz           = 3.
bmshelf_plateau  = 1.
bmshelf_abysses  = 10.
depth_talus      = -2000.
/                      
! 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
!___________________________________________________________ 

&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-5             ! 1.e-4  
tobmax          =       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 en Pa'
! toblim : (Pa) tes pour les iles '
!___________________________________________________________
&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

!____________________________________________________________
! loi de deformation 1          module deformation_mod_2lois
&loidef_1

exposant_1       = 3.
temp_trans_1     = -6.5
enhanc_fact_1    = 4.5       !3.0
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.0          
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.2
sf03           =        0.03
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

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

! kweert : coefficent, loigliss le type de loi
! coefbmax : facteur de normalisation pour influence eau
!___________________________________________________________ 
&clim_pert                             ! module climat_pertub_mois

coefT           =       0.
rappact         =       0.05                    ! 0.05
retroac         =       1
rapbmshelf      =       5.  
mincoefbmelt    =       0.
maxcoefbmelt    =       2.
filforc         = 'forcage-stationnaire_neem.dat'
/
!----------------------------------------------------------

&clim_init_mois                            ! module lect_clim_greeneem_mois_mod
filtr_t = "MAR_tt_70-79_45km.dat"
filtr_p = "MAR_precipes_70-79_45km.dat"

/
! un bloc de namelist par snapshot
!___________________________________________________________ 
&lapse_month                                          !  module lect_clim_greeneem_mois_mod
filtr_lr = "climato_actu_grad-temp_45km.dat"
/
! lapserate de la param gmint

!___________________________________________________________ 
& 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 T mensuelle                module ablation_month    

&ablation_month
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 

!----------------------------------------------------------