Changeset 244 for branches/iLoveclim/SOURCES/Ant16_files

Timestamp:

01/16/19 16:08:14 (5 years ago)

Author:

aquiquet

Message:

Grisli-iloveclim branch merged to trunk at revision 243

Location:

branches/iLoveclim

Files:

• . (modified) (1 prop)
• SOURCES/Ant16_files/bmelt-ant-regions-initmip_mod.f90 (modified) (2 diffs)
• SOURCES/Ant16_files/bmelt-seuil-profondeur_initmip_inv_mod.f90 (copied) (copied from trunk/SOURCES/Ant16_files/bmelt-seuil-profondeur_initmip_inv_mod.f90)
• SOURCES/Ant16_files/climat_InitMIP_years_perturb_mod.f90 (modified) (2 diffs)
• SOURCES/Ant16_files/module_choix-ant16km.f90 (modified) (7 diffs)
• SOURCES/Ant16_files/module_choix-ant16km.f90-iterbeta (modified) (5 diffs)
• SOURCES/Ant16_files/module_choix-ant16km.f90-spinupvitbil (modified) (2 diffs)
• SOURCES/Ant16_files/module_choix-ant16km.f90-std (modified) (9 diffs)
• SOURCES/Ant16_files/module_choix-ant16km.f90-std-abumip (copied) (copied from trunk/SOURCES/Ant16_files/module_choix-ant16km.f90-std-abumip)
• SOURCES/Ant16_files/module_choix-ant16km.f90-std-larmip (copied) (copied from trunk/SOURCES/Ant16_files/module_choix-ant16km.f90-std-larmip)
• SOURCES/Ant16_files/paradim-ant16_mod.f90 (modified) (1 diff)

branches/iLoveclim/SOURCES/Ant16_files/bmelt-ant-regions-initmip_mod.f90

@@ -64 +64 @@
     read(num_param,bmelt_ant_reg_initmip)
 
-    ! formats pour les ecritures dans 42
-428 format(A)
-
-    write(num_rep_42,428)'!___________________________________________________________' 
-    write(num_rep_42,428) '&bmelt-ant-regions-initmip                                 '
-    write(num_rep_42,*)
-    write(num_rep_42,'(A,18(f0.2,","))') 'bmelt_regions      = ', bmelt_regions(:)
-    write(num_rep_42,'(A,18(f0.2,","))') 'bmgrz_regions      = ', bmgrz_regions(:)
-    write(num_rep_42,*) 'bmelt_talus     = ', bmelt_talus 
-    write(num_rep_42,*) 'bmgrz_talus     = ', bmgrz_talus
-    write(num_rep_42,*) 'bmelt_coef      = ', bmelt_coef
-    write(num_rep_42,*) 'file_numer_shelves = ', file_number_shelves
-    write(num_rep_42,*) 'flag_dist = ', flag_dist
-    write(num_rep_42,*) 'file_dist_talus = ', file_dist_talus
-    write(num_rep_42,*) 'bmelt_dist0 = ', bmelt_dist0
-    write(num_rep_42,*)'/'                            
-    write(num_rep_42,428) '! bmelt_regions    :  fonte basale sous shelves pour les 18 regions initMIP'
-    write(num_rep_42,428) '! bmgrz_regions    :  fonte basale grounding zone pour les 18 regions initMIP'
-    write(num_rep_42,428) '! bmelt_talus & bmgrz_talus  :  fonte basale apres talus cont'
-    write(num_rep_42,428) '! bmelt_coef                 :  coef fonte (1 pour conserver val)'
-    write(num_rep_42,428) '! file_numer_ice-shelves     : fichier zones ice shelves'
-    write(num_rep_42,428) '! flag_dist : flag pour bmelt fnct distance talus'
-    write(num_rep_42,428) '! file_dist_talus : fichier de distance talus'
-    write(num_rep_42,428) '! bmelt_dist0 : coef de bmelt au talus'
-    write(num_rep_42,*)
+    write(num_rep_42,'(A)') '!  module bmelt_ant_regions_initmip '
+    write(num_rep_42,bmelt_ant_reg_initmip)                       
+    write(num_rep_42,'(A)') '! bmelt_regions    :  fonte basale sous shelves pour les 18 regions initMIP'
+    write(num_rep_42,'(A)') '! bmgrz_regions    :  fonte basale grounding zone pour les 18 regions initMIP'
+    write(num_rep_42,'(A)') '! bmelt_talus & bmgrz_talus  :  fonte basale apres talus cont'
+    write(num_rep_42,'(A)') '! bmelt_coef                 :  coef fonte (1 pour conserver val)'
+    write(num_rep_42,'(A)') '! file_numer_ice-shelves     : fichier zones ice shelves'
+    write(num_rep_42,'(A)') '! flag_dist : flag pour bmelt fnct distance talus'
+    write(num_rep_42,'(A)') '! file_dist_talus : fichier de distance talus'
+    write(num_rep_42,'(A)') '! bmelt_dist0 : coef de bmelt au talus'
+    write(num_rep_42,'(A)')'!_______________________________________________________________________'
 
 
@@ -138 +124 @@
   read(num_param,bmelt_anom_initMIP)
 
-  write(num_rep_42,428)'!_______________________________________________________________________' 
-  write(num_rep_42,428)'!  module bmelt_ant_regions_initmip_mod                                 '
+  write(num_rep_42,'(A)')'!  module bmelt_ant_regions_initmip                                     '
   write(num_rep_42,bmelt_anom_initMIP)
-  write(num_rep_42,428)'! file_bmelt_anom   = fichier anomalie bmelt                            '
-  write(num_rep_42,428)'! bmelt_time        = 0:fixe, 1:anomalies                               '
-  write(num_rep_42,428)'!_______________________________________________________________________' 
+  write(num_rep_42,'(A)')'! file_bmelt_anom   = fichier anomalie bmelt                            '
+  write(num_rep_42,'(A)')'! bmelt_time        = 0:fixe, 1:anomalies                               '
+  write(num_rep_42,'(A)')'!_______________________________________________________________________' 
 
   if ( bmelt_time .eq. 1 ) then  

branches/iLoveclim/SOURCES/Ant16_files/climat_InitMIP_years_perturb_mod.f90

@@ -16 +16 @@
 
 
-use module3d_phy,only: nx,ny,S,S0,H,Bsoc,acc,abl,BM,Tann,Tjuly,Ts,time,dt,num_param,num_rep_42,num_forc,dirforcage,dirnameinp,tafor,time,sealevel,coefbmshelf
+use module3d_phy,only: nx,ny,S,S0,H,Bsoc,acc,abl,BM,Tann,Tjuly,Ts,time,dt,num_param,num_rep_42,num_forc,dirforcage,dirnameinp,tafor,time,sealevel,sealevel_2d,coefbmshelf
 use netcdf
 use io_netcdf_grisli
@@ -285 +285 @@
 100 continue
 
+  sealevel_2d(:,:) = sealevel
+
   Tann (:,:) = Ta0 (:,:) + T_lapse_rate * (S(:,:)-S0(:,:)) +Tafor
   Ts(:,:)    = Tann(:,:)

branches/iLoveclim/SOURCES/Ant16_files/module_choix-ant16km.f90

@@ -12 +12 @@
 !! @note   - use lect_clim_act_anteis
 !! @note   - use climat_perturb_mod
-!! @note   - use no_lakes
+!! @note   - use no_rsl
 !! @note   - use isostasie_mod  
 !! @note   - use deformation_mod_2lois  
@@ -45 +45 @@
 use climat_Grice2sea_years_mod  ! climat force par fichier SMB directement (grice2sea)
 !use climat_Grice2sea_years_perturb_mod ! climat force par fichier SMB directement (grice2sea) + index temperature carotte de glace
-
+!use climat_InitMIP_years_perturb_mod ! climat pour experiences initMIP
 
 !use climat_forcage_mod
@@ -55 +55 @@
 use no_ablation  ! pas de calcul de l'ablation => lecture fichier SMB (necessaire avec climat_Grice2sea_years_mod)
 
-! pas de lacs proglaciaires
-use no_lakes
+! pas de differences locales de niveau marin
+use no_rsl
 
 ! suivi des traceurs ou pas ?
@@ -81 +81 @@
 
 
-!use dragging_prescr_beta 
 !use dragging_vitbil
 !use dragging_hwatstream
@@ -88 +87 @@
 !use dragging_hwat_contmaj ! frottement basal catherine grisliv8 ice stream peu actifs
 !use dragging_neff_slope ! dragging Aurelien et Christophe
-!use dragging_hwat_sedim version sediment HN
+! use dragging_hwat_sedim version sediment HN
+!use dragging_beta_iter_vitbil_mod ! dragging pour inversion beta
+use dragging_prescr_beta ! utilisation d'une carte de beta (issue de l'inversion)
 
 !------------ spinup -----------------------------------
-!use no_spinup           ! spinup=0
-use spinup_vitbil         ! plusieurs variantes de spinup : inclue le dragging
+use no_spinup           ! spinup=0
+!use spinup_vitbil         ! plusieurs variantes de spinup : inclue le dragging
                           ! enlever le use dragging dans ce cas là
 
@@ -102 +103 @@
 ! calving
 use calving_frange
+!use calving_frange_abuk ! pour les experience abuk ABUMIP
 
 ! choix resolution de la conservation de la masse
@@ -109 +111 @@
 
 !--------------Fusion basale------------------------
-use bmelt_ant_regions ! pour l'Antarctique avec régions
+use bmelt_ant_regions_larmip ! pour experience fusion basale LARMIP
+!use bmelt_ant_regions_initmip  ! pour initmip abumip melt et ctrl
+!use bmelt_seuil_prof  ! fusion basale 2 valeurs : abysses et talus et inversion du bmelt
+!use bmelt_ant_regions ! pour l'Antarctique avec régions
 !use bmelt_ant_regions_oce ! pour Ant region et temp ocean
 !use bmelt_nor_regions ! pour le nord avec régions

branches/iLoveclim/SOURCES/Ant16_files/module_choix-ant16km.f90-iterbeta

@@ -12 +12 @@
 !! @note   - use lect_clim_act_anteis
 !! @note   - use climat_perturb_mod
-!! @note   - use no_lakes
+!! @note   - use no_rsl
 !! @note   - use isostasie_mod  
 !! @note   - use deformation_mod_2lois  
@@ -45 +45 @@
 use climat_Grice2sea_years_mod  ! climat force par fichier SMB directement (grice2sea)
 !use climat_Grice2sea_years_perturb_mod ! climat force par fichier SMB directement (grice2sea) + index temperature carotte de glace
-
+!use climat_InitMIP_years_perturb_mod ! climat pour experiences initMIP
 
 !use climat_forcage_mod
@@ -55 +55 @@
 use no_ablation  ! pas de calcul de l'ablation => lecture fichier SMB (necessaire avec climat_Grice2sea_years_mod)
 
-! pas de lacs proglaciaires
-use no_lakes
+! pas de differences locales de niveau marin
+use no_rsl
 
 ! suivi des traceurs ou pas ?
@@ -88 +88 @@
 !use dragging_neff_slope ! dragging Aurelien et Christophe
 ! use dragging_hwat_sedim version sediment HN
-use beta_iter_vitbil_mod
+use dragging_beta_iter_vitbil_mod ! dragging pour inversion beta
+!use dragging_prescr_beta ! utilisation d'une carte de beta (issue de l'inversion)
 
 !------------ spinup -----------------------------------
@@ -109 +110 @@
 
 !--------------Fusion basale------------------------
-!use bmelt_seuil_prof  ! fusion basale 2 valeurs : abysses et talus
-use bmelt_ant_regions ! pour l'Antarctique avec régions
+!use bmelt_ant_regions_larmip ! pour experience fusion basale LARMIP
+!use bmelt_ant_regions_initmip
+use bmelt_seuil_prof  ! fusion basale 2 valeurs : abysses et talus et inversion du bmelt
+!use bmelt_ant_regions ! pour l'Antarctique avec régions
 !use bmelt_ant_regions_oce ! pour Ant region et temp ocean
 !use bmelt_nor_regions ! pour le nord avec régions

branches/iLoveclim/SOURCES/Ant16_files/module_choix-ant16km.f90-spinupvitbil

@@ -12 +12 @@
 !! @note   - use lect_clim_act_anteis
 !! @note   - use climat_perturb_mod
-!! @note   - use no_lakes
+!! @note   - use no_rsl
 !! @note   - use isostasie_mod  
 !! @note   - use deformation_mod_2lois  
@@ -55 +55 @@
 use no_ablation  ! pas de calcul de l'ablation => lecture fichier SMB (necessaire avec climat_Grice2sea_years_mod)
 
-! pas de lacs proglaciaires
-use no_lakes
+! pas de differences locales de niveau marin
+use no_rsl
 
 ! suivi des traceurs ou pas ?

branches/iLoveclim/SOURCES/Ant16_files/module_choix-ant16km.f90-std

@@ -1 @@
-!> \file module_choix-ant16km.f90
+!> \file module_choix-ant16.f90
 !! Ce module permet de choisir quels package on veut utiliser
 !! 
@@ -12 +12 @@
 !! @note   - use lect_clim_act_anteis
 !! @note   - use climat_perturb_mod
-!! @note   - use no_lakes
+!! @note   - use no_rsl
 !! @note   - use isostasie_mod  
 !! @note   - use deformation_mod_2lois  
@@ -34 +34 @@
 
 !--------------Lecture topo ------------------------
-!use lect_topo_heino      ! pour Heino
-use lect_topo_anteis    ! pour l'antarctique 40 km, 15km Lebrocq ou initMIP
+! use lect_topo_heino      ! pour Heino
+use lect_topo_anteis    ! pour l'antarctique 40 km
 !use lect_topo_nord40    ! pour l'hemisphere nord 40 km
 !use lect_topo_eurasie   ! pour l'eurasie  (toute resolutions)
-
-!--------------Lecture climat ref ------------------
-! Le climat de ref si necessaire est maintenant lu par le module de climat directement
-!use lect_clim_act_anteis
 
 
@@ -47 +43 @@
 !use climat_forcage_mois_mod ! forcage mensuel GCM 1 Snapshot Fev 2015
 !use climat_perturb_mod   ! climat perturbe a reverifier Dec 2015
-use climat_Grice2sea_years_mod  ! climat force par fichier SMB directement (grice2sea)
+!use climat_Grice2sea_years_mod  ! climat force par fichier SMB directement (grice2sea)
 !use climat_Grice2sea_years_perturb_mod ! climat force par fichier SMB directement (grice2sea) + index temperature carotte de glace
-
+use climat_InitMIP_years_perturb_mod ! climat pour experiences initMIP
 
 !use climat_forcage_mod
@@ -59 +55 @@
 use no_ablation  ! pas de calcul de l'ablation => lecture fichier SMB (necessaire avec climat_Grice2sea_years_mod)
 
-! pas de lacs proglaciaires
-use no_lakes
+! pas de differences locales de niveau marin
+use no_rsl
 
 ! suivi des traceurs ou pas ?
@@ -68 +64 @@
 !--------------Choix isostasie----------------------
 ! use isostasie_mod  ! module permettant de calculer la deflexion isostasique
- use noisostasie_mod ! module pour ne pas avoir d'isostasie
+use noisostasie_mod ! module pour ne pas avoir d'isostasie
 
 
@@ -89 +85 @@
 !use dragging_hwat_cont
 
-use dragging_prescr_beta  
-!use beta_iter_vitbil_mod ! dragging avec beta obtenu par inversion
 !use dragging_hwat_contmaj ! frottement basal catherine grisliv8 ice stream peu actifs
 !use dragging_neff_slope ! dragging Aurelien et Christophe
 ! use dragging_hwat_sedim version sediment HN
+!use beta_iter_vitbil_mod
+use dragging_prescr_beta
 
 !------------ spinup -----------------------------------
@@ -114 +110 @@
 
 !--------------Fusion basale------------------------
+use bmelt_ant_regions_initmip
+!use bmelt_seuil_prof  ! fusion basale 2 valeurs : abysses et talus
 !use bmelt_ant_regions ! pour l'Antarctique avec régions
 !use bmelt_ant_regions_oce ! pour Ant region et temp ocean
-use bmelt_seuil_prof ! bmelt fonction de la profondeur (2 valeurs)
 !use bmelt_nor_regions ! pour le nord avec régions
 !use bmelt_nor_depth   ! pour le nord avec profondeur d'eau
@@ -132 +129 @@
 !-----Suivi temporel
 ! use output_heminord_mod
- use output_antarcti_mod
+use output_antarcti_mod
 ! use output_heino
 !-----debugage...

branches/iLoveclim/SOURCES/Ant16_files/paradim-ant16_mod.f90

@@ -20 +20 @@
 ! character
 character(len=7),  parameter  ::  geoplace='ant16km'
-character(len=20), parameter  :: dirnameinp='../../INPUT/ANT16km/' !< input directory
-character(len=20), parameter  :: dirforcage='../../INPUT/Forcage/' !< input directory
+character(len=17), parameter  :: dirnameinp='../INPUT/ANT16km/' !< input directory
+character(len=17), parameter  :: dirforcage='../INPUT/Forcage/' !< input directory
 
 ! dimensionnement grilles