Re: give tick marks on a longitude-time plot

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
Received on Wed Nov 6 10:32:28 2013