Changeset 482 for branches/GRISLIv3/SOURCES/New-remplimat/diagno-L2_mod.f90

Timestamp:

02/09/24 17:08:19 (5 months ago)

Author:

aquiquet

Message:

Cleaning branch: some routines with explicit arguments in New-remplimat

File:

• branches/GRISLIv3/SOURCES/New-remplimat/diagno-L2_mod.f90 (modified) (12 diffs)

branches/GRISLIv3/SOURCES/New-remplimat/diagno-L2_mod.f90

@@ -16 +16 @@
 real, dimension(nx,ny) :: uyb1ramollo
 real, dimension(nx,ny) :: pvi_keep
-
-!cdc transfere dans module3d pour compatibilite avec furst_schoof_mod
-!cdc  integer, dimension(nx,ny) :: imx_diag
-!cdc  integer, dimension(nx,ny) :: imy_diag
 
 integer :: nneigh=max(int(400000./dx),1) !nb. points around the grline kept for back force comp. in case we reduce the extent
@@ -151 +147 @@
 !$OMP END WORKSHARE
 !$OMP END PARALLEL
-  ! avec couplage thermomecanique
-!  write(166,*) ' apres call calc_pvi',itour_pvi
 
 else
@@ -158 +152 @@
 end if
 
-  call imx_imy_nx_ny         ! pour rempli_L2 : calcule les masques imx et imy qui 
+  call imx_imy_nx_ny(imx_diag,imy_diag,flgzmx,flgzmy)         ! pour rempli_L2 : calcule les masques imx et imy qui 
 
   diagno_grline=sum(gr_line(2:(nx-1),2:(ny-1)))
   if ((Schoof.ge.1).and.(diagno_grline.gt.0)) then
-       call imx_imy_nx_ny_reduce(1) !afq: provides imx_reduce and imy_reduce
+       call imx_imy_nx_ny_reduce(1,flot,imx_diag,imy_diag,gr_line) !afq: provides imx_reduce and imy_reduce
   endif
-
-!cdc debug Schoof !!!!!!!!!!!!  
-!~   do j=1,ny
-!~              do i=1,nx
-!~                      write(578,*) uxbar(i,j)
-!~                      write(579,*) uybar(i,j)
-!~              enddo
-!~      enddo   
-  
-  !if (Schoof.eq.1.and.nt.GT.15000) then ! flux grounding line Schoof
-! afq -- below:  if (Schoof.eq.1) then ! flux grounding line Schoof
-! afq -- below:     call interpol_glflux ! calcul flux GL + interpolation sur voisins
-! afq -- below:  endif
-
-!~       do j=1,ny
-!~              do i=1,nx
-!~                      write(588,*) uxbar(i,j)
-!~                      write(589,*) uybar(i,j)
-!~              enddo
-!~      enddo   
-!~      print*,'ecriteure termineee !!!!!!'
-!~      read(*,*)
 
   ! donnent les cas de conditions aux limites
@@ -202 +174 @@
   nyd1=1
   nyd2=ny
-
-  !call rempli_L2(1,nx,1,ny,uxbar,uybar,uxb1,uyb1,imx_diag,imy_diag,ifail_diagno)    
-  !nxd1=15
-  !nxd2=19
-  !nyd1=30
-  !nyd2=34
-
-  !nxd1=35
-  !nxd2=60
-  !nyd1=35
-  !nyd2=60
-
        
-  !if (Schoof.eq.1.and.nt.GT.15000) then ! flux grounding line Schoof avec calcul de la back force par shelf ramollo
   if ((Schoof.ge.1).and.(diagno_grline.gt.0)) then ! flux grounding line Schoof avec calcul de la back force par shelf ramollo
 
@@ -224 +183 @@
     pvi_keep(:,:)=pvi(:,:)       
     where (flot(:,:).and.H(:,:).GT.2.)
-!       pvi(:,:)=1.e5    
        pvi(:,:)=pvimin
     endwhere
@@ -250 +208 @@
 
     if (ifail_diagno_ramollo.gt.0) then
-!       write(6,*) ' Probleme resolution systeme L2. ramollo ifail=',ifail_diagno_ramollo
-!       STOP
        write(*,*) ' Probleme resolution systeme L2. ramollo ifail=',ifail_diagno_ramollo
        write(*,*) '          ... we go on anyway!'
     endif
-!~   do j=1,ny
-!~              do i=1,nx
-!~                if (sqrt(uxb1(i,j)**2+ uyb1(i,j)**2).gt.0..and..not.flot(i,j)) then
-!~                              write(1034,*) sqrt(uxb1(i,j)**2+ uyb1(i,j)**2) / sqrt(uxb1ramollo(i,j)**2 + uyb1ramollo(i,j)**2)
-!~                else
-!~                              write(1034,*) 1.
-!~                      endif
-!~              enddo           
-!~      enddo
-
-!~   print*,'apres calcul rempli_L2'
-!~   read(*,*)
     
     call interpol_glflux ! calcul flux GL + interpolation sur voisins
@@ -503 +447 @@
 debug_3D(:,:,27)=pvi(:,:)
 !$OMP END WORKSHARE
-! attention run 35
-!-------------------- 
-!!$if (time.gt.10.) then
-!!$   where (flot(:,:))
-!!$      pvi(:,:)=pvimin
-!!$   end where
-!!$end if
 
 !  calcul de la viscosite integree au milieu des mailles (pvm)
@@ -528 +465 @@
 !------------------------------------------------------------------
 
-subroutine imx_imy_nx_ny
-
-use module3d_phy, only: imx_diag,imy_diag,flgzmx,flgzmy 
-         
+subroutine imx_imy_nx_ny(imx_diagl,imy_diagl,flgzmxl,flgzmyl)
+
 implicit none
+
+logical,dimension(:,:), intent(in)    :: flgzmxl,flgzmyl
+integer,dimension(:,:), intent(inout) :: imx_diagl,imy_diagl
 
 ! definition des masques
@@ -555 +493 @@
 !$OMP PARALLEL
 !$OMP WORKSHARE
-imx_diag(:,:)=0
-imy_diag(:,:)=0
+imx_diagl(:,:)=0
+imy_diagl(:,:)=0
 
 ! a l'interieur du domaine
 !-------------------------
 
-where (flgzmx(:,:))
-   imx_diag(:,:)=2             ! shelf ou stream        
+where (flgzmxl(:,:))
+   imx_diagl(:,:)=2             ! shelf ou stream        
 elsewhere
-   imx_diag(:,:)=1             ! vitesse imposee
+   imx_diagl(:,:)=1             ! vitesse imposee
 end where
 
-where (flgzmy(:,:))      
-   imy_diag(:,:)=2             ! shelf ou stream
+where (flgzmyl(:,:))      
+   imy_diagl(:,:)=2             ! shelf ou stream
 elsewhere
-   imy_diag(:,:)=1             ! vitesse imposee
+   imy_diagl(:,:)=1             ! vitesse imposee
 end where
 
 ! bord sud 
-imx_diag(:,1)=-1
-imy_diag(:,2)=-1
+imx_diagl(:,1)=-1
+imy_diagl(:,2)=-1
 
 ! bord nord
-imx_diag(:,ny)=-3
-imy_diag(:,ny)=-3
+imx_diagl(:,ny)=-3
+imy_diagl(:,ny)=-3
 
 ! bord Est
-imx_diag(1,:)=0    ! hors domaine a cause des mailles alternees
-imx_diag(2,:)=-4
-imy_diag(1,:)=-4
+imx_diagl(1,:)=0    ! hors domaine a cause des mailles alternees
+imx_diagl(2,:)=-4
+imy_diagl(1,:)=-4
 
 ! bord West
-imx_diag(nx,:)=-2
-imy_diag(nx,:)=-2
+imx_diagl(nx,:)=-2
+imy_diagl(nx,:)=-2
 
 ! Coins
-imx_diag(2,1)=-41       ! SW
-imy_diag(1,2)=-41
-
-imx_diag(nx,1)=-12      ! SE
-imy_diag(nx,2)=-12
-
-imx_diag(nx,ny)=-23     ! NE
-imy_diag(nx,ny)=-23
-
-imx_diag(2,ny)=-34      ! NW
-imy_diag(1,ny)=-34
+imx_diagl(2,1)=-41       ! SW
+imy_diagl(1,2)=-41
+
+imx_diagl(nx,1)=-12      ! SE
+imy_diagl(nx,2)=-12
+
+imx_diagl(nx,ny)=-23     ! NE
+imy_diagl(nx,ny)=-23
+
+imx_diagl(2,ny)=-34      ! NW
+imy_diagl(1,ny)=-34
 
 ! hors domaine
-imx_diag(1,:)=0     ! hors domaine a cause des mailles alternees
-imy_diag(:,1)=0     ! hors domaine a cause des mailles alternees
+imx_diagl(1,:)=0     ! hors domaine a cause des mailles alternees
+imy_diagl(:,1)=0     ! hors domaine a cause des mailles alternees
 !$OMP END WORKSHARE
 !$OMP END PARALLEL
@@ -612 +550 @@
 
 
-subroutine imx_imy_nx_ny_reduce(choix)
-
-use module3d_phy, only: flot,imx_diag,imy_diag,gr_line 
-         
+subroutine imx_imy_nx_ny_reduce(choix,flot,imx_diag,imy_diag,gr_line)
+
 implicit none
 
@@ -621 +557 @@
 !       We simply compute the velocities around (nneigh) the grounding line
 
-  integer, intent(in) :: choix
+  integer,                 intent(in)    :: choix
+  logical, dimension(:,:), intent(in)    :: flot
+  integer, dimension(:,:), intent(in)    :: imx_diag,imy_diag,gr_line
+
 
   integer i,j,nvx,nvy
@@ -717 +656 @@
 end subroutine mismip_boundary_cond
 
-!end module diagno_L2_mod
 end module diagno_mod