Changeset 3069 for branches/dev_2802_OBStools/NEMOGCM/TOOLS/OBSTOOLS

Timestamp:

2011-11-09T16:18:38+01:00 (13 years ago)

Author:

djlea

Message:

Update OBS and ASM documentation. Small updates and fixes to dataplot.

Location:

branches/dev_2802_OBStools/NEMOGCM/TOOLS/OBSTOOLS/dataplot

Files:

• prend2.pro (modified) (2 diffs)
• read_feedback.pro (modified) (2 diffs)
• read_sla.pro (modified) (4 diffs)

branches/dev_2802_OBStools/NEMOGCM/TOOLS/OBSTOOLS/dataplot/prend2.pro

@@ -10 +10 @@
 COMMON pr2, view_landscape
 
+; get filename of postscript device  
+
+r = fstat(!D.UNIT)
+filename = r.name
+
 device,/close          
 set_plot,'x'
@@ -15 +20 @@
 if keyword_set(view) then begin
 
-; get filename of postscript device  
-
-   r = fstat(!D.UNIT)
-   filename = r.name
-
    if (view_landscape) then begin
-      spawn,'gv -landscape -swap '+filename
+      spawn,'gv -orientation=landscape -swap '+filename
    endif else begin
       spawn,'gv '+filename 

branches/dev_2802_OBStools/NEMOGCM/TOOLS/OBSTOOLS/dataplot/read_feedback.pro

@@ -61 +61 @@
   if (NumObs gt 0) then begin  
     ByteNam = intarr(strnam, NumVars)
-    ByteEntries = intarr(strnam, NumEntries)
+    if (n_elements(NumEntries) eq 1) then ByteEntries = intarr(strnam, NumEntries)
     ByteId = intarr(strwmo, NumObs)
     ByteType = intarr(strtyp, NumObs)
     ByteJulDRef = intarr(strjuld)
     VarNames = strarr(NumVars)
-    Entries = strarr(NumEntries)
+    if (n_elements(NumEntries) eq 1) then Entries = strarr(NumEntries)
     Id = strarr(NumObs)
     Type = strarr(NumObs)
@@ -148 +148 @@
   
       varid = ncdf_varid(ncid, strtrim(VarNames(ivar))+'_Hx')
-      ncdf_varget, ncid, varid, tmp
-      Hx(*,*,ivar) = tmp
-      ncdf_attget, ncid, varid, '_Fillvalue', FillValHx        
+      if (varid gt -1) then begin
+         ncdf_varget, ncid, varid, tmp
+         Hx(*,*,ivar) = tmp
+         ncdf_attget, ncid, varid, '_Fillvalue', FillValHx
+      endif else begin
+         FillValHx=0
+      endelse 
       
       varid = ncdf_varid(ncid, strtrim(VarNames(ivar))+'_QC')

branches/dev_2802_OBStools/NEMOGCM/TOOLS/OBSTOOLS/dataplot/read_sla.pro

@@ -67 +67 @@
   if (varid eq -1) then begin 
    varid = ncdf_varid(ncid, 'Longitudes')
+        if (varid eq -1) then begin
+             varid = ncdf_varid(ncid, 'longitude')
+        endif
         ncdf_attget, ncid, varid, 'scale_factor', scale_factor
   endif
@@ -87 +90 @@
   if (varid eq -1) then begin
    varid = ncdf_varid(ncid, 'Latitudes')
+        if (varid eq -1) then begin
+             varid = ncdf_varid(ncid, 'latitude')
+        endif
         ncdf_attget, ncid, varid, 'scale_factor', scale_factor
   endif
@@ -103 +109 @@
     
   varid = ncdf_varid(ncid, 'JULD')
+  if (varid eq -1) then begin
+        varid = ncdf_varid(ncid, 'time')
+  endif
   if (varid ne -1) then begin 
    ncdf_varget, ncid, varid, Dats
@@ -188 +197 @@
   ncdf_attget, ncid, varid, '_FillValue', FillValue
   scale_factor=1.
-  if (ncinfo.ndims gt 2) then ncdf_attget, ncid, varid, 'scale_factor', scale_factor
+  add_offset=0.
+  ;if (ncinfo.ndims gt 2) then ncdf_attget, ncid, varid, 'scale_factor', scale_factor
+  ncviq=ncdf_varinq(ncid, varid)
+  for iatt=0,ncviq.natts-1 do begin
+    ncaiq=ncdf_attname(ncid, varid, iatt)
+    if ( ncaiq eq 'scale_factor') then ncdf_attget, ncid, varid, 'scale_factor', scale_factor
+    if ( ncaiq eq 'add_factor') then ncdf_attget, ncid, varid, 'add_factor', add_factor
+  endfor
+
   if not keyword_set(quiet) then print,'SLA scale factor: ',scale_factor
+  if not keyword_set(quiet) then print,'SLA add offset: ',add_offset
   if not keyword_set(quiet) then print,'SLA FillValue: ',FillValue
   pts = where(obsval eq FillValue, count)
 ;  if not keyword_set(quiet) then print,count
-  obsval=obsval*scale_factor
+  obsval=obsval*scale_factor+add_offset
   if count gt 0 then obsval(pts) = rmdi