Changeset 662 for trunk/AGRIF/AGRIF_FILES/modmask.F

Timestamp:

2007-05-25T17:58:52+02:00 (17 years ago)

Author:

opalod

Message:

RB: update Agrif internal routines with a new update scheme and performance improvment

File:

• trunk/AGRIF/AGRIF_FILES/modmask.F (modified) (10 diffs)

trunk/AGRIF/AGRIF_FILES/modmask.F

@@ -202 +202 @@
 C     Local scalar
       INTEGER                  :: i,ii,iii,jj,kk,ll,mm,nn 
-      INTEGER,DIMENSION(nbdim) :: imin,imax
+      INTEGER,DIMENSION(nbdim) :: imin,imax,idecal
       INTEGER                  :: Nbvals
       REAL                     :: Res
@@ -235 +235 @@
                   imin(iii) = max(indic(iii) - i,ppbtab(iii))
                   imax(iii) = min(indic(iii) + i,ppetab(iii))
+                  if (firsttest) then               
+                  if (indic(iii).GT.ppbtab(iii)) then
+
+                     idecal = indic
+                     idecal(iii) = idecal(iii)-1
+                     SELECT CASE(nbdim)
+                     CASE (1)
+                        if (tempP%var%array1(idecal(1)
+     &                            ) == Agrif_SpecialValue) then
+                           imin(iii) = imax(iii)
+                        endif                     
+                     CASE (2)
+                        if (tempP%var%array2(idecal(1),
+     &            idecal(2)) == Agrif_SpecialValue) then
+                           imin(iii) = imax(iii)
+                        endif
+                     CASE (3)
+                        if (tempP%var%array3(idecal(1),
+     &            idecal(2),idecal(3)) 
+     &                               == Agrif_SpecialValue) then
+                           imin(iii) = imax(iii)
+                        endif  
+                     CASE (4)
+                        if (tempP%var%array4(idecal(1),
+     &            idecal(2),idecal(3),idecal(4)) 
+     &                               == Agrif_SpecialValue) then
+                           imin(iii) = imax(iii)
+                        endif     
+                     CASE (5)
+                        if (tempP%var%array5(idecal(1),
+     &            idecal(2),idecal(3),idecal(4),idecal(5)) 
+     &                               == Agrif_SpecialValue) then
+                           imin(iii) = imax(iii)
+                        endif  
+                     CASE (6)
+                        if (tempP%var%array6(idecal(1),
+     &            idecal(2),idecal(3),idecal(4),idecal(5),idecal(6)) 
+     &                               == Agrif_SpecialValue) then
+                           imin(iii) = imax(iii)
+                        endif                                                                                           
+                     END SELECT
+                  endif
+                  endif
                endif            
             enddo
@@ -241 +284 @@
             Nbvals = 0
 C
-            if ( nbdim .EQ. 1 ) then
+            SELECT CASE(nbdim)
+            CASE (1)
                do ii = imin(1),imax(1)
                     ValParent = parent%var%array1(ii)
@@ -249 +293 @@
                     endif 
                enddo
-            endif
-C
-            if ( nbdim .EQ. 2 ) then
+C
+            CASE (2)
                do jj = imin(2),imax(2)
                do ii = imin(1),imax(1)
@@ -261 +304 @@
                enddo  
                enddo
-            endif
-C
-            if ( nbdim .EQ. 3 ) then
+               
+            CASE (3)
                do kk = imin(3),imax(3)
                do jj = imin(2),imax(2)
@@ -275 +317 @@
                   enddo  
                enddo
-            endif
-C
-            if ( nbdim .EQ. 4 ) then
+
+            CASE (4)
                do ll = imin(4),imax(4)
                do kk = imin(3),imax(3)
@@ -291 +332 @@
                   enddo  
                enddo
-            endif
-C
-            if ( nbdim .EQ. 5 ) then
+
+            CASE (5)
                do mm = imin(5),imax(5)
                do ll = imin(4),imax(4)
@@ -309 +349 @@
                   enddo  
                enddo
-            endif
-C
-            if ( nbdim .EQ. 6 ) then
+
+            CASE (6)
                do nn = imin(6),imax(6)
                do mm = imin(5),imax(5)
@@ -329 +368 @@
                   enddo  
                enddo
-            endif
+
+            END SELECT
 C
 C
@@ -336 +376 @@
               if (firsttest) then
                    firsttest = .FALSE.
+                   i=1
                    cycle
               endif
-              if ( nbdim .EQ. 1 ) tempP%var%array1(indic(1)) 
+            SELECT CASE(nbdim)
+            CASE (1)              
+              tempP%var%array1(indic(1)) 
      &           = Res/Nbvals
-              if ( nbdim .EQ. 2 ) tempP%var%array2(indic(1),
+            CASE (2)
+              tempP%var%array2(indic(1),
      &                            indic(2)) = Res/Nbvals
-              if ( nbdim .EQ. 3 ) tempP%var%array3(indic(1),
+            CASE (3)
+              tempP%var%array3(indic(1),
      &                            indic(2),indic(3)) = Res/Nbvals
-              if ( nbdim .EQ. 4 ) tempP%var%array4(indic(1),
+            CASE (4)
+              tempP%var%array4(indic(1),
      &                            indic(2),indic(3),indic(4))
      &                = Res/Nbvals
-              if ( nbdim .EQ. 5 ) tempP%var%array5(indic(1),
+            CASE (5)
+              tempP%var%array5(indic(1),
      &                            indic(2),indic(3),indic(4),
      &                   indic(5)) = Res/Nbvals
-              if ( nbdim .EQ. 6 ) tempP%var%array6(indic(1),
+            CASE (6)
+              tempP%var%array6(indic(1),
      &                            indic(2),indic(3),indic(4),
      &                           indic(5),indic(6)) = Res/Nbvals
+            END SELECT
               exit
             else