Hi Friends,
Thank you so much for your help. I tried that and it looks pretty good.
However, now another problem come out: two subtitles overlapped (see the
attachment). I think maybe plot the
*"FIELD NAME: nLw_412" *
**
*and *
**
*"UNITS: " + **data_pre@units* <data_pre@units>
(now the codes use *res@gsnLeftString* <res@gsnLeftString>* = "FIELD NAME:
nLw_412" and* *res@gsnRightString* <res@gsnRightString>* = "UNITS: " + **
data_pre@units* <data_pre@units>)
in two different rows probably can solve this problem. But I don't know how
to do this. Could you please help me see this? Thanks a lot.
Best,
Yi Wang
The hdf file:
ftp://ftp.hdfgroup.uiuc.edu/pub/outgoing/veer/MAP_DATA/files/OceanColor/OCTS/L2/O1996306152450.L2_GAC.hdf
The ncl file:
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 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("O1996306152450.L2_GAC.hdf", "r") ; Read file.
; Since the file is not HDF-EOS2, it is not necessary to append .he2
data_pre=eos_file->nLw_412 ; read data field
; Here, because the file is not HDF-EOS2, it is not necessary to append the
group name
data_pre@_FillValue=0h
data_pre@_FillValue=1h
theIntercept=data_pre@intercept ; intercept=0.0, given from the attributes
of nLw_412
theSlope=data_pre@slope; slope=0.0010, given from the attributes of nLw_412
data=theIntercept+data_pre*theSlope
;data@long_name = data_pre@long_name
lat=eos_file->latitude
lon=eos_file->longitude
dimsize = dimsizes(data)
numlat = dimsize(0)
numlon = dimsize(1)
; define new latitude
nlat=new( (/numlat,numlon/), float)
mlon=new( (/numlat,numlon/), float)
numCol = eos_file@Number_of_Pixel_Control_Points
valueInCell = eos_file->cntl_pt_cols
step1 = valueInCell(2) - valueInCell(1) ;step1=8
step2 = valueInCell(numCol-1) - valueInCell(numCol-2) ;step2=7
do i = 0, numlat-1
do j =0, numCol-1 ;numCol-1=50
;when j = 0
if(j.eq.0)
nlat(i,j) = lat(i,j)
mlon(i,j) = lon(i,j)
continue
end if
;when j = 1, ..., numCol-2 ;numCol-2=49
if((j.ge.1).and.(j.le.(numCol-2))) ;numCol-2=49
count=step1*(j-1)+1 ;step1=8
arr_fill=fspan(lat(i,(j-1)), lat(i, j), (step1+1)) ;step1=8
nlat(i, count:count+(step1-1)) = (/ arr_fill(1:(step1)) /) ;step1=8
arr_fill=fspan(lon(i,(j-1)), lon(i, j), (step1+1)) ;step1=8
mlon(i, count:count+(step1-1)) = (/ arr_fill(1:(step1)) /) ;step1=8
continue
end if
;when j = numCol-1
if(j.eq.(numCol-1)) ;numCol-1=50
count=step1*(j-1)+1 ;step1=8
array_fill=fspan(lat(i,(j-1)), lat(i, j), (step2+1)) ;step2=7
nlat(i, count:count+(step2-1)) = (/ array_fill(1:(step2)) /)
;step2=7
array_fill=fspan(lon(i,(j-1)), lon(i, j), (step2+1)) ;step2=7
mlon(i, count:count+(step2-1)) = (/ array_fill(1:(step2)) /)
;step2=7
continue
end if
end do
end do
data@lat2d=nlat
data@lon2d=mlon
xwks=gsn_open_wks("pdf","O1996306152450.L2_GAC_nLw_412_limitedArea") ; 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 <res@cnFillMode=%22RasterFill>" ; faster
res@lbOrientation="vertical <res@lbOrientation=%22vertical>" ; vertical
labels
res@cnMissingValFillPattern = 0 ; missing value pattern is set to
"SolidFill"
res@cnMissingValFillColor=0; white color for missing values
res@lbLabelAutoStride = True ; ensure no label overlap
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") ; define colormap
;You should also make sure you have a .hluresfile in your home directory
;so you don't get any PLCHHQ error messages.
;see: http://www.ncl.ucar.edu/Document/Graphics/hlures.shtml
*res@tiMainString* <res@tiMainString>* = "O1996306152450.L2_GAC.hdf"
**res@gsnLeftString* <res@gsnLeftString>* = "FIELD NAME: nLw_412"
**res@gsnRightString* <res@gsnRightString>* = "UNITS: " +
**data_pre@units*<data_pre@units>
plot=gsn_csm_contour_map_ce(xwks,data,res) ; plot on world map first
delete(plot) ; cleaning up resources used
delete(nlat)
delete(mlon)
delete(xwks)
delete(data)
delete(res)
delete(eos_file)
end
On Mon, Jun 7, 2010 at 11:19 AM, Mary Haley <haley@ucar.edu> wrote:
> There are a few ways you can do this.
>
> The "tiMainString" resource gives you a main title, and the three
> gsnLeftString, gsnCenterString, gsnRightString gives you three subtitles.
> You can use all four resources, or any combination of them.
>
> The gsnCenterString will end up directly under tiMainString.
>
> You can make these individual subtitles, or one long title.
>
> Here's one way you might do the subtitles with a main title:
>
> filename = *"A20021612002192.L3m_R32_NSST_4.hdf*<res@tiMainString=%22A20021612002192.L3m_R32_NSST_4.hdf>
> *"*
> res@gsnCenterString = "field name: l3m_data"
> res@gsnLeftString = "(units: " + eos_file@Units + ")"
> res@tiMainString = filename
>
> You can also combine them into one long title. For example.
>
> res@tiMainString = filename + ": field name: l3m_data (" +
> eos_file@Units + ")"
>
> You should also make sure you have a .hluresfile in your home directory so
> you don't get any PLCHHQ error messages.
>
> See:
>
> http://www.ncl.ucar.edu/Document/Graphics/hlures.shtml
>
> --Mary
>
>
> *;;;;;;;;;;;the following codes have
> problems;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;*
>
> *res@tiMainString="A20021612002192.L3m_R32_NSST_4.hdf*<res@tiMainString=%22A20021612002192.L3m_R32_NSST_4.hdf>
> *"
>
> **res@gsnCenterString="field* <res@gsnCenterString=%22field>* name:
> l3m_data" ; create center text*
> *res@txString="(units* <res@txString=%22(units>*:" + **eos_file@Units*<eos_file@Units>
> * + ")"; create title*
>
>
> On Jun 7, 2010, at 9:36 AM, Yi Wang wrote:
>
> Hi friends,
>
> I met some problems. I want to plot the title, the field name: l3m_data,
> and the units: deg-C above the picture. The title, and the field name can be
> hard codes, however, the unites need to be obtained from the attributes. I
> post my codes and errors as below. Could anyone please help me see the
> codes? You can download the corresponding hdf file from
> ftp://ftp.hdfgroup.uiuc.edu/pub/outgoing/veer/MAP_DATA/files/OceanColor/MODISA/L3SMI/A20021612002192.L3m_R32_NSST_4.hdf.
> I hope to see the plot with text (see the attachment) when open the pdf file
> which is created by the ncl codes. Thank you so much!
>
> Best,
> Yi Wang
>
>
> 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 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("A20021612002192.L3m_R32_NSST_4.hdf", "r") ; Read file.
> ;Not an HDFEOS file.
> ;For more information, consult section 4.3.2 of
> http://hdfeos.org/software/ncl.php.
>
> ;units=eos_file@Units
> data_unsigned=eos_file->l3m_data; read data field
> data_unsigned@_FillValue=65535h ; 'h' is appended to avoid type
> mismatching
> theIntercept=data_unsigned@Intercept ; Intercept=-2.0, given in the
> attributes of l3m_data
> theSlope=data_unsigned@Slope; Slope=7.17185E-4, given in the attributes of
> l3m_data
>
> data_signed=where(data_unsigned.lt.inttoshort(0), abs(data_unsigned)+32767,
> data_unsigned) ;
> ; Convert unsigned integer to signed integer for NCL to properly handle
>
> data=theIntercept+data_signed*theSlope
>
> data!0="lat" ; This file does not contain geolocation fields.
> data!1="lon" ; To properly display the data, the latitude/longitude must be
> remapped from its metadata.
> ; See http://www.ncl.ucar.edu/Applications/Scripts/regrid_10.ncl for more
> information
>
> nLat=eos_file@Number_of_Lines ;4320
> nLon=eos_file@Number_of_Columns ;8640
> latStep=eos_file@Latitude_Step ;0.041666668
> lonStep=eos_file@Longitude_Step ;0.041666668
> latStartpoint=eos_file@SW_Point_Latitude ;-89.979164
> lonStartpoint=eos_file@SW_Point_Longitude ;-179.97917
>
> lat=ispan(0,nLat-1,1)*latStep+latStartpoint
> lon=ispan(0,nLon-1,1)*lonStep+lonStartpoint
> ;about ispan, see
> http://www.ncl.ucar.edu/Document/Functions/Built-in/ispan.shtml
>
> lat@units = "degrees_north"
> ;lat@Latitude_Units="degrees_North" is wrong
> ;lat@units="degrees_North" is also wrong
>
> lon@units = "degrees_east"
> ;lon@Longitude_Units="degrees_East" is wrong
> ;lon@units="degrees_East" is also wrong
>
> lat=lat(::-1) ; use north-to-south latitude ordering
> data&lat=lat
> data&lon=lon
>
>
> xwks=gsn_open_wks("pdf","A20021612002192.L3m_R32_NSST_4_l3m_data") ; open
> workstation
> setvalues NhlGetWorkspaceObjectId() ; make maximum filesize larger
> "wsMaximumSize" : 200000000
> end setvalues
> res=True ; plot mods desired
> res@cnFillOn=True ; enable contour fill
> res@cnLinesOn=False ; turn off contour lines
> res@gsnSpreadColors=True ; use the entire color spectrum
> res@cnFillMode="RasterFill <res@cnFillMode=%22RasterFill>" ; faster
> res@lbLabelAutoStride=True ; ensure labels do not overlap
> res@lbOrientation = "vertical" ; vertical labels
> res@cnMissingValFillPattern = 0 ; missing value pattern is set to
> "SolidFill"
> res@cnMissingValFillColor=0; white color for missing values
> ;res@gsnLeftStringFontHeightF=10 ; make text smaller
> ;res@gsnCenterStringFontHeightF=10
> ;res@gsnRightStringFontHeightF=10
>
> gsn_define_colormap(xwks,"BlAqGrYeOrReVi200") ; define colormap
>
> *;;;;;;;;;;;the following codes have
> problems;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;*
>
> *res@tiMainString="A20021612002192.L3m_R32_NSST_4.hdf*<res@tiMainString=%22A20021612002192.L3m_R32_NSST_4.hdf>
> *"
>
> **res@gsnCenterString="field* <res@gsnCenterString=%22field>* name:
> l3m_data" ; create center text*
> *res@txString="(units* <res@txString=%22(units>*:" + **eos_file@Units*<eos_file@Units>
> * + ")"; create title*
>
> *
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> *
>
> ;print(eos_file@Units)
>
> plot=gsn_csm_contour_map_ce(xwks,data,res)
>
> delete(plot) ; cleaning up resources used
> delete(lat)
> delete(lon)
> delete(xwks)
> delete(data)
> delete(res)
> delete(eos_file)
> end
>
>
> *The errors are as below.*
> *PLCHHQ - CHARACTER NUMBER 16 (l) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 18 (m) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 19 (_) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 20 (d) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 21 (a) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 22 (t) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 23 (a) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 11 (d) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 12 (e) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 13 (g) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 14 (-) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 16 ()) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 11 (d) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 12 (e) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 13 (g) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 14 (-) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 16 ()) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 11 (d) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 12 (e) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 13 (g) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 14 (-) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 16 ()) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 11 (d) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 12 (e) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 13 (g) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 14 (-) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 16 ()) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 11 (d) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 12 (e) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 13 (g) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 14 (-) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 16 ()) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 11 (d) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 12 (e) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 13 (g) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 14 (-) IS NOT A LEGAL FUNCTION CODE
> PLCHHQ - CHARACTER NUMBER 16 ()) IS NOT A LEGAL FUNCTION CODE*
> <problems.png>_______________________________________________
> 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 Jun 7 16:16:15 2010
This archive was generated by hypermail 2.1.8 : Mon Jun 07 2010 - 16:48:44 MDT