Changeset 199


Ignore:
Timestamp:
07/04/18 17:06:06 (6 years ago)
Author:
aquiquet
Message:

A few additions for inputs / restart files / param_list templates

Location:
trunk
Files:
10 added
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/SOURCES/Fichiers-parametres/ant16km_param_list.dat

    r135 r199  
    33!___________________________________________________________ 
    44&runpar                       ! nom du bloc parametres du run 
    5  runname      =  "Ant16tut"   ! 8 caracteres  
    6  icompteur    =    1 
     5 runname      =  "3abumipc"   ! 8 caracteres  
     6 icompteur    =    1    
    77 iout         =    2           
    8  reprcptr     =  "../../../Restarts/Ant16rlz_2.nc" !Ant16-05_1.nc" !Ant16-03+k30.nc" 
     8 reprcptr     =  "../Fichier-CPTR/ABUMIP-cptr" !"../Ant16-temp-betacalc04-2_ispinup1/tp-c04-2+k30.nc" 
    99 itracebug    = 0 
    1010 num_tracebug = 163      
    11  comment_run  = "un commentaire" 
     11 comment_run  = "simule abuc abumip controle de 500 ans" 
    1212! runname    : nom de l experience (8 caracteres) 
    1313! icompteur  : reprise dans un fichier  0 -> non, 1 -> oui, 2 -> T et Hwat 
    … …  
    2424! igrdline :  1 ligne d echouage fixée, 3 recul forcé, 0 sans traitement specifique, 2 paleo gr line 
    2525! Schoof   :  0 pas de Schoof, 1 flux de Schoof 
    26 ! ibmelt_inv : 0 cas std, 1 inversion du bmelt (avec igrdline=1) 
     26 
    2727!___________________________________________________________ 
    2828&timesteps                ! bloc timestep 
    2929 
    30  tend      =   100. 
     30 tend      =   500. 
    3131 tbegin    =   0.         !1.e10      ! si tbegin > 1.e9 on prend le temps du fichier cptr 
    3232 dtmin     =   1.e-3 
    33  dtmax     =   0.5 !10. !1.   
    34  dtt       =   1. !10. !5.   
    35  testdiag  =   0.005      !0.005 cat 15km    
     33 dtmax     =   0.5 
     34 dtt         =   1.   
     35 testdiag  =   0.005      !0.005 cat 15km 
    3636 / 
    3737! tous les temps en annees. tbegin et tend : debut et fin du run  
    … …  
    4444!___________________________________________________________ 
    4545&topo_file 
    46  topo_ref     = "SHB_bedmap2_Ant16.nc"  
    47  topo_dep     = "SHB_bedmap2_Ant16.nc"  
     46 topo_ref     = "SHB_bedmap2_Ant16_relax20ans.nc"  
     47 topo_dep     = "SHB_bedmap2_Ant16_relax20ans.nc"  
    4848 grid_topo    = "coord-Ant-16km.dat" 
    49  ghf_fich     = "ghf-initMP-16km.nc" 
     49 ghf_fich     = "ghf_Ant_Lebrocq_16km.grd" 
    5050/ 
    5151! topo_ref= topo ref isostasie 
    … …  
    5757&mass_conserv            ! conservation de la masse 
    5858adv_frac        =  2.    ! 2-> advection seule 
    59 V_limit         = 3000.  ! maximum velocity for mass conservation : depend de la calotte 
     59V_limit         = 10000. !3000.  ! maximum velocity for mass conservation : depend de la calotte 
    6060/ 
    6161! Conservation de la masse avec equation advection-diffusion ' 
    … …  
    8383 hmax_till      =    20.00000     
    8484 poro_till      =   0.5000000     
    85  kond0          =   1.e-6   !1.000000E-06 
     85 kond0          =   6.14327918431256e-06   !1.000000E-06 
    8686  
    8787! hmax_till (m) : epaisseur max du sediment  
    … …  
    8989! conductivite du sediment :  kond0 (m/s) 
    9090 
     91!___________________________________________________________ 
     92&drag_param_beta                              ! drag_beta_param_mod 
     93 
     94beta_slope    = 0.00348349067112897    !1.e-3 
     95beta_expo     = 1.     !1.0 
     96betamax    = 5.e5     ! maximum value of beta 
     97betamin    = 1 !10.   ! minimum value of beta 
     98coef_ile   = 0.1 
     99/ 
     100 
    91101!____________________________________________________________ 
    92102&drag_neff_slope              ! nom du bloc dragging neff slope 
    93103 
    94 cf              =       20.e-5  !2.e-5             ! 1.e-4 
    95 betamax         =       1000. 
     104cf              =       PARAMCF  !2.e-5             ! 1.e-4 
     105betamax         =       BETAMAX 
    96106betamin         =       10. 
    97107toblim          =       0.7e5              ! 0.25e5 
    98 tostick         =       1.e5 
    99 seuil_neff      =       5.e5  !350.e5 
     108tostick         =       TOSTICK 
     109seuil_neff      =       PARAMSEUILNEFF  !350.e5 
    100110coef_gz         =       10. 
    101111coef_ile        =       0.1 
    102112slope_fich      = "ETOPO1_slope_R2min_Ant40km.nc" 
    103 expo_slope      =       1. 
     113expo_slope      =       EXPOSLOPE 
    104114pente_min      =        500. 
    105115pente_max      =        2000. 
    … …  
    173183&diagno_rheol              ! nom du bloc  diagno_rheol 
    174184 
    175 sf01           =        0.125   ! 0.125 
    176 sf03           =        0.125   ! 0.125 
     185sf01           =        0.8   ! 0.125 
     186sf03           =        0.8   ! 0.125 
    177187pvimin         =        1.5e3   ! 1.e3 
    178188/   
    … …  
    193203! gr_select  = ',gr_select 
    194204! gr_select = 1 : Tsai , 2 : Schoof' 
    195 ! 
     205!___________________________________________________________ 
     206 
     207! glissement                   module  sliding_Bindschadler  
     208&slid_bindsh                    ! nom du bloc 
     209  
     210kweert       =  PARAMKWEERT  !5.e-11 
     211loigliss     =  2 
     212coefbmax     =  10.  
     213 
     214! kweert : coefficent, loigliss le type de loi 
     215! coefbmax : facteur de normalisation pour influence eau 
    196216 
    197217!___________________________________________________________  
    … …  
    240260&clim_pert_massb                                   ! nom du bloc 
    241261 
    242 coefT           =       0. 
     262coefT           =       1. 
    243263filforc         =     'forcage-800k-2007-Jouzel-Bassinot.dat' !'LGM_permanent.dat' 
    244 pertsmb         =     0  !1 
     264pertsmb         =     1  !1 
    245265rapsmb          =     0.07 
    246266/ 
    247267! pertsmb = 1, smb pert via tpert / pertsmb = 0, smb = cst 
    248268!---------------------------------------------------------- 
    249 &meca_SIA_L1                        ! bloc resol_meca  
     269& meca_SIA_L1                        ! bloc resol_meca  
    250270 
    251271i_resolmeca     =        2 
    … …  
    285305!  module lect_climref_Ice2sea     
    286306&clim_smb_T_gen 
    287 smb_file         = 'SMB-MAR-ant35km-initMIP-16km.nc' ! smb m/an (en kg/m2/an) 
     307smb_file         = 'RACMO2.3p2_ANT27_smb_ltm_1979_2016_ant16.grd' ! smb m/an (en kg/m2/an) 
    288308coef_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 
    289 temp_annual_file = 'Tann-MAR-ant35km-initMIP-16km.nc'  ! annual surface temperature 
    290 / 
    291 !---------------------------------------------------------- 
    292 &smb_anom_initMIP 
    293 file_smb_anom     = 'smb_anomaly_16km.nc' 
    294 massb_time        = 0     ! 0=fixe, 1:anomalies 
     309temp_annual_file = 'RACMO2.3p2_ANT27_t2m_ltm_1979_2016_ant16.grd'  ! annual surface temperature 
    295310/ 
    296311!---------------------------------------------------------- 
    … …  
    340355bmgrz_talus     =       5.  
    341356bmelt_coef      =       1.  
    342 file_number_shelves = "numer-ice-shelves-juil07-Ant16km.nc" 
    343 / 
    344 !___________________________________________________________  
    345  
     357file_number_shelves = "numer-ice-shelves-sept16.nc" 
     358/ 
     359!___________________________________________________________ 
    346360&bmelt_ant_reg_initmip                    !module bmelt-ant-regions_mod 
    347 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 
    348 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 
     361bmelt_regions      = 0.01,0.01,0.01,0.01,0.2945439,0.01,0.01,0.01,0.01,2.329221,0.9675836,2.731526,0.01,0.2468391,0.01,0.01,1.272469,0.2656122 
     362bmgrz_regions      = 0.8843350,1.642938,1.219054,2.189940,1.221823,4.003392,15.97463,5.132639,0.01,3.371017,2.595406,4.477371,0.3003807,1.400497,2.378522,1.905212,9.594900,1.862373 
    349363bmelt_talus     =       5. 
    350 bmgrz_talus     =       5.  
    351 bmelt_coef      =       1.  
     364bmgrz_talus     =       5. 
     365bmelt_coef      =       1. 
    352366file_number_shelves = "shelves-from-initMIP.grd" 
    353 / 
     367flag_dist           = 0 
     368file_dist_talus     = "distance_talus_1500m.grd" 
     369bmelt_dist0         = 20. 
     370/ 
     371! 14&18=Ronnie-Filschner // 15&16=Ross // 17=PIG // 5=Amery 
    354372!---------------------------------------------------------- 
    355373! module bmelt-ant-regions_mod 
    … …  
    358376bmelt_time        = 0     ! 0=fixe, 1:anomalies 
    359377/ 
    360 !___________________________________________________________ 
    361378&beta_prescr                                   ! dragging_prescr_beta' 
    362 beta_c_file = 'beta_AI2Sn138-16km.grd' !'beta-run-inbeta01.grd'         ! beta cat Ant interpole sur 16km 
     379beta_c_file = 'ABUMIP-beta.grd'          !beta-GISCAM23_20ans.grd' 
    363380beta_limgz  = .5e6                             ! -1.e12 pour plastic ?    ! .5e6 en visqueux    
    364381beta_min    = 1 !10.                           ! for grounded ice (Pa) 
    … …  
    378395nb_iter_vitbil = 2                             ! nombre d'iterations par pas de temps (dtt) 
    379396coef_iter_vitbil = 1.                          ! coefficient pour rapport des vitesses <=1 
    380 Umag_bil_file  = 'vitbil-run-inbeta01.grd'  ! fichier des vitesses de bilan centrees 
     397Umag_bil_file  = 'vitbil-run-beta-ab4.grd'  ! fichier des vitesses de bilan centrees 
     398time_reiter    = 5000.                         ! temps entre 2 iterations 
    381399/ 
    382400! 
    … …  
    384402&drag_param_beta                              ! drag_beta_param_mod 
    385403 
    386 beta_slope    = 1.e-3    !1.e-3 
    387 beta_expo     = 1.     !1.0 
     404beta_intercept = PARAMBINTERCEPT !-0.08898 
     405beta_slope     = 0.00348349067112897     !0.57570 
    388406betamax    = 5.e5     ! maximum value of beta 
    389407betamin    = 1 !10.   ! minimum value of beta 
    … …  
    393411!___________________________________________________________ 
    394412&spinup                                ! warning : 2 different modules 
    395 ispinup      = 0 
    396 !type_vitbil  = 0 
     413ispinup = 0 
    397414/ 
    398415! with module no_spinup 
    … …  
    403420! ispinup = 2     conservation de la masse avec vitesses bilan ' 
    404421! ispinup = 3     equilibre temperature avec vitesses bilan' 
    405 ! type_vitbil type de lecture des vitesses de bilan 1: lect_vitbil, 2 : lect_vitbil_Lebrocq 
     422! 
    406423!___________________________________________________________ 
    407424&vitbil_upwind             ! nom du bloc vitbil calcule sur le mailles staggered 
    … …  
    417434! balance velocities on staggered grid Ux, Uy  
    418435!___________________________________________________________ 
    419 &vitbil_Lebrocq             ! nom du bloc vitbil obtenu par code Lebrocq 
    420 vit_balance_file = 'gwavel.mar.War10.grd' 
    421 flowdir_balance_file = 'gwaflowdir.mar.War10.grd' 
    422 / 
    423 ! vit_balance_file : balance velocity file 
    424 ! flowdir_balance_file : balance flow direction file 
    425 !___________________________________________________________ 
    426436&output_regions 
    427437nbregions       = 7 
Note: See TracChangeset for help on using the changeset viewer.