Re: Slow plot processing

Replace
    res@cnFillMode = "RasterFill"

with
    res@cnFillMode = "CellFill"

Also, a new plotting algorithim is being tested,
please ftp the file to

ftp ftp.cgd.ucar.edu
anonymous
email
cd incoming
put MYD09A1.A2007273.h03v07.005.2007285103507.hdf
quit

Send an email after the successful transfer.

On 12/14/12 6:33 AM, Filippo Locci wrote:
> Hi,
> I'm working with this example from
> http://www.hdfeos.org/zoo/index_openLPDAAC_Examples.php#MOD
> <http://www.hdfeos.org/zoo/index_openLPDAAC_Examples.php#MOD>
>
> ;This example code illustrates how to access and visualize LP_DAAC_MYD
> 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("MYD09A1.A2007273.h03v07.005.2007285103507.hdf",
> "r") ; Read file.
> he2_file = addfile("MYD09A1.A2007273.h03v07.005.2007285103507.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_b02
>
> ; print(data_raw)
>
> ; Apply scale factor according to the field attribute
> ; correct_data_value = (raw_data_value - offset) * scale_factor
> 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_MOD_Grid_500m_Surface_Reflectance
> data@lon2d=he2_file->GridLon_MOD_Grid_500m_Surface_Reflectance
>
> xwks=gsn_open_wks("pdf","MYD09A1.A2007273.h03v07.005.2007285103507_sur_refl_b02.ncl")
> ; open workstation
>
> res=True ; plot mods desired
> res@cnFillOn=True ; enable contour fill
> res@gsnMaximize=True; make plot large
> res@gsnPaperOrientation = "portrait" ; force portrait orientation
> res@cnLinesOn=False ; turn off contour lines
> res@cnLineLabelsOn = False; turn off contour line labels
> res@gsnSpreadColors = True ; use the entire color spectrum
> res@cnFillMode = "RasterFill" ; faster
> res@lbOrientation ="vertical" ;vertical labels
> res@cnMissingValFillPattern = 0 ; missing value pattern is set to "SolidFill"
> res@cnMissingValFillColor = 0; white color for missing values
>
> ; 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 = "MYD09A1.A2007273.h03v07.005.2007285103507.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(res)
> delete(xwks)
> delete(data)
> delete(data_raw)
> delete(eos_file)
> delete(he2_file)
>
> end
>
>
> The plotting elaboration is very slow, my PC seem not working;
>
> commenting the "zoom" part:
>
> ;res@mpMinLatF = min(data@lat2d)
> ;res@mpMaxLatF = max(data@lat2d)
> ;res@mpMinLonF = min(data@lon2d)
> ;res@mpMaxLonF = max(data@lon2d)
>
> plotting elaboration time is about 2 minutes.
>
> Are there errors in the script or my computer is a clunker?
>
> Thank you so much
>
> Filippo
>
>
>
>
> _______________________________________________
> 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 Fri Dec 14 07:44:34 2012