Changeset 135

Timestamp:

10/10/17 16:42:23 (7 years ago)

Author:

aquiquet

Message:

initMIP-Antarctica: setup for forward experiments

Location:

trunk

Files:

• INPUT/ANT16km/basal_melt_anomaly_16km.nc (added)
• INPUT/ANT16km/shelves-from-initMIP.grd (added)
• INPUT/ANT16km/smb_anomaly_16km.nc (added)
• SOURCES/Ant16_files/bmelt-ant-regions-initmip_mod.f90 (added)
• SOURCES/Ant16_files/climat_InitMIP_years_perturb_mod.f90 (added)
• SOURCES/Fichiers-parametres/ant16km_param_list.dat (modified) (10 diffs)
• SOURCES/Makefile.grisli.inc (modified) (2 diffs)

trunk/SOURCES/Fichiers-parametres/ant16km_param_list.dat

@@ -3 +3 @@
 !___________________________________________________________
 &runpar                       ! nom du bloc parametres du run
- runname      =  "inbeta02"   ! 8 caracteres 
- icompteur    =    2   
+ runname      =  "Ant16tut"   ! 8 caracteres 
+ icompteur    =    1
  iout         =    2          
- reprcptr     =  "../Ant16-03/Ant16-03+k30.nc"
+ reprcptr     =  "../../../Restarts/Ant16rlz_2.nc" !Ant16-05_1.nc" !Ant16-03+k30.nc"
  itracebug    = 0
  num_tracebug = 163     
- comment_run  = "Calcul beta"
+ comment_run  = "un commentaire"
 ! runname    : nom de l experience (8 caracteres)
 ! icompteur  : reprise dans un fichier  0 -> non, 1 -> oui, 2 -> T et Hwat
@@ -18 +18 @@
 &grdline                 ! bloc grounding line
 
- igrdline     =            1
+ igrdline     =            0
  Schoof       =            1
  ibmelt_inv   =            0
@@ -28 +28 @@
 &timesteps                ! bloc timestep
 
- tend      =   220
+ tend      =   100.
  tbegin    =   0.         !1.e10      ! si tbegin > 1.e9 on prend le temps du fichier cptr
  dtmin     =   1.e-3
- dtmax     =   0.5  
- dtt         =   1.  
+ dtmax     =   0.5 !10. !1.  
+ dtt       =   1. !10. !5.  
  testdiag  =   0.005      !0.005 cat 15km   
  /
@@ -83 +83 @@
  hmax_till      =    20.00000    
  poro_till      =   0.5000000    
- kond0          =   1.000000E-06   !1.000000E-06
+ kond0          =   1.e-6   !1.000000E-06
  
 ! hmax_till (m) : epaisseur max du sediment 
@@ -89 +89 @@
 ! conductivite du sediment :  kond0 (m/s)
 
-!___________________________________________________________
-&drag_param_beta                              ! drag_beta_param_mod
-
-beta_slope    = 0.00236050461661071    !1.e-3
-beta_expo     = 1.     !1.0
-betamax    = 5.e5     ! maximum value of beta
-betamin    = 1 !10.   ! minimum value of beta
-coef_ile   = 0.1
-/
-
 !____________________________________________________________
 &drag_neff_slope              ! nom du bloc dragging neff slope
 
-cf              =       PARAMCF  !2.e-5             ! 1.e-4
-betamax         =       BETAMAX
+cf              =       20.e-5  !2.e-5             ! 1.e-4
+betamax         =       1000.
 betamin         =       10.
 toblim          =       0.7e5              ! 0.25e5
-tostick         =       TOSTICK
-seuil_neff      =       PARAMSEUILNEFF  !350.e5
+tostick         =       1.e5
+seuil_neff      =       5.e5  !350.e5
 coef_gz         =       10.
 coef_ile        =       0.1
 slope_fich      = "ETOPO1_slope_R2min_Ant40km.nc"
-expo_slope      =       EXPOSLOPE
+expo_slope      =       1.
 pente_min      =        500.
 pente_max      =        2000.
@@ -250 +240 @@
 &clim_pert_massb                                   ! nom du bloc
 
-coefT           =       1.
+coefT           =       0.
 filforc         =     'forcage-800k-2007-Jouzel-Bassinot.dat' !'LGM_permanent.dat'
-pertsmb         =     1  !1
+pertsmb         =     0  !1
 rapsmb          =     0.07
 /
 ! pertsmb = 1, smb pert via tpert / pertsmb = 0, smb = cst
 !----------------------------------------------------------
-& meca_SIA_L1                        ! bloc resol_meca 
+&meca_SIA_L1                        ! bloc resol_meca 
 
 i_resolmeca     =        2
@@ -298 +288 @@
 coef_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
 temp_annual_file = 'Tann-MAR-ant35km-initMIP-16km.nc'  ! annual surface temperature
+/
+!----------------------------------------------------------
+&smb_anom_initMIP
+file_smb_anom     = 'smb_anomaly_16km.nc'
+massb_time        = 0     ! 0=fixe, 1:anomalies
 /
 !----------------------------------------------------------
@@ -347 +342 @@
 file_number_shelves = "numer-ice-shelves-juil07-Ant16km.nc"
 /
-
+!___________________________________________________________ 
+
+&bmelt_ant_reg_initmip                    !module bmelt-ant-regions_mod
+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
+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
+bmelt_talus     =       5.
+bmgrz_talus     =       5. 
+bmelt_coef      =       1. 
+file_number_shelves = "shelves-from-initMIP.grd"
+/
+!----------------------------------------------------------
+! module bmelt-ant-regions_mod
+&bmelt_anom_initMIP
+file_bmelt_anom     = 'basal_melt_anomaly_16km.nc'
+bmelt_time        = 0     ! 0=fixe, 1:anomalies
+/
+!___________________________________________________________
 &beta_prescr                                   ! dragging_prescr_beta'
-beta_c_file = 'beta-run-inbeta01.grd'         ! beta cat Ant interpole sur 16km
+beta_c_file = 'beta_AI2Sn138-16km.grd' !'beta-run-inbeta01.grd'         ! beta cat Ant interpole sur 16km
 beta_limgz  = .5e6                             ! -1.e12 pour plastic ?    ! .5e6 en visqueux   
 beta_min    = 1 !10.                           ! for grounded ice (Pa)
@@ -373 +384 @@
 &drag_param_beta                              ! drag_beta_param_mod
 
-beta_intercept = PARAMBINTERCEPT !-0.08898
-beta_slope     = 0.00236050461661071     !0.57570
+beta_slope    = 1.e-3    !1.e-3
+beta_expo     = 1.     !1.0
 betamax    = 5.e5     ! maximum value of beta
 betamin    = 1 !10.   ! minimum value of beta
@@ -383 +394 @@
 &spinup                                ! warning : 2 different modules
 ispinup      = 0
+!type_vitbil  = 0
 /
 ! with module no_spinup

trunk/SOURCES/Makefile.grisli.inc

@@ -43 +43 @@
         climat-forcage-insolation_mod.o climat_GrIce2sea_years_mod.o \
         climat_GrIce2sea_years_perturb_mod.o \
+        climat_InitMIP_years_perturb_mod.o \
         climat-perturb_mod-0.4.o \
         ablation_mod.o no_ablation_mod.o 
@@ -195 +196 @@
                 lect-anteis_mod.o \
                 bmelt-ant-regions_mod.o \
+                bmelt-ant-regions-initmip_mod.o \
                 fake-routines-ant_mod.o \
                 beta_iter_vitbil_mod.o \