Changeset 2671 for vendors/AGRIF/current/AGRIF_FILES/modupdatebasic.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/modupdatebasic.F (modified) (3 diffs)

vendors/AGRIF/current/AGRIF_FILES/modupdatebasic.F

@@ -571 +571 @@
       INTEGER :: i1,i2
       REAL :: invsumweight
-      REAL :: weights(-(coeffraf-1):coeffraf-1)
+      REAL :: weights(-(coeffraf):coeffraf)
       
 C
@@ -594 +594 @@
        ENDIF
 C
-             
-       it1 = -(coeffraf-1)
-       i1 = -(coeffraf-1)+locind_child_left+diffmod
-       i2 = 2*coeffraf - 2
        
-      invsumweight=1./coeffraf**2
-      do i=-(coeffraf-1),0
-        weights(i) = invsumweight*(coeffraf + i)
-      enddo
-      do i=1,coeffraf-1
-        weights(i) = invsumweight*(coeffraf - i)
-      enddo
+      if (diffmod == 1) THEN
+        invsumweight=1./(2.*coeffraf**2)
+        do i=-coeffraf,-1
+          weights(i) = invsumweight*(2*(coeffraf+i)+1)
+        enddo
+        do i=0,coeffraf-1
+          weights(i)=weights(-(i+1))
+        enddo
+        it1 = -coeffraf
+        i1 = -(coeffraf-1)+locind_child_left
+        i2 = 2*coeffraf - 1
+      else
+        invsumweight=1./coeffraf**2
+        do i=-(coeffraf-1),0
+          weights(i) = invsumweight*(coeffraf + i)
+        enddo
+        do i=1,coeffraf-1
+          weights(i) = invsumweight*(coeffraf - i)
+        enddo
+        it1 = -(coeffraf-1)
+        i1 = -(coeffraf-1)+locind_child_left
+        i2 = 2*coeffraf - 2
+      endif
 
       sumweight = 0                    
@@ -638 +650 @@
       enddo    
 C
+
       Return 
 C