Changeset 501 for trunk/SRC/Grid/ncdf_meshroms.pro

Timestamp:

09/20/16 20:15:00 (8 years ago)

Author:

smasson

Message:

set of bugfixes...

File:

• trunk/SRC/Grid/ncdf_meshroms.pro (modified) (4 diffs)

trunk/SRC/Grid/ncdf_meshroms.pro

@@ -273 +273 @@
     nothing = execute(command)
   ENDFOR
-  d1t = 1.e3*(shift(d1t, -1, 0) - d1t)
-  d2t = 1.e3*(shift(d2t, 0, -1) - d2t)
+  d1t = shift(d1t, -1, 0) - d1t
+  d2t = shift(d2t, 0, -1) - d2t
   for i = 0, n_elements(namebase2)-1 do begin
     command = namebase2[i]+'t = '+namebase2[i]+'t[0:jpi-2, 0:jpj-2]'
@@ -306 +306 @@
     nothing = execute(command)
   ENDFOR
-  tmpsave = 2. * 1.e3 * d1u[0, 0:jpj-2]
-  d1u = 1.e3*(shift(d1u, -1, 0) - d1u)
-  d2u = 1.e3*(shift(d2u, 0, -1) - d2u)
+  tmpsave = 2. * d1u[0, 0:jpj-2]
+  d1u = shift(d1u, -1, 0) - d1u
+  d2u = shift(d2u, 0, -1) - d2u
   for i = 0, n_elements(namebase2)-1 do begin
     command = namebase2[i]+'u = '+namebase2[i]+'u[*, 0:jpj-2]'
@@ -342 +342 @@
     nothing = execute(command)
   ENDFOR
-  d1v = 1.e3*(shift(d1v, -1, 0) - d1v)
-  tmpsave = 2. * 1.e3 * d2v[0:jpi-2, 0]
-  d2v = 1.e3*(shift(d2v, 0, -1) - d2v)
+  d1v = shift(d1v, -1, 0) - d1v
+  tmpsave = 2. * d2v[0:jpi-2, 0]
+  d2v = shift(d2v, 0, -1) - d2v
   for i = 0, n_elements(namebase2)-1 do begin
     command = namebase2[i]+'v = '+namebase2[i]+'v[0:jpi-2, *]'
@@ -372 +372 @@
     nothing = execute(command)
   ENDFOR
-  tmpsave1 = 2. * 1.e3 * d1f[0, *]
-  d1f = 1.e3*(shift(d1f, -1, 0) - d1f)
-  tmpsave2 = 2. * 1.e3 * d2f[*, 0]
-  d2f = 1.e3*(shift(d2f, 0, -1) - d2f)
+  tmpsave1 = 2. * d1f[0, *]
+  d1f = shift(d1f, -1, 0) - d1f
+  tmpsave2 = 2. * d2f[*, 0]
+  d2f = shift(d2f, 0, -1) - d2f
   fmaskredy = byte(maskf[jpi-1, *])
   IF jpk GT 1 THEN fmaskredy = reform(fmaskredy[*]#replicate(1b, jpk), 1, jpj, jpk, /overwrite)