Changeset 572

Timestamp:

2006-11-13T15:16:21+01:00 (17 years ago)

Author:

opalod

Message:

RB: ensure correct vectorization in modbc.F in Agrif routines

File:

• trunk/AGRIF/AGRIF_FILES/modbc.F (modified) (3 diffs)

trunk/AGRIF/AGRIF_FILES/modbc.F

@@ -746 +746 @@
       INTEGER                      :: nb,ndir,n,sizetab
       REAL, DIMENSION(:), Allocatable :: tab ! Array used for the interpolation
+      INTEGER, DIMENSION(nbdim)  :: ztab ! Array used for the interpolation
       REAL    :: c1t,c2t               ! Coefficients for the time interpolation
                                        !    (c2t=1-c1t) 
@@ -928 +929 @@
 C      
               sizetab = 1
-C
-              do i = 1,nbdim
+              ztab(1) = ptres2(1,2,ndir,nb)-ptres2(1,1,ndir,nb)+1
+C
+              do i = 2,nbdim
 C          
-                sizetab = sizetab 
+                ztab(i) = ztab(i-1) 
      &              * (ptres2(i,2,ndir,nb)-ptres2(i,1,ndir,nb)+1)
 C      
               enddo
-C              
+              sizetab=ztab(nbdim)  
               allocate(tab(sizetab))
 C
@@ -991 +993 @@
      &        .AND. loctab_child(nb) /= -3) then
 C      
-              sizetab=1   
-C
-              do i = 1,nbdim
+              sizetab=1  
+              ztab(1) = ptres2(1,2,ndir,nb)-ptres2(1,1,ndir,nb)+1
+              do i = 2,nbdim
 C          
-                sizetab = sizetab
+                ztab(i) = ztab(i-1)
      &             * (ptres2(i,2,ndir,nb)-ptres2(i,1,ndir,nb)+1)
 C      
               enddo
 C              
+              sizetab = ztab(nbdim) 
               allocate(tab(sizetab))
 C