Changeset 731 in lmdz_wrf

Timestamp:

Apr 27, 2016, 10:28:47 AM (9 years ago)

Author:

lfita

Message:

Adding `iptpos' as control of syncronizing file in 'EntireGlobalMap?'

File:

• trunk/tools/nc_var_tools.py (modified) (5 diffs)

trunk/tools/nc_var_tools.py

@@ -18773 +18773 @@
         mindiffLl = np.zeros((Nallpts), dtype=np.float)
 
+    iptpos=0
     for iv in range(Nallpts):
         lonv = longvs[iv]
@@ -18793 +18794 @@
         ilatlonfrac = index_mat(difffraclonlat, mindifffracLl)
 
-#    print 'mindifffracLl:', mindifffracLl, 'ilatlonfrac:', ilatlonfrac
-#    print 'frac lon, lat:', fraclon[ilatlonfrac[0],ilatlonfrac[1]], ',',             \
-#      fraclat[ilatlonfrac[0],ilatlonfrac[1]]
-#    print 'values lon, lat:', lonv, latv
+#        print 'values lon, lat:', lonv, latv
+#        print 'mindifffracLl:', mindifffracLl, 'ilatlonfrac:', ilatlonfrac
+#        print 'frac lon, lat:', fraclon[ilatlonfrac[0],ilatlonfrac[1]], ',',         \
+#          fraclat[ilatlonfrac[0],ilatlonfrac[1]]
       
 # Providing fraction range
@@ -18883 +18884 @@
 #        iend = [iyend, ixend]
 #        print 'find window; beginning', ibeg, 'end:', iend
-        lon = ilon[iybeg:iyend,ixbeg:ixend]
-        lat = ilat[iybeg:iyend,ixbeg:ixend]
+#        lon = ilon[iybeg:iyend,ixbeg:ixend]
+#        lat = ilat[iybeg:iyend,ixbeg:ixend]
+#        print 'Looking within _______'
+#        print 'lon:',lon[0,0],',',lon[-1,-1]
+#        print 'lat:',lat[0,0],',',lat[-1,-1]
     
 #        print 'lon _______'
@@ -18908 +18912 @@
 #        quit()
 
-        if mindiffLl == 0. and type(newvar[ilatlon[0],ilatlon[1]]) == type(amsk): 
-            percendone(iv,Ninpts,0.5,'done:')
+#        if mindiffLl == 0. and type(newvar[ilatlon[0],ilatlon[1]]) == type(amsk): 
+#            percendone(iv,Ninpts,0.5,'done:')
+        if mindiffLl == 0.: 
+            iptpos = iptpos + 1
 # Speeding it up!
 #            if mindiffLl > mindiff:
@@ -18942 +18948 @@
 #                  'ilatlon:', ilatlon[1]
 #                quit(-1)
-
-            if np.mod(iv,fracs) == 0: 
+            if np.mod(iptpos,fracs) == 0: 
+                print 'Lluis syncronizing!'
                 onc.sync()
     if out: