Changeset 769 in lmdz_wrf

Timestamp:

May 10, 2016, 3:50:40 PM (9 years ago)

Author:

lfita

Message:

Fortran routines provide values from 1,dimx; in 'ilonlatv', they should be 0,dimx-1!!!

File:

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

trunk/tools/nc_var_tools.py

@@ -14355 +14355 @@
         else:
             newvar = newnc.createVariable('lon','f8',('lon'))
-        basicvardef(newvar, 'lon', 'Longitudes','degrees_East')
+        basicvardef(newvar, 'longitude', 'Longitude','degrees_East')
         newvar[:] = lonmap
         newvar.setncattr('axis', 'X')
@@ -14364 +14364 @@
         else:
             newvar = newnc.createVariable('lat','f8',('lat'))
-        basicvardef(newvar, 'lat', 'Latitudes','degrees_North')
+        basicvardef(newvar, 'latitude', 'Latitude','degrees_North')
         newvar[:] = latmap
         newvar.setncattr('axis', 'Y')
@@ -14374 +14374 @@
         newvarindiff = newnc.createVariable('locindiff','f4',('inpts'))
         basicvardef(newvarindiff, 'locindiff', 'distance between input point and its final location','degree')
+        set_attributek(newvarindiff,'authorized_minimum_distance',mindiff,'R')
 
 # map variable
@@ -14391 +14392 @@
 
         basicvardef(newvarin, 'inpts', 'Equivalent point from the input source', '-')
-        newvar.setncattr('coordinates', 'lon lat')
+        newvarin.setncattr('coordinates', 'lon lat')
 
     else:
@@ -14425 +14426 @@
     Nptsf = len(ptsf)
     print Ninpts,'Npoints to find:', len(ptsf), ptsf[0:10], newvarindiff[ptsf[0:10]]
+    ipoint = ptsf[0]
     print '  ' + fname + ': from:', Ninpts,'re-locating:',Nptsf,'points starting at',\
       ipoint,'...'
@@ -14446 +14448 @@
             idiff, ilonlatv = fin.module_forinterpolate.interpolate(projlon, projlat,\
               lonvss, latvss, np.float64(mindiff), inptss)
+# Fortran routine provides ranges from 1,dimx!
+            ilonlatv = ilonlatv - 1
             for i in range(np.min([len(idiff),fracs])):
                 newvar[ilonlatv[i,0],ilonlatv[i,1]] = ovar[iri + i]