source: trunk/SOURCES/Fichiers-parametres/ant16km_param_list.dat @ 135

!  Parametres du run : Ant-16

!___________________________________________________________
&runpar                       ! nom du bloc parametres du run
 runname      =  "Ant16tut"   ! 8 caracteres 
 icompteur    =    1
 iout         =    2          
 reprcptr     =  "../../../Restarts/Ant16rlz_2.nc" !Ant16-05_1.nc" !Ant16-03+k30.nc"
 itracebug    = 0
 num_tracebug = 163     
 comment_run  = "un commentaire"
! 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
! ibmelt_inv : 0 cas std, 1 inversion du bmelt (avec igrdline=1)
!___________________________________________________________
&timesteps                ! bloc timestep

 tend      =   100.
 tbegin    =   0.         !1.e10      ! si tbegin > 1.e9 on prend le temps du fichier cptr
 dtmin     =   1.e-3
 dtmax     =   0.5 !10. !1.  
 dtt       =   1. !10. !5.  
 testdiag  =   0.005      !0.005 cat 15km   
 /
! 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_bedmap2_Ant16.nc" 
 topo_dep     = "SHB_bedmap2_Ant16.nc" 
 grid_topo    = "coord-Ant-16km.dat"
 ghf_fich     = "ghf-initMP-16km.nc"
/
! topo_ref= topo ref isostasie
! topo_dep= topo de depart
! grid_topo : fichier i,j,x,y,lon,lat
! ghf_fich  : fichier flux geothermique

!___________________________________________________________
&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          =   1.e-6   !1.000000E-06
 
! hmax_till (m) : epaisseur max du sediment 
! poro_till : porosite du sediment 
! conductivite du sediment :  kond0 (m/s)

!____________________________________________________________
&drag_neff_slope              ! nom du bloc dragging neff slope

cf              =       20.e-5  !2.e-5             ! 1.e-4
betamax         =       1000.
betamin         =       10.
toblim          =       0.7e5              ! 0.25e5
tostick         =       1.e5
seuil_neff      =       5.e5  !350.e5
coef_gz         =       10.
coef_ile        =       0.1
slope_fich      = "ETOPO1_slope_R2min_Ant40km.nc"
expo_slope      =       1.
pente_min      =        500.
pente_max      =        2000.
/
! cf coefficient de la loi de frottement fonction Neff'
! betamax : (Pa) frottement maxi '
! betamin : (Pa) frottement mini '
! toblim : (Pa) tes pour les iles '
! tostick : (Pa)  pour les points non flgzmx '
! seuil_neff (Pa) : seuil sur la pression effective pour avoir glissement'
! coef_gz : coef frottement zones stream std
! coef_ile : coef frottement zones iles

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

Hcoup_plateau   =      200     !   tres petit quand shelves fixes sinon 250
Hcoup_abysses   =       200     !   tres petit quand shelves fixes sinon 250
prof_plateau   =       50     !   tres petit quand shelves fixes sinon 250
prof_abysses   =       1000     !   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    = 1.  ! 3.52 ! 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 = 2
/
! frot_coef  = ',frot_coef
! gr_select  = ',gr_select
! gr_select = 1 : Tsai , 2 : Schoof'
!

!___________________________________________________________ 
&snap_forcage_mois                            ! module climat_forcage_mois_mod'

filtr_t1        = 'PnGrmin360_SE_2431_2440_1M_histmthij-t2m.dat'    
filtr_p1        = 'PnGrmin360_SE_2431_2440_1M_histmthij-precip.dat'
/
! fichiers temperature et precip : 12 mois

!___________________________________________________________ 
&lapse_month                                  !  module climat_forcage_mois_mod'

lect_lapse      = 6.0, 6.0, 6.0, 6.0, 6.0, 6.0, 6.0, 6.0, 6.0, 6.0, 6.0, 6.0  ! 12 mois
/
! laspe rates janvier -> decembre en deg/km

!___________________________________________________________ 
&ablation                                     !  module ablation_mod'

pdd_type        = 1 
annual          = .false.
Cice            = 0.008
Csnow           = 0.003
csi             = 0.6
sigma_ice       = 5.0
/
! 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-800k-2007-Jouzel-Bassinot.dat' !'LGM_permanent.dat'
/
!___________________________________________________________ 
&clim_pert_massb                                   ! nom du bloc

coefT           =       0.
filforc         =     'forcage-800k-2007-Jouzel-Bassinot.dat' !'LGM_permanent.dat'
pertsmb         =     0  !1
rapsmb          =     0.07
/
! 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 

!----------------------------------------------------------
! 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_green_gen
&topo_groen_gen
topo_surf      = 'Surface_mean_5km_Tamsin_ieq-param.grd'           !  surface
correc_surf    = 'no'                                              ! ice-real correction
topo_bed       = 'Bedrock_mean_5km_Tamsin.grd'                     ! bedrock
mask_grounded  = 'Mask_5km_checked_Tamsin-Fjords.grd'              ! mask
longitude      = 'long_ZBL_CISM_5km.grd'                           ! longitude ! a garder en .dat
latitude       = 'lat_ZBL_CISM_5km.grd'                            ! latitude  ! car donne xmin ymin 
heatflux       = 'ghf-cism.grd'                                    ! geothermal heat flux
/
! remarques : l'epaisseur n'est pas lue mais calculee a partir du masque
! 
!----------------------------------------------------------
!  module lect_climref_Ice2sea    
&clim_smb_T_gen
smb_file         = 'SMB-MAR-ant35km-initMIP-16km.nc' ! smb m/an (en kg/m2/an)
coef_smb_unit    = 1.    ! 1.0989e-3          ! 1.0989e-3  pour1/910  1.0893e-3 pour 1/918 (valeur habituelle) kg/m2/an -> m(glace)/an
temp_annual_file = 'Tann-MAR-ant35km-initMIP-16km.nc'  ! annual surface temperature
/
!----------------------------------------------------------
&smb_anom_initMIP
file_smb_anom     = 'smb_anomaly_16km.nc'
massb_time        = 0     ! 0=fixe, 1:anomalies
/
!----------------------------------------------------------
! module climat_Grice2sea_mod
&lapse_rates
T_lapse_rate     =  0. !-6.309e-3    ! celui pour la T annuelle de Fausto
/
!----------------------------------------------------------
&clim_snap
nb_snap           = 200
time_depart_snaps = 2000000
ecart_snap        = 1.
file_smb_snap     = 'smb_ECHAM5_A1B_MAR_highres_all_time.nc'
massb_time        = 0     ! 0=fixe, 1:interpolation snapshots
/
!----------------------------------------------------------
&grad_smb                        !routine init_grad_smb 
file_grad_smb     = 'GRISLI_gradients_MAR_00.txt'
/
!---------------------------------------------------------
&bmelt_seuil                     ! module  bmelt_seuil_prof
bm_grz           = 3.      !3
bmshelf_plateau  = 1.        !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
!___________________________________________________________ 

&bmelt_ant_reg                    !module bmelt-ant-regions_mod
bmelt_Ross      =       0.21 !0.10 !0.25 !0.8 !0.5 !0.20
bmgrz_Ross      =       0.18 !0.80 !0.15 !1.0 !0.6 !0.25
bmelt_FRis      =       0.54 !0.50 !0.5 !0.2 !0.5 !0.20 !1.
bmgrz_FRis      =       0.27 !1.0  !0.65 !0.4 !0.6 !0.25 !1.3
bmelt_Amery     =       0.49 !0.8 !1.
bmgrz_Amery     =       1.97 !10. !12.
bmelt_PIG       =       2.31 !15. !0.7 !0.2 !1.
bmgrz_PIG       =       1.79 !15. !0.4 !2.
bmelt_Pen       =       0.13 !0.6 !0.6 !1.
bmgrz_Pen       =       1.08 !1.0 !0.5 !2.
bmelt_other     =       0.50 !10. !0.6 !1.
bmgrz_other     =       2.02 !15. !0.4 !2.
bmelt_talus     =       5.
bmgrz_talus     =       5. 
bmelt_coef      =       1. 
file_number_shelves = "numer-ice-shelves-juil07-Ant16km.nc"
/
!___________________________________________________________ 

&bmelt_ant_reg_initmip                    !module bmelt-ant-regions_mod
bmelt_regions      =       1, 1, 1, 1, 1,  1, 1, 1, 1, 1,  0.2, 0.2, 0.2, 0.2, 0.2,  0.2, 0.2, 0.2
bmgrz_regions      =       2, 2, 2, 2, 2,  2, 2, 2, 2, 2,  0.5, 0.5, 0.5, 0.5, 0.5,  0.5, 0.5, 0.5
bmelt_talus     =       5.
bmgrz_talus     =       5. 
bmelt_coef      =       1. 
file_number_shelves = "shelves-from-initMIP.grd"
/
!----------------------------------------------------------
! module bmelt-ant-regions_mod
&bmelt_anom_initMIP
file_bmelt_anom     = 'basal_melt_anomaly_16km.nc'
bmelt_time        = 0     ! 0=fixe, 1:anomalies
/
!___________________________________________________________
&beta_prescr                                   ! dragging_prescr_beta'
beta_c_file = 'beta_AI2Sn138-16km.grd' !'beta-run-inbeta01.grd'         ! beta cat Ant interpole sur 16km
beta_limgz  = .5e6                             ! -1.e12 pour plastic ?    ! .5e6 en visqueux   
beta_min    = 1 !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_iter_vitbil                              ! beta_iter_vitbil

time_iter      = 5                             ! temps de demarrage des iterations
nb_iter_vitbil = 2                             ! nombre d'iterations par pas de temps (dtt)
coef_iter_vitbil = 1.                          ! coefficient pour rapport des vitesses <=1
Umag_bil_file  = 'vitbil-run-inbeta01.grd'  ! fichier des vitesses de bilan centrees
/
!
!___________________________________________________________
&drag_param_beta                              ! drag_beta_param_mod

beta_slope    = 1.e-3    !1.e-3
beta_expo     = 1.     !1.0
betamax    = 5.e5     ! maximum value of beta
betamin    = 1 !10.   ! minimum value of beta
coef_ile   = 0.1
/
!
!___________________________________________________________
&spinup                                ! warning : 2 different modules
ispinup      = 0
!type_vitbil  = 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'
! type_vitbil type de lecture des vitesses de bilan 1: lect_vitbil, 2 : lect_vitbil_Lebrocq
!___________________________________________________________
&vitbil_upwind             ! nom du bloc vitbil calcule sur le mailles staggered
balance_Ux_file = 'Uy_moy_stag_5km.nc'
balance_Uy_file = 'Ux_moy_stag_5km.nc'
/
!balance_Ux_file = 'Ux_moy_stag_5km.nc'
!balance_Uy_file = 'Uy_moy_stag_5km.nc'
!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 
!___________________________________________________________
&vitbil_Lebrocq             ! nom du bloc vitbil obtenu par code Lebrocq
vit_balance_file = 'gwavel.mar.War10.grd'
flowdir_balance_file = 'gwaflowdir.mar.War10.grd'
/
! vit_balance_file : balance velocity file
! flowdir_balance_file : balance flow direction file
!___________________________________________________________
&output_regions
nbregions       = 7
region_file     = 'Basins_Greenland_aurel.grd'
/
! nbregions   = nombre de regions
! region_file = carte des regions
!_______________________________________________________________________________