source: trunk/lauching_area/degla_PMIP_euras40_param_list.dat @ 334

Last change on this file since 334 was 313, checked in by dumas, 4 years ago

Add degla_PMIP param_list files for euras40, euras10 and hemin40

File size: 15.0 KB
Line 
1!  Parametres du run : geoplace GI2S_5T Groenland Ice2sea 5 km cut Tamsin
2
3!___________________________________________________________
4&runpar                       ! nom du bloc parametres du run
5 runname      =  "RUNNAME"   ! 8 caracteres
6 icompteur    =    0
7 iout         =    2 
8 reprcptr     =  "../../../Fichier-CPTR/LOV21aT0+k50.nc"
9 itracebug    = 0
10 num_tracebug = 166     
11 comment_run  = "Eurasie Degla IPSL clim_forcage"
12/
13! runname    : nom de l experience (8 caracteres)
14! icompteur  : reprise dans un fichier  0 -> non, 1 -> oui, 2 -> T et Hwat
15!              3-> T seulement
16! iout       : 1-> sortie cptr pour reprise,2 -> sortie nc pour reprise
17! reprcptr   : nom du fichier
18!___________________________________________________________
19&grdline                 ! bloc grounding line
20
21 igrdline     =            0
22 Schoof       =            1
23 ibmelt_inv   =            0
24 /
25! igrdline :  1 ligne d echouage fixée, 3 recul forcé, 0 sans traitement specifique, 2 paleo gr line
26! Schoof   :  0 pas de Schoof, 1 flux de Schoof
27
28!___________________________________________________________
29&timesteps                ! bloc timestep
30
31 tend      =   RUNTEND
32 tbegin    =   RUNTBEGIN
33 dtmin     =   2.e-3      ! 2.e-3 40km
34 dtmax     =   1.         ! 1.    40km
35 dtt         =   1.         ! 1.    40km
36 testdiag  =   0.025      !0.025  40km      !0.016     
37 /
38! tous les temps en annees. tbegin et tend : debut et fin du run
39! pour equation masse, pas de temps mini -> dtmin, maxi -> dtmax
40! dtt : pas de temps long
41! testdiag, pour gerer le pas de temps dynamique dt
42! ordres de grandeur (a moduler selon dx) :
43! 40 km dtmin=2.e-3, dtmax=1., dtt=5., tesdiag=0.02
44
45!___________________________________________________________
46&topo_file
47 topo_ref        = "SHB_etopo1_euras40.nc"
48 topo_dep        = "SHB_etopo1_euras40.nc"
49 grid_topo       = "coord-euras40.dat"
50 ghf_fich        = "ghf_euras40.nc"
51/
52! topo_ref= topo ref isostasie
53! topo_dep= topo de depart
54! grid_topo : fichier i,j,x,y,lon,lat
55! ghf_fich  : fichier flux geothermique
56!___________________________________________________________
57&mass_conserv            ! conservation de la masse
58adv_frac        =  2.    ! 2-> advection seule
59V_limit         = 5000.  ! 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        =   1.e-3   !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          =   0.000193342957286257   !10.000000E-05
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. !50     !   tres petit quand shelves fixes sinon 250
95Hcoup_abysses   =       250. !400     !   tres petit quand shelves fixes sinon 250
96prof_plateau   =       50. !100     !   tres petit quand shelves fixes sinon 250
97prof_abysses   =       1000. !800     !   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    = 1.82595387172575          ! 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 
147&diagno_rheol              ! nom du bloc  diagno_rheol
148
149sf01           =        0.125   ! 0.125
150sf03           =        0.125   ! 0.125
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! pvimin : valeur de pvi pour les noeuds fictifs ~ 1.e3
157! tres petit par rapport aux valeurs standards ~ 1.e10
158!___________________________________________________________
159
160! furst_schoof
161&furst_schoof
162
163frot_coef = 0.6   !0.035
164gr_select = 1
165/
166! frot_coef  = ',frot_coef
167! gr_select  = ',gr_select
168! gr_select = 1 : Tsai , 2 : Schoof'
169!___________________________________________________________
170&lapse_month                                  !  module climat_forcage_mois_mod'
171
172lect_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
173/
174! laspe rates janvier -> decembre en deg/km
175!___________________________________________________________
176&clim_coupl_atm                                ! nom du bloc
177
178mincoefbmelt  =  0.
179maxcoefbmelt  =  2.
180filforc       = "forcage-429k-quiquet-et-al_CP2013_glacial-index_Waelbroeck_AMOC.dat" !"signal-cycle-hemin.dat" !""40k.dat" !"778k.dat"
181snow          =       .true.
182! mincoefbmelt : borne mini de coefbmshelf
183! maxcoefbmelt : borne maxi de coefbmshelf
184! filforc : fichier de forcage nivx des mers
185! snow   : .true. si neige, .false. si precip
186/
187!___________________________________________________________
188&ablation                                     !  module ablation_mod'
189
190pdd_type        = RUNPDDTYPE
191annual          = .false.
192Cice            = 0.008
193Csnow           = 0.003
194csi             = 0.6
195sigma_ice       = 5.5
196/
197! pdd_type : 0 reeh, 1 Fausto, 2 Tarasov
198! annual : T = annuel, F = mensuel
199! Cice and Csnow, melting factors for ice and snow
200! sigma variabilite Tday
201! csi proportion of melted water that can refreeze
202
203!___________________________________________________________
204&clim_pert                                    ! nom du bloc
205
206coefT           =       1.
207rappact         =       0.07                    ! 0.05->0.75  0.07 -> 0.5
208retroac         =       1
209rapbmshelf      =       5. 
210mincoefbmelt    =       0.
211maxcoefbmelt    =       2.
212filforc         =     'forcage-900k-2007-Parrenin-Bassinot.dat' !'LGM_permanent.dat'
213/
214!___________________________________________________________
215&clim_pert_massb                                   ! nom du bloc
216
217coefT           =     0. ! 1.
218filforc         =   'forcage-429k-2007-Jouzel-Waelbroeck_AMOC.dat' ! 'quiquet-et-al_CP2013_forcage.dat' !'forcage-900k-2007-Parrenin-Bassinot.dat' !'LGM_permanent.dat'
219pertsmb         =     0  !1  !1
220rapsmb          =     0. !0.07 !7
221pertbmb         =     1
222coefbmb   =     10.
223/
224! pertsmb = 1, smb pert via tpert / pertsmb = 0, smb = cst
225!----------------------------------------------------------
226&meca_SIA_L1                        ! bloc resol_meca
227
228i_resolmeca     =        2
229/
230! i_resolmeca type d association entre SIA et L1'
231! i_resolmeca=0  chacun dans sa zone'
232! i_resolmeca=1  dans les zones stream, addition si uxdef > uxL1 (MIS11 Cairns)
233! i_resolmeca=2  addition systematique dans les zones stream
234
235!----------------------------------------------------------
236! pour le module tracer (suivi d'isotopes)
237&tracer
238file_tr_dat     = "../INPUT/tracer.dat"
239file_tr_out     = "../RESULTATS/tracer.out"
240file_tr_dep     = "../RESULTATS/deparchive.out"
241lambdab         = 0.0001
242agemax          = 1000000
243coefT_tra       = 1.
244rappact_tra     = 0.07
245filforc         = "quiquet-et-al_CP2013_forcage.dat" !"forcage-900k-2007-Parrenin-Bassinot.dat"
246pert_type       = 0
247/
248!----------------------------------------------------------
249!  module lect_climref_Ice2sea   
250&clim_smb_T_gen
251smb_file         = 'smb_MARv3.5.2-25km_1958_1989mean_hemin40.nc' !'smb_ERAI_MAR_highres_1989-2008mean_hemin40.nc'   ! smb (en kg/m2/an
252coef_smb_unit    = 1.0893e-3    ! 1.0989e-3          ! 1.0989e-3  pour1/910  1.0893e-3 pour 1/918 (valeur habituelle)
253temp_annual_file ='Tann-cism-hemin40.nc'                            ! annual surface temperature
254/
255!----------------------------------------------------------
256! module climat_Grice2sea_mod
257&lapse_rates
258T_lapse_rate     =  0. !-6.309e-3    ! celui pour la T annuelle de Fausto
259/
260!----------------------------------------------------------
261! module climat_forcage
262&clim_forcage
263clim_ref_file = 'climato_monthly_t2m_precip_orog_ERA5-1979-2008-euras40.nc'
264ttr = -21000, 0.
265forcage_file1 = 'IPSL-CM5A-LR-lgm-euras40.nc' !'MIROC-ESM-lgm-hemin40.nc' !'MRI-CGCM3-lgm-hemin40.nc' !'MPI-ESM-P-lgm-hemin40.nc' !'FGOALS-g2-lgm-hemin40.nc' !'IPSL-CM5A-LR-lgm-hemin40.nc'
266forcage_file2 = 'IPSL-CM5A-LR-piControl-euras40.nc' !'MIROC-ESM-piControl-hemin40.nc' !'MRI-CGCM3-piControl-hemin40.nc' !'MPI-ESM-P-piControl-hemin40.nc' !'FGOALS-g2-piControl-hemin40.nc' !'IPSL-CM5A-LR-piControl-hemin40.nc'
267typerun = 0
268lapserate = 0.007
269rappact = 0.11
270filforc = 'forcage-429k-quiquet-et-al_CP2013_glacial-index_Waelbroeck_AMOC.dat' !'forcage-429k-2007-Jouzel-Waelbroeck_AMOC.dat'
271pertbmb = 1
272coefT   = 1.
273coefbmb = 10.
274r_atmvar = RUNRATMVAR
275/
276!----------------------------------------------------------
277&clim_snap
278nb_snap           = 200
279time_depart_snaps = 2000000
280ecart_snap        = 1.
281file_smb_snap     = 'smb_ECHAM5_A1B_MAR_highres_all_time.nc'
282massb_time        = 0     ! 0=fixe, 1:interpolation snapshots
283/
284!----------------------------------------------------------
285&grad_smb                        !routine init_grad_smb
286file_grad_smb     = 'GRISLI_gradients_MAR_00.txt'
287/
288!---------------------------------------------------------
289&imposed_regions_rsl               ! module  imposed_regions_rsl
290nb_reg           = 1
291filreg           = Hemin40_Hudson.grd
292filforc          = rsl_hudson_forc.txt
293/
294!
295!___________________________________________________________
296&bmelt_seuil                     ! module  bmelt_seuil_prof
297bm_grz           = 0.3      !3
298bmshelf_plateau  = 0.2        !0.4
299bmshelf_abysses  = 10.
300depth_talus      = -1000.
301/
302! Pour l actuel : bm_grz a la grounding line
303!                 bmshelf_plateau sur le plateau continental
304!                 bmshelf_abysses pour les grandes profondeurs
305!       depth_talus, negative, separation entre les 2 domaines
306!___________________________________________________________
307&beta_prescr                                   ! dragging_prescr_beta'
308beta_c_file = 'beta-GISCAM23_20ans.grd'        ! 'beta_Fab_10avril12_bord.grd'    ! 'beta_fab_comb_alpha_v4.grd'
309beta_limgz  = .5e6                             ! -1.e12 pour plastic ?    ! .5e6 en visqueux   
310beta_min    = 1 !10.                           ! for grounded ice (Pa)
311beta_mult   = 1.                               ! coefficient multiplicateur   
312/
313! read beta on centered grid
314! beta_file : nom des fichiers qui contiennent les betamx et betamy
315! beta-estime-run07.dat calcul direct beta sans modif
316! beta-staggered-15km-opt1.dat premier test en utilisant beta calcule  avec iterations
317! above beta_limgz, gzmx is false if negative : no sliding
318! betamx_file = not read anymore
319! betamy_file = not read anymore
320!___________________________________________________________
321&beta_iter_vitbil                              ! beta_iter_vitbil
322
323time_iter      = 5                             ! temps de demarrage des iterations
324nb_iter_vitbil = 2                             ! nombre d'iterations par pas de temps (dtt)
325coef_iter_vitbil = 1.                          ! coefficient pour rapport des vitesses <=1
326Umag_bil_file  = 'ebauche_vitbil_GrIce2sea.grd'  ! fichier des vitesses de bilan centrees
327time_reiter    = 200.                          ! temps entre 2 iterations
328/
329!
330!___________________________________________________________
331&drag_param_beta                              ! drag_beta_param_mod
332
333beta_slope     = 0.001584297807517     !1.e-3
334beta_expo      = 1.      !1.0
335betamax    = 5.e5     ! maximum value of beta
336betamin    = 100. !10.   ! minimum value of beta
337coef_ile   = 1.
338/
339!
340!___________________________________________________________
341&drag_param_beta_sedim                       ! drag_beta_param_sedim_mod
342
343beta_slope    = 0.001584297807517    !1.e-3
344beta_expo     = 1.     !1.0
345betamax    = 5.e5     ! maximum value of beta
346betamin    = 100. !10.   ! minimum value of beta
347coef_ile   = 1
348file_sedim  = 'sedmap_euras40.nc'
349seuil_sedim = 200.
350coef_sedim  = 1.e-1
351/
352!
353!___________________________________________________________
354&spinup                                ! warning : 2 different modules
355ispinup = 0
356/
357! with module no_spinup
358! ispinup = 0     run standard ou calcul du beta
359! ispinup = 1     temperature equilibrium with grisli velocities'
360!
361! with module spinup_vitbil
362! ispinup = 2     conservation de la masse avec vitesses bilan '
363! ispinup = 3     equilibre temperature avec vitesses bilan'
364!
365!___________________________________________________________
366&vitbil_upwind             ! nom du bloc vitbil calcule sur le mailles staggered
367balance_Ux_file = 'Uy_moy_stag_5km.nc'
368balance_Uy_file = 'Ux_moy_stag_5km.nc'
369/
370!balance_Ux_file = 'Ux_moy_stag_5km.nc'
371!balance_Uy_file = 'Uy_moy_stag_5km.nc'
372!balance_vel_file = 'Umoy-for-calc-beta_LBq15-06_ZBL_15km.dat'
373! balance_vel_file = 'Uslid-for-calc-beta_LBq15-06_ZBL_15km.dat'
374! balance_vel_file = 'Umoy-for-calc-beta_LBq15-06_ZBL_15km.dat'
375! balance_vel_file = 'Vitbil-stag-pour-spinup_merge_runs_Lbq_ZBL_15km.dat'
376! balance velocities on staggered grid Ux, Uy
377!___________________________________________________________
378&output_regions
379nbregions       = 7
380region_file     = 'Basins_Greenland_aurel.grd'
381/
382! nbregions   = nombre de regions
383! region_file = carte des regions
384!_______________________________________________________________________________
Note: See TracBrowser for help on using the repository browser.