source: trunk/SOURCES/Fichiers-parametres/A-LBq15_param_list-standard.dat @ 10

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

initial import GRISLI trunk

File size: 13.4 KB
Line 
1!  Parametres du run : Antar  A_LBq15
2
3
4!___________________________________________________________
5&runpar                  ! nom du bloc parametres du run
6 runname      =  "AntI2S11"   ! 8 caracteres
7 icompteur    =    2   
8 iout         =    2         
9 reprcptr     =   "../Fichier-CPTR/LB15C-33_2.nc" ! "Fichier-CPTR/LGMAn024_2k.nc" !"../Fichier-CPTR/LGMAn004_80k.nc" ! "../Fichier-CPTR/LGMAn022_7k.nc" 
10 itracebug    = 0
11 num_tracebug = 168     
12 comment_run  = "ice2sea idem run 5 sorties differentes"
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     =            1
22 /
23! igrdline :  1 ligne d echouage fixée, sinon 0
24! paleo grounding line -> 2
25!___________________________________________________________
26&timesteps                ! bloc timestep
27
28 tend      =   1.
29 tbegin    =   0.    !1.e10      ! si tbegin > 1.e9 on prend le temps du fichier cptr
30 dtmin     =   2.e-3
31 dtmax     =   1.         ! 1.        ! 10.   ! 0.1
32 dtt       =   1.   ! 5
33 testdiag  =   0.016      ! 0.025      ! 0.016
34 /
35! tous les temps en annees. tbegin et tend : debut et fin du run
36! pour equation masse, pas de temps mini -> dtmin, maxi -> dtmax
37! dtt : pas de temps long
38! testdiag, pour gerer le pas de temps dynamique dt
39! ordres de grandeur (a moduler selon dx) :
40! 40 km dtmin=2.e-2, dtmax=1., dtt=5., tesdiag=0.02
41 
42!___________________________________________________________
43&eaubasale1              ! nom du premier bloc eau basale
44 
45 ecoulement_eau =  T
46 hwatermax      =    5000.000   
47 infiltr        =   1.0000001E-03
48 /
49! ecoulement eau : .false. -> modele bucket, sinon equ. diffusion
50! hwatermax :  hauteur d eau basale maximum dans le sediment (m)
51! infiltr est la quantite d eau qui peut s infiltrer dans le sol (m/an)
52
53!___________________________________________________________
54&param_hydr             ! nom du  bloc parametres hydrauliques
55 
56 hmax_till      =    20.00000   
57 poro_till      =   0.5000000   
58 kond0          =   1.000000E-06
59 
60! hmax_till (m) : epaisseur max du sediment
61! poro_till : porosite du sediment
62! conductivite du sediment :  kond0 (m/s)
63
64!____________________________________________________________
65&drag_vit_bil_CISM_gen             ! nom du bloc dragging_vit_bil_CISM_gen
66
67hwatstream      =       50. 
68cf              =       1.e-5              ! 1.e-5 
69tobmax          =       20000.           
70toblim          =       0.7e5              ! 0.25e5
71seuil_vel       =       200.
72balance_vel_file = 'balance-vel_ZBL_15km.dat'
73/
74! hwatstream (m) :  critere de passage en stream en partant de la cote'
75! si hwater > hwatstream '
76! cf coefficient de la loi de frottement fonction Neff'
77! tobmax : (Pa) frottement maxi en Pa'
78! toblim : (Pa) tes pour les iles '
79
80!___________________________________________________________
81&drag_LGM                                       !  module dragging_LGM
82
83Tstream_lim     =       -1.                      ! temperature limite 
84cf              =       1.e-5                    ! 1.e-5
85betamax         =       1.e6
86beta_prescrit   =       1000.                    ! valeur du dragging
87masque_stream   = 'masque_ice_stream_LGM.grd'
88/
89! Tstream_lim :  critere de passage en stream en partant de la cote'
90! autorisé si T-Tpmp > Tstream_lim
91! cf coefficient de la loi de frottement fonction Neff'
92! betamax : (Pa) frottement maxi en Pa dans les streams'
93! toblim : (Pa)  pour les iles '
94!___________________________________________________________
95&drag_plastic_LGM                                !  module dragging_plastic_LGM
96
97Tstream_lim     =       -5.                      ! temperature limite 
98betamax         =       1.e6
99tob_prescrit    =       1.e5                     ! valeur du dragging
100masque_stream   =      'mask_U_crxULGM_50km.grd' !'masque-mixte-bed-fleuves_fin.grd'  ! 'masque_ice_stream_LGM.grd' ! 'mask_U_crxULGM_50km.grd' 'masque_outcrops_rough0.5.grd'
101/
102! Tstream_lim :  critere de passage en stream en partant de la cote'
103! autorisé si T-Tpmp > Tstream_lim
104! cf coefficient de la loi de frottement fonction Neff'q:qq
105! betamax : (Pa) frottement maxi en Pa dans les streams'
106! toblim : (Pa)  pour les iles '
107
108!___________________________________________________________
109&calving            !   nom du bloc calving méthode Vincent
110 
111Hcoup           =       1.2     !   tres petit quand shelves fixes sinon 250
112ifrange         =       4
113meth_Hcoup      =       0
114/
115
116! Hcoup epaisseur de coupure
117! ifrange=0 -> pas de traitement particulier sur les bords'
118! ifrange=1 -> traitement de Vincent avec ice shelves frangeants'
119! ifrange=2 -> ice shelves frangeant seulement si bm-bmelt positif
120!____________________________________________________________
121! loi de deformation 1          module deformation_mod_2lois
122&loidef_1
123
124exposant_1       = 3.
125temp_trans_1     = -6.5
126enhanc_fact_1    = 3.         
127coef_cold_1      = 1.660E-16
128Q_cold_1         = 7.820E+04
129coef_warm_1      = 2.000E-16
130Q_warm_1         = 9.545E+04
131/
132! exposant (glen), temperature de transition (ttrans)
133! enhancement factor (sf)
134! pour les temperatures inf. a Temp_trans :
135!            coef_cold (Bat1) et Q_cold (Q1)
136! pour les temperatures sup. a Temp_trans :
137!            coef_warm (Bat2) et Q_warm (Q2)
138!________________________________________________________
139! loi de deformation 2          module deformation_mod_2lois
140&loidef_2
141
142exposant_2       = 1.
143temp_trans_2     = -10.
144enhanc_fact_2    = 3.         
145coef_cold_2      = 8.313E-08
146Q_cold_2         = 4.000E+04
147coef_warm_2      = 8.313E-08
148Q_warm_2         = 6.000E+04
149/
150! exposant (glen), temperature de transition (ttrans)
151! enhancement factor (sf)
152! pour les temperatures inf. a Temp_trans :
153!            coef_cold (Bat1) et Q_cold (Q1)
154! pour les temperatures sup. a Temp_trans :
155!            coef_warm (Bat2) et Q_warm (Q2)
156!___________________________________________________________
157 
158&diagno_rheol              ! nom du bloc  diagno_rheol
159
160sf01           =        0.125  ! 1/8
161sf03           =        0.125
162pvimin         =        1.e4   ! 1.e3
163
164! coefficients par rapport a la loi glace posee                       
165! sf01 : coefficient viscosite loi lineaire
166! sf03 : coefficient viscosite loi n=3
167! pvimin : valeur de pvi pour les noeuds fictifs ~ 1.e3
168! tres petit par rapport aux valeurs standards ~ 1.e10
169!___________________________________________________________
170
171! glissement                   module  sliding_Bindschadler
172&slid_bindsh                    ! nom du bloc
173 
174kweert       =  0.             ! 5.e-11
175loigliss     =  2
176coefbmax     =  10.
177
178! kweert : coefficent, loigliss le type de loi
179! coefbmax : facteur de normalisation pour influence eau
180!___________________________________________________________
181&clim_pert                                    ! nom du bloc
182
183coefT           =       1.
184rappact         =       0.07                    ! 0.05->0.75  0.07 -> 0.5
185retroac         =       1
186rapbmshelf      =       5. 
187mincoefbmelt    =       0.
188maxcoefbmelt    =       2.
189filforc         =     'forcage-900k-2007-Parrenin-Bassinot.dat' !'LGM_permanent.dat'
190/
191!----------------------------------------------------------
192& meca_SIA_L1                        ! bloc resol_meca
193
194i_resolmeca     =        2
195/                           
196! i_resolmeca type d association entre SIA et L1'
197! i_resolmeca=0  chacun dans sa zone'
198! i_resolmeca=1  dans les zones stream, addition si uxdef > uxL1 (MIS11 Cairns)
199! i_resolmeca=2  addition systematique dans les zones stream
200
201!----------------------------------------------------------
202! PDD base Tann et Tjuly                module ablation_ann   
203
204&ablation_ann
205Cice          = 0.008
206Csnow         = 0.005
207Sigma         = 5.
208csi           = 0.6
209/
210! Cice and Csnow, melting factors for ice and snow
211! sigma variabilite Tday
212! csi proportion of melted water that can refreeze
213
214!----------------------------------------------------------
215!  module  lect_topo_ant_CISM    ATTENTION CE N'EST PLUS CELUI LA   
216& topo_CISM_gen
217topo_surf      = 'usrf-ice-equ-fin_ZBL_15km.dat'                   ! surface
218correc_surf    = 'no'                                              ! ice-real correction
219topo_thick     = 'thick2-ice-equ_ZBL_15km.dat'                     ! thickness
220topo_bed       = 'topg2_ZBL_15km.dat'                              ! bedrock
221mask_grounded  = 'newmask2grisli-15km.dat'                         ! mask
222longitude      = 'long_ZBL_LBq_15km.dat'                           ! longitude
223latitude       = 'lat_LBq_15km.grd'                        ! latitude
224heatflux       = 'ghfsr_15km.grd'                          ! geothermal heat flux en mW/m2
225/
226! here LeBrocq files
227! remark : in the mask, island causing trouble have been removed
228!----------------------------------------------------------
229!  module  lect_topo_ant_gen
230& topo_ant_gen
231topo_surf      = 'usrf-ice-equ_aout2010_15km.grd'                  !  surface
232correc_surf    = 'no'                                              ! ice-real correction
233topo_thick     = thick2-ice-equ_aout2010_15km.grd                  ! thickness
234topo_bed       = 'topg2_aout2010_15km.grd'                         ! bedrock
235mask_grounded  = 'newmask_aout2010_15km.grd'                       ! mask
236longitude      = 'long_ZBL_LBq_15km.dat'                           ! longitude a garder en .dat
237latitude       = 'lat_ZBL_LBq_15km.dat'                            ! latitude
238heatflux       = 'ghfsr_15km.grd'                                  ! geothermal heat flux
239/
240! here LeBrocq files
241! remark : in the mask, island causing trouble have been removed
242!----------------------------------------------------------
243!  module  lect_topo_ant_gen      lecture d une nouvelle topo de depart
244& topo_ant_startingpoint
245surf_start     = 'S_LGM_mini.grd'                  !  surface
246thick_start    = 'H_LGM_mini.grd'                  ! thickness
247bed_start      = 'topg2_aout2010_15km.grd'         ! bedrock
248mask_start     = 'grounded-LGM-15.grd'             ! mask
249sealevel       = -120
250/
251! grounding line from Anderson
252!----------------------------------------------------------
253!  module  lect_clim_acc_T_ant_gen 
254&climat_acc_T_gen
255
256precip_file    = 'acca_15km.grd'                                   ! precipitation
257coef_dens      = 1.                                                ! direct. en m glace
258temp_annual_file ='temp_15km.grd'                      ! annual surface temperature
259/
260!----------------------------------------------------------
261&bmelt_seuil                     ! module  bmelt_seuil_prof
262bm_grz           = 3.
263bmshelf_plateau  = .4
264bmshelf_abysses  = 10.
265depth_talus      = -3000.
266/                     
267! Pour l actuel : bm_grz a la grounding line
268!                 bmshelf_plateau sur le plateau continental
269!                 bmshelf_abysses pour les grandes profondeurs
270!       depth_talus, negative, separation entre les 2 domaines
271!___________________________________________________________
272&vitbil_upwind             ! nom du bloc vitbil calcule sur le mailles staggered
273balance_Ux_file = 'Ux-balvel-01.grd'
274balance_Uy_file = 'Uy-balvel-01.grd'
275/
276!balance_vel_file = 'Umoy-for-calc-beta_LBq15-06_ZBL_15km.dat'
277! balance_vel_file = 'Uslid-for-calc-beta_LBq15-06_ZBL_15km.dat'
278! balance_vel_file = 'Umoy-for-calc-beta_LBq15-06_ZBL_15km.dat'
279! balance_vel_file = 'Vitbil-stag-pour-spinup_merge_runs_Lbq_ZBL_15km.dat'
280! balance velocities on staggered grid Ux, Uy
281!___________________________________________________________
282
283! module dragging_vit_bil_LBq_gen
284 &drag_vit_bil_LBq_gen
285 HWATSTREAM              =   50.00000    ,
286 CF                      =  9.9999997E-06,
287 TOBMAX                  =   5000.000    ,
288 TOBLIM                  =   70000.00    ,
289 SEUIL_VEL               =   200.0000    ,
290 BALANCE_VEL_FILE        = 'gwavel-LeBrocq-grounded_ZBL_15km.dat'
291 /
292! hwatstream (m) :  critere de passage en stream vitesses de bilan
293!  si hwater > hwatstream
294! cf coefficient de la loi de frottement fonction Neff
295! seulement pour les points cotiers
296! tobmax : (Pa) frottement maxi sous les streams
297! toblim : (Pa)  pour les iles
298! seuil_vel (m/an) : seuil sur les vitesses pour definir le masque stream
299! balance_vel_file : nom du fichier qui contient les vitesse de bilan
300!___________________________________________________________
301&beta_prescr                                   ! dragging_prescr_beta'
302beta_c_file = 'beta-driving-over-sliding.grd'  !'beta_driving-over-velmax.grd' ! 'beta-from-vel.grd'   !beta_c-LB15C-29.nc'    !
303beta_limgz = .5e5                              ! -1.e12 pour plastic ?    ! .5e6 en visqueux   
304beta_min = 10.                                 ! for grounded ice (Pa)
305beta_mult = 1.                                 ! coefficient multiplicateur   
306/
307! read beta on centered grid
308! beta_file : nom des fichiers qui contiennent les betamx et betamy
309! beta-estime-run07.dat calcul direct beta sans modif
310! beta-staggered-15km-opt1.dat premier test en utilisant beta calcule  avec iterations
311! above beta_limgz, gzmx is false if negative : no sliding
312! betamx_file = not read anymore
313! betamy_file = not read anymore
314!___________________________________________________________
315&beta_stag                           !read beta on staggered grid'
316betamx_file = ''
317betamy_file = ''
318beta_limgz = 2.e3       
319/
320
321! above beta_limgz, gzmx is false
322!___________________________________________________________
323&spinup                                ! warning : 2 different modules
324ispinup = 0
325/
326! with module no_spinup
327! ispinup = 0     run standard ou calcul du beta
328! ispinup = 1     temperature equilibrium with grisli velocities'
329!
330! with module spinup_vitbil
331! ispinup = 2     conservation de la masse avec vitesses bilan '
332! ispinup = 3     equilibre temperature avec vitesses bilan'
333!
334!___________________________________________________________
Note: See TracBrowser for help on using the repository browser.