Changeset 2301 in lmdz_wrf for trunk/tools

Timestamp:

Jan 30, 2019, 2:26:02 AM (7 years ago)

Author:

lfita

Message:

Getting there ?

File:

• trunk/tools/module_scientific.f90 (modified) (3 diffs)

trunk/tools/module_scientific.f90

@@ -6255 +6255 @@
   END SUBROUTINE coincident_list_2Dcoords
 
-  SUBROUTINE coincident_gridsin2D(dxA, dyA, dxyA, NpointsA, pointsA, dxB, dyB, dxyB, NpointsB,        &
+  SUBROUTINE coincident_gridsin2D_old(dxA, dyA, dxyA, NpointsA, pointsA, dxB, dyB, dxyB, NpointsB,        &
     pointsB, Npoints, points, inpointsA, inpointsB)
   ! Subroutine to determine which lists of 2D gridsin points of an A list are also found in a B list
@@ -6292 +6292 @@
 ! inpointsB: coincident points list B
 
-    fname = 'coincident_gridsin2D'
+    fname = 'coincident_gridsin2D_old'
 
     Npoints = 0
@@ -6327 +6327 @@
     END DO
 
+  END SUBROUTINE coincident_gridsin2D_old
+
+  SUBROUTINE coincident_gridsin2D(dxA, dyA, dxyA, NpointsA, pointsA, dxB, dyB, dxyB, NpointsB,        &
+    pointsB, Npoints, points, inpointsA, inpointsB)
+  ! Subroutine to determine which lists of 2D gridsin points of an A list are also found in a B list
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: dxA, dyA, dxyA
+    INTEGER, INTENT(in)                                  :: dxB, dyB, dxyB
+    INTEGER, DIMENSION(dxA, dyA), INTENT(in)             :: NpointsA
+    INTEGER, DIMENSION(dxB, dyB), INTENT(in)             :: NpointsB
+    INTEGER, DIMENSION(dxA, dyA, dxyA, 2), INTENT(in)    :: pointsA
+    INTEGER, DIMENSION(dxB, dyB, dxyB, 2), INTENT(in)    :: pointsB
+    INTEGER, DIMENSION(dxA, dyA, dxB, dyB), INTENT(out)  :: Npoints
+    INTEGER, DIMENSION(dxA, dyA, dxB, dyB, dxyA, 2),                                                  &
+      INTENT(out)                                        :: points
+    INTEGER, DIMENSION(dxA, dyA, dxyA), INTENT(out)      :: inpointsA
+    INTEGER, DIMENSION(dxB, dyB, dxyA), INTENT(out)      :: inpointsB
+
+    ! Local
+    INTEGER                                              :: ixA, iyA, ixB, iyB, iv, iv1, iv2
+    INTEGER                                              :: NA, NB
+    INTEGER, DIMENSION(dxyA)                             :: ptsA, ptsB
+    INTEGER, DIMENSION(dxyA, 2)                          :: pts
+
+
+!!!!!!! Variables
+! dxA, dyA: 2D shape of the list A
+! NpointsA: 2D Number of points of the list A
+! pointsA: points of the list A
+! dxB, dyB: 2D shape of the list B
+! NpointsB: 2D Number of points of the list B
+! pointsB: points of the list B
+! Npoints: Number of coincident points
+! points: coincident points
+! inpointsA: coincident points list A
+! inpointsB: coincident points list B
+
+    fname = 'coincident_gridsin2D'
+
+    Npoints = 0
+    points = 0
+    inpointsA = 0
+    inpointsB = 0
+
+    DO ixA=1, dxA
+      DO iyA=1, dyA
+        NA = NpointsA(ixA,iyA)
+        DO ixB=1, dxB
+          DO iyB=1, dyB
+            NB = NpointsB(ixB,iyB)
+            iv = 0
+            DO iv1=1, NA
+              DO iv2=1, NB
+                IF ( (pointsA(ixA,iyA,iv1,1) == pointsB(ixB,iyB,iv2,1)) .AND.                         &
+                  (pointsA(ixA,iyA,iv1,2) == pointsB(ixB,iyB,iv2,2)) ) THEN
+                  iv = iv + 1
+                  points(ixA,iyA,ixB,iyB,iv,1) = pointsA(ixA,iyA,iv1,1)
+                  points(ixA,iyA,ixB,iyB,iv,2) = pointsA(ixA,iyA,iv1,2)
+                  inpointsA(ixA,iyA,iv) = iv1
+                  inpointsB(ixB,iyB,iv) = iv2
+                END IF
+              END DO
+            END DO
+            Npoints(ixA,iyA,ixB,iyB) = iv
+          END DO
+        END DO
+      END DO
+    END DO
+
   END SUBROUTINE coincident_gridsin2D