Changeset 626 in lmdz_wrf

Timestamp:

Sep 14, 2015, 5:40:15 PM (9 years ago)

Author:

lfita

Message:

Testing the radial cross section program

File:

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

trunk/tools/nc_var_tools.py

@@ -15251 +15251 @@
           'wrf': WRF kind
         [Nangle]: Number of angles
-        [radii]: length in grtid points of the angle
+        [radii]: length in grid points of the radius
       ncfile= netCDF file to use
       varn= ',' list of variables' name ('all', for all variables)
@@ -15343 +15343 @@
 
     trajpos = np.zeros(Ttraj,Nangle,raddi)
+
+# Grid position at each time-step and angle
+    varvalst = np.ones((Traj, Nangle, radii, 2), dtype=np.float)*fillValue
+    varvalstSQ = np.zeros((Traj, radii, radii), dtype=bool)
 
     it = 0
@@ -15382 +15386 @@
 # Assuming time as the first dimension in a C/python like way (t=0)
 
-            varvalst = np.ones((Nangle, radii), dtype=np.float)*fillValue
             posangle = np.zeros((radii,2), dtype=np.float)
 
@@ -16001 +16004 @@
         ivar = ivar + 1
 
-
 # global attributes
     objofile.setncattr('author', 'L. Fita')
@@ -16022 +16024 @@
 
     return 
+
+#        [trajfile]: ASCII file with the trajectory ('#' not readed)
+#          [time] [xpoint] [ypoint]
+#        [Tbeg]: equivalent first time-step of the trajectory within the netCDF file
+#        [lonname],[latname],[zname],[timename]: longitude, latitude, z and time variables names
+#        [timekind]: kind of time
+#          'cf': cf-compilant
+#          'wrf': WRF kind
+#        [Nangle]: Number of angles
+#        [radii]: length in grid points of the radius
+
+compute_tevolboxtraj('/bdd/PCER/workspace/lfita/etudes/FF/medic051215_2km/run/trajectory.dat@6,XLONG,XLAT,ZNU,Times,wrf,4,10',       \
+  '001/full_concatenated.nc', 'U10')
+
+ompute_tevolboxtraj_radialsec(values, ncfile, varn)
+quit()
+
 
 def DatesFiles(values, ncfile, varn):