Changeset 1622 in lmdz_wrf

Timestamp:

Sep 7, 2017, 7:39:11 PM (8 years ago)

Author:

lfita

Message:

Fixing various issues

File:

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

trunk/tools/module_scientific.f90

@@ -103 +103 @@
   meanwctrs = fillval64
   areas = fillval64
-  nvals = -1
+  nvals = fillvalI
   xvals = fillval64
   mvals = fillval64
@@ -109 +109 @@
   stdvals = fillval64
   quants = fillval64
-  nvcoords = -1
-  xvcoords = -1
+  nvcoords = fillvalI
+  xvcoords = fillvalI
   meanvnctrs = fillval64
   meanvxctrs = fillval64
@@ -255 +255 @@
     DO j=1,dy
       IF (poly(i,j)) THEN
-        polypts(ip,:) = (/ i,j /)
+        polypts(ip,1) = i
+        polypts(ip,2) = j
         polyvals(ip) = values(i,j)
         SELECT CASE (TRIM(projS))
@@ -291 +292 @@
 
   IF (dbg) THEN
-    PRINT *,'  Values: ', polyvals
+    PRINT *,'  grdi_coord lon lat value _______ '
+    DO ip=1, Npolypts
+      PRINT *, polypts(ip,:), ';', lon(polypts(ip,1),polypts(ip,2)), lat(polypts(ip,1),polypts(ip,2)),&
+         ':', polyvals(ip)
+    END DO
   END IF
 
@@ -336 +341 @@
         meanvxctr(2) = meanvxctr(2)+lat(i,j)*distvx(ip)*yres*ycorr(i,j)
 
-        sumDISTnlon = sumDISTnlon + lon(i,j)*distvn(ip)*xres*xcorr(i,j)
-        sumDISTnlat = sumDISTnlat + lat(i,j)*distvn(ip)*yres*ycorr(i,j)
-        sumDISTxlon = sumDISTxlon + lon(i,j)*distvx(ip)*xres*xcorr(i,j)
-        sumDISTxlat = sumDISTxlat + lat(i,j)*distvx(ip)*yres*ycorr(i,j)
+        sumDISTnlon = sumDISTnlon + distvn(ip)*xres*xcorr(i,j)
+        sumDISTnlat = sumDISTnlat + distvn(ip)*yres*ycorr(i,j)
+        sumDISTxlon = sumDISTxlon + distvx(ip)*xres*xcorr(i,j)
+        sumDISTxlat = sumDISTxlat + distvx(ip)*yres*ycorr(i,j)
 
         ip = ip + 1