Re: about ocean figure error from ORCA grid

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Mon Dec 16 2013 - 09:22:25 MST

Dear Moon-Hyun Kim,

As far as I can tell, you are doing everything correctly. It looks like your "nav_lat" and "nav_lon" values may be problematic.

Did you look at the lat/lon values to makes sure they are okay?

print("min/max lat = " + min(lat) + "/" + max(lat))
print("min/max lon = " + min(lon) + "/" + max(lon))

Can you provide us with the data file? You can use our ftp account:

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

I don't think this will fix your problem, but since your array is a bit large, you might try this setting for faster plotting:

  res@cnFillMode = "RasterFill"

--Mary

On Dec 15, 2013, at 10:21 PM, 김문현 <corps798@gmail.com> wrote:

> Dear, ncl talkers.
> >
> >
> >
> > Hello,
> >
> > I would like to use ORCA grid of ocean model to analyze the air-sea interaction,
> >
> > but, have a problem such as attached figure file.
> >
> >
> >
> > My ORCA grid point is 1021 x 1442.
> >
> > I guess that I may need some infomation such as axis or bounds info, to use the sfXArray and sfYArray of each 2 dimensional x, y axis.
> >
> >
> >
> > How can I set up it? Please, help me.
> >
> >
> >
> > from Moon-Hyun Kim
> >
> >
> >
> > Ps: Attached files are err figure using below and my nc data head info.
> >
> > Below contents are my just to draw ORCA grid.
> >
> > -------------------------------------------
> >
> > load "$NCARG_ROOT/lib/ncarg/ncls/csm/gsn_code.ncl"
> > load "$NCARG_ROOT/lib/ncarg/ncls/csm/gsn_csm.ncl"
> >
> > begin
> >
> > dir = "/crlnas/vol31/gpc/GS5Archive/prodh/2013nc/onm/T/"
> > srcFileName = "prodh_rs_sfc_onm_T_13_20090509_200911_001.nc"
> >
> > orca = addfile(dir+srcFileName,"r")
> > sst = orca->votemper(0,0,:1019,1:) ; Read data
> > lon = orca->nav_lon
> > lat = orca->nav_lat
> >
> > ;---- PLOT
> > wks = gsn_open_wks("ps","orca")
> > gsn_define_colormap(wks,"gui_default")
> >
> > res = True
> > res@gsnMaximize = True ; Maximize plot in frame
> > res@trGridType = "TriangularMesh" ; Define the mesh type.
> >
> > res@sfXArray = lon
> > res@sfYArray = lat
> > res@sfMissingValueV = sst@_FillValue
> >
> > res@cnFillOn = True
> > res@cnFillColors = (/ 2, 3, 4, 5, 7, 8, 9, 11, 12, 13, 15, 16, 17, \
> > 19, 20, 21, 23/)
> >
> > res@cnLinesOn = False ; Turn lines off
> > res@cnLineLabelsOn = False ; Turn labels off
> >
> > res@mpProjection = "Orthographic"
> > res@mpCenterLatF = 50
> > res@mpGridLineColor = -1 ; Turn off grid lat/lon lines
> > res@mpGridAndLimbOn = True
> >
> > res@cnFillMode = "AreaFill" ; The default.
> > res@tiMainString = "Orca Grid - Default fill mode"
> >
> > map = gsn_csm_contour_map(wks,sst,res)
> >
> > end
> <info_gs5_ocean.txt><orca.gif>_______________________________________________
> 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 Dec 16 09:22:36 2013

This archive was generated by hypermail 2.1.8 : Wed Dec 18 2013 - 14:29:51 MST