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>
> Send ncl-talk mailing list submissions to
> 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
>
> You can reach the person managing the list at
> 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>
> Subject: XY Plots over a Map
> To: Ncl Talk <ncl-talk@ucar.edu>
> Message-ID: <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>
> Subject: Re: XY Plots over a Map
> To: ncl-talk@ucar.edu
> Message-ID: <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
> 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>
> Subject: about Taylor diagram ...
> To: NCL USERS <ncl-talk@ucar.edu>
> Message-ID: <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>
> Subject: Re: Track one storm when there are also two other
> storms in the domain
> To: Yongxin Zhang <yongxin.fred@gmail.com>
> Cc: ncl-talk <ncl-talk@ucar.edu>, wrfhelp <wrfhelp@ucar.edu>
> Message-ID: <50748F12.20107@ucar.edu>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> WRF questions should be sent to 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>
> Subject: Re: Track one storm when there are also two other
> storms in the domain
> To: Yongxin Zhang <yongxin.fred@gmail.com>
> Cc: ncl-talk <ncl-talk@ucar.edu>, wrfhelp <wrfhelp@ucar.edu>
> Message-ID: <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
> >
> > ---
> >
> > 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>
> Subject: EOF question
> To: ncl-talk <ncl-talk@ucar.edu>
> Message-ID:
> <
> 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>
> Subject: A question about resource tmXBMode
> To: ncl-talk@ucar.edu
> Message-ID:
> <CAA=sSiyXVwVtT=O1NPeFUW=o77zVL+RNxNQmtJ=A=
> 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>
> Subject: Re: cmip5 ?lon2d,lat2d?data
> To: "Mary Haley" <haley@ucar.edu>, <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
> 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>
> Subject: register
> To: ncl-talk@ucar.edu
> Message-ID:
> <CAFxT1sbE_iO-_nCJ8=Niff5552WzcB2wo=
> 4t7PvWNBF8gtWeVQ@mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> i want to register to the the 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>
> Subject: Change color within text string?
> To: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID:
> <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>
> Subject: stat_dispersion missing values
> To: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID:
> <
> 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>
> Subject: Re: about Taylor diagram ...
> To: "Ufuk Utku Turuncoglu (BE)" <u.utku.turuncoglu@be.itu.edu.tr>
> Cc: NCL USERS <ncl-talk@ucar.edu>
> Message-ID: <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>
> Subject: Re: EOF question
> To: Shinn Wong <shinnshinnwong@gmail.com>
> Cc: ncl-talk <ncl-talk@ucar.edu>
> Message-ID: <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>
> Subject: Re: stat_dispersion missing values
> To: "Schiffer, Nicole June" <nschiff2@illinois.edu>
> Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID: <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>
> Subject: Re: A question about resource tmXBMode
> To: Zheng Lu <cfluzheng@gmail.com>
> Cc: ncl-talk@ucar.edu
> Message-ID: <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>
> Subject: Re: cmip5 ?lon2d,lat2d?data
> To: ?yan <sunyan_10@hotmail.com>
> Cc: ncl-talk@ucar.edu
> Message-ID: <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
> > 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>
> Subject: Re: Change color within text string?
> To: Alan Crouse <ac_soaring@yahoo.com>
> Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Message-ID: <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>
> Subject: Re: Coordinate subscript type mismatch.
> To: H.Dang <danghy@gmail.com>
> Cc: "ncl-talk@ucar.edu USERS" <ncl-talk@ucar.edu>
> Message-ID: <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>
> Subject: date
> To: ncl-talk <ncl-talk@ucar.edu>
> Message-ID:
> <CABtrY=
> pNCHsrUzJs1PNKtcYKcfrrXxmnv2XcouM+C-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>
> Subject: Re: cmip5 ?lon2d,lat2d?data
> To: "Mary Haley" <haley@ucar.edu>, <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
> 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
> > 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>
> Subject: Re: cmip5 ?lon2d,lat2d?data
> To: ?yan <sunyan_10@hotmail.com>
> Cc: ncl-talk@ucar.edu
> Message-ID: <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
> > 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
> >> 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>
> Subject: Re: date
> To: virginie hergault <virginie.hergault@gmail.com>
> Cc: ncl-talk <ncl-talk@ucar.edu>
> Message-ID: <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>
> Subject: Re: date
> To: virginie hergault <virginie.hergault@gmail.com>
> Cc: ncl-talk <ncl-talk@ucar.edu>
> Message-ID: <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>
> Subject: Re: Coordinate subscript type mismatch.
> To: "ncl-talk@ucar.edu USERS" <ncl-talk@ucar.edu>
> Message-ID: <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
> 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
Received on Thu Oct 11 03:52:19 2012
This archive was generated by hypermail 2.1.8 : Fri Oct 12 2012 - 15:38:19 MDT