res@cnFillMode = "CellFill" no pdf or ps output

Hi,

res@cnFillMode = "CellFill" working well with "x11" option in xwks =
gsn_open_wks("x11","pippo"),
but void pdf or ps file has been produced, why?
Thank you for your support
Filippo

;This example code illustrates how to access and visualize LP_DAAC_MOD Grid
file in NCL.
;If you have any questions, suggestions, comments on this example, please
use the HDF-EOS Forum (http://hdfeos.org/forums).
;If you would like to see an example of any other NASA HDF/HDF-EOS data
product that is not listed in the HDF-EOS Comprehensive Examples page (
http://hdfeos.org/zoo),
;feel free to contact us at eoshelp@hdfgroup.org or post it at the HDF-EOS
Forum (http://hdfeos.org/forums).

load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

; This is an example of a 2-D grid file data field.
; It is assumed users know how to obtain information such as _FillValue
from HDFView.
; For information about HDFView, visit
http://www.hdfgroup.org/hdf-java-html/hdfview/.

begin

  eos_file = addfile("MOD09GA.A2007268.h10v08.005.2007272184810.hdf", "r")
; Read file for the first time.
  he2_file = addfile("MOD09GA.A2007268.h10v08.005.2007272184810.hdf.he2",
"r") ; Read file for the second time to efficiently obtain lat and lon.

; List available fields and attributes.
; print(eos_file)

; read data field
  data_raw = eos_file->sur_refl_b01_1

; print(data_raw)

; Apply scale factor according to the field attribute
; correct_data_value = (raw_data_value - offset) / scale
  data = (data_raw - data_raw@add_offset) / data_raw@scale_factor

  data@_FillValue = (data_raw@_FillValue - data_raw@add_offset) /
data_raw@scale_factor

  data@long_name = data_raw@long_name
  data@units = data_raw@units

; Associate longitude and latitude
  data@lat2d = he2_file->GridLat_MODIS_Grid_500m_2D
  data@lon2d = he2_file->GridLon_MODIS_Grid_500m_2D

; open workstation
  xwks = gsn_open_wks("x11","pippo")

   res = True
   res@cnFillOn = True ; turn on color
   res@cnFillMode = "CellFill"
   res@cnLinesOn = False ; Turn off contour lines

   res@trGridType = "TriangularMesh" ; Necessary b/c lat, lon
   res@gsnAddCyclic = False ; data not global

   res@gsnPaperOrientation = "portrait"
          ; arrays have missing values.

   res@gsnSpreadColors = True ; use full colormap
   res@gsnAddCyclic = False ; data not cyclic
   res@gsnMaximize = True ; ps, pdf, pdf
   res@lbLabelAutoStride = True ; let NCL decide spacing
   res@lbOrientation = "Vertical" ; just for something different
   res@pmTickMarkDisplayMode = "Always" ; use NCL default

  res@pmTickMarkDisplayMode = "Always"

  res@sfXArray = data@lon2d
  res@sfYArray = data@lat2d
;
; Use the "CellFill" method to fill the contours, and
; outline the cell edges.
;

  res@gsnLeftStringFontHeightF = 16 ; make font smaller
  res@gsnRightStringFontHeightF = 16 ; make font smaller

; the following 5 sentences are used to create a zoomed image
 res@mpLimitMode = "LatLon"
 res@mpMinLatF = min(data@lat2d) ; Set limits of map, based on the min/max
of the dataset latitude/longitude
 res@mpMaxLatF = max(data@lat2d)
 res@mpMinLonF = min(data@lon2d)
 res@mpMaxLonF = max(data@lon2d)

  gsn_define_colormap(xwks,"BlAqGrYeOrReVi200") ; choose colormap
  res@tiMainString = "MOD09GA.A2007268.h10v08.005.2007272184810.hdf" ;
create title
  res@gsnLeftString = "Field name: " + data@long_name
  res@gsnRightString = "Units: " + data@units

  plot=gsn_csm_contour_map_ce(xwks,data,res) ; create plot

; cleaning up used resources
  delete(plot)
  delete(xwks)
  delete(data_raw)
  delete(data)
  delete(res)
  delete(eos_file)
  delete(he2_file)

end

-- 
*Filippo Locci*
*PhD Student*
International PhD in
Environmental Science and Engineering
University of Cagliari
Tel. +390706757701
Mob. +393471372292

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Dec 17 08:53:11 2012