Changeset 98 for trunk/SOURCES

Timestamp:

01/19/17 15:21:53 (7 years ago)

Author:

dumas

Message:

Schoof update : Schoof only if nt > 15000

File:

• trunk/SOURCES/New-remplimat/diagno-L2_mod.f90 (modified) (5 diffs)

trunk/SOURCES/New-remplimat/diagno-L2_mod.f90

@@ -28 +28 @@
 
 integer :: ifail_diagno  ! pour recuperation d'erreur
+integer :: ifail_diagno_ramollo  ! pour recuperation d'erreur shelf ramollo
 integer :: iplus1,jplus1
 integer ::  ctvisco,iumax,jumax
@@ -140 +141 @@
 !~      enddo   
   
-  if (Schoof.eq.1) then ! flux grounding line Schoof
+  if (Schoof.eq.1.and.nt.GT.15000) then ! flux grounding line Schoof
                 call interpol_glflux ! calcul flux GL + interpolation sur voisins
         endif   
@@ -180 +181 @@
   !nyd2=60
 
-
   call rempli_L2(nxd1,nxd2,nyd1,nyd2,uxbar(nxd1:nxd2,nyd1:nyd2),uybar(nxd1:nxd2,nyd1:nyd2), & 
        uxb1(nxd1:nxd2,nyd1:nyd2),uyb1(nxd1:nxd2,nyd1:nyd2),  &
        imx_diag(nxd1:nxd2,nyd1:nyd2),imy_diag(nxd1:nxd2,nyd1:nyd2),ifail_diagno)
        
-  if (Schoof.eq.1) then ! flux grounding line Schoof avec calcul de la back force par shelf ramollo
+  if (Schoof.eq.1.and.nt.GT.15000) then ! flux grounding line Schoof avec calcul de la back force par shelf ramollo
     pvi_keep(:,:)=pvi(:,:)       
-                where (flot(:,:))
-                                pvi(:,:)=pvimin
-                endwhere                        
-            
+                where (flot(:,:).and.H(:,:).GT.2.)
+                                pvi(:,:)=1.e5    
+!                               pvi(:,:)=pvimin
+                endwhere
+                                        
                 call rempli_L2(nxd1,nxd2,nyd1,nyd2,uxbar(nxd1:nxd2,nyd1:nyd2),uybar(nxd1:nxd2,nyd1:nyd2), & 
        uxb1ramollo(nxd1:nxd2,nyd1:nyd2),uyb1ramollo(nxd1:nxd2,nyd1:nyd2),  &
-       imx_diag(nxd1:nxd2,nyd1:nyd2),imy_diag(nxd1:nxd2,nyd1:nyd2),ifail_diagno)
+       imx_diag(nxd1:nxd2,nyd1:nyd2),imy_diag(nxd1:nxd2,nyd1:nyd2),ifail_diagno_ramollo)
+       
     pvi(:,:)=pvi_keep(:,:) 
+    
     where ((uxb1ramollo(:,:)**2 + uyb1ramollo(:,:)**2).GT.0.)
                         back_force(:,:) = sqrt(uxb1(:,:)**2+ uyb1(:,:)**2) / sqrt(uxb1ramollo(:,:)**2 + uyb1ramollo(:,:)**2)
@@ -200 +203 @@
                         back_force(:,:)=1.
                 endwhere
-  endif
-  
+                
+                if (ifail_diagno_ramollo.gt.0) then
+                        write(6,*) ' Probleme resolution systeme L2. ramollo ifail=',ifail_diagno_ramollo
+                        STOP
+                endif
 !~   do j=1,ny
 !~              do i=1,nx
@@ -214 +220 @@
 !~   print*,'apres calcul rempli_L2'
 !~   read(*,*)      
+  endif
+  
 
   ! Dans rempli_L2