Changeset 41

Timestamp:

02/03/16 09:38:41 (8 years ago)

Author:

dumas

Message:

Nouveau module de dragging : dragging_neff_contmaj_mod.f90 | correction d'un bug sur les indices dans relaxation_water_diffusion.f90 | suppression du coef 0.1 dans le calcul de betamx et betamy pour les points cotiers dans dragging_hwat_sedim_mod.f90 pour eviter le recul rapide de la grounding line en Antarctique

Location:

trunk/SOURCES

Files:

• Ant40_files/module_choix-antar40-0.4.f90 (modified) (1 diff)
• Fichiers-parametres/anteis1_TEMPS-NETCDF.dat (modified) (1 diff)
• Makefile.grisli-gfortran.inc (modified) (2 diffs)
• Makefile.grisli.inc (modified) (2 diffs)
• dragging_hwat_sedim_mod.f90 (modified) (2 diffs)
• dragging_neff_contmaj_mod.f90 (added)
• relaxation_water_diffusion.f90 (modified) (1 diff)

trunk/SOURCES/Ant40_files/module_choix-antar40-0.4.f90

@@ -87 +87 @@
 
 use dragging_hwat_contmaj ! frottement basal catherine grisliv8 ice stream peu actifs
+! use dragging_neff_contmaj ! dragging Aurelien et Christophe Fev 2016
 ! use dragging_hwat_sedim version sediment HN
 

trunk/SOURCES/Fichiers-parametres/anteis1_TEMPS-NETCDF.dat

@@ -20 +20 @@
 1.e10
 -------------------------------------------------------------------------------------------
-15      npredeft        lecture de predef_tsort
+8       npredeft        lecture de predef_tsort
 10
 100
 200
 500
-100010
-100100
-200010
-200100
--20990
--20900
--20500
--19500
--121990
--125990
--125900
--419990
--419500
--419000
-
-
+1000
+2000
+5000
+7000
 -------------------------------------------------------------------------------------------

trunk/SOURCES/Makefile.grisli-gfortran.inc

@@ -72 +72 @@
         dragging_hwatermax_0.2_mod.o dragging_calc_beta_mod.o  \
         dragging_hwat-contigu_mod.o dragging_hwat_contmaj_mod.o \
-        dragging_hwat_sedim_mod.o \
+        dragging_hwat_sedim_mod.o dragging_neff_contmaj_mod.o \
         calving_frange.o no_calving.o no_lakes.o \
         out_profile_mod.o printtable_mod.o mix-SIA-L1_mod.o \
@@ -700 +700 @@
 dragging_hwat_sedim_mod.o : dragging_hwat_sedim_mod.f90
         $(FT) dragging_hwat_sedim_mod.f90
+
+dragging_neff_contmaj_mod.o : dragging_neff_contmaj_mod.f90
+        $(FT) dragging_neff_contmaj_mod.f90
 
 dragging_calc_beta_mod.o : Antarctique_general_files/dragging_calc_beta_mod.f90

trunk/SOURCES/Makefile.grisli.inc

@@ -72 +72 @@
         dragging_hwatermax_0.2_mod.o dragging_calc_beta_mod.o  \
         dragging_hwat-contigu_mod.o dragging_hwat_contmaj_mod.o \
-        dragging_hwat_sedim_mod.o \
+        dragging_hwat_sedim_mod.o dragging_neff_contmaj_mod.o \
         calving_frange.o no_calving.o no_lakes.o \
         out_profile_mod.o printtable_mod.o mix-SIA-L1_mod.o \
@@ -704 +704 @@
 dragging_hwat_sedim_mod.o : dragging_hwat_sedim_mod.f90
         $(FT) dragging_hwat_sedim_mod.f90
+
+dragging_neff_contmaj_mod.o : dragging_neff_contmaj_mod.f90
+        $(FT) dragging_neff_contmaj_mod.f90
 
 dragging_calc_beta_mod.o : Antarctique_general_files/dragging_calc_beta_mod.f90

trunk/SOURCES/dragging_hwat_sedim_mod.f90

@@ -325 +325 @@
 
       if (cotemx(i,j)) then                        ! point cotier
-         betamx(i,j)=cf*neffmx(i,j)*0.1 
+         betamx(i,j)=cf*neffmx(i,j) 
 ! cat         betamx(i,j)=min(betamx(i,j),tobmax)
 ! version Jorge
@@ -379 +379 @@
 
       if (cotemy(i,j)) then                        ! point cotier
-         betamy(i,j)=cf*neffmy(i,j)*0.1 
+         betamy(i,j)=cf*neffmy(i,j) 
 ! cat         betamy(i,j)=min(betamy(i,j),tobmax)
 ! version Jorge

trunk/SOURCES/relaxation_water_diffusion.f90

@@ -128 +128 @@
           erelax(i,j)=-kmy(i,j+1)*dtwdx2      ! drelax : diagonale i,j+1 
 
-          crelax(i,j)=1.+((kmx(i,j)+kmx(i+1,j))+(kmy(i,j+1)+kmy(i,j+1)))*dtwdx2 
+          crelax(i,j)=1.+((kmx(i,j)+kmx(i+1,j))+(kmy(i,j)+kmy(i,j+1)))*dtwdx2 
                                               !crelax : diagonale i,j