There is no explicit function. People usually
write their own code. Did you try?
-- Also, for people to respond, more information is required. Is this one station? A grid? What form are the units of time? Ascii, text, binary? -- Also, in the future please make the Subject line specific as requested on the digest > When replying, please edit your Subject line so it is more specific > than "Re: Contents of ncl-talk digest..." People generally ignore question with no specific Subject. On 10/11/12 3:51 AM, Tidiane Diedhiou wrote: > hi, > I ask for a script that converts hourly data (06 hours) as a daily > average. I have hourly data 1989-2012 > > 2012/10/10 <ncl-talk-request@ucar.edu <mailto:ncl-talk-request@ucar.edu>> > > Send ncl-talk mailing list submissions to > ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> > > To subscribe or unsubscribe via the World Wide Web, visit > http://mailman.ucar.edu/mailman/listinfo/ncl-talk > or, via email, send a message with subject or body 'help' to > ncl-talk-request@ucar.edu <mailto:ncl-talk-request@ucar.edu> > > You can reach the person managing the list at > ncl-talk-owner@ucar.edu <mailto:ncl-talk-owner@ucar.edu> > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of ncl-talk digest..." > > > Today's Topics: > > 1. XY Plots over a Map (Soumik Basu) > 2. Re: XY Plots over a Map (Adam Phillips) > 3. about Taylor diagram ... (Ufuk Utku Turuncoglu (BE)) > 4. Re: Track one storm when there are also two other storms in > the domain (Dennis Shea) > 5. Re: Track one storm when there are also two other storms in > the domain (Dennis Shea) > 6. EOF question (Shinn Wong) > 7. A question about resource tmXBMode (Zheng Lu) > 8. Re: cmip5 ?lon2d,lat2d?data (?yan) > 9. register (Tidiane Diedhiou) > 10. Change color within text string? (Alan Crouse) > 11. stat_dispersion missing values (Schiffer, Nicole June) > 12. Re: about Taylor diagram ... (Dennis Shea) > 13. Re: EOF question (Dennis Shea) > 14. Re: stat_dispersion missing values (Dennis Shea) > 15. Re: A question about resource tmXBMode (Mary Haley) > 16. Re: cmip5 ?lon2d,lat2d?data (Mary Haley) > 17. Re: Change color within text string? (Mary Haley) > 18. Re: Coordinate subscript type mismatch. (Mary Haley) > 19. date (virginie hergault) > 20. Re: cmip5 ?lon2d,lat2d?data (?yan) > 21. Re: cmip5 ?lon2d,lat2d?data (Mary Haley) > 22. Re: date (Mary Haley) > 23. Re: date (Dennis Shea) > 24. Re: Coordinate subscript type mismatch. (David Brown) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 09 Oct 2012 10:17:47 -0800 > From: Soumik Basu <sbasu@iarc.uaf.edu <mailto:sbasu@iarc.uaf.edu>> > Subject: XY Plots over a Map > To: Ncl Talk <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: <0da992a4af35ab43bc9922081c245e34@arthur.iarc.uaf.edu > <mailto:0da992a4af35ab43bc9922081c245e34@arthur.iarc.uaf.edu>> > Content-Type: text/plain; charset=UTF-8; format=flowed > > Hi, > > I have XY plots of PDF of 6 variable for four different regions of > North America. I want to put these XY plots on the top of a map over > their corresponding regions ( six XY plots per region). Is it possible > to do so in NCL? > > Thanks, > Soumik > > -- > "Numbing the pain for a while will make it worse when you finally feel > it. > ~J.K. Rowling" > ***************************************** > Soumik Basu > Graduate Student, Research Assistant > International Arctic Research Center > University of Alaska Fairbanks > Fairbanks,Alaska,USA > ***************************************** > > > ------------------------------ > > Message: 2 > Date: Tue, 09 Oct 2012 12:55:54 -0600 > From: Adam Phillips <asphilli@ucar.edu <mailto:asphilli@ucar.edu>> > Subject: Re: XY Plots over a Map > To: ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> > Message-ID: <5074733A.20400@ucar.edu <mailto:5074733A.20400@ucar.edu>> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Hi Soumik, > Yes. You can do this by manually overlaying, setting gsnFrame = False, > and by using the various vp resources. > Something like this: > > wks = gsn_open_wks("ps","test") > > res = True ; spatial plot resource list > res@mpMinLatF = 10. > res@mpMaxLatF = 75. > res@mpMinLonF = -170 > res@mpMaxLonF = -45 > res@mpCenterLonF = (/res@mpMinLonF+res@mpMaxLonF)/2 > ..... > res@gsnFrame = False ; do not advance the frame, as we wish to > manually overlay xy plots onto page > plot = gsn_csm_contour_map(wks,data,res) > > res2 = True ; xy resource list > res2@vpXF = 0.2 ; placement of left border of xy plot along x-axis > of page in NDC units (0->1) > res2@vpYF = 0.8 ; placement of top border of xy plot along y-axis of > page in NDC units (0->1) > res2@vpWidthF = 0.2 ; width of XY plot in NDC units (0->1) > res2@vpHeightF = 0.2 ; height of XY plot in NDC units (0->1) > .... > res2@gsnFrame = False ; do not advance the frame > plot2 = gsn_csm_xy(wks,val,pdf0,res2) > > res2@vpXF = 0.35 > res2@vpYF = 0.6 > plot2 = gsn_csm_xy(wks,val,pdf1,res2) > drawNDCGrid(wks) ; draw the NDC grid in the background so you can > better choose your settings for vpXF/vpYF > frame(wks) ; advance the frame > > See the vp resources here: > http://www.ncl.ucar.edu/Document/Graphics/Resources/vp.shtml > and the drawNDCGrid documentation here: > http://www.ncl.ucar.edu/Document/Functions/Shea_util/drawNDCGrid.shtml > > Good luck, > Adam > > On 10/09/2012 12:17 PM, Soumik Basu wrote: > > Hi, > > > > I have XY plots of PDF of 6 variable for four different regions of > > North America. I want to put these XY plots on the top of a map over > > their corresponding regions ( six XY plots per region). Is it > possible > > to do so in NCL? > > > > Thanks, > > Soumik > > > > -- > ______________________________________________________________ > Adam Phillips asphilli@ucar.edu <mailto:asphilli@ucar.edu> > NCAR/Climate and Global Dynamics Division (303) 497-1726 > P.O. Box 3000 > Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli > > > > ------------------------------ > > Message: 3 > Date: Tue, 09 Oct 2012 23:44:38 +0300 > From: "Ufuk Utku Turuncoglu (BE)" <u.utku.turuncoglu@be.itu.edu.tr > <mailto:u.utku.turuncoglu@be.itu.edu.tr>> > Subject: about Taylor diagram ... > To: NCL USERS <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: <50748CB6.5020508@be.itu.edu.tr > <mailto:50748CB6.5020508@be.itu.edu.tr>> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Hi, > > I am trying to create a Taylor Diagram with modified correlation axis. > In my case correlation axis range is between 0.6 and 1.0. So, i want to > zoom in to higher correlation values. At this point, there is no any > problem and i could create the zoomed Taylor Diagram successfully. I > just wonder that if i change the limits of the correlation axis do i > need to modify the code which creates the RMS (root-mean-square) error > circles? In the taylor_diagram.ncl (Part 3), there is a comment that > says it is in the test mode. How can i be sure that the RMS circle > placed in the correct position? > > Best Regards, > > --ufuk > > > ------------------------------ > > Message: 4 > Date: Tue, 09 Oct 2012 14:54:42 -0600 > From: Dennis Shea <shea@ucar.edu <mailto:shea@ucar.edu>> > Subject: Re: Track one storm when there are also two other > storms in the domain > To: Yongxin Zhang <yongxin.fred@gmail.com > <mailto:yongxin.fred@gmail.com>> > Cc: ncl-talk <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>>, wrfhelp > <wrfhelp@ucar.edu <mailto:wrfhelp@ucar.edu>> > Message-ID: <50748F12.20107@ucar.edu <mailto:50748F12.20107@ucar.edu>> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > WRF questions should be sent to wrfhelp@ucar.edu > <mailto:wrfhelp@ucar.edu> > > --- > > This support question is being forwarded. > > > On 02/01/2012 06:31 AM, Yongxin Zhang wrote: > > My WRF domain covers the Atlantic basin and I will need to track only > > one storm while there are also two other storms in the domain. I used > > So-Young's script to plot the track but it did not work out well > since > > her script finds the minimum sea level pressure in the entire domain > > which in my case corresponds to one storm at one time and another > storm > > at another time. I am wondering if any of you worked on similar cases > > before and have a script or scripts to share. Any suggestions are > also > > highly appreciated. > > > > Thanks, > > Yongxin > > > > > > _______________________________________________ > > ncl-talk mailing list > > List instructions, subscriber options, unsubscribe: > > http://mailman.ucar.edu/mailman/listinfo/ncl-talk > > > > > > ------------------------------ > > Message: 5 > Date: Tue, 09 Oct 2012 14:55:45 -0600 > From: Dennis Shea <shea@ucar.edu <mailto:shea@ucar.edu>> > Subject: Re: Track one storm when there are also two other > storms in the domain > To: Yongxin Zhang <yongxin.fred@gmail.com > <mailto:yongxin.fred@gmail.com>> > Cc: ncl-talk <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>>, wrfhelp > <wrfhelp@ucar.edu <mailto:wrfhelp@ucar.edu>> > Message-ID: <50748F51.60002@ucar.edu <mailto:50748F51.60002@ucar.edu>> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > All ... SORRY ... something went wrong ... ignore lastr email. > > > On 10/09/2012 02:54 PM, Dennis Shea wrote: > > WRF questions should be sent to wrfhelp@ucar.edu > <mailto:wrfhelp@ucar.edu> > > > > --- > > > > This support question is being forwarded. > > > > > > On 02/01/2012 06:31 AM, Yongxin Zhang wrote: > >> My WRF domain covers the Atlantic basin and I will need to track > only > >> one storm while there are also two other storms in the domain. I > used > >> So-Young's script to plot the track but it did not work out well > since > >> her script finds the minimum sea level pressure in the entire domain > >> which in my case corresponds to one storm at one time and > another storm > >> at another time. I am wondering if any of you worked on similar > cases > >> before and have a script or scripts to share. Any suggestions > are also > >> highly appreciated. > >> > >> Thanks, > >> Yongxin > >> > >> > >> _______________________________________________ > >> 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 > > > > > > ------------------------------ > > Message: 6 > Date: Wed, 10 Oct 2012 10:12:39 +0800 > From: Shinn Wong <shinnshinnwong@gmail.com > <mailto:shinnshinnwong@gmail.com>> > Subject: EOF question > To: ncl-talk <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: > > <CAFvJNsQi4kTDsehz6OzGVx7rbzNaP1MiG_DWKZ1fF3uB5cF1xw@mail.gmail.com > <mailto:CAFvJNsQi4kTDsehz6OzGVx7rbzNaP1MiG_DWKZ1fF3uB5cF1xw@mail.gmail.com>> > Content-Type: text/plain; charset="iso-8859-1" > > Dear all, > I would like to perform a seaon-reliant EOF (S-EOF; Wang and An > 2005) based > on SST. Does anyone know how to do this task using > ncl? Thanks in advance for any helps and suggestions. > Regards, > Shinn > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: > http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20121010/460fa44b/attachment.html > > ------------------------------ > > Message: 7 > Date: Tue, 9 Oct 2012 23:19:30 -0400 > From: Zheng Lu <cfluzheng@gmail.com <mailto:cfluzheng@gmail.com>> > Subject: A question about resource tmXBMode > To: ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> > Message-ID: > > <CAA=sSiyXVwVtT=O1NPeFUW=o77zVL+RNxNQmtJ=A=4GGnEDcjg@mail.gmail.com > <mailto:4GGnEDcjg@mail.gmail.com>> > Content-Type: text/plain; charset="iso-8859-1" > > Hi, all, > > I have a question about resources of tmXBMode and tmYLMode. I tried > to plot > the contour of a 31*32 array using gsn_csm_contour with explicit > tick mark > mode. However, the figure looks wrong. First, the X-axis and Y-axis > labels > do not show. Also the location of colorbar is not right also. Could > someone > help me with this question? Thanks much! > > Zheng > > P.S. The code: > > 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" > load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl" > > > > ;********************************* > begin > ;********************************* > > > ss=asciiread("ss.txt",(/31,32/),"float") > > > wks = gsn_open_wks("x11","stat") ; open a ncgm file > gsn_define_colormap(wks,"precip_diff_12lev") ; choose a colormap > > res = True ; plot mods desired > res@cnFillMode = "RasterFill" > res@cnFillOn = True ; turn on color > res@cnLinesOn = False > > res@gsnDraw =False > res@gsnFrame=False > > plevel=(/-90,-60,-40,-30,-20,-10,10,20,30,40,60,90/) > pcolor= (/2,3,4,5,6,7,8,9,10,11,12,13,14 /) > res@cnLevelSelectionMode = "ExplicitLevels" > res@cnLevels = plevel > res@cnFillColors = pcolor > > res@tmXBMode="Explicit" > > res@tmXBValues=(/0,5,10,15.5,21,26,31/) > res@tmXBLables=(/"-1.5","-1.0","-0.5","0","0.5","1.0","1.5"/) > res@tmXBLabelFontHeightF=0.01 > > res@tmYLMode="Explicit" > res@tmYLValues=(/0,5,10,15,20,25,30/) > res@tmYLLables=(/"0","1","2","3","4","5","6"/) > res@tmYLLabelFontHeightF=0.01 > > res@vpWidthF = 0.5 > res@vpHeightF = 0.3 > > plot = gsn_csm_contour(wks,ss,res) > res_lines=True > > xx=(/15.5,15.5/) > yy=(/0,30/) > > line = gsn_add_polyline(wks,plot,xx,yy,res_lines) > draw(plot) > frame(wks) > > > end > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: > http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20121009/716d07ee/attachment.html > > ------------------------------ > > Message: 8 > Date: Wed, 10 Oct 2012 15:24:28 +0800 > From: ?yan <sunyan_10@hotmail.com <mailto:sunyan_10@hotmail.com>> > Subject: Re: cmip5 ?lon2d,lat2d?data > To: "Mary Haley" <haley@ucar.edu <mailto:haley@ucar.edu>>, > <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: <SNT124-DS12D74DB85B9521A89858328F8E0@phx.gbl> > Content-Type: text/plain; format=flowed; charset="gb2312"; > reply-type=original > > thank you very much,Mary Haley . your suggest is right. > I see example ESMF_regrid_6.ncl. in ESMF_regrid_6.ncl,the time > dimension is > only one,it works well.when I use the whole time dimensions,new problem > appears: > > my ncl script is: > 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" > load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl" > > begin > ;---Interpolation methods > methods = (/"patch"/) > > ;---Input file > srcFileName = > "tos_Omon_MPI-ESM-P_historical_r1i1p1_185001-200512.nc" > wgtFile = "CMIP5_2_World_" + methods + ".nc" > > ;---Get data and lat/lon grid from CMIP5 Grid > sfile = addfile(srcFileName,"r") > thetao = sfile->tos(0,:,:) > thetao@lat2d = sfile->lat > thetao@lon2d = sfile->lon > printVarSummary(thetao) > > Opt = True > Opt@SrcFileName = "CMIP5_SCRIP.nc" ; source file name > Opt@DstFileName = "World1deg_SCRIP.nc" ; destination > file name > Opt@ForceOverwrite = True > > Opt@SrcGridCornerLat = sfile->lat_vertices ; corners are > necessary > Opt@SrcGridCornerLon = sfile->lon_vertices ; for "conserve" > method > Opt@SrcMask2D = where(.not.ismissing(thetao),1,0) > Opt@DstGridType = "1x1" ; Destination grid > Opt@DstTitle = "World Grid 1-degree Resolution" > Opt@DstLLCorner = (/-89.75d, 0.00d /) > Opt@DstURCorner = (/ 89.75d, 359.75d /) > > print("Generating interpolation weights from CMIP5 to") > print("World 1 degree grid using the " + methods + " method.") > > Opt@WgtFileName = wgtFile > Opt@InterpMethod = methods > > ;---------------------------------------------------------------------- > ; Interpolate data from CMIP5 to World 1-degree grid. > ;---------------------------------------------------------------------- > thetao_regrid = ESMF_regrid(thetao,Opt) > printVarSummary(thetao_regrid) > > end > > > fatal:Eq: Dimension size, for dimension number 0, of operands does not > match, can't continue > > fatal:["Execute.c":7556]:Execute: Error occurred at or near line 1071 in > file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl > > fatal:["Execute.c":7556]:Execute: Error occurred at or near line 2810 in > file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl > > fatal:["Execute.c":7556]:Execute: Error occurred at or near line 2921 in > file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl > > I don't know what the problem is. thank you. > ------------------------------------------------------------------- > -----????----- > From: Mary Haley > Sent: Monday, October 01, 2012 1:49 AM > To: ?yan > Cc: ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> > Subject: Re: cmip5 ?lon2d,lat2d?data > > > On Sep 27, 2012, at 6:28 PM, ?yan wrote: > > > hello, > > (1) I use cmip5 data, many datas are ?lon2d,lat2d?data.when I use > > ?rcm2rgrid? function to regrid,some region is FillValue.why? > > (2) how to deal with ?lon2d,lat2d?data? > > thank you > > > > There's not enough information to help you. It doesn't surprise me > that you > would have missing values, especially if the rectilinear grid goes > outside > the range of your curvilinear grid. > > Did you print out any values, or try plotting the regridded data to > see if > it looks correct? > > You can also visit our ESMF regridding page. See example > ESMF_regrid_6.ncl: > > http://www.ncl.ucar.edu/Applications/ESMF.shtml#ex6 > > --Mary > > > > > ------------------------------------------------------------------ > > > > National Key Laboratory of Atmospheric Sciences and > > Geophysical Fluid Dynamics (LASG), > > Institute of Atmospheric Physics (IAP), > > Chinese Academy of Sciences (CAS), > > > > ------------------------------------------------------------------- > > _______________________________________________ > > ncl-talk mailing list > > List instructions, subscriber options, unsubscribe: > > http://mailman.ucar.edu/mailman/listinfo/ncl-talk > > > > ------------------------------ > > Message: 9 > Date: Wed, 10 Oct 2012 08:33:02 +0000 > From: Tidiane Diedhiou <tige83diedhiou@gmail.com > <mailto:tige83diedhiou@gmail.com>> > Subject: register > To: ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> > Message-ID: > > <CAFxT1sbE_iO-_nCJ8=Niff5552WzcB2wo=4t7PvWNBF8gtWeVQ@mail.gmail.com > <mailto:4t7PvWNBF8gtWeVQ@mail.gmail.com>> > Content-Type: text/plain; charset="iso-8859-1" > > i want to register to the the ncl-talk@ucar.edu > <mailto:ncl-talk@ucar.edu> mailing list > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: > http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20121010/0bf84ced/attachment.html > > ------------------------------ > > Message: 10 > Date: Wed, 10 Oct 2012 06:22:19 -0700 (PDT) > From: Alan Crouse <ac_soaring@yahoo.com <mailto:ac_soaring@yahoo.com>> > Subject: Change color within text string? > To: "ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>" > <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: > <1349875339.15890.YahooMailNeo@web120504.mail.ne1.yahoo.com > <mailto:1349875339.15890.YahooMailNeo@web120504.mail.ne1.yahoo.com>> > Content-Type: text/plain; charset="iso-8859-1" > > I have located the functions to change font characteristics within a > text string > (http://www.ncl.ucar.edu/Document/Graphics/function_code.shtml), but > changing color within a string does not seem to be an option.? Am I > missing something? > > I am using ncl 6.1.0-beta. > > Thank you, > > Alan > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: > http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20121010/07346d2a/attachment.html > > ------------------------------ > > Message: 11 > Date: Wed, 10 Oct 2012 15:15:30 +0000 > From: "Schiffer, Nicole June" <nschiff2@illinois.edu > <mailto:nschiff2@illinois.edu>> > Subject: stat_dispersion missing values > To: "ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>" > <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: > > <73548662BBDE994AA86657740D530A6410D36FF1@CITESMBX5.ad.uillinois.edu > <mailto:73548662BBDE994AA86657740D530A6410D36FF1@CITESMBX5.ad.uillinois.edu>> > Content-Type: text/plain; charset="us-ascii" > > Stat_disperion is working fine for the first 22 stats it provides > (through > % missing), but the last 8 are missing. Any idea why? The program > does not > throw any errors. > > Thanks, > Nicole > > ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ > * Nicole Schiffer > * Science Writing Intern (NCSA) > * Graduate Research Fellow (Dept. of Energy) > * Department of Atmospheric Sciences > * University of Illinois, Urbana-Champaign > * Email: nschiff2 [at] illinois [dot] edu > ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ > > > > > > ------------------------------ > > Message: 12 > Date: Wed, 10 Oct 2012 09:20:03 -0600 > From: Dennis Shea <shea@ucar.edu <mailto:shea@ucar.edu>> > Subject: Re: about Taylor diagram ... > To: "Ufuk Utku Turuncoglu (BE)" <u.utku.turuncoglu@be.itu.edu.tr > <mailto:u.utku.turuncoglu@be.itu.edu.tr>> > Cc: NCL USERS <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: <50759223.9050003@ucar.edu > <mailto:50759223.9050003@ucar.edu>> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Hi Ufuk > > On 10/9/12 2:44 PM, Ufuk Utku Turuncoglu (BE) wrote: > > Hi, > > > > I am trying to create a Taylor Diagram with modified correlation > axis. > > In my case correlation axis range is between 0.6 and 1.0. So, i > want to > > zoom in to higher correlation values. At this point, there is no any > > problem and i could create the zoomed Taylor Diagram successfully. I > > just wonder that if i change the limits of the correlation axis do i > > need to modify the code which creates the RMS (root-mean-square) > error > > circles? > > I think the answer is "yes" > > >In the taylor_diagram.ncl (Part 3), there is a comment that > > says it is in the test mode. How can i be sure that the RMS circle > > placed in the correct position? > > The Taylor Diagram code is not in test mode. > That comment has been removed. > > Good luck > > > > > Best Regards, > > > > --ufuk > > _______________________________________________ > > ncl-talk mailing list > > List instructions, subscriber options, unsubscribe: > > http://mailman.ucar.edu/mailman/listinfo/ncl-talk > > > > > ------------------------------ > > Message: 13 > Date: Wed, 10 Oct 2012 09:21:53 -0600 > From: Dennis Shea <shea@ucar.edu <mailto:shea@ucar.edu>> > Subject: Re: EOF question > To: Shinn Wong <shinnshinnwong@gmail.com > <mailto:shinnshinnwong@gmail.com>> > Cc: ncl-talk <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: <50759291.4030103@ucar.edu > <mailto:50759291.4030103@ucar.edu>> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > I suggest you contact the author(s). > > A person here at NCAR said this sounds like extended EOF analysis. > > > http://dx.doi.org/10.1029/2005GL022709 > > Wang, B., and S.-I. An (2005), A method for detecting season-dependent > modes of climate variability: S-EOF analysis, Geophys. Res. Lett., 32, > L15710, doi:10.1029/2005GL022709. > > Good Luck > > On 10/9/12 8:12 PM, Shinn Wong wrote: > > Dear all, > > I would like to perform a seaon-reliant EOF (S-EOF; Wang and An > > 2005) based on SST. Does anyone know how to do this task using > > ncl? Thanks in advance for any helps and suggestions. > > Regards, > > Shinn > > > > > > _______________________________________________ > > ncl-talk mailing list > > List instructions, subscriber options, unsubscribe: > > http://mailman.ucar.edu/mailman/listinfo/ncl-talk > > > > > ------------------------------ > > Message: 14 > Date: Wed, 10 Oct 2012 10:07:33 -0600 > From: Dennis Shea <shea@ucar.edu <mailto:shea@ucar.edu>> > Subject: Re: stat_dispersion missing values > To: "Schiffer, Nicole June" <nschiff2@illinois.edu > <mailto:nschiff2@illinois.edu>> > Cc: "ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>" > <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: <50759D45.1000604@ucar.edu > <mailto:50759D45.1000604@ucar.edu>> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > The documentation states > > "The results are returned in a one dimensional array of size [30]. > Double if x is type "double"; otherwise, "float" > > Up thru 6.0.0 the first 22 locations are filled with assorted > statistics. The last 8 were reserved for future use. > > You can examine the actual code via (search for stat_dispersion') > > %> less $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl > > === > Note: the current documentation for 6.1.0 (soon to be released) > > http://www.ncl.ucar.edu/Document/Functions/Contributed/stat_dispersion.shtml > > shows all 30 locations now contain statistics. > This was done at a users request. > > > > On 10/10/12 9:15 AM, Schiffer, Nicole June wrote: > > Stat_disperion is working fine for the first 22 stats it provides > (through > > % missing), but the last 8 are missing. Any idea why? The program > does not > > throw any errors. > > > > Thanks, > > Nicole > > > > ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ > > * Nicole Schiffer > > * Science Writing Intern (NCSA) > > * Graduate Research Fellow (Dept. of Energy) > > * Department of Atmospheric Sciences > > * University of Illinois, Urbana-Champaign > > * Email: nschiff2 [at] illinois [dot] edu > > ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ > > > > > > > > _______________________________________________ > > ncl-talk mailing list > > List instructions, subscriber options, unsubscribe: > > http://mailman.ucar.edu/mailman/listinfo/ncl-talk > > > > > ------------------------------ > > Message: 15 > Date: Wed, 10 Oct 2012 10:19:04 -0600 > From: Mary Haley <haley@ucar.edu <mailto:haley@ucar.edu>> > Subject: Re: A question about resource tmXBMode > To: Zheng Lu <cfluzheng@gmail.com <mailto:cfluzheng@gmail.com>> > Cc: ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> > Message-ID: <AACD7B7C-E9D1-435B-966A-C4D02C1BCB16@ucar.edu > <mailto:AACD7B7C-E9D1-435B-966A-C4D02C1BCB16@ucar.edu>> > Content-Type: text/plain; charset=iso-8859-1 > > Hi Zheng, > > You have two resources misspelled. You probably got a warning > message that looks something like this: > > "Warning: tmXBLables is not a resource in XXXXX at this time" > > Both "tmXBLabels" and "tmYLLabels" need to be spelled with "Labels", > and not "Lables". > > If you still have problems with the labelbar, then please include > the image. > > --Mary > > On Oct 9, 2012, at 9:19 PM, Zheng Lu wrote: > > > Hi, all, > > > > I have a question about resources of tmXBMode and tmYLMode. I > tried to plot the contour of a 31*32 array using gsn_csm_contour > with explicit tick mark mode. However, the figure looks wrong. > First, the X-axis and Y-axis labels do not show. Also the location > of colorbar is not right also. Could someone help me with this > question? Thanks much! > > > > Zheng > > > > P.S. The code: > > > > 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" > > load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl" > > > > > > > > ;********************************* > > begin > > ;********************************* > > > > > > ss=asciiread("ss.txt",(/31,32/),"float") > > > > > > wks = gsn_open_wks("x11","stat") ; open a ncgm file > > gsn_define_colormap(wks,"precip_diff_12lev") ; choose a colormap > > > > res = True ; plot mods desired > > res@cnFillMode = "RasterFill" > > res@cnFillOn = True ; turn on color > > res@cnLinesOn = False > > > > res@gsnDraw =False > > res@gsnFrame=False > > > > plevel=(/-90,-60,-40,-30,-20,-10,10,20,30,40,60,90/) > > pcolor= (/2,3,4,5,6,7,8,9,10,11,12,13,14 /) > > res@cnLevelSelectionMode = "ExplicitLevels" > > res@cnLevels = plevel > > res@cnFillColors = pcolor > > > > res@tmXBMode="Explicit" > > > > res@tmXBValues=(/0,5,10,15.5,21,26,31/) > > res@tmXBLables=(/"-1.5","-1.0","-0.5","0","0.5","1.0","1.5"/) > > res@tmXBLabelFontHeightF=0.01 > > > > res@tmYLMode="Explicit" > > res@tmYLValues=(/0,5,10,15,20,25,30/) > > res@tmYLLables=(/"0","1","2","3","4","5","6"/) > > res@tmYLLabelFontHeightF=0.01 > > > > res@vpWidthF = 0.5 > > res@vpHeightF = 0.3 > > > > plot = gsn_csm_contour(wks,ss,res) > > res_lines=True > > > > xx=(/15.5,15.5/) > > yy=(/0,30/) > > > > line = gsn_add_polyline(wks,plot,xx,yy,res_lines) > > draw(plot) > > frame(wks) > > > > > > end > > _______________________________________________ > > ncl-talk mailing list > > List instructions, subscriber options, unsubscribe: > > http://mailman.ucar.edu/mailman/listinfo/ncl-talk > > > > ------------------------------ > > Message: 16 > Date: Wed, 10 Oct 2012 10:21:35 -0600 > From: Mary Haley <haley@ucar.edu <mailto:haley@ucar.edu>> > Subject: Re: cmip5 ?lon2d,lat2d?data > To: ?yan <sunyan_10@hotmail.com <mailto:sunyan_10@hotmail.com>> > Cc: ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> > Message-ID: <5BE3E87E-30DA-418E-9513-F7CC682A7C51@ucar.edu > <mailto:5BE3E87E-30DA-418E-9513-F7CC682A7C51@ucar.edu>> > Content-Type: text/plain; charset=GB2312 > > Hi, > > Did you download a newer version of the "ESMF_regridding.ncl" script > from: > > http://www.ncl.ucar.edu/Applications/Files/ESMF_regridding.ncl > > Once you do this, you either need to copy it to > $NCARG_ROOT/lib/ncarg/nclscripts/esmf/, > or, if you don't have permission to do this, then put it in the same > directory as > your script, and edit your script to load this script instead: > > ; load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl" > load "./ESMF_regridding.ncl" > > --Mary > > > On Oct 10, 2012, at 1:24 AM, ?yan wrote: > > > thank you very much,Mary Haley . your suggest is right. > > I see example ESMF_regrid_6.ncl. in ESMF_regrid_6.ncl,the time > dimension is > > only one,it works well.when I use the whole time dimensions,new > problem > > appears: > > > > my ncl script is: > > 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" > > load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl" > > > > begin > > ;---Interpolation methods > > methods = (/"patch"/) > > > > ;---Input file > > srcFileName = > "tos_Omon_MPI-ESM-P_historical_r1i1p1_185001-200512.nc" > > wgtFile = "CMIP5_2_World_" + methods + ".nc" > > > > ;---Get data and lat/lon grid from CMIP5 Grid > > sfile = addfile(srcFileName,"r") > > thetao = sfile->tos(0,:,:) > > thetao@lat2d = sfile->lat > > thetao@lon2d = sfile->lon > > printVarSummary(thetao) > > > > Opt = True > > Opt@SrcFileName = "CMIP5_SCRIP.nc" ; source file name > > Opt@DstFileName = "World1deg_SCRIP.nc" ; destination > file name > > Opt@ForceOverwrite = True > > > > Opt@SrcGridCornerLat = sfile->lat_vertices ; corners are > necessary > > Opt@SrcGridCornerLon = sfile->lon_vertices ; for "conserve" > method > > Opt@SrcMask2D = where(.not.ismissing(thetao),1,0) > > Opt@DstGridType = "1x1" ; Destination grid > > Opt@DstTitle = "World Grid 1-degree Resolution" > > Opt@DstLLCorner = (/-89.75d, 0.00d /) > > Opt@DstURCorner = (/ 89.75d, 359.75d /) > > > > print("Generating interpolation weights from CMIP5 to") > > print("World 1 degree grid using the " + methods + " method.") > > > > Opt@WgtFileName = wgtFile > > Opt@InterpMethod = methods > > > > > ;---------------------------------------------------------------------- > > ; Interpolate data from CMIP5 to World 1-degree grid. > > > ;---------------------------------------------------------------------- > > thetao_regrid = ESMF_regrid(thetao,Opt) > > printVarSummary(thetao_regrid) > > > > end > > > > > > fatal:Eq: Dimension size, for dimension number 0, of operands > does not > > match, can't continue > > > > fatal:["Execute.c":7556]:Execute: Error occurred at or near line > 1071 in > > file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl > > > > fatal:["Execute.c":7556]:Execute: Error occurred at or near line > 2810 in > > file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl > > > > fatal:["Execute.c":7556]:Execute: Error occurred at or near line > 2921 in > > file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl > > > > I don't know what the problem is. thank you. > > ------------------------------------------------------------------- > > -----????----- > > From: Mary Haley > > Sent: Monday, October 01, 2012 1:49 AM > > To: ?yan > > Cc: ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> > > Subject: Re: cmip5 ?lon2d,lat2d?data > > > > > > On Sep 27, 2012, at 6:28 PM, ?yan wrote: > > > >> hello, > >> (1) I use cmip5 data, many datas are ?lon2d,lat2d?data.when I use > >> ?rcm2rgrid? function to regrid,some region is FillValue.why? > >> (2) how to deal with ?lon2d,lat2d?data? > >> thank you > >> > > > > There's not enough information to help you. It doesn't surprise > me that you > > would have missing values, especially if the rectilinear grid > goes outside > > the range of your curvilinear grid. > > > > Did you print out any values, or try plotting the regridded data > to see if > > it looks correct? > > > > You can also visit our ESMF regridding page. See example > ESMF_regrid_6.ncl: > > > > http://www.ncl.ucar.edu/Applications/ESMF.shtml#ex6 > > > > --Mary > > > >> > >> ------------------------------------------------------------------ > >> > >> National Key Laboratory of Atmospheric Sciences and > >> Geophysical Fluid Dynamics (LASG), > >> Institute of Atmospheric Physics (IAP), > >> Chinese Academy of Sciences (CAS), > >> > >> ------------------------------------------------------------------- > >> _______________________________________________ > >> 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 > > > > ------------------------------ > > Message: 17 > Date: Wed, 10 Oct 2012 10:24:25 -0600 > From: Mary Haley <haley@ucar.edu <mailto:haley@ucar.edu>> > Subject: Re: Change color within text string? > To: Alan Crouse <ac_soaring@yahoo.com <mailto:ac_soaring@yahoo.com>> > Cc: "ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>" > <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: <370D1C78-3669-4FFD-8E36-AB0D99A4BB91@ucar.edu > <mailto:370D1C78-3669-4FFD-8E36-AB0D99A4BB91@ucar.edu>> > Content-Type: text/plain; charset=iso-8859-1 > > Hi Alan, > > I'm afraid there's no way to do this. We have a ticket open on this > request specifically (NCL-220 for your reference). > There hasn't been a huge request for this feature, so it hasn't been > high on our priority list. If more people request > is, then we'll bump it up. I will be sure to add your email to the > ticket. > > > --Mary > > On Oct 10, 2012, at 7:22 AM, Alan Crouse wrote: > > > I have located the functions to change font characteristics > within a text string > (http://www.ncl.ucar.edu/Document/Graphics/function_code.shtml), but > changing color within a string does not seem to be an option. Am I > missing something? > > > > I am using ncl 6.1.0-beta. > > > > Thank you, > > > > Alan > > _______________________________________________ > > ncl-talk mailing list > > List instructions, subscriber options, unsubscribe: > > http://mailman.ucar.edu/mailman/listinfo/ncl-talk > > > > ------------------------------ > > Message: 18 > Date: Wed, 10 Oct 2012 10:30:09 -0600 > From: Mary Haley <haley@ucar.edu <mailto:haley@ucar.edu>> > Subject: Re: Coordinate subscript type mismatch. > To: H.Dang <danghy@gmail.com <mailto:danghy@gmail.com>> > Cc: "ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> USERS" > <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: <F46B75D7-CB96-43B2-8C8F-B8D406744260@ucar.edu > <mailto:F46B75D7-CB96-43B2-8C8F-B8D406744260@ucar.edu>> > Content-Type: text/plain; charset=iso-8859-1 > > Hongyan, > > I don't think I've see this type of error before, so I don't fully > understand it. (Dave, can you help out here?) > > Line 168 is > > > 168 nox_modelEY_ave_40_60NS(ipressure) = > avg(vNOx({pre_MIPAS},{40:60},:)) > > I need to see what those variables look like: > > printVarSummary(nox_modelEY_ave_40_60NS) > printVarSummary(vNOx) > > My feeling is that there's an issue with the way you are > subscripting vNOx, so > I need to see what its coordinate arrays look like. > > --Mary > > On Oct 8, 2012, at 8:37 PM, H.Dang wrote: > > > Dear all, > > > > I got the following error information in my NCL code: > > > > fatal:Coordinate subscript type mismatch. Subscript (0) can not > be coerced to type of coordinate variable > > fatal:["Execute.c":7556]:Execute: Error occurred at or near line > 168 in file obs_mdl.ncl > > > > > > 161 do ipressure = 0,28 - 1 > > 162 pre_MIPAS = plev_MIPAS(ipressure) > > 163 pre_MIPAS = dble2flt(plev_MIPAS(ipressure)) > > 164 > > 165 printVarSummary(pre_MIPAS) > > 166 printVarSummary(ple40) > > 167 ; > > 168 nox_modelEY_ave_40_60NS(ipressure) = > avg(vNOx({pre_MIPAS},{40:60},:)) > > 179 end do > > > > The summary of the variables are: > > > > Variable: pre_MIPAS > > Type: double > > Total Size: 8 bytes > > 1 values > > Number of Dimensions: 1 > > Dimensions and sizes: [1] > > Coordinates: > > Number Of Attributes: 7 > > typeConversion_op_ncl : double converted to float > > units : hPa > > long_name : pressure > > axis : Z > > positive : down > > standard_name : air_pressure > > plev : 300 > > > > Variable: ple40 > > Type: float > > Total Size: 160 bytes > > 40 values > > Number of Dimensions: 1 > > Dimensions and sizes: [ple | 40] > > Coordinates: > > ple: [964..0.1] > > Number Of Attributes: 1 > > units : mb > > > > I though the error is caused by the variable type, but when I > used "dble2flt" to convert the type from double to float, the error > message is till the same. So is my idea of converting correct or > not? How should I fix this problem? Thank you! > > > > Hongyan > > > > _______________________________________________ > > ncl-talk mailing list > > List instructions, subscriber options, unsubscribe: > > http://mailman.ucar.edu/mailman/listinfo/ncl-talk > > > > ------------------------------ > > Message: 19 > Date: Wed, 10 Oct 2012 18:40:27 +0200 > From: virginie hergault <virginie.hergault@gmail.com > <mailto:virginie.hergault@gmail.com>> > Subject: date > To: ncl-talk <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: > > <CABtrY=pNCHsrUzJs1PNKtcYKcfrrXxmnv2XcouM+C-T5FUGZ5g@mail.gmail.com > <mailto:pNCHsrUzJs1PNKtcYKcfrrXxmnv2XcouM%2BC-T5FUGZ5g@mail.gmail.com>> > Content-Type: text/plain; charset=ISO-8859-1 > > Hi users, > > I would like to know how to convert a date YYYYMMDD to the day of the > year (between 1 and 365) ? > Is there a function available in ncl for this? > > Thank you, > Virginie > > > ------------------------------ > > Message: 20 > Date: Thu, 11 Oct 2012 00:41:44 +0800 > From: ?yan <sunyan_10@hotmail.com <mailto:sunyan_10@hotmail.com>> > Subject: Re: cmip5 ?lon2d,lat2d?data > To: "Mary Haley" <haley@ucar.edu <mailto:haley@ucar.edu>>, > <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: <SNT124-DS1966A785D6C7B9902809AF8F8E0@phx.gbl> > Content-Type: text/plain; format=flowed; charset="gb2312"; > reply-type=original > > thank you,Mary Haley. I had downloaded the newer version of the > "ESMF_regridding.ncl" script,I am using the newer version, this problem > exists under the newer version "ESMF_regridding.ncl" script. > I don't know why,thank you very much! > > -----????----- > From: Mary Haley > Sent: Thursday, October 11, 2012 12:21 AM > To: ?yan > Cc: ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> > Subject: Re: cmip5 ?lon2d,lat2d?data > > Hi, > > Did you download a newer version of the "ESMF_regridding.ncl" script > from: > > http://www.ncl.ucar.edu/Applications/Files/ESMF_regridding.ncl > > Once you do this, you either need to copy it to > $NCARG_ROOT/lib/ncarg/nclscripts/esmf/, > or, if you don't have permission to do this, then put it in the same > directory as > your script, and edit your script to load this script instead: > > ; load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl" > load "./ESMF_regridding.ncl" > > --Mary > > > On Oct 10, 2012, at 1:24 AM, ?yan wrote: > > > thank you very much,Mary Haley . your suggest is right. > > I see example ESMF_regrid_6.ncl. in ESMF_regrid_6.ncl,the time > dimension > > is > > only one,it works well.when I use the whole time dimensions,new > problem > > appears: > > > > my ncl script is: > > 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" > > load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl" > > > > begin > > ;---Interpolation methods > > methods = (/"patch"/) > > > > ;---Input file > > srcFileName = > "tos_Omon_MPI-ESM-P_historical_r1i1p1_185001-200512.nc" > > wgtFile = "CMIP5_2_World_" + methods + ".nc" > > > > ;---Get data and lat/lon grid from CMIP5 Grid > > sfile = addfile(srcFileName,"r") > > thetao = sfile->tos(0,:,:) > > thetao@lat2d = sfile->lat > > thetao@lon2d = sfile->lon > > printVarSummary(thetao) > > > > Opt = True > > Opt@SrcFileName = "CMIP5_SCRIP.nc" ; source file name > > Opt@DstFileName = "World1deg_SCRIP.nc" ; destination > file name > > Opt@ForceOverwrite = True > > > > Opt@SrcGridCornerLat = sfile->lat_vertices ; corners are > necessary > > Opt@SrcGridCornerLon = sfile->lon_vertices ; for "conserve" > method > > Opt@SrcMask2D = where(.not.ismissing(thetao),1,0) > > Opt@DstGridType = "1x1" ; Destination grid > > Opt@DstTitle = "World Grid 1-degree Resolution" > > Opt@DstLLCorner = (/-89.75d, 0.00d /) > > Opt@DstURCorner = (/ 89.75d, 359.75d /) > > > > print("Generating interpolation weights from CMIP5 to") > > print("World 1 degree grid using the " + methods + " method.") > > > > Opt@WgtFileName = wgtFile > > Opt@InterpMethod = methods > > > > > ;---------------------------------------------------------------------- > > ; Interpolate data from CMIP5 to World 1-degree grid. > > > ;---------------------------------------------------------------------- > > thetao_regrid = ESMF_regrid(thetao,Opt) > > printVarSummary(thetao_regrid) > > > > end > > > > > > fatal:Eq: Dimension size, for dimension number 0, of operands > does not > > match, can't continue > > > > fatal:["Execute.c":7556]:Execute: Error occurred at or near line > 1071 in > > file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl > > > > fatal:["Execute.c":7556]:Execute: Error occurred at or near line > 2810 in > > file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl > > > > fatal:["Execute.c":7556]:Execute: Error occurred at or near line > 2921 in > > file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl > > > > I don't know what the problem is. thank you. > > ------------------------------------------------------------------- > > -----????----- > > From: Mary Haley > > Sent: Monday, October 01, 2012 1:49 AM > > To: ?yan > > Cc: ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> > > Subject: Re: cmip5 ?lon2d,lat2d?data > > > > > > On Sep 27, 2012, at 6:28 PM, ?yan wrote: > > > >> hello, > >> (1) I use cmip5 data, many datas are ?lon2d,lat2d?data.when I use > >> ?rcm2rgrid? function to regrid,some region is FillValue.why? > >> (2) how to deal with ?lon2d,lat2d?data? > >> thank you > >> > > > > There's not enough information to help you. It doesn't surprise > me that > > you > > would have missing values, especially if the rectilinear grid > goes outside > > the range of your curvilinear grid. > > > > Did you print out any values, or try plotting the regridded data > to see if > > it looks correct? > > > > You can also visit our ESMF regridding page. See example > > ESMF_regrid_6.ncl: > > > > http://www.ncl.ucar.edu/Applications/ESMF.shtml#ex6 > > > > --Mary > > > >> > >> ------------------------------------------------------------------ > >> > >> National Key Laboratory of Atmospheric Sciences and > >> Geophysical Fluid Dynamics (LASG), > >> Institute of Atmospheric Physics (IAP), > >> Chinese Academy of Sciences (CAS), > >> > >> ------------------------------------------------------------------- > >> _______________________________________________ > >> 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 > > > > ------------------------------ > > Message: 21 > Date: Wed, 10 Oct 2012 10:49:06 -0600 > From: Mary Haley <haley@ucar.edu <mailto:haley@ucar.edu>> > Subject: Re: cmip5 ?lon2d,lat2d?data > To: ?yan <sunyan_10@hotmail.com <mailto:sunyan_10@hotmail.com>> > Cc: ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> > Message-ID: <8E9F57F2-8D73-47C0-8628-2CE96298C008@ucar.edu > <mailto:8E9F57F2-8D73-47C0-8628-2CE96298C008@ucar.edu>> > Content-Type: text/plain; charset=GB2312 > > Can you send me the error messages you are getting from using the > new file, so I can see what > line numbers the errors occur on. > > Or, if you can provide your dataset, that would help. > > http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP > > Thanks, > > --Mary > > On Oct 10, 2012, at 10:41 AM, ?yan wrote: > > > thank you,Mary Haley. I had downloaded the newer version of the > > "ESMF_regridding.ncl" script,I am using the newer version, this > problem > > exists under the newer version "ESMF_regridding.ncl" script. > > I don't know why,thank you very much! > > > > -----????----- > > From: Mary Haley > > Sent: Thursday, October 11, 2012 12:21 AM > > To: ?yan > > Cc: ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> > > Subject: Re: cmip5 ?lon2d,lat2d?data > > > > Hi, > > > > Did you download a newer version of the "ESMF_regridding.ncl" > script from: > > > > http://www.ncl.ucar.edu/Applications/Files/ESMF_regridding.ncl > > > > Once you do this, you either need to copy it to > > $NCARG_ROOT/lib/ncarg/nclscripts/esmf/, > > or, if you don't have permission to do this, then put it in the same > > directory as > > your script, and edit your script to load this script instead: > > > > ; load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl" > > load "./ESMF_regridding.ncl" > > > > --Mary > > > > > > On Oct 10, 2012, at 1:24 AM, ?yan wrote: > > > >> thank you very much,Mary Haley . your suggest is right. > >> I see example ESMF_regrid_6.ncl. in ESMF_regrid_6.ncl,the time > dimension > >> is > >> only one,it works well.when I use the whole time dimensions,new > problem > >> appears: > >> > >> my ncl script is: > >> 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" > >> load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl" > >> > >> begin > >> ;---Interpolation methods > >> methods = (/"patch"/) > >> > >> ;---Input file > >> srcFileName = > "tos_Omon_MPI-ESM-P_historical_r1i1p1_185001-200512.nc" > >> wgtFile = "CMIP5_2_World_" + methods + ".nc" > >> > >> ;---Get data and lat/lon grid from CMIP5 Grid > >> sfile = addfile(srcFileName,"r") > >> thetao = sfile->tos(0,:,:) > >> thetao@lat2d = sfile->lat > >> thetao@lon2d = sfile->lon > >> printVarSummary(thetao) > >> > >> Opt = True > >> Opt@SrcFileName = "CMIP5_SCRIP.nc" ; source file name > >> Opt@DstFileName = "World1deg_SCRIP.nc" ; destination > file name > >> Opt@ForceOverwrite = True > >> > >> Opt@SrcGridCornerLat = sfile->lat_vertices ; corners are > necessary > >> Opt@SrcGridCornerLon = sfile->lon_vertices ; for "conserve" > method > >> Opt@SrcMask2D = where(.not.ismissing(thetao),1,0) > >> Opt@DstGridType = "1x1" ; Destination grid > >> Opt@DstTitle = "World Grid 1-degree Resolution" > >> Opt@DstLLCorner = (/-89.75d, 0.00d /) > >> Opt@DstURCorner = (/ 89.75d, 359.75d /) > >> > >> print("Generating interpolation weights from CMIP5 to") > >> print("World 1 degree grid using the " + methods + " method.") > >> > >> Opt@WgtFileName = wgtFile > >> Opt@InterpMethod = methods > >> > >> > ;---------------------------------------------------------------------- > >> ; Interpolate data from CMIP5 to World 1-degree grid. > >> > ;---------------------------------------------------------------------- > >> thetao_regrid = ESMF_regrid(thetao,Opt) > >> printVarSummary(thetao_regrid) > >> > >> end > >> > >> > >> fatal:Eq: Dimension size, for dimension number 0, of operands > does not > >> match, can't continue > >> > >> fatal:["Execute.c":7556]:Execute: Error occurred at or near line > 1071 in > >> file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl > >> > >> fatal:["Execute.c":7556]:Execute: Error occurred at or near line > 2810 in > >> file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl > >> > >> fatal:["Execute.c":7556]:Execute: Error occurred at or near line > 2921 in > >> file $NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl > >> > >> I don't know what the problem is. thank you. > >> ------------------------------------------------------------------- > >> -----????----- > >> From: Mary Haley > >> Sent: Monday, October 01, 2012 1:49 AM > >> To: ?yan > >> Cc: ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> > >> Subject: Re: cmip5 ?lon2d,lat2d?data > >> > >> > >> On Sep 27, 2012, at 6:28 PM, ?yan wrote: > >> > >>> hello, > >>> (1) I use cmip5 data, many datas are ?lon2d,lat2d?data.when > I use > >>> ?rcm2rgrid? function to regrid,some region is FillValue.why? > >>> (2) how to deal with ?lon2d,lat2d?data? > >>> thank you > >>> > >> > >> There's not enough information to help you. It doesn't surprise > me that > >> you > >> would have missing values, especially if the rectilinear grid > goes outside > >> the range of your curvilinear grid. > >> > >> Did you print out any values, or try plotting the regridded data > to see if > >> it looks correct? > >> > >> You can also visit our ESMF regridding page. See example > >> ESMF_regrid_6.ncl: > >> > >> http://www.ncl.ucar.edu/Applications/ESMF.shtml#ex6 > >> > >> --Mary > >> > >>> > >>> ------------------------------------------------------------------ > >>> > >>> National Key Laboratory of Atmospheric Sciences and > >>> Geophysical Fluid Dynamics (LASG), > >>> Institute of Atmospheric Physics (IAP), > >>> Chinese Academy of Sciences (CAS), > >>> > >>> ------------------------------------------------------------------- > >>> _______________________________________________ > >>> 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 > > > > ------------------------------ > > Message: 22 > Date: Wed, 10 Oct 2012 10:59:22 -0600 > From: Mary Haley <haley@ucar.edu <mailto:haley@ucar.edu>> > Subject: Re: date > To: virginie hergault <virginie.hergault@gmail.com > <mailto:virginie.hergault@gmail.com>> > Cc: ncl-talk <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: <FE2BEF1D-D6A0-4E07-BC1F-5A57F2C489B3@ucar.edu > <mailto:FE2BEF1D-D6A0-4E07-BC1F-5A57F2C489B3@ucar.edu>> > Content-Type: text/plain; charset="us-ascii" > > You can use "day_of_year", but you need to split the value into > separate year, month, and day integers. > > http://www.ncl.ucar.edu/Document/Functions/Built-in/day_of_year.shtml > > Here's some sample code: > > d = 20120305 ; assuming "d" is an int or a long > > year = toint(d/10000) > month = toint((d-(year*10000))/100) > day = toint(d-(year*10000)-(month*100)) > > print("year = " + year) > print("month = " + month) > print("day = " + day) > > doy = day_of_year(year,month,day) > > print("doy = " + doy) > > > If YYYYMMDD is a string, then convert it to an integer first using > "toint". > > --Mary > > On Oct 10, 2012, at 10:40 AM, virginie hergault wrote: > > > Hi users, > > > > I would like to know how to convert a date YYYYMMDD to the day of the > > year (between 1 and 365) ? > > Is there a function available in ncl for this? > > > > Thank you, > > Virginie > > _______________________________________________ > > ncl-talk mailing list > > List instructions, subscriber options, unsubscribe: > > http://mailman.ucar.edu/mailman/listinfo/ncl-talk > > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: > http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20121010/eb2adf14/attachment.html > > ------------------------------ > > Message: 23 > Date: Wed, 10 Oct 2012 11:13:12 -0600 > From: Dennis Shea <shea@ucar.edu <mailto:shea@ucar.edu>> > Subject: Re: date > To: virginie hergault <virginie.hergault@gmail.com > <mailto:virginie.hergault@gmail.com>> > Cc: ncl-talk <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: <5075ACA8.5050001@ucar.edu > <mailto:5075ACA8.5050001@ucar.edu>> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > I styerted to write this before MH's response > > ==== > Please look at the NCL function list. > > All of NCL's date function are summarized under the "Date" category. > > http://www.ncl.ucar.edu/ > > (1) Click "Functions" and 'Category Listing' > (2) Click "Date" > > There is no 'yyyymmdd_to_day_of_year'. > There can not be a function for everything! > However, it is trivial to create your own. > > undef ("ymd2doy") > function ymd2doy(yyyymmdd:integer) > local yyyy, mm, dd > begin > yyyy = yyyymmdd/10000 > mm = (yyyymmdd - (yyyy*10000))/100 > dd = yyyymmdd%100 > return( day_of_year(yyyy,mm,dd) ) > end > > On 10/10/12 10:59 AM, Mary Haley wrote: > > You can use "day_of_year", but you need to split the value into > separate > > year, month, and day integers. > > > > http://www.ncl.ucar.edu/Document/Functions/Built-in/day_of_year.shtml > > > > Here's some sample code: > > > > d = 20120305 ; assuming "d" is an int or a long > > > > year = toint(d/10000) > > month = toint((d-(year*10000))/100) > > day = toint(d-(year*10000)-(month*100)) > > > > print("year = " + year) > > print("month = " + month) > > print("day = " + day) > > > > doy = day_of_year(year,month,day) > > > > print("doy = " + doy) > > > > > > If YYYYMMDD is a string, then convert it to an integer first > using "toint". > > > > --Mary > > > > On Oct 10, 2012, at 10:40 AM, virginie hergault wrote: > > > >> Hi users, > >> > >> I would like to know how to convert a date YYYYMMDD to the day > of the > >> year (between 1 and 365) ? > >> Is there a function available in ncl for this? > >> > >> Thank you, > >> Virginie > >> _______________________________________________ > >> 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 > > > > > ------------------------------ > > Message: 24 > Date: Wed, 10 Oct 2012 11:15:41 -0600 > From: David Brown <dbrown@ucar.edu <mailto:dbrown@ucar.edu>> > Subject: Re: Coordinate subscript type mismatch. > To: "ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> USERS" > <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>> > Message-ID: <F4C83608-E6D6-4895-9A99-EFA631DABA0D@ucar.edu > <mailto:F4C83608-E6D6-4895-9A99-EFA631DABA0D@ucar.edu>> > Content-Type: text/plain; charset=us-ascii > > I think the problem is that you initially create the variable > 'pre_MIPAS' as a double with the statement: > > pre_MIPAS = plev_MIPAS(ipressure) > > In NCL once a variable has been created its type cannot be changed, > so the next statement: > > pre_MIPAS = dble2flt(plev_MIPAS(ipressure)) > > does indeed create a float value on the right side of the > assignment, but then it is converted using normal type coercion back > into a double for > assignment to a pre-existing double variable. (See > http://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclDataTypes.shtml#Coercion.) > > The cure, I think, is simply to eliminate the first statement: > > pre_MIPAS = plev_MIPAS(ipressure) > > so that you are creating 'pre_MIPAS' as a float in the first place. > -dave > > > On Oct 10, 2012, at 10:30 AM, Mary Haley wrote: > > > Hongyan, > > > > I don't think I've see this type of error before, so I don't > fully understand it. (Dave, can you help out here?) > > > > Line 168 is > > > >> 168 nox_modelEY_ave_40_60NS(ipressure) = > avg(vNOx({pre_MIPAS},{40:60},:)) > > > > I need to see what those variables look like: > > > > printVarSummary(nox_modelEY_ave_40_60NS) > > printVarSummary(vNOx) > > > > My feeling is that there's an issue with the way you are > subscripting vNOx, so > > I need to see what its coordinate arrays look like. > > > > --Mary > > > > On Oct 8, 2012, at 8:37 PM, H.Dang wrote: > > > >> Dear all, > >> > >> I got the following error information in my NCL code: > >> > >> fatal:Coordinate subscript type mismatch. Subscript (0) can not > be coerced to type of coordinate variable > >> fatal:["Execute.c":7556]:Execute: Error occurred at or near line > 168 in file obs_mdl.ncl > >> > >> > >> 161 do ipressure = 0,28 - 1 > >> 162 pre_MIPAS = plev_MIPAS(ipressure) > >> 163 pre_MIPAS = dble2flt(plev_MIPAS(ipressure)) > >> 164 > >> 165 printVarSummary(pre_MIPAS) > >> 166 printVarSummary(ple40) > >> 167 ; > >> 168 nox_modelEY_ave_40_60NS(ipressure) = > avg(vNOx({pre_MIPAS},{40:60},:)) > >> 179 end do > >> > >> The summary of the variables are: > >> > >> Variable: pre_MIPAS > >> Type: double > >> Total Size: 8 bytes > >> 1 values > >> Number of Dimensions: 1 > >> Dimensions and sizes: [1] > >> Coordinates: > >> Number Of Attributes: 7 > >> typeConversion_op_ncl : double converted to float > >> units : hPa > >> long_name : pressure > >> axis : Z > >> positive : down > >> standard_name : air_pressure > >> plev : 300 > >> > >> Variable: ple40 > >> Type: float > >> Total Size: 160 bytes > >> 40 values > >> Number of Dimensions: 1 > >> Dimensions and sizes: [ple | 40] > >> Coordinates: > >> ple: [964..0.1] > >> Number Of Attributes: 1 > >> units : mb > >> > >> I though the error is caused by the variable type, but when I > used "dble2flt" to convert the type from double to float, the error > message is till the same. So is my idea of converting correct or > not? How should I fix this problem? Thank you! > >> > >> Hongyan > >> > >> _______________________________________________ > >> 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 > ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk > > > End of ncl-talk Digest, Vol 107, Issue 13 > ***************************************** > > > > > _______________________________________________ > 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-talkReceived on Thu Oct 11 07:21:00 2012
This archive was generated by hypermail 2.1.8 : Fri Oct 12 2012 - 15:38:19 MDT