Changeset 454 for trunk/SRC

Timestamp:

04/07/11 10:03:36 (13 years ago)

Author:

smasson

Message:

bugfix to keep double precision of vertical scale factors in all cases

Location:

trunk/SRC

Files:

• Computation/e3t_3d.pro (modified) (1 diff)
• Computation/e3u_3d.pro (modified) (1 diff)
• Computation/e3v_3d.pro (modified) (1 diff)
• Computation/e3w_3d.pro (modified) (1 diff)
• Grid/ncdf_meshread.pro (modified) (1 diff)

trunk/SRC/Computation/e3t_3d.pro

@@ -44 +44 @@
     keyword_set(e1) AND keyword_set(e2):arr2d = e1t[fstx:lstx, fsty:lsty] * e2t[fstx:lstx, fsty:lsty]
     keyword_set(e1)                    :arr2d = e1t[fstx:lstx, fsty:lsty]
-    keyword_set(e2):arr2d =                             e2t[fstx:lstx, fsty:lsty]
-    ELSE:arr2d = replicate(1., nx*ny)
+                        keyword_set(e2):arr2d =                             e2t[fstx:lstx, fsty:lsty]
+    ELSE:arr2d = replicate(1.d, nx*ny)
   ENDCASE
   e3t3d = arr2d[*] # e3t[firstzt:lastzt]

trunk/SRC/Computation/e3u_3d.pro

@@ -70 +70 @@
   ENDCASE
 ;
-  IF n_elements(arr2d) NE 0 THEN e3u3d = temporary(e3u3d) * ( arr2d[*] # replicate(1., nzt) )
+  IF n_elements(arr2d) NE 0 THEN e3u3d = temporary(e3u3d) * ( arr2d[*] # replicate(1.d, nzt) )
 ;
   return, e3u3d

trunk/SRC/Computation/e3v_3d.pro

@@ -59 +59 @@
   ENDCASE
 ;
-  IF n_elements(arr2d) NE 0 THEN e3v3d = temporary(e3v3d) * ( arr2d[*] # replicate(1., nzt) )
+  IF n_elements(arr2d) NE 0 THEN e3v3d = temporary(e3v3d) * ( arr2d[*] # replicate(1.d, nzt) )
 ;
   return, e3v3d

trunk/SRC/Computation/e3w_3d.pro

@@ -39 +39 @@
     keyword_set(e1):arr2d = e1t[firstxt:lastxt, firstyt:lastyt]
     keyword_set(e2):arr2d = e2t[firstxt:lastxt, firstyt:lastyt]
-    ELSE:arr2d = replicate(1., nxt*nyt)
+    ELSE:arr2d = replicate(1.d, nxt*nyt)
   ENDCASE
   e3w_3d = arr2d[*] # e3w[firstzw:lastzw]

trunk/SRC/Grid/ncdf_meshread.pro

@@ -553 +553 @@
     ENDCASE
     nothing = execute(command)
-    command = namevar[i]+'=float('+namevar[i]+')'
-    nothing = execute(command)
+    IF strmid(fnamevar[i], 0, 2) NE  'e3' THEN BEGIN
+      command = namevar[i]+'=float('+namevar[i]+')'    
+      nothing = execute(command)
+    ENDIF
     command = 'if size('+namevar[i]+', /n_dimension) gt 0 then '+namevar[i]+' = reform('+namevar[i]+', /over)'
     nothing = execute(command)