Changeset 1747 in lmdz_wrf for trunk/tools

Timestamp:

Dec 17, 2017, 4:08:32 PM (8 years ago)

Author:

lfita

Message:

Adding:

• `fill3DI_2Dvec': Subroutine to fill a 3D integer matrix from a series of indices from a 2D matrix
• `fill3DR_2Dvec': Subroutine to fill a 3D float matrix from a series of indices from a 2D matrix

File:

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

trunk/tools/module_scientific.f90

@@ -15 +15 @@
 ! clean_polygons: Subroutine to clean polygons from non-used paths, polygons only left as path since they are inner path of a hole
 ! FindMinimumR_K*: Function returns the location of the minimum in the section between Start and End.
+! fill3DI_2Dvec: Subroutine to fill a 3D integer matrix from a series of indices from a 2D matrix
+! fill3DR_2Dvec: Subroutine to fill a 3D float matrix from a series of indices from a 2D matrix
 ! gridpoints_InsidePolygon: Subroutine to determine if a series of grid points are inside a polygon 
 !   following ray casting algorithm
@@ -53 +55 @@
 
   CONTAINS
+
+  SUBROUTINE fill3DI_2Dvec(matind, inmat, id1, id2, od1, od2, od3, outmat)
+! Subroutine to fill a 3D integer matrix from a series of indices from a 2D matrix
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: id1, id2, od1, od2, od3
+    INTEGER, DIMENSION(id1,id2), INTENT(in)              :: inmat
+    INTEGER, DIMENSION(od1,od2), INTENT(in)              :: matind
+    INTEGER, DIMENSION(od1,od2,od3), INTENT(out)         :: outmat
+
+! Local
+    INTEGER                                              :: i, j
+
+!!!!!!! Variables
+! matind: equivalence on the 2D space of the location in inmat
+! inmat: values of the input matrix (1Dvec, time)
+! id1/2: dimensions of the input matrix
+! od1/3: dimensions of the output matrix
+! outmat: output matrix
+! NOTE: id2 == od3
+
+    fname = 'fill3DR_2Dvec'
+
+    DO i=1, od1 
+      DO j=1, od2
+        IF (matind(i,j) /= -1) THEN
+          outmat(i,j,:) = inmat(matind(i,j),:)
+        ELSE
+          outmat(i,j,:) = fillvalI
+        END IF
+      END DO
+    END DO
+
+  END SUBROUTINE fill3DI_2Dvec
+
+  SUBROUTINE fill3DR_2Dvec(matind, inmat, id1, id2, od1, od2, od3, outmat)
+! Subroutine to fill a 3D float matrix from a series of indices from a 2D matrix
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: id1, id2, od1, od2, od3
+    REAL(r_k), DIMENSION(id1,id2), INTENT(in)            :: inmat
+    INTEGER, DIMENSION(od1,od2), INTENT(in)              :: matind
+    REAL(r_k), DIMENSION(od1,od2,od3), INTENT(out)       :: outmat
+
+! Local
+    INTEGER                                              :: i, j
+
+!!!!!!! Variables
+! matind: equivalence on the 2D space of the location in inmat
+! inmat: values of the input matrix (1Dvec, time)
+! id1/2: dimensions of the input matrix
+! od1/3: dimensions of the output matrix
+! outmat: output matrix
+! NOTE: id2 == od3
+
+    fname = 'fill3DR_2Dvec'
+
+    DO i=1, od1 
+      DO j=1, od2
+        IF (matind(i,j) /= -1) THEN
+          outmat(i,j,:) = inmat(matind(i,j),:)
+        ELSE
+          outmat(i,j,:) = fillval64
+        END IF
+      END DO
+    END DO
+
+  END SUBROUTINE fill3DR_2Dvec
 
   SUBROUTINE reconstruct_matrix(vectorXpos, vectorYpos, dvec, Xmin, Xmax, Ymin, Ymax, dmatx, dmaty,   &