source: trunk/SOURCES/Fichiers-parametres/A-LBq15_param_list_Schoof.dat @ 242

Last change on this file since 242 was 93, checked in by dumas, 8 years ago

First version with Schoof flux parameterisation at the grounding line. | New module furst_schoof_mod.f90 | New flag Schoof in grdline namelist (see in SOURCES/Fichiers-parametres/A-LBq15_param_list_Schoof.dat)
File size: 16.8 KB
Line 
1!  Parametres du run : Ant-15
2
3!___________________________________________________________
4&runpar                       ! nom du bloc parametres du run
5 runname      =  "Schoof04"   ! 8 caracteres
6 icompteur    =    1   
7 iout         =    2         
8 reprcptr     =  "../Schoof01/Schoof01_100ans.nc" !  "../../Fichier-CPTR/RGL-0194_an2040.nc"  !restart Cat
9 itracebug    = 0
10 num_tracebug = 163     
11 comment_run  = "Equilibre avec BM Lebrocq Athern et al"
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     =            1
21 Schoof       =            1
22 /
23! igrdline :  1 ligne d echouage fixée, 3 recul forcé, 0 sans traitement specifique
24! Schoof   :  0 pas de Schoof, 1 flux de Schoof
25! paleo grounding line -> 2
26!___________________________________________________________
27&timesteps                ! bloc timestep
28
29 tend      =   100.
30 tbegin    =   0.         !1.e10      ! si tbegin > 1.e9 on prend le temps du fichier cptr
31 dtmin     =   1.e-3
32 dtmax     =   1. 
33 dtt       =     0.5 
34 testdiag  =   1.e-2      !0.016     
35 /
36! tous les temps en annees. tbegin et tend : debut et fin du run
37! pour equation masse, pas de temps mini -> dtmin, maxi -> dtmax
38! dtt : pas de temps long
39! testdiag, pour gerer le pas de temps dynamique dt
40! ordres de grandeur (a moduler selon dx) :
41! 40 km dtmin=2.e-3, dtmax=1., dtt=5., tesdiag=0.02
42
43!___________________________________________________________
44&topo_file
45 topo_ref     = "SHB_Ant_Lebrocq_15km.nc"
46 topo_dep     = "SHB_Ant_Lebrocq_15km.nc"
47 grid_topo    = "coord-Ant-15km.dat"
48 ghf_fich     = "ghf_Ant_Lebrocq_15km.nc"
49/
50! topo_ref= topo ref isostasie
51! topo_dep= topo de depart
52! grid_topo : fichier i,j,x,y,lon,lat
53! ghf_fich  : fichier flux geothermique
54
55! here LeBrocq files
56! remark : in the mask, island causing trouble have been removed
57!----------------------------------------------------------
58!  module  lect_topo_ant_gen
59& topo_ant_gen
60topo_surf      = 'usrf-ice-equ-final_15km_i2s.grd'          !'usrf-ice-equ_aout2010_15km.grd'   !  surface
61correc_surf    = 'no'                                       ! ice-real correction
62topo_thick     = 'thick2-ice-equ_15km_i2s.grd'              !'thick2-ice-equ_aout2010_15km.grd' ! thickness
63topo_bed       = 'topg2_15km_i2s.grd'                       ! topg2_aout2010_15km.grd'          ! bedrock
64mask_grounded  = 'new_mask_avril2012_i2s.grd'               ! newmask_aout2010_15km.grd'        ! mask
65longitude      = 'long_LBq_15km.grd'                        ! long_ZBL_LBq_15km.dat'      ! longitude a garder en .dat
66latitude       = 'lat_LBq_15km.grd'                     !lat_ZBL_LBq_15km.dat'        ! latitude
67heatflux       = 'ghfsr_15km.grd'                           ! geothermal heat flux
68/
69
70!___________________________________________________________
71!  module  lect_clim_acc_T_ant_gen                         
72 &CLIMAT_ACC_T_GEN
73 PRECIP_FILE     = acca_15km.grd                                                                                       ,
74 COEF_DENS       =   1.000000    ,
75 TEMP_ANNUAL_FILE        = temp_15km.grd                                                                                       
76 /
77
78!___________________________________________________________
79&clim_pert                                    ! nom du bloc
80
81coefT           =   1.
82rappact               =   0.07                    ! 0.05->0.75  0.07 -> 0.5
83retroac               =         1
84rapbmshelf      =   5. 
85mincoefbmelt      =     0.
86maxcoefbmelt      =     2.
87filforc         =   'forcage-900k-2007-Parrenin-Bassinot.dat' !'LGM_permanent.dat'
88/
89
90!___________________________________________________________
91&mass_conserv            ! conservation de la masse
92adv_frac        =  2.    ! 2-> advection seule
93V_limit         = 3000.  ! maximum velocity for mass conservation : depend de la calotte
94/
95! Conservation de la masse avec equation advection-diffusion '
96! la repartition depend de adv_frac'
97! >1  -> advection seule'
98! 0   -> diffusion seule'
99! 0<*<1   -> fraction de l advection'
100! -1 -> zones diffusion + zones advecttion'
101! V_limit depend de la calotte : typiquement 3000 en Antarctique, 10000 au Groenland
102
103!___________________________________________________________
104!  read parameters of retreat
105 &RETREAT_ICE2SEA
106 REGION_FILE     = mask_regions_15km-0.grd                                                                             ,
107 FILE_EXP        = antarctic_stability_1000_samples-V3.dat                                                             
108 /
109!  file_exp est le fichier qui contient toutes les experiences
110! le numero de l experience correspond au numero du job
111
112!--------------------------------------------------------------------
113&eaubasale1              ! nom du premier bloc eau basale
114 
115 ecoulement_eau =  T
116 hwatermax      =    5000.000   
117 infiltr        =    1.e-3   !1.0000001E-03
118 /
119! ecoulement eau : .false. -> modele bucket, sinon equ. diffusion
120! hwatermax :  hauteur d eau basale maximum dans le sediment (m)
121! infiltr est la quantite d eau qui peut s infiltrer dans le sol (m/an)
122
123!___________________________________________________________
124&param_hydr             ! nom du  bloc parametres hydrauliques
125 
126 hmax_till      =    20.00000   
127 poro_till      =   0.5000000   
128 kond0          =   1.e-6   !1.000000E-06
129 
130! hmax_till (m) : epaisseur max du sediment
131! poro_till : porosite du sediment
132! conductivite du sediment :  kond0 (m/s)
133
134!____________________________________________________________
135&drag_neff_slope              ! nom du bloc dragging neff slope
136
137cf              =       0.000336909217323254  !2.e-5             ! 1.e-4
138betamax         =       1000.
139betamin         =       10.
140toblim          =       0.7e5              ! 0.25e5
141tostick         =       1.e5
142seuil_neff      =       5.e5   !350.e5
143coef_gz         =       10.
144coef_ile        =       0.1
145slope_fich      = "ETOPO1_slope_R2min_Ant15km.nc"
146expo_slope      =       9.83623988446314
147pente_min      =        500.
148pente_max      =        2000.
149/
150! cf coefficient de la loi de frottement fonction Neff'
151! betamax : (Pa) frottement maxi '
152! betamin : (Pa) frottement mini '
153! toblim : (Pa) tes pour les iles '
154! tostick : (Pa)  pour les points non flgzmx '
155! seuil_neff (Pa) : seuil sur la pression effective pour avoir glissement'
156! coef_gz : coef frottement zones stream std
157! coef_ile : coef frottement zones iles
158
159!____________________________________________________________
160&calving                        !   nom du bloc calving méthode Vincent
161 
162Hcoup             =          1.2     !   tres petit quand shelves fixes sinon 250
163ifrange           =          4
164meth_Hcoup      =          0
165/
166! Hcoup epaisseur de coupure
167! ifrange=0 -> pas de traitement particulier sur les bords'
168! ifrange=1 -> traitement de Vincent avec ice shelves frangeants'
169! ifrange=2 -> ice shelves frangeant seulement si bm-bmelt positif
170!____________________________________________________________
171! loi de deformation 1          module deformation_mod_2lois
172&loidef_1
173
174exposant_1       = 3.
175temp_trans_1     = -6.5
176enhanc_fact_1    = 1.8  ! 3.52 ! Grice2sea 3.         
177coef_cold_1      = 1.660E-16
178Q_cold_1         = 7.820E+04
179coef_warm_1      = 2.000E-16
180Q_warm_1         = 9.545E+04
181/
182! exposant (glen), temperature de transition (ttrans)
183! enhancement factor (sf)
184! pour les temperatures inf. a Temp_trans :
185!            coef_cold (Bat1) et Q_cold (Q1)
186! pour les temperatures sup. a Temp_trans :
187!            coef_warm (Bat2) et Q_warm (Q2)
188!________________________________________________________
189! loi de deformation 2          module deformation_mod_2lois
190&loidef_2
191
192exposant_2       = 1.
193temp_trans_2     = -10.
194enhanc_fact_2    = 1.8            ! Grice2sea 3.         
195coef_cold_2      = 8.313E-08
196Q_cold_2         = 4.000E+04
197coef_warm_2      = 8.313E-08
198Q_warm_2         = 6.000E+04
199/
200! exposant (glen), temperature de transition (ttrans)
201! enhancement factor (sf)
202! pour les temperatures inf. a Temp_trans :
203!            coef_cold (Bat1) et Q_cold (Q1)
204! pour les temperatures sup. a Temp_trans :
205!            coef_warm (Bat2) et Q_warm (Q2)
206!___________________________________________________________
207 
208&diagno_rheol              ! nom du bloc  diagno_rheol
209
210sf01           =        0.2083300   ! 0.125
211sf03           =        0.2083300   ! 0.125
212pvimin         =  1.e4   ! 1.e3
213/ 
214! coefficients par rapport a la loi glace posee                       
215! sf01 : coefficient viscosite loi lineaire
216! sf03 : coefficient viscosite loi n=3
217! pvimin : valeur de pvi pour les noeuds fictifs ~ 1.e3
218! tres petit par rapport aux valeurs standards ~ 1.e10
219!___________________________________________________________
220
221! furst_schoof
222&furst_schoof
223
224frot_coef = 0.6   !0.035
225gr_select = 1
226/
227! frot_coef  = ',frot_coef
228! gr_select  = ',gr_select
229! gr_select = 1 : Tsai , 2 : Schoof'
230!___________________________________________________________
231
232! glissement                   module  sliding_Bindschadler
233&slid_bindsh                    ! nom du bloc
234 
235kweert       =  PARAMKWEERT  !5.e-11
236loigliss     =  2
237coefbmax     =  10.
238
239! kweert : coefficent, loigliss le type de loi
240! coefbmax : facteur de normalisation pour influence eau
241
242!___________________________________________________________
243&snap_forcage_mois                            ! module climat_forcage_mois_mod'
244
245filtr_t1        = 'PnGrmin360_SE_2431_2440_1M_histmthij-t2m.dat'   
246filtr_p1        = 'PnGrmin360_SE_2431_2440_1M_histmthij-precip.dat'
247/
248! fichiers temperature et precip : 12 mois
249
250!___________________________________________________________
251&lapse_month                                  !  module climat_forcage_mois_mod'
252
253lect_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
254/
255! laspe rates janvier -> decembre en deg/km
256
257!___________________________________________________________
258&ablation                                     !  module ablation_mod'
259
260pdd_type        = 0
261annual          = .false.
262Cice            = 0.008
263Csnow           = 0.005
264csi             = 0.6
265sigma_ice       = 5.0
266/
267! pdd_type : 0 reeh, 1 Fausto, 2 Tarasov
268! annual : T = annuel, F = mensuel
269! Cice and Csnow, melting factors for ice and snow
270! sigma variabilite Tday
271! csi proportion of melted water that can refreeze
272
273
274!___________________________________________________________
275&clim_pert_massb                                   ! nom du bloc
276
277coefT           =       1.
278filforc         =     'forcage-900k-2007-Parrenin-Bassinot.dat' !'LGM_permanent.dat'
279pertsmb               =     1  !1
280rapsmb                =     0.07
281/
282! pertsmb = 1, smb pert via tpert / pertsmb = 0, smb = cst
283!----------------------------------------------------------
284& meca_SIA_L1                        ! bloc resol_meca
285
286i_resolmeca     =        2
287/                           
288! i_resolmeca type d association entre SIA et L1'
289! i_resolmeca=0  chacun dans sa zone'
290! i_resolmeca=1  dans les zones stream, addition si uxdef > uxL1 (MIS11 Cairns)
291! i_resolmeca=2  addition systematique dans les zones stream
292
293!----------------------------------------------------------
294! PDD base Tann et Tjuly                module ablation_ann   
295
296&ablation_ann
297Cice          = 0.008
298Csnow         = 0.005
299Sigma         = 5.
300csi           = 0.6
301/
302! Cice and Csnow, melting factors for ice and snow
303! sigma variabilite Tday
304! csi proportion of melted water that can refreeze
305
306!----------------------------------------------------------
307!  module  lect_topo_green_gen
308&topo_groen_gen
309topo_surf      = 'Surface_mean_5km_Tamsin_ieq-param.grd'           !  surface
310correc_surf    = 'no'                                              ! ice-real correction
311topo_bed       = 'Bedrock_mean_5km_Tamsin.grd'                     ! bedrock
312mask_grounded  = 'Mask_5km_checked_Tamsin-Fjords.grd'              ! mask
313longitude      = 'long_ZBL_CISM_5km.grd'                           ! longitude ! a garder en .dat
314latitude       = 'lat_ZBL_CISM_5km.grd'                            ! latitude  ! car donne xmin ymin
315heatflux       = 'ghf-cism.grd'                                    ! geothermal heat flux
316/
317! remarques : l'epaisseur n'est pas lue mais calculee a partir du masque
318!
319!----------------------------------------------------------
320!  module lect_climref_Ice2sea   
321&clim_smb_T_gen
322smb_file         = 'acc_Ant_Lebrocq_15km.nc' ! smb m/an (en kg/m2/an)
323coef_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
324temp_annual_file = 'TANN_Ant_Lebrocq_15km.nc'  ! annual surface temperature
325/
326!----------------------------------------------------------
327! module climat_Grice2sea_mod
328&lapse_rates
329T_lapse_rate     =  0. !-6.309e-3    ! celui pour la T annuelle de Fausto
330/
331!----------------------------------------------------------
332&clim_snap
333nb_snap           = 200
334time_depart_snaps = 2000000
335ecart_snap        = 1.
336file_smb_snap     = 'smb_ECHAM5_A1B_MAR_highres_all_time.nc'
337massb_time        = 0     ! 0=fixe, 1:interpolation snapshots
338/
339!----------------------------------------------------------
340&grad_smb                        !routine init_grad_smb
341file_grad_smb     = 'GRISLI_gradients_MAR_00.txt'
342/
343!---------------------------------------------------------
344&bmelt_seuil                     ! module  bmelt_seuil_prof
345bm_grz           = 3.      !3
346bmshelf_plateau  = 0.4        !0.4
347bmshelf_abysses  = 10.
348depth_talus      = -3000.
349/                     
350! Pour l actuel : bm_grz a la grounding line
351!                 bmshelf_plateau sur le plateau continental
352!                 bmshelf_abysses pour les grandes profondeurs
353!       depth_talus, negative, separation entre les 2 domaines
354!___________________________________________________________
355
356&bmelt_ant_reg                    !module bmelt-ant-regions_mod
357bmelt_Ross      =       0.25
358bmgrz_Ross      =       0.15
359bmelt_FRis      =       0.5
360bmgrz_FRis      =       0.65
361bmelt_Amery     =       0.8
362bmgrz_Amery     =       10.
363bmelt_PIG       =       0.8
364bmgrz_PIG       =       0.5
365bmelt_Pen       =       0.7
366bmgrz_Pen       =       0.6
367bmelt_other     =       0.7
368bmgrz_other     =       0.6
369bmelt_talus     =       5.
370bmgrz_talus     =       5.
371bmelt_coef      =       1.36138040668331
372file_number_shelves = "numer-ice-shelves-juil07-Ant15km.nc"
373/
374
375&beta_prescr                                   ! dragging_prescr_beta'
376beta_c_file = 'beta_AI2Sn138.grd'              !'beta-GISCAM23_20ans.grd' 'beta_Fab_10avril12_bord.grd'    ! 'beta_fab_comb_alpha_v4.grd'
377beta_limgz  = 5000000.                         ! -1.e12 pour plastic ?    ! .5e6 en visqueux   
378beta_min    = 1.                               ! for grounded ice (Pa)  ! 10.
379beta_mult   = 1.                               ! coefficient multiplicateur   
380/
381! read beta on centered grid
382! beta_file : nom des fichiers qui contiennent les betamx et betamy
383! beta-estime-run07.dat calcul direct beta sans modif
384! beta-staggered-15km-opt1.dat premier test en utilisant beta calcule  avec iterations
385! above beta_limgz, gzmx is false if negative : no sliding
386! betamx_file = not read anymore
387! betamy_file = not read anymore
388!___________________________________________________________
389&beta_iter_vitbil                              ! beta_iter_vitbil
390
391time_iter      = 5                             ! temps de demarrage des iterations
392nb_iter_vitbil = 2                             ! nombre d'iterations par pas de temps (dtt)
393coef_iter_vitbil = 1.                          ! coefficient pour rapport des vitesses <=1
394Umag_bil_file  = 'ebauche_vitbil_GrIce2sea.grd'  ! fichier des vitesses de bilan centrees
395/
396!
397!___________________________________________________________
398&spinup                                ! warning : 2 different modules
399ispinup = 0
400/
401! with module no_spinup
402! ispinup = 0     run standard ou calcul du beta
403! ispinup = 1     temperature equilibrium with grisli velocities'
404!
405! with module spinup_vitbil
406! ispinup = 2     conservation de la masse avec vitesses bilan '
407! ispinup = 3     equilibre temperature avec vitesses bilan'
408!
409!___________________________________________________________
410&vitbil_upwind             ! nom du bloc vitbil calcule sur le mailles staggered
411balance_Ux_file = 'Uy_moy_stag_5km.nc'
412balance_Uy_file = 'Ux_moy_stag_5km.nc'
413/
414!balance_Ux_file = 'Ux_moy_stag_5km.nc'
415!balance_Uy_file = 'Uy_moy_stag_5km.nc'
416!balance_vel_file = 'Umoy-for-calc-beta_LBq15-06_ZBL_15km.dat'
417! balance_vel_file = 'Uslid-for-calc-beta_LBq15-06_ZBL_15km.dat'
418! balance_vel_file = 'Umoy-for-calc-beta_LBq15-06_ZBL_15km.dat'
419! balance_vel_file = 'Vitbil-stag-pour-spinup_merge_runs_Lbq_ZBL_15km.dat'
420! balance velocities on staggered grid Ux, Uy
421!___________________________________________________________
422&output_regions
423nbregions       = 14
424region_file     = 'mask_regions_15km-0.grd'
425/
426! nbregions   = nombre de regions
427! region_file = carte des regions
428!_______________________________________________________________________________
Note: See TracBrowser for help on using the repository browser.