Changeset 462

Timestamp:

01/12/24 16:28:09 (6 months ago)

Author:

aquiquet

Message:

Cleaning branch: change the formulation of nefmin so that it follows PISM with a fraction of driving stress, at present: nefmin=0 still for compatibility with older simulations

File:

• branches/GRISLIv3/SOURCES/neffect-0.4.f90 (modified) (3 diffs)

branches/GRISLIv3/SOURCES/neffect-0.4.f90

@@ -27 +27 @@
   integer :: i,j
   !cdc      real :: Nefmin=1.e5                  ! Pression effective minimum (~ 10 m de glace)
-  real :: Nefmin=0.  !cdc 22/01/2019 on autorise une pression effective nulle
+  !afq  - no longer used - real :: Nefmin=0.  !cdc 22/01/2019 on autorise une pression effective nulle
+  real :: nefminfrac = 0. ! put 0.04 to have something similar to PISM
 
   if (itracebug.eq.1)  call tracebug(' Entree dans routine neffect')
@@ -66 +67 @@
 
            !             la pression effective mini est Nefmin (en Pa) pour la glace posée 
-           Neffmy(i,j)=max(Neffmy(i,j),Nefmin)
+! afq --           Neffmy(i,j)=max(Neffmy(i,j),Nefmin)
+           Neffmy(i,j)=max(Neffmy(i,j),nefminfrac*RO*G*HMY(I,J)) !PISM: nefminfrac=4%
 
            !     coefmybmelt est utilise dans la vitesse de glissement 
@@ -108 +110 @@
 
            !             la pression effective mini est 1 Pascal pour la glace posée 
-           Neffmx(i,j)=max(Neffmx(i,j),Nefmin)
+! afq --           Neffmx(i,j)=max(Neffmx(i,j),Nefmin)
+           Neffmx(i,j)=max(Neffmx(i,j),nefminfrac*RO*G*HMX(I,J)) !PISM: nefminfrac=4%
 
            !     coefmybmelt est utilise dans la vitesse de glissement