source: branches/GRISLIv3/launching_area/Master-files/Param/anteis1_param_list_lhs-nolin.dat @ 489

Last change on this file since 489 was 489, checked in by aquiquet, 4 months ago

Cleaning branch: adding necessary files to run LHS simulations / module choix and param has to be adapted though
File size: 13.0 KB
Line 
1!  Parametres du run : Ant-40
2
3!___________________________________________________________
4&runpar                       ! nom du bloc parametres du run
5 runname      =  "RUNNAME"   ! 8 caracteres
6 icompteur    =    2
7 iout         =    2         
8 reprcptr     =  "../../../Fichier-CPTR/AN40T039+k100.nc"  ! AN40K048+k100.nc !40temp03+k60.nc
9 itracebug    = 0
10 num_tracebug = 163     
11 comment_run  = "Equilibre avec BM RACMO"
12! runname    : nom de l experience (8 caracteres)
13! icompteur  : reprise dans un fichier  0 -> non, 1 -> oui, 2 -> T et Hwat
14!              3-> T seulement
15! iout       : 1-> sortie cptr pour reprise,2 -> sortie nc pour reprise
16! reprcptr   : nom du fichier
17!___________________________________________________________
18&grdline                 ! bloc grounding line
19
20 igrdline     =            0
21 Schoof       =            1
22 ibmelt_inv   =            0
23 /
24! igrdline :  1 ligne d echouage fixée, 3 recul forcé, 0 sans traitement specifique, 2 paleo gr line
25! Schoof   :  0 pas de Schoof, 1 flux de Schoof
26
27!___________________________________________________________
28&timesteps                ! bloc timestep
29
30 tend      =   RUNTEND
31 tbegin    =   RUNTBEGIN         !1.e10      ! si tbegin > 1.e9 on prend le temps du fichier cptr
32 dtmin     =   2.e-3
33 dtmax     =   1. 
34 dtt       =   5. 
35 testdiag  =   0.016 !0.025      !0.016     
36 /
37! tous les temps en annees. tbegin et tend : debut et fin du run
38! pour equation masse, pas de temps mini -> dtmin, maxi -> dtmax
39! dtt : pas de temps long
40! testdiag, pour gerer le pas de temps dynamique dt
41! ordres de grandeur (a moduler selon dx) :
42! 40 km dtmin=2.e-3, dtmax=1., dtt=5., tesdiag=0.02
43
44!___________________________________________________________
45&topo_file
46 topo_ref     = "BedMachine-Antv3_Ant40.nc" !"SHB_bedmap2_Ant40_relax20ans.nc" !"SHB_bedmap2_Ant40.nc" !"SHB_Ant_Lebrocq_40km.nc"
47 topo_dep     = "BedMachine-Antv3_Ant40.nc" !"SHB_bedmap2_Ant40_relax20ans.nc" !"SHB_bedmap2_Ant40.nc" !"SHB_Ant_Lebrocq_40km.nc"
48 grid_topo    = "coord-Ant-40km.dat"
49 ghf_fich     = "ghf_Haeger-et-al_2022_ant40.nc"
50/
51! topo_ref= topo ref isostasie
52! topo_dep= topo de depart
53! grid_topo : fichier i,j,x,y,lon,lat
54! ghf_fich  : fichier flux geothermique
55
56!___________________________________________________________
57&mass_conserv            ! conservation de la masse
58adv_frac        =  2.    ! 2-> advection seule
59V_limit         = 3000.  ! maximum velocity for mass conservation : depend de la calotte
60/
61! Conservation de la masse avec equation advection-diffusion '
62! la repartition depend de adv_frac'
63! >1  -> advection seule'
64! 0   -> diffusion seule'
65! 0<*<1   -> fraction de l advection'
66! -1 -> zones diffusion + zones advecttion'
67! V_limit depend de la calotte : typiquement 3000 en Antarctique, 10000 au Groenland
68
69!--------------------------------------------------------------------
70&eaubasale1              ! nom du premier bloc eau basale
71 
72 ecoulement_eau =  T
73 hwatermax      =    5000.000   
74 infiltr        =   PARAMINFILTR   !1.0000001E-03
75 /
76! ecoulement eau : .false. -> modele bucket, sinon equ. diffusion
77! hwatermax :  hauteur d eau basale maximum dans le sediment (m)
78! infiltr est la quantite d eau qui peut s infiltrer dans le sol (m/an)
79
80!___________________________________________________________
81&param_hydr             ! nom du  bloc parametres hydrauliques
82 
83 hmax_till      =    20.00000   
84 poro_till      =   0.5000000   
85 kond0          =   PARAMKOND0   !1.000000E-06
86 
87! hmax_till (m) : epaisseur max du sediment
88! poro_till : porosite du sediment
89! conductivite du sediment :  kond0 (m/s)
90
91!____________________________________________________________
92&calving                        !   nom du bloc calving méthode Vincent
93
94Hcoup_plateau   =      250     !   tres petit quand shelves fixes sinon 250
95Hcoup_abysses   =       250     !   tres petit quand shelves fixes sinon 250
96prof_plateau   =       50     !   tres petit quand shelves fixes sinon 250
97prof_abysses   =       1000     !   tres petit quand shelves fixes sinon 250
98ifrange         =       4
99meth_Hcoup      =       0
100/
101
102! Hcoup epaisseurs de coupure pour les zones peu prodondes et profondes
103! Hcoup_plateau<Hcoup_abysses && prof_plateau<prof_abysse
104! prof profondeur delimitant les zones peu prodondes et profondes
105! ifrange=0 -> pas de traitement particulier sur les bords
106! ifrange=1 -> traitement de Vincent avec ice shelves frangeants
107! ifrange=2 -> ice shelves frangeant seulement si bm-bmelt positif
108! meth_hcoup pour faire eventuellement varier Hcoup avec le climat
109!____________________________________________________________
110! loi de deformation 1          module deformation_mod_2lois
111&loidef_1
112
113exposant_1       = 3.
114temp_trans_1     = -6.5
115enhanc_fact_1    = PARAMENHFACT  ! 3.52 ! Grice2sea 3.         
116coef_cold_1      = 1.660E-16
117Q_cold_1         = 7.820E+04
118coef_warm_1      = 2.000E-16
119Q_warm_1         = 9.545E+04
120/
121! exposant (glen), temperature de transition (ttrans)
122! enhancement factor (sf)
123! pour les temperatures inf. a Temp_trans :
124!            coef_cold (Bat1) et Q_cold (Q1)
125! pour les temperatures sup. a Temp_trans :
126!            coef_warm (Bat2) et Q_warm (Q2)
127!________________________________________________________
128! loi de deformation 2          module deformation_mod_2lois
129&loidef_2
130
131exposant_2       = 1.
132temp_trans_2     = -10.
133enhanc_fact_2    = 1.            ! Grice2sea 3.         
134coef_cold_2      = 8.313E-08
135Q_cold_2         = 4.000E+04
136coef_warm_2      = 8.313E-08
137Q_warm_2         = 6.000E+04
138/
139! exposant (glen), temperature de transition (ttrans)
140! enhancement factor (sf)
141! pour les temperatures inf. a Temp_trans :
142!            coef_cold (Bat1) et Q_cold (Q1)
143! pour les temperatures sup. a Temp_trans :
144!            coef_warm (Bat2) et Q_warm (Q2)
145!___________________________________________________________
146&diagno_rheol              ! nom du bloc  diagno_rheol
147
148sf01           =        0.125
149sf03           =        0.125
150sffunc         =        1
151pvimin         =        1.5e3   ! 1.e3
152/ 
153! coefficients par rapport a la loi glace posee                       
154! sf01 : coefficient viscosite loi lineaire
155! sf03 : coefficient viscosite loi n=3
156! sffunc : 1 for visco coef that depends on enh factor
157! pvimin : valeur de pvi pour les noeuds fictifs ~ 1.e3
158! tres petit par rapport aux valeurs standards ~ 1.e10
159!___________________________________________________________
160
161! furst_schoof
162&furst_schoof
163
164frot_coef = 0.6   !0.035
165gr_select = 1
166/
167! frot_coef  = ',frot_coef
168! gr_select  = ',gr_select
169! gr_select = 1 : Tsai , 2 : Schoof'
170!___________________________________________________________
171&ablation                                     !  module ablation_mod'
172
173pdd_type        = 1
174annual          = .false.
175Cice            = 0.008
176Csnow           = 0.003
177csi             = 0.6
178sigma_ice       = 5.0
179/
180! pdd_type : 0 reeh, 1 Fausto, 2 Tarasov
181! annual : T = annuel, F = mensuel
182! Cice and Csnow, melting factors for ice and snow
183! sigma variabilite Tday
184! csi proportion of melted water that can refreeze
185
186!___________________________________________________________
187&clim_pert_massb                                   ! nom du bloc
188
189coefT           =       1.
190filforc         =     'forcage-429k-2007-Jouzel-Waelbroeck_AMOC.dat' !'LGM_permanent.dat'
191pertsmb         =     0  !1
192rapsmb          =     0.07
193pertbmb         =     0
194coefbmb   =     10.
195/
196! pertsmb = 1, smb pert via tpert / pertsmb = 0, smb = cst
197! pertbmb = 1, bmelt pert via derniere colonne du forcage, avec coef rapbmb
198!----------------------------------------------------------
199&meca_SIA_L1                        ! bloc resol_meca
200
201i_resolmeca     =        2
202/                           
203! i_resolmeca type d association entre SIA et L1'
204! i_resolmeca=0  chacun dans sa zone'
205! i_resolmeca=1  dans les zones stream, addition si uxdef > uxL1 (MIS11 Cairns)
206! i_resolmeca=2  addition systematique dans les zones stream
207
208!----------------------------------------------------------
209!  module lect_climref_Ice2sea   
210&clim_smb_T_gen
211smb_file         = 'RACMO2.3p2_ANT27_smb_ltm_1979_2016_ant40.grd' !'SMB-MAR-ant35km-grisli-40km.grd' !'acc_Ant_Lebrocq_40km.nc' ! smb m/an (en kg/m2/an)
212coef_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
213temp_annual_file = 'RACMO2.3p2_ANT27_t2m_ltm_1979_2016_ant40.grd' !'Tann-MAR-ant35km-grisli-40km.grd' !'TANN_Ant_Lebrocq_40km.nc'  ! annual surface temperature
214/
215!----------------------------------------------------------
216! module climat_Grice2sea_mod
217&lapse_rates
218T_lapse_rate     =  0.
219/
220!---------------------------------------------------------
221! module bmelt_beckmann_gcm
222&bmelt_beckmann_gcm_mod
223K_t = PARAMKT
224TO_file ='TO-SO_clim-from-obs_1995-2017_ant40.nc'
225SO_file ='TO-SO_clim-from-obs_1995-2017_ant40.nc'
226Bassin_file='shelves-from-initMIP-40km-forBeckmann.nc'
227nbassins     = 18
228deltaT_bool = 1
229deltaT_file = 'coeff_gamma0_DeltaT_quadratic_non_local_median_ant40.nc'
230bmelt_empty = 10
231nonlocal_bool = 1
232grlmelt_bool = 0
233/
234!---------------------------------------------------------
235&bmelt_ismip6_param                    !module bmel-ismip6-param_mod
236file_TF            = "obs_thermal_forcing_1995-2017_8km_x_60m_ant40.nc"
237file_basinNumbers  = "imbie2_basin_numbers_8km_ant40.nc"
238file_coef          = "coeff_gamma0_DeltaT_quadratic_non_local_median_ant40.nc"
239gamma0             = 14477.3367602277
240! file_TF: 3d thermal forcing
241! file_basinNumbers: identifier for the Imbie2 basins
242! file_coef: deltaT
243! gamma0: value associated with the deltaT file
244/
245!___________________________________________________________
246! module bmelt_ismip6_param
247&bmelt_anom_initMIP
248file_bmelt_anom     = 'basal_melt_anomaly_16km.nc' ! 'basal_melt_anomaly_16km.nc'
249nb_snap             = 106
250time_depart_snaps   = 0.
251bmelt_time          = 0     ! 0=fixe, 1:anomalies 2:interp snapshots
252/
253!___________________________________________________________
254&beta_prescr                                   ! dragging_prescr_beta'
255beta_c_file = 'beta-run-beta-9-4.grd' !beta-run-beta10-5.grd' !'beta_ef-01-04.grd' !'beta-GISCAM23_20ans.grd'        ! 'beta_Fab_10avril12_bord.grd'    ! 'beta_fab_comb_alpha_v4.grd'
256beta_limgz  = .5e6                             ! -1.e12 pour plastic ?    ! .5e6 en visqueux   
257beta_min    = 1. !10.                           ! for grounded ice (Pa)
258beta_mult   = 1.                               ! coefficient multiplicateur   
259/
260! read beta on centered grid
261! beta_file : nom des fichiers qui contiennent les betamx et betamy
262! beta-estime-run07.dat calcul direct beta sans modif
263! beta-staggered-15km-opt1.dat premier test en utilisant beta calcule  avec iterations
264! above beta_limgz, gzmx is false if negative : no sliding
265! betamx_file = not read anymore
266! betamy_file = not read anymore
267!___________________________________________________________
268&beta_iter_vitbil                              ! beta_iter_vitbil
269
270time_iter      = 5                             ! temps de demarrage des iterations
271nb_iter_vitbil = 2                             ! nombre d'iterations par pas de temps (dtt)
272coef_iter_vitbil = 1.                          ! coefficient pour rapport des vitesses <=1
273Umag_bil_file  = 'ebauche_vitbil_GrIce2sea.grd'  ! fichier des vitesses de bilan centrees
274/
275!
276!___________________________________________________________
277&drag_coulomb_friction                        ! drag_coulomb_friction
278
279cf            = PARAMBSLOPE    !1.e-3
280m_nolin       = PARAMMNOLIN
281niter_nolin   = PARAMNBITERNOLIN
282betamax    = 5.e5     ! maximum value of beta
283betamin    = 1.   ! minimum value of beta
284bool_sedim = .F.
285file_sedim  = 'Ant40_hsedim.nc'
286seuil_sedim = 200.
287coef_sedim  = 0.5e-1
288/
289! cf: a friction coefficient (to be tuned)
290! m_nolin: non-linear exponent, from 1 to infinity (put -1 for infinity)
291!          m_nolin=1/q in Pattyn TC 2017, q in [0:1]
292! niter_nolin: number of iterations to solve the non-linearity (expensive!)
293!
294!___________________________________________________________
295&drag_coulomb_friction_simplhydro           ! drag_coulomb_friction_simplhydro
296
297till_c0       = 0.
298TFA_min     = 5.
299TFA_max     = 30.
300bed_crit     = -1000.
301Ob_coeff      = PARAMOB !1. !0.15843    !1.e-3
302m_nolin       = PARAMMNOLIN
303niter_nolin   = PARAMNBITERNOLIN
304betamax    = 5.e5     ! maximum value of beta
305betamin    = 1.   ! minimum value of beta
306bool_sedim = .F.
307file_sedim  = 'Ant40_hsedim.nc'
308seuil_sedim = 200.
309coef_sedim  = 0.5e-1
310/
311! till_c0: sediment cohesion (usually = 0)
312! TFA_coeff: till friction angle
313! Ob_coeff: a friction coefficient
314! m_nolin: non-linear exponent, from 1 to infinity (put -1 for infinity)
315!          m_nolin=1/q in Pattyn TC 2017, q in [0:1]
316! niter_nolin: number of iterations to solve the non-linearity (expensive!)
317!
318!___________________________________________________________
319&spinup                                ! warning : 2 different modules
320ispinup = 0
321/
322! with module no_spinup
323! ispinup = 0     run standard ou calcul du beta
324! ispinup = 1     temperature equilibrium with grisli velocities'
325!
326! with module spinup_vitbil
327! ispinup = 2     conservation de la masse avec vitesses bilan '
328! ispinup = 3     equilibre temperature avec vitesses bilan'
329!
330!___________________________________________________________
Note: See TracBrowser for help on using the repository browser.