Changeset 70

Timestamp:

05/16/06 14:41:08 (18 years ago)

Author:

pinsard

Message:

update according to idl-NetCDF-v04.tar.gz available via ​http://www.ittvis.com/codebank/search.asp?FID=418; side effects are not tested

Location:

trunk/ReadWrite/idl-NetCDF

Files:

• ncdf_listdims.pro (modified) (1 diff)
• ncdf_quickread/README_ncdf_quickread.txt (modified) (1 diff)
• ncdf_quickread/ncdf_quickread.pro (modified) (1 diff)
• ncdf_quickread/ncdf_quickread_helper.pro (modified) (4 diffs)
• ncdf_quickwrite/README_ncdf_quickwrite.txt (modified) (2 diffs)
• ncdf_quickwrite/ncdf_quickwrite_verbose.pro (modified) (1 diff)
• ncdf_read.pro (modified) (2 diffs)

trunk/ReadWrite/idl-NetCDF/ncdf_listdims.pro

@@ -1 @@
-FUNCTION ncdf_listdims, ncid, sizes = sizes
+FUNCTION ncdf_listdims,ncid
 
-  n = (ncdf_inquire(ncid)).ndims
+n=(ncdf_inquire(ncid)).ndims
 
-  sizes = lonarr(n)
-  names = strarr(n)
+names=strarr(n)
 
-  for i = 0, n-1 do begin
-    ncdf_diminq, ncid, i, name, size
-    names[i] = name
-    sizes[i] = size
-  endfor
+for i=0,n-1 do begin
+    ncdf_diminq,ncid,i,name,size
+    names[i]=name
+endfor
 
-  return, names
+return,names
 
 end

trunk/ReadWrite/idl-NetCDF/ncdf_quickread/README_ncdf_quickread.txt

@@ -28 +28 @@
 
          ncverbose : set equal to 1 if you want to include attributes
+
+         ncfields : Set to a string which is a space-separated list of
+                    variables to read.  Only the variables specified will be
+                    read in (and where ncverbose is set, only the attributes
+                    associated with these variables and the global attributes).
+                    If unset or the null string, then reads all fields.
 
 4) Type:  @ncdf_quickread

trunk/ReadWrite/idl-NetCDF/ncdf_quickread/ncdf_quickread.pro

@@ -5 +5 @@
 
 __commands= $
-  ncdf_quickread_helper(ncfile,verbose=ncverbose,prefix=ncprefix,/reform)
+           ncdf_quickread_helper(ncfile,verbose=ncverbose,prefix=ncprefix, $
+                                 fields=ncfields,/reform)
 
 for __loop=0,n_elements(__commands)-1 do $

trunk/ReadWrite/idl-NetCDF/ncdf_quickread/ncdf_quickread_helper.pro

@@ -34 +34 @@
 ;                           is used to remove dimensions of size 1.
 ;  02 Oct 03  A.M.Iwi       Change into helper routine for ncdf_quickread
+;  11 Aug 04  A.M.Iwi       Add "fields" option to read only certain fields.
+;                           Also, only stringify attributes of type CHAR.
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 
@@ -139 +141 @@
 
 function ncdf_quickread_helper, infile, verbose=verbose,  $
-                                prefix=prefix, reform=reform
+                                prefix=prefix, fields=fields, $
+                                reform=reform
 on_error,2
 compile_opt hidden
@@ -169 +172 @@
 commands="__ncid = NCDF_OPEN('"+infile+"')"
 
+subset=0
+if n_elements(fields) ne 0 then begin
+  if (fields ne '') then begin
+      subset=1
+      subfields=strsplit(fields,/extract)
+  endif
+endif
+
 ; Place the desired variables in local arrays.
 for i=0, ncidinfo.Nvars-1 do begin 
    vardata = NCDF_VARINQ(ncid, i)
-   varname = ncdf_quickread_validatename(vardata.Name) 
-   if keyword_set(prefix) then varname=prefix+varname
-   commands=$
-     [commands,"NCDF_VARGET, __ncid, "+strcompress(string(i))+", "+varname]
-   if keyword_set(reform) and vardata.ndims ge 2 $
-     then commands=[commands,varname+'=reform('+varname+')']
-   if (keyword_set(verbose)) then begin
-      for j=0, vardata.Natts-1 do begin
-         att = NCDF_ATTNAME(ncid, i, j)
-         attname = strcompress(varname+"_"+strcompress(att,/REMOVE_ALL))
-         commands=$
-           [commands,"NCDF_ATTGET, __ncid, "+strcompress(string(i))+$
-            ", '"+att+"', "+attname]
-         commands=$
-           [commands,attname+" = STRING("+attname+")"]
-      endfor
+   if not subset then begin
+       wanted=1
+   endif else begin
+       match = where(subfields eq vardata.Name, nmatch)
+       wanted=(nmatch ne 0)
+   endelse
+   if wanted then begin
+       varname = ncdf_quickread_validatename(vardata.Name) 
+       if keyword_set(prefix) then varname=prefix+varname
+       commands=$
+                [commands,"NCDF_VARGET, __ncid, "+strcompress(string(i))+", "+varname]
+       if keyword_set(reform) and vardata.ndims ge 2 $
+         then commands=[commands,varname+'=reform('+varname+')']
+       if (keyword_set(verbose)) then begin
+           for j=0, vardata.Natts-1 do begin
+               att = NCDF_ATTNAME(ncid, i, j)
+               attname = strcompress(varname+"_"+strcompress(att,/REMOVE_ALL))
+               commands=$
+                        [commands,"NCDF_ATTGET, __ncid, "+strcompress(string(i))+$
+                         ", '"+att+"', "+attname]
+               attinfo = ncdf_attinq(ncid, i, att)
+               if attinfo.datatype eq 'CHAR' then $
+                 commands=[commands,attname+" = STRING("+attname+")"]
+           endfor
+       endif
    endif
 endfor
@@ -198 +218 @@
      commands=$
        [commands,"NCDF_ATTGET, __ncid, /GLOBAL, '"+name+"', "+attname]
-     commands=$
-       [commands,attname+" = STRING("+attname+")"]
+     attinfo = ncdf_attinq(ncid, /global, name)
+     if attinfo.datatype eq 'CHAR' then $
+       commands=[commands,attname+" = STRING("+attname+")"]
   endfor
 endif
 
+ncdf_close,ncid
 commands=[commands,"NCDF_CLOSE, __ncid"]
 

trunk/ReadWrite/idl-NetCDF/ncdf_quickwrite/README_ncdf_quickwrite.txt

@@ -188 +188 @@
                    time = UNLIMITED ; // (100 currently)
            variables:
-                   float pressure(time, latitude, longitude) ;
+                   double pressure(time, latitude, longitude) ;
                            pressure:units = "pascals" ;
                            pressure:missing_value = 1.e+10f ;
@@ -195 +195 @@
                    float latitude(latitude) ;
                            latitude:units = "degrees" ;
-                   float ubar(time, latitude) ;
+                   double ubar(time, latitude) ;
                            ubar:units = "m s-1" ;
                    int year(time) ;

trunk/ReadWrite/idl-NetCDF/ncdf_quickwrite/ncdf_quickwrite_verbose.pro

@@ -16 +16 @@
 print,*__ncdf.commands,format='(a)'
 
-;for __loop=0,__ncdf.ncommands-1 do $
-;  if not execute((*__ncdf.commands)[__loop]) then stop
-
-for __loop=0,__ncdf.ncommands-1 do begin & print,(*__ncdf.commands)[__loop] & if not execute((*__ncdf.commands)[__loop]) then stop & endfor
+for __loop=0,__ncdf.ncommands-1 do $
+  if not execute((*__ncdf.commands)[__loop]) then stop
 
 ncdf_quickwrite_helper3,__ncdf

trunk/ReadWrite/idl-NetCDF/ncdf_read.pro

@@ -87 +87 @@
 ;; global attributes
 
-gatts=replicate({name:'', $
-                 datatype:'', $
-                 length:0L, $
-                 values:ptr_new()}, $
-                info.ngatts)
-for iatt=0,info.ngatts-1 do begin
-    name=ncdf_attname(id,iatt,/global)
-    gatts[iatt].name=name
-    att=ncdf_attinq(id,name,/global)
-    gatts[iatt].length=att.length
-    gatts[iatt].datatype=att.datatype
-    ncdf_attget,id,name,vals,/global
-    gatts[iatt].values=ptr_new(vals)
-endfor
+if info.ngatts gt 0 then begin
+    gatts=replicate({name:'', $
+                     datatype:'', $
+                     length:0L, $
+                     values:ptr_new()}, $
+                    info.ngatts)
+    for iatt=0,info.ngatts-1 do begin
+        name=ncdf_attname(id,iatt,/global)
+        gatts[iatt].name=name
+        att=ncdf_attinq(id,name,/global)
+        gatts[iatt].length=att.length
+        gatts[iatt].datatype=att.datatype
+        ncdf_attget,id,name,vals,/global
+        gatts[iatt].values=ptr_new(vals)
+    endfor
+endif else begin
+    ;; arbitary scalar value
+    ;; an empty list would be sensible but IDL doesn't support it
+    gatts=-1
+endelse
 
 
@@ -108 +114 @@
 vatts=replicate(ptr_new(),info.nvars)
 for ivar=0,info.nvars-1 do begin
-    vatts[ivar]=ptr_new(replicate({name:'', $
-                                      datatype:'', $
-                                      length:0L, $
-                                      values:ptr_new()}, $
+    if vinfo[ivar].natts gt 0 then begin
+        vatts[ivar]=ptr_new(replicate({name:'', $
+                                       datatype:'', $
+                                       length:0L, $
+                                       values:ptr_new()}, $
                                       vinfo[ivar].natts))
-    for iatt=0,vinfo[ivar].natts-1 do begin
-        name=ncdf_attname(id,ivar,iatt)
-        (*vatts[ivar])[iatt].name=name
-        att=ncdf_attinq(id,ivar,name)
-        (*vatts[ivar])[iatt].length=att.length
-        (*vatts[ivar])[iatt].datatype=att.datatype
-        ncdf_attget,id,ivar,name,vals
-        (*vatts[ivar])[iatt].values=ptr_new(vals)
-    endfor
+        for iatt=0,vinfo[ivar].natts-1 do begin
+            name=ncdf_attname(id,ivar,iatt)
+            (*vatts[ivar])[iatt].name=name
+            att=ncdf_attinq(id,ivar,name)
+            (*vatts[ivar])[iatt].length=att.length
+            (*vatts[ivar])[iatt].datatype=att.datatype
+            ncdf_attget,id,ivar,name,vals
+            (*vatts[ivar])[iatt].values=ptr_new(vals)
+        endfor
+    endif else begin
+        vatts[ivar]=ptr_new(-1)
+        ;; Pointer to arbitrary scalar -- analogous to case of lack of
+        ;; global attributes above.  We could put a <NullPointer> here
+        ;; instead, but try to be friendlier to code that might try
+        ;; to dereference it.
+    endelse
 endfor