Changeset 92 for trunk/procs/nc_put.pro

Timestamp:

04/30/08 22:21:22 (16 years ago)

Author:

ericg

Message:

Misc EG updates following merge with new nc_read.pro (r86)

File:

• trunk/procs/nc_put.pro (modified) (6 diffs)

trunk/procs/nc_put.pro

@@ -59 +59 @@
    tid = ncdf_dimdef(cdfid, 'time', /UNLIMITED)
 
+   IF debug_w THEN print, '    dim check ', dim_check
    
    IF dim_check NE dim(0) THEN BEGIN
@@ -69 +70 @@
 
    ; horizontal zonal axis
+   IF debug_w THEN print, '    building zonal axis, xid: ', xid
 ;   IF strpos(type, 'x') NE -1 AND keyword_set(X) THEN BEGIN 
       IF strpos (type, 'xy') NE -1 THEN BEGIN 
@@ -81 +83 @@
 
    ; horizontal meridional axis
-;   IF strpos(type, 'y') NE -1 AND keyword_set(Y) THEN BEGIN 
-      IF strpos (type, 'xy') NE -1 THEN BEGIN 
-         latid = ncdf_vardef(cdfid,'lat', [xid, yid], /FLOAT)
-      ENDIF ELSE latid = ncdf_vardef(cdfid,'nav_lat', [yid], /FLOAT)
-      NCDF_ATTPUT, cdfid, latid, 'units', y.units
-      NCDF_ATTPUT, cdfid, latid, 'long_name', y.long_name
-      NCDF_ATTPUT, cdfid, latid, 'valid_min', min(y.data), /float
-      NCDF_ATTPUT, cdfid, latid, 'valid_max', max(y.data), /float
-      swy = 1
-;   ENDIF 
+   IF debug_w THEN print, '    building Y axis, yid: ', yid
+   IF strpos (type, 'xy') NE -1 THEN BEGIN 
+      latid = ncdf_vardef(cdfid,'lat', [xid, yid], /FLOAT)
+   ENDIF ELSE latid = ncdf_vardef(cdfid,'nav_lat', [yid], /FLOAT)
+   NCDF_ATTPUT, cdfid, latid, 'units', y.units
+   NCDF_ATTPUT, cdfid, latid, 'long_name', y.long_name
+   NCDF_ATTPUT, cdfid, latid, 'valid_min', min(y.data), /float
+   NCDF_ATTPUT, cdfid, latid, 'valid_max', max(y.data), /float
+   swy = 1
 
    ; vertical axis
-      IF zpos NE -1 AND keyword_set(Z) THEN BEGIN 
-         vertid = ncdf_vardef(cdfid,z.name, [zid], /FLOAT)
-         NCDF_ATTPUT, cdfid, vertid, 'units', z.units
-         NCDF_ATTPUT, cdfid, vertid, 'long_name', z.long_name
-         NCDF_ATTPUT, cdfid, vertid, 'valid_min', min(z.data), /float
-         NCDF_ATTPUT, cdfid, vertid, 'valid_max', max(z.data), /float
-         swz = 1
-      ENDIF 
+   IF debug_w THEN print, '    building Z axis, zid: ', zid
+   IF zpos NE -1 AND keyword_set(Z) THEN BEGIN 
+      vertid = ncdf_vardef(cdfid,z.name, [zid], /FLOAT)
+      NCDF_ATTPUT, cdfid, vertid, 'units', z.units
+      NCDF_ATTPUT, cdfid, vertid, 'long_name', z.long_name
+      NCDF_ATTPUT, cdfid, vertid, 'valid_min', min(z.data), /float
+      NCDF_ATTPUT, cdfid, vertid, 'valid_max', max(z.data), /float
+      swz = 1
+   ENDIF 
    
    ; time axis
+   IF debug_w THEN print, '    building T axis, tid: ', tid
    IF tpos NE -1 AND keyword_set(T) THEN BEGIN 
       mid = NCDF_VARDEF(cdfid, 'time', [tid], /long)
@@ -136 +139 @@
 ; field attributes
 
+   IF debug_w THEN print, '    type, size fld.data ', type, size(fld.data)
+   
+   varn = strcompress(fld.short_name, /remove_all)
+
    CASE type OF
-      'x': fldid = NCDF_VARDEF(cdfid,fld.short_name, [xid], /FLOAT)
-      'y': fldid = NCDF_VARDEF(cdfid,fld.short_name, [yid], /FLOAT)
-      'z': fldid = NCDF_VARDEF(cdfid,fld.short_name, [zid], /FLOAT)
-      't': fldid = NCDF_VARDEF(cdfid,fld.short_name, [tid], /FLOAT)
-      'xy': fldid = NCDF_VARDEF(cdfid,fld.short_name, [xid, yid], /FLOAT)
-      'xz': fldid = NCDF_VARDEF(cdfid,fld.short_name, [xid, zid], /FLOAT)
-      'yz': fldid = NCDF_VARDEF(cdfid,fld.short_name, [yid, zid], /FLOAT)
-      'xyz': fldid = NCDF_VARDEF(cdfid,fld.short_name, [xid, yid, zid], /FLOAT)
-      'xyt': fldid = NCDF_VARDEF(cdfid,fld.short_name, [xid, yid, tid], /FLOAT)
-      'xzt': fldid = NCDF_VARDEF(cdfid,fld.short_name, [xid, zid, tid], /FLOAT)
-      'yzt':  fldid = NCDF_VARDEF(cdfid,fld.short_name, [xid, yid, zid, tid], /FLOAT)
-      'xyzt': fldid = NCDF_VARDEF(cdfid,fld.short_name, [xid, yid, zid, tid], /FLOAT)
+      'x': fldid = NCDF_VARDEF(cdfid,varn, [xid], /FLOAT)
+      'y': fldid = NCDF_VARDEF(cdfid,varn, [yid], /FLOAT)
+      'z': fldid = NCDF_VARDEF(cdfid,varn, [zid], /FLOAT)
+      't': fldid = NCDF_VARDEF(cdfid,varn, [tid], /FLOAT)
+      'xy': fldid = NCDF_VARDEF(cdfid,varn, [xid, yid], /FLOAT)
+      'xz': fldid = NCDF_VARDEF(cdfid,varn, [xid, zid], /FLOAT)
+      'yz': fldid = NCDF_VARDEF(cdfid,varn, [yid, zid], /FLOAT)
+      'xyz': fldid = NCDF_VARDEF(cdfid,varn, [xid, yid, zid], /FLOAT)
+      'xyt': fldid = NCDF_VARDEF(cdfid,varn, [xid, yid, tid], /FLOAT)
+      'xzt': fldid = NCDF_VARDEF(cdfid,varn, [xid, zid, tid], /FLOAT)
+      'yzt':  fldid = NCDF_VARDEF(cdfid,varn, [xid, yid, zid, tid], /FLOAT)
+      'xyzt': fldid = NCDF_VARDEF(cdfid,varn, [xid, yid, zid, tid], /FLOAT)
       ELSE:
    ENDCASE 
@@ -157 +164 @@
    NCDF_ATTPUT, cdfid, fldid, 'units', fld.units
    NCDF_ATTPUT, cdfid, fldid, 'long_name', fld.long_name
-   NCDF_ATTPUT, cdfid, fldid, 'short_name', fld.short_name
+   NCDF_ATTPUT, cdfid, fldid, 'short_name', varn
    NCDF_ATTPUT, cdfid, fldid, 'missing_value', fld.missing_value
    NCDF_ATTPUT, cdfid, fldid, 'valid_min', min(fld.data(where (fld.data NE fld.missing_value))), /float
@@ -178 +185 @@
 ;
    print, ' '
-   print, '    Writing '+fld.short_name+' to netCDF file ',ncdf_db+file_name
+   print, '    Writing '+varn+' to netCDF file ',ncdf_db+file_name