Changeset 2318 in lmdz_wrf

Timestamp:

Feb 5, 2019, 4:18:20 PM (6 years ago)

Author:

lfita

Message:

Adding which dimension coincide with slice in `multi_spaceweightstats_in1DRKno_slc3v3'

File:

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

trunk/tools/module_scientific.f90

@@ -5642 +5642 @@
   END SUBROUTINE spaceweightstats
 
-  SUBROUTINE multi_spaceweightstats_in1DRKno_slc3v3(varin, Ngridsin, gridsin, percentages, varout,     &
-    di1, ds1, ds2, ds3, maxNgridsin)
+  SUBROUTINE multi_spaceweightstats_in1DRKno_slc3v3(varin, idv, Ngridsin, gridsin, percentages,       &
+    varout, di1, ds1, ds2, ds3, maxNgridsin)
   ! Subroutine to compute an spatial statistics value from a 1D RK matrix without running one into a 
   !   matrix of 3-variables slices of rank 3 using spatial weights
@@ -5649 +5649 @@
     IMPLICIT NONE
 
-    INTEGER, INTENT(in)                                  :: di1, ds1, ds2, ds3
+    INTEGER, INTENT(in)                                  :: di1, idv, ds1, ds2, ds3
     INTEGER, INTENT(in)                                  :: maxNgridsin
     INTEGER, DIMENSION(ds1,ds2,ds3), INTENT(in)          :: Ngridsin
@@ -5676 +5676 @@
 ! maxNgridsin: maximum number of grid points from the 3D matrix in any slice
 ! varin: 1D RK variable to be used
+! idv: which dimension of the sliced grids coincide with the dimension of 1D varin
 ! Ngridsin: number of grids from 3D RK matrix for each slice
 ! gridsin: coordinates of grids of the 3D RK matrix B to matrix of slices
@@ -5720 +5721 @@
             ! Getting the values
             DO iv=1, Nin
-              i1 = gin(iv,1)
+              i1 = gin(iv,idv)
               vin(iv) = varin(i1)
             END DO
@@ -5747 +5748 @@
             varout(s1,s2,s3,7) = Nin*1.
           ELSE
-            i1 = gridsin(s1,s2,s3,1,1)
+            i1 = gridsin(s1,s2,s3,1,idv)
             varout(s1,s2,s3,1) = varin(i1)
             varout(s1,s2,s3,2) = varin(i1)