Re: give tick marks on a longitude-time plot

You can have both, something like ...

  yrStrt = 1981
  yrLast = 1989
  year = ispan(yrStrt, yrLast,1) ; integer
  yearDir= year + "/" ; string
  print(year)
  print(yearDir)

On 11/7/13 4:08 PM, Yi-Chih Huang wrote:
> Mary,
>
> Because year is the name of directories in my current working
> environment, "/" and string setting are necessary to work on directory in
> the script as follows.
>
> dir1= "/fs3/SysE_DB/saji/OISST/Monthly/Mean/"
> dir2=
> (/"1981/","1982/","1983/","1984/","1985/","1986/","1987/","1988/","1989/",\
>
> "1990/","1991/","1992/","1993/","1994/","1995/","1996/","1997/","1998/","1999/","2000/",\
> "2001/","2002/","2003/","2004/","2005/","2006/","2007/","2008/"/)
> dir = str_concat(dir1+dir2(0))
>
> As for the ind statement below, I don't understand why it is "yStrt =
> (/"1981/","1985/","1996/"/)", rather than "yStrt =
> (/"1982/","1986/","1997/"/)" to select the year 1982-83, 1986-87, and
> 1997-98. What do you think of this?
>
> 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 Fri, Nov 8, 2013 at 12:17 AM, Mary Haley <haley@ucar.edu> wrote:
>
>> 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
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Nov 7 16:21:09 2013