Re: about ocean figure error from ORCA grid

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Wed Dec 18 2013 - 14:28:26 MST

This file indicates why it is so important to look at your data!

Upon quick glance, it looks like your lat/lon values got divided by 10.0 before they were written to the file.
I use a "cnFillMode" of "RasterFill" as this is slightly faster.

According to the information on the file:

      float nav_lon ( y, x )
         standard_name : longitude
         units : degrees_east
         valid_min : -180
         valid_max : 180

      float nav_lat ( y, x )
         standard_name : latitude
         units : degrees_north
         valid_min : -77.01048
         valid_max : 89.94787

Note the values for "valid_min" and "valid_max".

They indicate that "nav_lon" is supposed to be -180 to 180, but your "print" call indicates it is -18 to 18.

"nav_lat" is supposed to be -77.01 to 89.95, but your "print" call indicates it is -7.70105 to 8.99479

If I simply multiply your lat/lon values by 10, then I get what looks like a correct plot.

See the attached.

--Mary

On Dec 17, 2013, at 9:34 PM, 김문현 <corps798@gmail.com> wrote:

>
> > Dear Mary Haley,
> >
> >
> >
> > I appreciate about your response with some coments.
> >
> >
> >
> > My data indicates lat/lon info like below;
> >
> > (0) min/max lat=-7.70105/8.99479
> > (0) min/max lon=-18/18
> >
> > Absolutely, It has some problems.
> >
> > But, when I showed global distrubution of sst and salinity using ncview,
> >
> > It seems like correct.
> >
> >
> >
> > I put this file at "incoming" directory in your server, ftp.cgd.ucar.edu.
> >
> > You just use gunzip for this gzip file.
> >
> >
> >
> > My info is below;
> >
> > (1) file name is "prodh_rs_sfc_onm_T_13_20090509_200911_002.nc.gz"
> >
> > I guess that *003.nc.gz file may transfer incorrectly.
> >
> > (2) NCL version is 6.1.2
> >
> > (3) System type is "Linux cirrus2 2.6.32.59-0.3-default #1 SMP 2012-04-27 11:14:44 +0200 x86_64 x86_64 x86_64 GNU/Linux".
> >
> > (4) There is no error message in my results, but some ploblem is explained previous mail.
> >
> >
> >
> > Thank you.
> >
> >
> >
> > From, Moon-Hyun
> >
> >
> >
> >
> >
> >
> >
> >
> >>
> >> ------------ 원본 메일 내용 ------------
> >> 보낸 사람 : "김문현"<corps798@gmail.com>
> >> 받는 사람 : "mhkim77@korea.kr" <mhkim77@korea.kr>
> >> 메일 제목 : Fwd: Re: [ncl-talk] about ocean figure error from ORCA grid
> >> 보낸 날짜 : 2013/12/18 12:08
> >>
> >> ---------- 전달된 메일 ----------
> >> 보낸사람: "Mary Haley" <haley@ucar.edu>
> >> 날짜: 2013. 12. 17. 오전 1:22
> >> 제목: Re: [ncl-talk] about ocean figure error from ORCA grid
> >> 받는사람: "김문현" <corps798@gmail.com>
> >> 참조: <ncl-talk@ucar.edu>
> >>
> >> 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

orca_RasterFill.png
Received on Wed Dec 18 14:28:44 2013

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