Changeset 2671 for vendors/AGRIF/current/AGRIF_FILES/modinterpbasic.F

Timestamp:

2011-03-08T15:08:49+01:00 (13 years ago)

Author:

rblod

Message:

Load working_directory into vendors/AGRIF/current.

File:

• vendors/AGRIF/current/AGRIF_FILES/modinterpbasic.F (modified) (11 diffs)

vendors/AGRIF/current/AGRIF_FILES/modinterpbasic.F

@@ -1 +1 @@
 !
-! $Id: modinterpbasic.F 898 2008-04-22 15:35:20Z rblod $
+! $Id: modinterpbasic.F 779 2007-12-22 17:04:17Z rblod $
 !
 C     AGRIF (Adaptive Grid Refinement In Fortran)
@@ -113 +113 @@
            locind_parent_left = locind_parent_left + 1.
            globind_parent_right = globind_parent_right + 1.
+           ypos2 = ypos*invds+(i-1)*invds2           
         endif
         
         diff=(globind_parent_right - ypos2)
-        
         y(i) = (diff*x(locind_parent_left)
      &        + (1.-diff)*x(locind_parent_left+1))
@@ -170 +170 @@
       Real, Dimension(:,:), Allocatable :: coeffparent_tmp
       REAL    :: ypos,globind_parent_left,globind_parent_right
-      REAL    :: invds, invds2
+      REAL    :: invds, invds2, invds3
       REAL :: ypos2,diff
 C
@@ -189 +189 @@
       invds = 1./ds_parent
       invds2 = ds_child/ds_parent
+      invds3 = 0.5/real(coeffraf)
 
       ypos2 = ypos*invds
@@ -216 +217 @@
            locind_parent_left = locind_parent_left + 1
            globind_parent_right = globind_parent_right + 1.
+           ypos2 = ypos*invds+(i-1)*invds2
         endif
-        
+
         diff=(globind_parent_right - ypos2)
+        diff = invds3*nint(2*coeffraf*diff)
         indparent(i,dir) = locind_parent_left
+
         coeffparent(i,dir) = diff
         
@@ -238 +242 @@
 C      
        indparent(nc,dir) = locind_parent_left
-       
-         coeffparent(nc,dir) = (globind_parent_left + ds_parent - ypos)
+       diff = (globind_parent_left + ds_parent - ypos)
      &       * invds
-      endif                                           
+        diff = invds3*nint(2*coeffraf*diff)
+         coeffparent(nc,dir) = diff
+        endif                                           
 
       do i=2, np2
@@ -280 +285 @@
       INTEGER :: i,coeffraf,locind_parent_left
       REAL    :: ypos,globind_parent_left,globind_parent_right
-      REAL    :: invds, invds2
+      REAL    :: invds, invds2, invds3
       REAL :: ypos2,diff
 C
@@ -305 +310 @@
       invds = 1./ds_parent
       invds2 = ds_child/ds_parent
+      invds3 = 0.5/real(coeffraf)      
       
       ypos2 = ypos*invds
@@ -323 +329 @@
            locind_parent_left = locind_parent_left + 1
            globind_parent_right = globind_parent_right + 1.
+           ypos2 = ypos*invds+(i-1)*invds2           
         endif
         
         diff=(globind_parent_right - ypos2)
+        
+        diff = invds3*nint(2*coeffraf*diff)
+                
         indparent(i,1) = locind_parent_left
+
         coeffparent(i,1) = diff
         ypos2 = ypos2 + invds2
@@ -344 +355 @@
 C      
        indparent(nc,1) = locind_parent_left
-       
-         coeffparent(nc,1) = (globind_parent_left + ds_parent - ypos)
+ 
+       diff = (globind_parent_left + ds_parent - ypos)
      &       * invds
-      endif                                          
+        diff = invds3*nint(2*coeffraf*diff)
+         coeffparent(nc,1) = diff
+         endif                                          
 C           
       Return
@@ -462 +475 @@
      &                        + (locind_parent_left - 1)*ds_parent 
      
-        deltax = invdsparent*(ypos-globind_parent_left)
+C        deltax = invdsparent*(ypos-globind_parent_left)
+        deltax = nint(coeffraf*deltax)/real(coeffraf)
+        
         ypos = ypos + ds_child
          if (abs(deltax).LE.0.0001) then