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
-- *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 Fri Dec 14 06:33:53 2012
This archive was generated by hypermail 2.1.8 : Fri Dec 21 2012 - 10:43:23 MST