Re: Slow plot processing

An alternative to Dennis' suggestion might be to set the resource
res@trGridType = "TriangularMesh"
This can often speed up the rasterfill plot style significantly when 2D coordinates are used. But we would still like to see the data per Dennis' instructions.
 -dave

On Dec 14, 2012, at 7:44 AM, Dennis Shea wrote:

> 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Dec 14 11:36:00 2012