! Parametres du run :CY-S-ANT anteis1 ! sources dans ../../SOURCES !___________________________________________________________ &runpar ! nom du bloc parametres du run runname = "RUNNAME" icompteur = 1 iout = 2 reprcptr = "../../../Fichier-CPTR/TB40S252-k126-cptr.nc" itracebug = 0 num_tracebug = 163 comment_run = "Cycle Antarctique Schoof" / ! runname : nom de l experience (8 caracteres) ! icompteur : reprise dans un fichier 0 -> non, 1 -> oui, 2 -> T et Hw ! icompteur : 3-> T seulement ! iout : 1-> sortie cptr pour reprise,2 -> sortie nc pour reprise ! reprcptr : nom du fichier restart ! itracebug: 1-> ecriture de traces au debut des routines ! num_tracebug numero d unite ecriture traces ! comment_run: commentaire court sur le run !___________________________________________________________ &grdline ! bloc grounding line igrdline = 0 Schoof = 1 ibmelt_inv = 0 / ! igrdline : 1 ligne d echouage fixée, sinon 0 ! Schoof : 0 pas de Schoof, 1 flux de Schoof ! ibmelt_inv : 0 cas std, 1 inversion du bmelt (avec igrdline=1) !___________________________________________________________ ×teps ! bloc time steps tend = RUNTEND tbegin = RUNTBEGIN dtmin = 2.0000001E-03 dtmax = 1.000000 dtt = 5.000000 testdiag = 1.6000001E-02 / ! tous les temps en annees. tbegin et tend : debut et fin du run ! pour equation masse, pas de temps mini -> dtmin, maxi -> dtmax ! dtt : pas de temps long ! testdiag, pour gérer le pas de temps dynamique dt ! ordres de grandeur (a moduler selon dx) : ! 40 km dtmin=2.e-2, dtmax=1., dtt=5., tesdiag=0.02 !___________________________________________________________ &topo_file ! input_topo topo_ref = SHB_bedmap2_Ant40_relax20ans.nc topo_dep = SHB_bedmap2_Ant40_relax20ans.nc grid_topo =coord-Ant-40km.dat ghf_fich = ghf_Ant_Lebrocq_40km.nc / ! topo_ref= topo ref isostasie ! topo_dep= topo de depart ! grid_topo : fichier i,j,x,y,lon,lat ! ghf_fich : fichier flux geothermique ! module climat_Grice2sea_years_perturb_mod &LAPSE_RATES T_LAPSE_RATE = -8.0000004E-03 / !T_lapse_rate = lapse rate temp annuelle !________________________________________________________________ ! module climat_Grice2sea_years_perturb_mod &CLIM_PERT_MASSB COEFT = 1.000000 , FILFORC = forcage-429k-2007-Jouzel-Waelbroeck_AMOC.dat , PERTSMB = 1, RAPSMB = 7.0000000E-02, PERTBMB = 1, COEFBMB = 10.00000 / ! coefT : coef amplitude perturbation T ! filforc : fichier de forcage temporel ! pertsmb : 0 SMB fixe, 1 SMB exp temperature ! rapsmb : coef variation SMB ! pertbmb : 0 bmelt fixe, 1 bmelt variable ! coefbmb : coef variation bmelt !________________________________________________________________ ! module climat_Grice2sea_years_perturb_mod &CLIM_SMB_T_GEN SMB_FILE = RACMO2.3p2_ANT27_smb_ltm_1979_2016_ant40.grd , COEF_SMB_UNIT = 1.000000 , TEMP_ANNUAL_FILE = RACMO2.3p2_ANT27_t2m_ltm_1979_2016_ant40.grd / ! smb_file = fichier SMB (kg/m2/an) ! coef_smb_unit = coef passage m glace/an (1/910 ou 1/918) ! temp_annual_file = Temp moy annuelle (°C) !________________________________________________________________ ! module climat_Grice2sea_years_perturb_mod &CLIM_SNAP NB_SNAP = 200, TIME_DEPART_SNAPS = 2000000. , ECART_SNAP = 1.000000 , FILE_SMB_SNAP = smb_ECHAM5_A1B_MAR_highres_all_time.nc , MASSB_TIME = 0 / ! nb_snap = nombre de snapshots ! time_depart_snaps = debut du forçage ! ecart_snap = ecart entre les snapshots ! file_smb_snap = fichier serie temp anomalie SMB de GCM ! massb_time = 0:fixe, 1:interp snapshots, 2:snapsh+interp vert !_______________________________________________________________________ !___________________________________________________________ &eaubasale1 ! nom du premier bloc eau basale ecoulement_eau = T hwatermax = 5000.000 infiltr = 1.0000000E-03 / ! ecoulement eau : .false. -> modele bucket, sinon equ. diffusion ! hwatermax : hauteur d eau basale maximum dans le sediment (m) ! infiltr est la quantite d eau qui peut s infiltrer dans le sol (m/an) !___________________________________________________________ ¶m_hydr ! nom du bloc parametres hydrauliques hmax_till = 20.00000 poro_till = 0.5000000 kond0 = 1.7748050E-04 / ! hmax_till (m) : epaisseur max du sediment ! poro_till : porosite du sediment ! conductivite du sediment : kond0 (m/s) ! module bmelt_ant_regions_initmip &BMELT_ANT_REG_INITMIP BMELT_REGIONS = 2*0.5000000 , 0.1000000 , 1.000000 , 0.3000000 , 2*1.000000 , 0.5000000 , 1.000000 , 2.000000 , 1.000000 , 0.8000000 , 0.9000000 , 0.2200000 , 0.1200000 , 0.2000000 , 2*0.5000000 , BMGRZ_REGIONS = 6.000000 , 2*1.500000 , 6.000000 , 3.000000 , 6.500000 , 7.000000 , 4.000000 , 6.000000 , 5.000000 , 6.000000 , 8.000000 , 3.500000 , 0.5400000 , 0.4500000 , 0.2500000 , 4.700000 , 0.1000000 , BMELT_TALUS = 20.00000 , BMGRZ_TALUS = 20.00000 , BMELT_COEF = 0.9688509 , FILE_NUMBER_SHELVES = shelves-from-initMIP-40km.nc , FLAG_DIST = 1, FILE_DIST_TALUS = distance_talus_1500m.grd , BMELT_DIST0 = 20.00000 / ! bmelt_regions : fonte basale sous shelves pour les 18 regions initMIP ! bmgrz_regions : fonte basale grounding zone pour les 18 regions initMIP ! bmelt_talus & bmgrz_talus : fonte basale apres talus cont ! bmelt_coef : coef fonte (1 pour conserver val) ! file_numer_ice-shelves : fichier zones ice shelves ! flag_dist : flag pour bmelt fnct distance talus ! file_dist_talus : fichier de distance talus ! bmelt_dist0 : coef de bmelt au talus !_______________________________________________________________________ ! module bmelt_ant_regions_initmip &BMELT_ANOM_INITMIP FILE_BMELT_ANOM = basal_melt_anomaly_16km.nc , BMELT_TIME = 0 / ! file_bmelt_anom = fichier anomalie bmelt ! bmelt_time = 0:fixe, 1:anomalies !_______________________________________________________________________ &LOIDEF_1 EXPOSANT_1 = 3.000000 , TEMP_TRANS_1 = -6.500000 , ENHANC_FACT_1 = 2.732374 , COEF_COLD_1 = 1.6599999E-16, Q_COLD_1 = 78200.00 , COEF_WARM_1 = 2.0000000E-16, Q_WARM_1 = 95450.00 / &LOIDEF_2 EXPOSANT_2 = 1.000000 , TEMP_TRANS_2 = -10.00000 , ENHANC_FACT_2 = 1.000000 , COEF_COLD_2 = 8.3129997E-08, Q_COLD_2 = 40000.00 , COEF_WARM_2 = 8.3129997E-08, Q_WARM_2 = 60000.00 / !___________________________________________________________ ! loi de deformation module deformation_mod_2lois ! exposant (glen), temperature de transition (ttrans) ! enhancement factor (sf) ! pour les temperatures inf. a Temp_trans : ! coef_cold (Bat1) et Q_cold (Q1) ! pour les temperatures sup. a Temp_trans : ! coef_warm (Bat2) et Q_warm (Q2) !___________________________________________________________ !___________________________________________________________ &meca_SIA_L1 ! bloc resol_meca i_resolmeca = 2 / ! i_resolmeca type d association entre SIA et L1 ! i_resolmeca=0 chacun dans sa zone ! i_resolmeca=1 dans les zones stream, addition si uxdef > uxL1 (MIS11 Cairns) ! i_resolmeca=2 addition systematique dans les zones stream !___________________________________________________________ &furst_schoof ! nom du bloc frot_coef = 0.6000000 gr_select = 2 / ! frot_coef : solid friction law coef 0.6 in GMD 2018 ! gr_select = 1 : Tsai , 2 : Schoof !___________________________________________________________ &diagno_rheol ! nom du bloc diagno_rheol sf01 = 0.1250000 sf03 = 0.1250000 pvimin = 1500.000 / ! coefficients par rapport a la loi glace posee ! sf01 : coefficient viscosite loi lineaire ! sf03 : coefficient viscosite loi n=3 ! pvimin : valeur de pvi pour les noeuds fictifs ~ 1.e3 ! tres petit par rapport aux valeurs standards ~ 1.e10 !___________________________________________________________ ! spinup module no_spinup &SPINUP ISPINUP = 0 / ! ispinup = 0 run standard ! ispinup = 1 equilibre temperature avec vitesses grisli ! ispinup = 2 conservation de la masse avec vitesses bilan ! ispinup = 3 equilibre temperature avec vitesses bilan !___________________________________________________________ &drag_param_beta ! nom du bloc dragging param beta beta_slope = 1.1707746E-03 beta_expo = 1.000000 betamax = 500000.0 betamin = 1.000000 / ! Slope & expo of beta = - slope x Neff ** expo &MASS_CONSERV ADV_FRAC = 2.000000 , V_LIMIT = 3000.000 / ! Conservation de la masse avec equation advection-diffusion ! la repartition depend de adv_frac ! >1 -> advection seule ! 0 -> diffusion seule ! 0<*<1 -> fraction de l advection ! -1 -> zones diffusion + zones advection ! V_limit depend de la calotte : ! typiquement 3000 m/an en Antarctique, 10000 m/an au Groenland !___________________________________________________________ !___________________________________________________________ &calving ! nom du bloc calving méthode Vincent Hcoup_plateau = 250.0000 Hcoup_abysses = 250.0000 prof_plateau = 50.00000 prof_abysses = 1000.000 ifrange = 4 meth_hcoup = 0 / ! Hcoup epaisseurs de coupure pour les zones peu prodondes et profondes ! Hcoup_plateau pas de traitement particulier sur les bords ! ifrange=1 -> traitement de Vincent avec ice shelves frangeants partout ! ifrange=2 -> ice shelves frangeants seulement si bm-bmelt positif ! meth_hcoup pour faire eventuellement varier Hcoup avec le climat