Changeset 221 for trunk/SOURCES/furst_schoof_mod.f90

Timestamp:

12/17/18 15:13:04 (6 years ago)

Author:

aquiquet

Message:

Bug correction (minor impact) in furst_schoof + velocity limitation based on what we declare in fichier param

File:

• trunk/SOURCES/furst_schoof_mod.f90 (modified) (2 diffs)

trunk/SOURCES/furst_schoof_mod.f90

@@ -464 +464 @@
   ! afq -- if we discard the points with multiple fluxes coming, uncomment:
   where (countx(:,:).ge.2)
-     uxbartemp(:,:)=uxbar(i,j)*countx(:,:)
+     uxbartemp(:,:)=uxbar(:,:)*countx(:,:)
      imx_diag(:,:) = imx_diag_in(:,:)
   end where
@@ -473 +473 @@
   ! afq --
   
-  where (countx(:,:).ne.0) uxbar(:,:)= max(min(uxbartemp(:,:)/countx(:,:),1000.),-1000.)
-  where (county(:,:).ne.0) uybar(:,:)= max(min(uybartemp(:,:)/county(:,:),1000.),-1000.)
+  where (countx(:,:).ne.0) uxbar(:,:)= max(min(uxbartemp(:,:)/countx(:,:),V_limit),-V_limit)
+  where (county(:,:).ne.0) uybar(:,:)= max(min(uybartemp(:,:)/county(:,:),V_limit),-V_limit)
            
   do j= 3,ny-3