Re: Execute: Error occurred

[1]
xy_5.ncl is a plot of a time series index: time[:], dsoi[:]

dsoi is an index that has a range of -3 to +3

[2]
You are trying to plot:
      time[*], u(:,{20:50},{30:120})
the *actual* precipitation values. There are *no* negative values.

You have daily values from one year:
    time = UNLIMITED ; // (366 currently)

Where have you calculated anomalies?
Further, anomalies from what reference?

[3]
It might help if your read the Mini-Language and Graphics manuals at:

        http://test.www.ncl.ucar.edu/Document/Manuals/

Good luck

Kim wrote:
> Hi,
>
> Thank you for your reply but I face another problem given below. Please
> could you guide me how to fix it? Actually I want to plot it as like
> example (xy_5.ncl
> <http://www.ncl.ucar.edu/Applications/Scripts/xy_5.ncl>) & my code is
> also attached, any help will be appreciated.
>
> Thanks again
>
> ‘’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’
>
> fatal:Number of subscripts do not match number of dimensions of
> variable,(2) Subscripts used, (3) Subscripts expected
>
> fatal:Execute: Error occurred at or near line 3922 in file
> $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl
>
>
>
> fatal:Execute: Error occurred at or near line 38 in file x1.ncl
>
> ‘’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’’
>
> netcdf precip.mon.mean {
>
> dimensions:
>
> lat = 72 ;
>
> lon = 144 ;
>
> time = UNLIMITED ; // (366 currently)
>
> variables:
>
> float lat(lat) ;
>
> lat:units = "degrees_north" ;
>
> lat:actual_range = 88.75f, -88.75f ;
>
> lat:long_name = "Latitude" ;
>
> float lon(lon) ;
>
> lon:units = "degrees_east" ;
>
> lon:long_name = "Longitude" ;
>
> lon:actual_range = 1.25f, 358.75f ;
>
> float precip(time, lat, lon) ;
>
> precip:long_name = "Average Monthly Rate of
> Precipitation" ;
>
> precip:valid_range = 0.f, 100.f ;
>
> precip:units = "mm/day" ;
>
> precip:add_offset = 0.f ;
>
> precip:scale_factor = 1.f ;
>
> precip:actual_range = 0.f, 85.967f ;
>
> precip:missing_value = -9.96921e+36f ;
>
> precip:precision = 32767s ;
>
> precip:least_significant_digit = 2s ;
>
> precip:var_desc = "Precipitation" ;
>
> precip:dataset = "GPCP Version 2.1 Combined
> Precipitation Dataset" ;
>
> precip:level_desc = "Surface" ;
>
> precip:statistic = "Mean" ;
>
> precip:parent_stat = "Mean" ;
>
> double time(time) ;
>
> time:units = "days since 1700-1-1 00:00:0.0" ;
>
> time:long_name = "Time" ;
>
> time:delta_t = "0000-01-00 00:00:00" ;
>
> time:actual_range = 101902., 113011. ;
>
> time:avg_period = "0000-01-00 00:00:00" ;
>
>
>
> // global attributes:
>
> :Conventions = "COARDS" ;
>
> :title = "GPCP Version 2.1 Combined Precipitation
> Dataset" ;
>
> :history = "created jul 2009 by CAS at NOAA/ESRL PSD" ;
>
> :platform = "Observation" ;
>
> :source = "GPCP Polar Satellite Precipitation Data
> Centre - Emission (SSM/I emission estimates).\n",
>
>
> --- On *Sun, 1/31/10, Dennis Shea /<shea@ucar.edu>/* wrote:
>
>
> From: Dennis Shea <shea@ucar.edu>
> Subject: Re: [ncl-talk] Execute: Error occurred
> To: "Kim" <r4rid@yahoo.com>
> Cc: ncl-talk@ucar.edu
> Date: Sunday, January 31, 2010, 6:15 AM
>
> You do not indicate the units of "time".
> You are treating this as "yyyymm". I speculate this is not correct.
>
> The 1st rule of data processing is "look at the data".
> If you use 'ncdump -h' or 'ncl_filedum' on the file, I speculate
> you will see:
>
> double time(time) ;
> time:units = "hours since 1-1-1 00:00:0.0" ;
>
> You could use
> time = f->time
> yyyymm = ut_calendar(time, 1)
> yrfrac = yyyymm_to_yyyyfrac(yyyymm, 0.0) ; x (time) axis
>
> Please look at:
>
> http://www.ncl.ucar.edu/Document/Functions/index.shtml
>
> the "Date" category
>
> Good luck
>
> Kim wrote:
> > Dear ncl users, anyone could help me to fix the problem given below.
> >
> > Thanks in advance,
> >
>[snip]
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Feb 1 07:24:51 2010