Changeset 1178 in lmdz_wrf for trunk/tools/module_ForInterpolate.F90

Timestamp:

Oct 11, 2016, 6:26:01 PM (8 years ago)

Author:

lfita

Message:

Testing and fixing `npp' wrong minimum values for mindifflonlat

File:

• trunk/tools/module_ForInterpolate.F90 (modified) (4 diffs)

trunk/tools/module_ForInterpolate.F90

@@ -739 +739 @@
   REAL(r_k)                                              :: mindiffLl, dist
   REAL(r_k), DIMENSION(idimx,idimy)                      :: difflonlat
-  REAL(r_k), DIMENSION(2)                                :: extremelon, extremelat
+  REAL(r_k), DIMENSION(2)                                :: extremelon, extremelat, ipos
   INTEGER, DIMENSION(2)                                  :: iLl
   CHARACTER(LEN=50)                                      :: fname
@@ -797 +797 @@
             ix = INT(outLlw(1,iv,i,j))
             iy = INT(outLlw(2,iv,i,j))
-            dist = SQRT( (projlon(i,j)-inlonv(ix,iy))**2 + (projlat(i,j)-inlatv(ix,iy))**2 )
+            dist = SQRT( (projlon(i,j)-inlonv(ix,iy))**2. + (projlat(i,j)-inlatv(ix,iy))**2. )
             IF ( dist /= 0.) THEN
               outLlw(3,iv,i,j) = 1./dist
@@ -814 +814 @@
           difflonlat = SQRT((projlon(i,j)-inlonv)**2. + (projlat(i,j)-inlatv)**2.)
           mindiffLl = MINVAL(difflonlat)
-          outLlw(1:2,1,i,j) = index2DArrayR(difflonlat, idimx, idimy, mindiffLl)
+          ipos = index2DArrayR(difflonlat, idimx, idimy, mindiffLl)
+          outLlw(1:2,1,i,j) = ipos
           outLlw(3,1,i,j) = 1.
+          ix = outLlw(1,1,i,j)
+          iy = outLlw(2,1,i,j)
+          PRINT *,i,j,projlon(i,j),projlat(i,j),':',mindiffLl,'!',outLlw(1:2,1,i,j),'|',inlonv(ix,iy), &
+            inlatv(ix,iy), '<>', MINVAL(inlonv), MINVAL(inlatv)
+          IF (i > 1) THEN
+            PRINT *,inlonv(1,1),inlatv(1,1),':',index2DArrayR(inlonv, idimx, idimy, MINVAL(inlonv))
+            STOP
+          END IF
         END DO
       END DO
@@ -902 +911 @@
   REAL(r_k), DIMENSION(idimx,idimy), INTENT(in)          :: inlonv, inlatv
   CHARACTER(LEN=50), INTENT(in)                          :: intkind
-  REAL(r_k), DIMENSION(pdimx,pdimy,d3), INTENT(in)       :: var3Din
+  REAL(r_k), DIMENSION(idimx,idimy,d3), INTENT(in)       :: var3Din
   REAL(r_k), DIMENSION(pdimx,pdimy,d3), INTENT(out)      :: varout