Changeset 418 for branches/GRISLIv3/SOURCES

Timestamp:

04/24/23 11:21:51 (15 months ago)

Author:

dumas

Message:

use only in neffect

File:

• branches/GRISLIv3/SOURCES/neffect-0.4.f90 (modified) (1 diff)

branches/GRISLIv3/SOURCES/neffect-0.4.f90

@@ -12 +12 @@
 !! @note    - use printtable
 !<
-      subroutine Neffect()
+subroutine Neffect
 
-!     ===============================================================
-!     *    Pression effective                 Cat   18 janvier 2005
-!     ===============================================================
+  !     ===============================================================
+  !     *    Pression effective                 Cat   18 janvier 2005
+  !     ===============================================================
 
-       USE module3D_phy
-       use module_choix
-       use printtable
+  use param_phy_mod, only: ro,row,g
+  use runparam, only: itracebug
+  use module3D_phy, only: nx,ny,flotmx,flotmy,coefmxbmelt,coefmybmelt,hwater,neffmx,neffmy,hmx,hmy,hwatermax
+  use module_choix, only: eaubasale
 
-      implicit none
+  implicit none
 
-!      real,dimension(nx,ny) :: table_interm
-!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
+  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
 
-if (itracebug.eq.1)  call tracebug(' Entree dans routine neffect')
+  if (itracebug.eq.1)  call tracebug(' Entree dans routine neffect')
 
 
 
 
-! le calcul de pente est fait dans diffusiv
-
-    
-! calcul de la hauteur d'eau dans eaubasale
+  ! le calcul de pente est fait dans diffusiv
 
 
-      call  EAUBASALE
+  ! calcul de la hauteur d'eau dans eaubasale
 
-if (itracebug.eq.1)  call tracebug(' Apres eau basale')
-! calcul de la presion effective Neffmy
 
-        do J=2,NY
-          do I=2,NX-1
-pose_y :    if (.not.flotmy(i,j)) then
+  call eaubasale
 
-!     coefmybmelt est la hauteur d'eau au point considere
-!     moyenne pondérée des points voisins.
-!     coefmybmelt est utilise dans la vitesse de glissement 
+  if (itracebug.eq.1)  call tracebug(' Apres eau basale')
+  ! calcul de la presion effective Neffmy
 
-!    ancienne version : moyenne large
-!             COEFMYBMELT(I,J)=(HWATER(I,J)+HWATER(I,J-1))/4.+
-!     &               ((HWATER(I+1,J)+HWATER(I+1,J-1))+
-!     &                (HWATER(I-1,J)+HWATER(I-1,J-1)))/8.
+  do j=2,ny
+     do i=2,nx-1
+        pose_y :    if (.not.flotmy(i,j)) then
 
-!    nouvelle version, moyenne sur les deux plus proches voisins pour mieux 
-!    localiser les fleuves de glace
+           !     coefmybmelt est la hauteur d'eau au point considere
+           !     moyenne pondérée des points voisins.
+           !     coefmybmelt est utilise dans la vitesse de glissement 
 
-             COEFMYBMELT(I,J)=(HWATER(I,J)+HWATER(I,J-1))*0.5
-             Neffmy(i,j)=RO*G*HMY(I,J)-ROW*G*COEFMYBMELT(I,J)
+           !    ancienne version : moyenne large
+           !             COEFMYBMELT(I,J)=(HWATER(I,J)+HWATER(I,J-1))/4.+
+           !     &               ((HWATER(I+1,J)+HWATER(I+1,J-1))+
+           !     &                (HWATER(I-1,J)+HWATER(I-1,J-1)))/8.
+
+           !    nouvelle version, moyenne sur les deux plus proches voisins pour mieux 
+           !    localiser les fleuves de glace
+
+           coefmybmelt(i,j)=(hwater(i,j)+HWATER(i,j-1))*0.5
+           Neffmy(i,j)=ro*g*hmy(i,j)-row*g*coefmybmelt(i,j)
 !!!!!        Neffmy(i,j)=RO*G*HMY(I,J)-pwater(i,j)
 
-!             la pression effective mini est Nefmin (en Pa) pour la glace posée 
-             Neffmy(i,j)=max(Neffmy(i,j),Nefmin)
+           !             la pression effective mini est Nefmin (en Pa) pour la glace posée 
+           Neffmy(i,j)=max(Neffmy(i,j),Nefmin)
 
-!     coefmybmelt est utilise dans la vitesse de glissement 
-!     et pour cela est limité par hwatermax
-             COEFMYBMELT(i,j)=min(hwatermax,COEFMYBMELT(i,j))
+           !     coefmybmelt est utilise dans la vitesse de glissement 
+           !     et pour cela est limité par hwatermax
+           coefmybmelt(i,j)=min(hwatermax,coefmybmelt(i,j))
 
 
 
-! on calcule tobmy dans diagnoshelf apres etre passe dans remplimat
-! pour utiliser les nouvelles vitesses.
+           ! on calcule tobmy dans diagnoshelf apres etre passe dans remplimat
+           ! pour utiliser les nouvelles vitesses.
 
-          else
-!           points completement flottants
-              Neffmy(i,j)=0.
+        else
+           !           points completement flottants
+           Neffmy(i,j)=0.
 
-           endif  pose_y
-        end do
-      end do
+        endif  pose_y
+     end do
+  end do
 
-! calcul de la presion effective Neffmx
+  ! calcul de la presion effective Neffmx
 
-        do J=2,NY-1
-          do I=2,NX
-pose_x :       if (.not.flotmx(i,j)) then
+  do j=2,Ny-1
+     do i=2,nx
+        pose_x :       if (.not.flotmx(i,j)) then
 
-!     coefmybmelt est la hauteur d'eau au point considere
-!     moyenne pondérée des points voisins.
-!     coefmybmelt est utilise dans la vitesse de glissement 
+           !     coefmybmelt est la hauteur d'eau au point considere
+           !     moyenne pondérée des points voisins.
+           !     coefmybmelt est utilise dans la vitesse de glissement 
 
-!    ancienne version : moyenne large
-!             COEFMYBMELT(I,J)=(HWATER(I,J)+HWATER(I-1,J))/4.+
-!     &               ((HWATER(I,J+1)+HWATER(I-1,J+1))+
-!     &                (HWATER(I,J-1)+HWATER(I-1,J-1)))/8.
+           !    ancienne version : moyenne large
+           !             COEFMYBMELT(I,J)=(HWATER(I,J)+HWATER(I-1,J))/4.+
+           !     &               ((HWATER(I,J+1)+HWATER(I-1,J+1))+
+           !     &                (HWATER(I,J-1)+HWATER(I-1,J-1)))/8.
 
-!    nouvelle version, moyenne sur les deux plus proches voisins pour mieux 
-!    localiser les fleuves de glace
+           !    nouvelle version, moyenne sur les deux plus proches voisins pour mieux 
+           !    localiser les fleuves de glace
 
-             COEFMXBMELT(I,J)=(HWATER(I,J)+HWATER(I-1,J))*0.5
-             Neffmx(i,j)=RO*G*HMX(I,J)-ROW*G*COEFMXBMELT(I,J)
+           coefmxbmelt(i,j)=(hwater(i,j)+hwater(i-1,j))*0.5
+           Neffmx(i,j)=ro*g*hmx(I,J)-row*g*coefmxbmelt(i,j)
 !!!!!        Neffmx(i,j)=RO*G*HMX(I,J)-pwater(i,j)
 
-!             la pression effective mini est 1 Pascal pour la glace posée 
-             Neffmx(i,j)=max(Neffmx(i,j),Nefmin)
+           !             la pression effective mini est 1 Pascal pour la glace posée 
+           Neffmx(i,j)=max(Neffmx(i,j),Nefmin)
 
-!     coefmybmelt est utilise dans la vitesse de glissement 
-!     et pour cela est limité par hwatermax
-             COEFMXBMELT(i,j)=min(hwatermax,COEFMXBMELT(i,j))
+           !     coefmybmelt est utilise dans la vitesse de glissement 
+           !     et pour cela est limité par hwatermax
+           coefmxbmelt(i,j)=min(hwatermax,coefmxbmelt(i,j))
 
 
 
 
-! on calcule tobmy dans diagnoshelf apres etre passe dans remplimat
-! pour utiliser les nouvelles vitesses.
+           ! on calcule tobmy dans diagnoshelf apres etre passe dans remplimat
+           ! pour utiliser les nouvelles vitesses.
 
-           else
-!           points completement flottants
-              Neffmx(i,j)=0.
- 
-           endif  pose_x 
-        end do
-      end do
+        else
+           !           points completement flottants
+           Neffmx(i,j)=0.
 
-! nom_table='pwater'
-! table_interm=ROW*G*COEFMXBMELT/10e5      
-! call printtable_r(table_interm,nom_table)         
-   
- if (itracebug.eq.1)  call tracebug(' Sortie routine neffect')
-      return
-    end subroutine Neffect
+        endif  pose_x
+     end do
+  end do
+
+  ! nom_table='pwater'
+  ! table_interm=ROW*G*COEFMXBMELT/10e5      
+  ! call printtable_r(table_interm,nom_table)         
+
+  if (itracebug.eq.1)  call tracebug(' Sortie routine neffect')
+  return
+end subroutine Neffect