From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>

Date: Thu Nov 07 2013 - 08:17:44 MST

Date: Thu Nov 07 2013 - 08:17:44 MST

Yi-Chih,

These two lines really don't make sense:

*> yStrt = (/"1981/","1985/","1996/"/)
*

*> yLast = (/"1983/","1987/","1998/"/)
*

Why do you have these as strings and not integers? Also, why do you have the extra "/" after the year values?

I think these two lines should be integers, because they eventually are being compared against float values, not string values.

*> yStrt = (/1981,1985,1996/)
*

*> yLast = (/1983,1987,1998/)
*

--Mary

On Nov 6, 2013, at 6:43 PM, Yi-Chih Huang <dscpln@gmail.com> wrote:

*> Mary,
*

*>
*

*> I used pcps!0= "year" to remove the warning. Besides, I want to choose 1982-83, 1986-87, and 1997-98. At first I believed the command yStrt = (/"1982/","1986/","1997/"/) would have chosen the years 1982-83, 1986-87, and 1997-98. However, it is yStrt = (/"1981/","1985/","1996/"/) which was working. Could you explain why yStrt = (/"1982/","1986/","1997/"/) only chose the year 1983?
*

*>
*

*> yStrt = (/"1981/","1985/","1996/"/)
*

*> yLast = (/"1983/","1987/","1998/"/)
*

*> iYr = ind(YMDHMS(:,0).ge.yStrt(ny) .and. YMDHMS(:,0).le.yLast(ny))
*

*>
*

*> Thanks much,
*

*>
*

*> Yi-Chih
*

*>
*

*>
*

*> On Thu, Nov 7, 2013 at 2:32 AM, Mary Haley <haley@ucar.edu> wrote:
*

*> You are getting these warnings because you are using a _Wrap function on "pcps", which has no metadata.
*

*>
*

*> The _Wrap functions are supposed to be used when you have metadata that you want to retain.
*

*>
*

*> I believe when you do something like this:
*

*>
*

*> > pcps(ny,:,:)= dim_avg_n_Wrap( preci,1 ) ; (time,lat)
*

*>
*

*> since the variable "pcps" already exists, and since you are using it in a subscripting fashion, NCL will not copy over any metadata.
*

*>
*

*> The warning is probably harmless, but you can also just use "dim_avg_n", and then use copy_VarMeta to copy metadata to pcps later.
*

*>
*

*> --Mary
*

*>
*

*> On Nov 5, 2013, at 1:01 AM, Yi-Chih Huang <dscpln@gmail.com> wrote:
*

*>
*

*> > Hello,
*

*> >
*

*> > I modified Dennis' script to calculate anomies. I don't understand why there is "Coordinate variables" problem, especially I copy variable coordinate to the variable in the Hovmoller diagram and didn't change the commands about time axis. In addition, I executed "dimsizes(pcps)". Why does the warning "warning:Dimension (0) has not been defined" still show up?
*

*> >
*

*> > Thanks,
*

*> >
*

*> > Yi-Chih
*

*> >
*

*> > ########
*

*> > yhuang@arc-env:/fs3/yhuang/scripts$ ncl pcpSstTime.ncl
*

*> > Copyright (C) 1995-2013 - All Rights Reserved
*

*> > University Corporation for Atmospheric Research
*

*> > NCAR Command Language Version 6.1.2
*

*> > The use of this software is governed by a License Agreement.
*

*> > See http://www.ncl.ucar.edu/ for more details.
*

*> > (0) a
*

*> > warning:Dimension (0) has not been defined
*

*> > (0) b
*

*> > fatal:Coordinate variables must be the same dimension as their dimension
*

*> > fatal:No coordinate variable exists for dimension (time) in variable (dpcp)
*

*> > fatal:["Execute.c":8128]:Execute: Error occurred at or near line 89 in file pcpSstTime.ncl
*

*> >
*

*> > ###
*

*> > load "$SysE/lib/ncl/helper_libs.ncl"
*

*> >
*

*> > begin
*

*> > latS = -5.0
*

*> > latN = 5.0
*

*> > lonE = 120.0
*

*> > lonW = 280.0
*

*> >
*

*> > yrStrt = 1979
*

*> > yrLast = 2009
*

*> >
*

*> > in = addfile("/fs3/yhuang/precip.mon.mean.nc","r") ; 197901-201306
*

*> > YMDHMS = cd_calendar(in->time, 0) ; all times
*

*> > ; index for desired time
*

*> > iYear = ind(YMDHMS(:,0).ge.yrStrt .and. YMDHMS(:,0).le.yrLast)
*

*> > ; read time and region
*

*> > precip = in->precip(iYear,{latS:latN},{lonE:lonW})
*

*> > pcpClm = clmMonTLL( precip )
*

*> >
*

*> > xt= precip(0:23,{latS:latN},{lonE:lonW})
*

*> > ds= dimsizes(xt)
*

*> > nyr= yrLast -yrStrt +1
*

*> > pcpm= new((/ds(0),ds(2)/),float)
*

*> >
*

*> > yyyymmM = cd_calendar(precip&time,-1) ; verify times are correct
*

*> >
*

*> > dimp = dimsizes(precip)
*

*> > ntim = dimp(0)
*

*> > nlat = dimp(1)
*

*> > mlon = dimp(2)
*

*> >
*

*> > pcpt = dim_avg_n_Wrap( pcpClm,1 ) ; (time,lat)
*

*> > do nm= 0,11
*

*> > pcpm(nm,:) = pcpt(nm,:)
*

*> > pcpm(nm+12,:)= pcpt(nm,:)
*

*> > end do
*

*> >
*

*> > yStrt = (/"1981/","1985/","1996/"/)
*

*> > yLast = (/"1983/","1987/","1998/"/)
*

*> > nyrs= dimsizes(yStrt)
*

*> >
*

*> > in= addfile("/fs3/yhuang/precip.mon.mean.nc","r") ; 197901-201306
*

*> > prec= in->precip
*

*> > xt= prec(0:23,{latS:latN},{lonE:lonW})
*

*> > ds= dimsizes(xt)
*

*> > pcps= new((/nyrs,ds(0),ds(2)/),float)
*

*> >
*

*> > do ny= 0,nyrs-1 ; precipitation in ENSO years
*

*> > iYr = ind(YMDHMS(:,0).ge.yStrt(ny) .and. YMDHMS(:,0).le.yLast(ny))
*

*> > ; read time and region
*

*> > preci = in->precip(iYr,{latS:latN},{lonE:lonW})
*

*> >
*

*> > yyyymmE = cd_calendar(preci&time,-1) ; verify times are correct
*

*> >
*

*> > dimp = dimsizes(preci)
*

*> > ntim = dimp(0)
*

*> > nlat = dimp(1)
*

*> > mlon = dimp(2)
*

*> >
*

*> > pcps(ny,:,:)= dim_avg_n_Wrap( preci,1 ) ; (time,lat)
*

*> > end do
*

*> >
*

*> > dimp = dimsizes(pcps)
*

*> > ntim = dimp(0)
*

*> > print("a")
*

*> > pcpE= dim_avg_n_Wrap( pcps, 0 )
*

*> > print("b")
*

*> > dpcp= pcpE - pcpm
*

*> > copy_VarCoords(pcpE,dpcp)
*

*> >
*

*> > wks= gsn_open_wks("X11","pcpa")
*

*> > gsn_define_colormap(wks,"BlueRed")
*

*> >
*

*> > res = True
*

*> > res@gsnAddCyclic = False
*

*> >
*

*> > res@cnFillOn = True
*

*> > res@cnLinesOn = False ; turn off the contour lines
*

*> > res@cnLineLabelsOn = False
*

*> > res@cnSpanFillPalette = True
*

*> >
*

*> > res@gsnSpreadColors = True ; use full colormap
*

*> > res@gsnSpreadColorStart = 10 ; for BlueRed
*

*> > res@gsnSpreadColorEnd = 250 ; for BlueRed
*

*> >
*

*> > yr= ispan(1,ntim,1)
*

*> > yr@long_name = "months"
*

*> > delete(dpcp&time)
*

*> > dpcp&time = yr
*

*> > plot = gsn_csm_hov(wks,dpcp,res)
*

*> > end
*

*> > _______________________________________________
*

*> > 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

Received on Thu Nov 7 08:17:52 2013

*
This archive was generated by hypermail 2.1.8
: Mon Nov 11 2013 - 09:45:34 MST
*