Changeset 2356 in lmdz_wrf

Timestamp:

Feb 20, 2019, 2:46:48 PM (6 years ago)

Author:

lfita

Message:

Adding:

• `from_ptlist_2DRKmatrix': Subroutine to construct a 2D RK matrix from a list of values accompaigned by a list of coordinates
• module_ForGen.so': python binding to the funcionts in module_generic.f90'

Location:

trunk/tools

Files:

• Makefile.trueno (modified) (5 diffs)
• module_generic.f90 (modified) (2 diffs)

trunk/tools/Makefile.trueno

@@ -54 +54 @@
 intsrcs = $(srcs) module_ForInterpolate.f90
 scisrcs = $(srcs)
+gensrcs = module_definitions.f90 module_basic.f90 module_generic.f90
 defsrcs = module_definitions.f90
 
@@ -92 +93 @@
         pyintmods.o \
         pyscimods.o \
+        pygenmods.o \
         pydefmods.o \
         trajectories_overlap.o
@@ -101 +103 @@
         pyintmods.o \
         pyscimods.o \
+        pygenmods.o \
         pydefmods.o
 
@@ -116 +119 @@
 
 clean :
-        $(RM) *.mod *.o interpolate module_ForDiag.so module_ForInt.so module_ForDistriCorrect.so module_ForSci.so module_ForDef.so
+        $(RM) *.mod *.o interpolate module_ForDiag.so module_ForInt.so module_ForDistriCorrect.so module_ForSci.so module_ForGen.so module_ForDef.so
 
 ########    #######
@@ -170 +173 @@
         $(F2PY) -c $(F2PYF) -I$(NCINCFOLD) -m module_ForSci $(scisrcs) -L$(NCLIBFOLD)
 
+pygenmods.o:
+        $(F2PY) -c $(F2PYF) -I$(NCINCFOLD) -m module_ForGen $(gensrcs) -L$(NCLIBFOLD)
+
 pydefmods.o:
         $(F2PY) -c $(F2PYF) -I$(NCINCFOLD) -m module_ForDef $(defsrcs) -L$(NCLIBFOLD)

trunk/tools/module_generic.f90

@@ -5 +5 @@
 ! continguos_homogene_zones: Subroutine to look for contiguous zones by looking by continuous grid points
 ! freeunit: provides the number of a free unit in which open a file
+! from_ptlist_2DRKmatrix: Subroutine to construct a 2D RK matrix from a list of values accompaigned 
+!   by a list of coordinates
 ! GetInNamelist: Subroutine to get a paramter from a namelistfile
 ! get_xyconlimits: Subroutine for getting the limits of contiguous values from a given point in a 2D matrix
@@ -985 +987 @@
   END SUBROUTINE get_xyconlimits
 
+  SUBROUTINE from_ptlist_2DRKmatrix(Npts, pts, vals, dx, dy, NOval, matrix)
+  ! Subroutine to construct a 2D RK matrix from a list of values accompaigned by a list of coordinates
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: Npts, dx, dy
+    REAL(r_k), INTENT(in)                                :: NOval
+    INTEGER, DIMENSION(Npts,2), INTENT(in)               :: pts
+    REAL(r_k), DIMENSION(Npts), INTENT(in)               :: vals
+    REAL(r_k), DIMENSION(dx,dy), INTENT(out)             :: matrix
+
+! Local
+    INTEGER                                              :: iv, i, j
+
+!!!!!!! Variables
+! Npts: Number of values of the list
+! pts: 2D matrix coordinates of the values
+! vals: list of values correspondant to the coordinates
+! dx, dy: shape of the 2D matrix
+! NOval: Value to assign when there is no coordinate
+! matrix: resultant matrix
+
+    fname = 'from_ptlist_2DRKmatrix'
+
+    matrix = NOval
+
+    DO iv=1, Npts
+      i = pts(iv,1)
+      j = pts(iv,2)
+      matrix(i,j) = vals(iv)
+    END DO
+
+    RETURN
+
+  END SUBROUTINE from_ptlist_2DRKmatrix
+
 END MODULE module_generic