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

From: Wei Huang <huangwei_at_nyahnyahspammersnyahnyah>
Date: Mon Dec 17 2012 - 09:59:30 MST

Filippo,

Have you tried edit your script and change "x11" to "pdf" or "ps"?
Now your script has x11, so you get a x11 term.
If you change x11 to pdf, then you should get a pdf file.
And change to ps will let you get a ps file.
You won't be able to get all three at the same time.

Wei

huangwei@ucar.edu
VETS/CISL
National Center for Atmospheric Research
P.O. Box 3000 (1850 Table Mesa Dr.)
Boulder, CO 80307-3000 USA
(303) 497-8924

On Dec 17, 2012, at 8:52 AM, Filippo Locci wrote:

> 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_at_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

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

This archive was generated by hypermail 2.1.8 : Fri Dec 21 2012 - 10:43:23 MST