<div>Dear friends,</div>
<div> </div>
<div>I met a problem when I used NCL. The warning is "<font size="3" face="Calibri">ContourPlotInitialize: 0.0 not currently supported as a missing value; expect inaccurate plot". Is there any kind friend can help me handle this problem? Thanks a lot. I attach the codes as below.</font></div>
<div><font size="3" face="Calibri"></font> </div>
<div><font size="3" face="Calibri">Best,</font></div>
<div><font size="3" face="Calibri">Yi</font></div>
<div><font size="3" face="Calibri"></font> </div>
<div><font size="3" face="Calibri"></font> </div>
<div><font size="3" face="Calibri">load "$NCARG_ROOT/lib/ncarg/nclex/gsun/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</font></div>
<div><font size="3" face="Calibri">; This is an example of a 2-D data field.<br>; It is assumed users know how to obtain information such as _FillValue from HDFView.<br>; For information about HDFView, visit <a href="http://www.hdfgroup.org/hdf-java-html/hdfview/">http://www.hdfgroup.org/hdf-java-html/hdfview/</a>.</font></div>
<div><font size="3" face="Calibri">begin<br>eos_file=addfile("C1978303124834.L2_MLAC.hdf", "r") ; Read file. <br>; Since the file is not HDF-EOS2, it is not necessary to append .he2</font></div>
<div><font size="3" face="Calibri">data=eos_file->nLw_443 ; read data field<br>; Here, because the file is not HDF-EOS2, it is not necessary to append the group name</font></div>
<div><font size="3" face="Calibri"><a href="mailto:data@_FillValue=0b">data@_FillValue=0b</a></font></div>
<div><font size="3" face="Calibri">lat=eos_file->latitude<br>lon=eos_file->longitude</font></div>
<div><font size="3" face="Calibri">dimsize = dimsizes(data)<br>numlat = dimsize(0)<br>numlon = dimsize(1)</font></div>
<div><font size="3" face="Calibri">; define new latitude<br>nlat=new( (/numlat,numlon/), float)<br>mlon=new( (/numlat,numlon/), float)</font></div><font size="3" face="Calibri">
<div><br>do i = 0, numlat-1<br> do j =0, 246</div>
<div> ;when j = 0<br> if(j.eq.0)<br> nlat(i,j) = lat(i,j)<br> mlon(i,j) = lon(i,j)<br> continue<br> end if</div>
<div> ;when j = 1, ..., 245<br> if((j.ge.1).and.(j.le.245))<br> count=8*(j-1)+1<br> arr_fill=fspan(lat(i,(j-1)), lat(i, j), 9)<br> nlat(i, count:count+7) = (/ arr_fill(1:8) /)<br> arr_fill=fspan(lon(i,(j-1)), lon(i, j), 9)<br>
mlon(i, count:count+7) = (/ arr_fill(1:8) /)<br> continue<br> end if</div>
<div> ;when j = 246<br> if(j.eq.246)<br> count=8*(j-1)+1<br> array_fill=fspan(lat(i,(j-1)), lat(i, j), 8)<br> nlat(i, count:count+6) = (/ array_fill(1:7) /) <br> array_fill=fspan(lon(i,(j-1)), lon(i, j), 8)<br>
mlon(i, count:count+6) = (/ array_fill(1:7) /) <br> continue<br> end if<br> <br> end do<br>end do</div>
<div> </div>
<div><a href="mailto:data@lat2d=nlat">data@lat2d=nlat</a><br><a href="mailto:data@lon2d=mlon">data@lon2d=mlon</a></div>
<div> </div>
<div>xwks=gsn_open_wks("pdf","C1978303124834.L2_MLAC_nLw_443") ; open workstation</div>
<div>res=True ; plot mods desired<br><a href="mailto:res@cnFillOn=True">res@cnFillOn=True</a> ; enable contour fill<br><a href="mailto:res@gsnMaximize=True">res@gsnMaximize=True</a>; make plot large<br><a href="mailto:res@gsnPaperOrientation">res@gsnPaperOrientation</a> = "portrait" ; force portrait orientation<br>
<a href="mailto:res@cnLinesOn=False">res@cnLinesOn=False</a> ; turn off contour lines<br><a href="mailto:res@cnLineLabelsOn">res@cnLineLabelsOn</a> = False; turn off contour line labels<br><a href="mailto:res@gsnSpreadColors=True">res@gsnSpreadColors=True</a> ; use the entire color spectrum<br>
<a href="mailto:res@cnFillMode=%22RasterFill">res@cnFillMode="RasterFill</a>" ; faster<br><a href="mailto:res@lbOrientation=%22vertical">res@lbOrientation="vertical</a>" ; vertical labels<br><a href="mailto:res@cnMissingValFillPattern">res@cnMissingValFillPattern</a> = 0 ; missing value pattern is set to "SolidFill"<br>
<a href="mailto:res@cnMissingValFillColor=0">res@cnMissingValFillColor=0</a>; white color for missing values<br><a href="mailto:res@lbLabelAutoStride">res@lbLabelAutoStride</a> = True ; ensure no label overlap</div>
<div> </div>
<div><br>gsn_define_colormap(xwks,"BlAqGrYeOrReVi200") ; define colormap</div>
<div><br><a href="mailto:res@tiMainString=%22C1978303124834.L2_MLAC.hdf">res@tiMainString="C1978303124834.L2_MLAC.hdf</a>" ; create title<br><a href="mailto:res@gsnCenterString=%22nLw_443">res@gsnCenterString="nLw_443</a>" ; create center text<br>
plot=gsn_csm_contour_map_ce(xwks,data,res) ; plot on world map first</div>
<div> </div>
<div> </div>
<div> </div>
<div>delete(plot) ; cleaning up resources used<br>delete(nlat)<br>delete(mlon)<br>delete(xwks)<br>delete(data)<br>delete(res)<br>delete(eos_file)<br>end</div>
<div> </div>
<div> </div>
<div></div></font>