Re: plotting

Hi,

There's no way to guess at the lat/lon values, but there might be some information about them on the file. If there isn't, and your data is already projected onto a map, then you will need to get the map projection exactly right using NCL resources. It looks like you have started to do this, but as you said, the plot doesn't look quite right.

What happens if you add:

   res@gsnAddCyclic = False

?

Finally, can you provide us with the file? If so, you can use our ftp site:

http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP

Otherwise, can you email ncl-talk with the full output from:

    ncl_filedump 1992-2002MDOT060401.gx1v6.nc

Thanks,

--Mary

On Dec 26, 2012, at 2:37 PM, Alexander Semenov wrote:

> Hey
>
> I have a question regarding plotting the data which is supposed to be only over the oceans. For some reason my plot gets messed up (you can look at it). I was thinking of assigning the variable (SSH) lats and lons, but I don't know how to do that given the data has no lats and lons as variables (only as dimensions). below you find the file info, script and plot. Thank you!
>
> ncdump - h:
>
> netcdf \1992-2002MDOT060401.gx1v6 {
> dimensions:
> nlon = 320 ;
> nlat = 384 ;
> variables:
> float SSH(nlat, nlon) ;
> SSH:missing_value = 1.e+30f ;
> SSH:long_name = "Maximenko/Niiler Sea Surface Height" ;
>
> the script:
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>
> begin
>
> g = addfile("/import/c/w/asemenov/CESM_plots/obs/1992-2002MDOT060401.gx1v6.nc","r")
>
> pwv = g->SSH ; sea-ice fraction
>
> ;************************************************
>
> wks = gsn_open_wks("x11","Spring")
> gsn_define_colormap(wks,"BlWhRe")
>
>
> res = True
> res@gsnFrame = True ; So we can draw markers
> res@gsnMaximize = True
> res@gsnSpreadColors = True
> res@mpGridAndLimbOn = True
> res@mpGridLineDashPattern = 2
> res@mpGridLineDashPattern = 1
> res@mpGeophysicalLineColor = "Black"
> res@mpGeophysicalLineThicknessF = 1.5
> res@mpNationalLineColor = "Black"
> res@mpNationalLineThicknessF = 1.5
> res@mpUSStateLineColor = "Black"
> res@mpUSStateLineThicknessF = 1.5
> res@mpPerimLineColor = "Black"
> res@mpDataBaseVersion = "LowRes"
> res@cnFillOn = True
> res@mpFillColors = (/-1,-1,0,-1/)
> res@mpGridLatSpacingF = 30.
> res@mpGridLonSpacingF = 30.
> res@mpOutlineOn = True
> res@tiMainOn = True
> res@tiMainPosition = "Center"
> res@tiMainFontHeightF = 0.025
> res@lbLabelFontHeightF = 0.02
> res@gsnPolarLabelDistance = 1.08
> res@gsnPolarLabelFontHeightF= 0.02
> res@gsnLeftString = " "
> res@gsnRightString = " "
> res@pmLabelBarDisplayMode = "Always"
> res@lbLabelAutoStride = True
> res@lbOrientation = "Vertical"
> res@lbLabelFontHeightF = 0.025
> res@lbPerimOn = False
> res@lbAutoManage = False
> res@pmTickMarkDisplayMode = "Always"
> res@tmXTOn = False
> res@tmXBOn = True
> res@tmYLOn = False
> res@tmYROn = False
> res@tmXTLabelFontHeightF = 0.025
> res@tmXBLabelFontHeightF = 0.025
> res@tmYLLabelFontHeightF = 0.025
> res@tmYRLabelFontHeightF = 0.025
> ; tmres@tmYLValues = (/ 60,70,80 /)
> res@tmXBValues = (/ 150., -180., -150. /)
> res@tmXBLabelDeltaF = -0.8
> res@tmYLLabelDeltaF = -0.8
> res@cnLinesOn = False
> res@mpProjection = "Stereographic"
> res@mpLimitMode = "Corners"
> res@mpLeftCornerLatF = 41.0
> res@mpLeftCornerLonF = 140.8014
> res@mpRightCornerLatF = 44.0
> res@mpRightCornerLonF = -50.0
> res@mpCenterLatF = 90
> res@mpCenterLonF = -180
> res@mpFillOn = True
> res@cnLevelSelectionMode = "ManualLevels"
> res@cnMinLevelValF = -80
> res@cnMaxLevelValF = 80
> res@cnLevelSpacingF = 10
>
> res@tiMainString = "Sea-surface height, cm~C~Model-Obs difference"
> printVarSummary(pwv)
> map = gsn_csm_contour_map(wks,pwv(:,:),res)
>
> end
>
>
> <Screen Shot 2012-12-26 at 12.28.58 PM.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 Thu Dec 27 08:57:58 2012