Changeset 1178 in lmdz_wrf for trunk/tools

Timestamp:

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

Author:

lfita

Message:

Testing and fixing `npp' wrong minimum values for mindifflonlat

Location:

trunk/tools

Files:

• module_ForInterpolate.F90 (modified) (4 diffs)
• nc_var_tools.py (modified) (1 diff)

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
 

trunk/tools/nc_var_tools.py

@@ -15900 +15900 @@
             quit(-1)
     
+    # Checking right lon, lat ranges
+    iminlon = np.min(ilonvals)
+    imaxlon = np.max(ilonvals)
+    rminlon = np.min(rlonvals)
+    rmaxlon = np.max(rlonvals)
+
+    if np.abs(imaxlon - rmaxlon) > 100.: 
+        print warnmsg
+        print '  ' + fname + ': reshaping longitudes!'
+        print '    input maximum lon:', imaxlon   
+        print '    target maximum lon:', rmaxlon   
+        print '    shifting input longitudes...'
+        if rmaxlon > 180.:
+            ilonvals = np.where(ilonvals < 0., ilonvals + 360., ilonvals)
+        else:
+            ilonvals = np.where(ilonvals > 180., ilonvals - 360., ilonvals)
+
     idx = ilonvals.shape[1]
     idy = ilonvals.shape[0]