Changeset 100

Timestamp:

04/12/17 13:57:18 (7 years ago)

Author:

aquiquet

Message:

Calving_frange: hcoup as a function of bathymetry. Updated hemin40_param_list.dat commited as an example for other geometries.

Location:

trunk/SOURCES

Files:

• Fichiers-parametres/hemin40_param_list.dat (modified) (4 diffs)
• calving_frange.f90 (modified) (5 diffs)

trunk/SOURCES/Fichiers-parametres/hemin40_param_list.dat

@@ -20 +20 @@
 
  igrdline     =            0
+ Schoof       =            0
  /
 ! igrdline :  1 ligne d echouage fixée, sinon 0
@@ -167 +168 @@
 &calving                        !   nom du bloc calving méthode Vincent 
  
-Hcoup           =       200     !   tres petit quand shelves fixes sinon 250
-ifrange         =       4 
-meth_Hcoup      =       0 
-/
-
-! Hcoup epaisseur de coupure 
-! ifrange=0 -> pas de traitement particulier sur les bords'
-! ifrange=1 -> traitement de Vincent avec ice shelves frangeants'
+Hcoup_plateau   =       50     !   tres petit quand shelves fixes sinon 250
+Hcoup_abysses   =       300     !   tres petit quand shelves fixes sinon 250
+prof_plateau   =       50     !   tres petit quand shelves fixes sinon 250
+prof_abysses   =       1000     !   tres petit quand shelves fixes sinon 250
+ifrange         =       4
+meth_Hcoup      =       0
+/
+
+! Hcoup epaisseurs de coupure pour les zones peu prodondes et profondes
+! Hcoup_plateau<Hcoup_abysses && prof_plateau<prof_abysse
+! prof profondeur delimitant les zones peu prodondes et profondes
+! ifrange=0 -> pas de traitement particulier sur les bords
+! ifrange=1 -> traitement de Vincent avec ice shelves frangeants
 ! ifrange=2 -> ice shelves frangeant seulement si bm-bmelt positif
+! meth_hcoup pour faire eventuellement varier Hcoup avec le climat
 !____________________________________________________________
 ! loi de deformation 1          module deformation_mod_2lois
@@ -226 +233 @@
 ! tres petit par rapport aux valeurs standards ~ 1.e10
 !___________________________________________________________ 
+
+! furst_schoof
+&furst_schoof
+
+frot_coef = 0.6   !0.035
+gr_select = 1
+/
+! frot_coef  = ',frot_coef
+! gr_select  = ',gr_select
+! gr_select = 1 : Tsai , 2 : Schoof'
+!___________________________________________________________
 
 ! glissement                   module  sliding_Bindschadler 
@@ -407 +425 @@
 !
 !___________________________________________________________
+&drag_param_beta                              ! drag_beta_param_mod
+
+beta_intercept = -2. !-0.08898
+beta_slope     = 0.95     !0.57570
+betamax    = 5.e5     ! maximum value of beta
+betamin    = 1 !10.   ! minimum value of beta
+coef_ile   = 0.1
+/
+!
+!___________________________________________________________
 &spinup                                ! warning : 2 different modules
 ispinup = 0

trunk/SOURCES/calving_frange.f90

@@ -21 +21 @@
 real, dimension (nx,ny) :: hmhc  ! hauteur au dessus de la coupure
 real, dimension (nx,ny) :: bil_tot ! bilan surface et fond (bm-bmelt)
-real :: hcoup ! epaiseur de coupure au temps time
-real :: hcoup_max  ! epaisseur max
-real :: hcoup_min=100.  ! epaisseur min
+real, dimension (nx,ny) :: hcoup ! epaiseur de coupure au temps time
+real :: hcoup_plateau  ! coupure points peu profonds
+real :: hcoup_abysses  ! coupure points ocean profond
+real :: prof_plateau  ! profondeur max des points peu profonds
+real :: prof_abysses  ! profondeur min des points ocean profond
 integer :: meth_hcoup  ! pour avoir hcoup dépendant du coefbmshelf
 integer :: ifrange        !  0 pas de traitement particulier pres du bord, 1 -> franges
@@ -36 +38 @@
 subroutine init_calving
 
-namelist/calving/Hcoup,ifrange,meth_hcoup
+namelist/calving/Hcoup_plateau,Hcoup_abysses,prof_plateau,prof_abysses,ifrange,meth_hcoup
 
 
@@ -55 +57 @@
 write(num_rep_42,428) '&calving              ! nom du bloc calving méthode Vincent '
 write(num_rep_42,*)
-write(num_rep_42,*)   'Hcoup      = ', Hcoup
-write(num_rep_42,*)   'ifrange    = ',ifrange
+write(num_rep_42,*)   'Hcoup_plateau  = ', Hcoup_plateau
+write(num_rep_42,*)   'Hcoup_abysses  = ', Hcoup_abysses
+write(num_rep_42,*)   'prof_plateau  = ', prof_plateau
+write(num_rep_42,*)   'prof_abysses  = ', prof_abysses
+write(num_rep_42,*)   'ifrange    = ', ifrange
 write(num_rep_42,*)   'meth_hcoup = ', meth_hcoup
 write(num_rep_42,*)'/'
                             
-write(num_rep_42,428) '! Hcoup epaisseur de coupure, le max si attache a coefbmshelf '
+write(num_rep_42,428) '! Hcoup epaisseurs de coupure pour les zones peu prodondes et profondes'
+write(num_rep_42,428) '! Hcoup_plateau<Hcoup_abysses && prof_plateau<prof_abysses'
+write(num_rep_42,428) '! prof profondeur delimitant les zones peu prodondes et profondes'
 write(num_rep_42,428) '! ifrange=0 -> pas de traitement particulier sur les bords'
 write(num_rep_42,428) '! ifrange=1 -> traitement de Vincent avec ice shelves frangeants partout'
 write(num_rep_42,428) '! ifrange=2 -> ice shelves frangeants seulement si bm-bmelt positif'
 write(num_rep_42,*)   '! meth_hcoup pour faire eventuellement varier Hcoup avec le climat'
-write(num_rep_42,*)   '! Hcoup_min=',Hcoup_min
+write(num_rep_42,*)   '! Hcoup_plateau=',Hcoup_plateau
+write(num_rep_42,*)   '! Hcoup_abysses=',Hcoup_abysses
+write(num_rep_42,*)   '! prof_plateau=',prof_plateau
+write(num_rep_42,*)   '! prof_abysses=',prof_abysses
 write(num_rep_42,*)
 
-Hcoup_max=Hcoup
+
+! afq -- coupure depend de la profondeur:
+!
+!  hcoup                       prof_abysses
+!   ^                               v 
+!   |                                _______ hcoup_abysses
+!   |                               /
+!   |                              /
+!   |                             /
+!   |                            /
+!   |      hcoup_plateau _______/
+!                               ^
+!                          prof_plateau     
+
+Hcoup(:,:) = min ( max(                                         &
+     (-(Bsoc0(:,:)-sealevel) - prof_plateau)/(prof_abysses-prof_plateau)    &
+     *(hcoup_abysses-hcoup_plateau)+hcoup_plateau               &
+     , hcoup_plateau), hcoup_abysses )
+
+if (meth_hcoup.eq.1) then
+   Hcoup(:,:)=coefbmshelf*Hcoup(:,:)
+   Hcoup(:,:)=min( max(Hcoup (:,:),Hcoup_plateau),Hcoup_abysses)
+else if (meth_hcoup.eq.2) then
+   Hcoup(:,:)=coefbmshelf*Hcoup(:,:)
+   Hcoup(:,:)=max(Hcoup(:,:),Hcoup_plateau)
+endif
+
 
 end subroutine init_calving
@@ -83 +119 @@
 bil_tot(:,:)=0.5*(bm(:,:)-bmelt(:,:))
 
-! coupure
-if (meth_hcoup.eq.0) then
-   Hcoup=Hcoup_max
-else if (meth_hcoup.eq.1) then
-   Hcoup=coefbmshelf*Hcoup_max
-   Hcoup=min(Hcoup,Hcoup_max)
-   Hcoup=max(Hcoup,Hcoup_min)
-
+Hcoup(:,:) = min ( max(                                         &
+     (-(Bsoc(:,:)-sealevel) - prof_plateau)/(prof_abysses-prof_plateau)    &
+     *(hcoup_abysses-hcoup_plateau)+hcoup_plateau               &
+     , hcoup_plateau), hcoup_abysses )
+
+if (meth_hcoup.eq.1) then
+   Hcoup(:,:)=coefbmshelf*Hcoup(:,:)
+   Hcoup(:,:)=min( max(Hcoup (:,:),Hcoup_plateau),Hcoup_abysses)
 else if (meth_hcoup.eq.2) then
-   Hcoup=coefbmshelf*Hcoup_max
-   Hcoup=max(Hcoup,Hcoup_min)
+   Hcoup(:,:)=coefbmshelf*Hcoup(:,:)
+   Hcoup(:,:)=max(Hcoup(:,:),Hcoup_plateau)
 endif
 
 ! hauteur au dessus de la coupure
-hmhc(:,:)=H(:,:)-Hcoup
+hmhc(:,:)=H(:,:)-Hcoup(:,:)
 
 ! coupure de l'ice shelf
@@ -108 +144 @@
 
 
-ifext:  if((flot(i,j)).and.(h(i,j).le.hcoup))  then
+ifext:  if((flot(i,j)).and.(h(i,j).le.hcoup(i,j)))  then
 ! ifext: pour les noeuds flottants avec h < hcoup