Re: help

[1]
The time varable must have a "units" attribute.

  float time(time) ;
  time:long_name = "days since 199601" ;

You provided a long_name but not a units attribute.

Make it: time@units = "days since 199601"

[2]
SDF is a GrADS related issue not NCL.

Good luck

On 10/25/10 1:00 AM, Mohamad Ibrahim OMAR wrote:
> Dear Shea
> I try your suggestion but the same problem didn't solve I try to remove
> the function dbl2flt from my script the interpolation netcdf file
> created but when I opened it by grads gave me the following:
> SDF file has no discernable time coordinate -- using default values.
> gadsdf: SDF file does not have any non-coordinate variables.
> and the ncdump -h to this file gived:
>
> netcdf GPCPprecip_199601-200412TEST {
> dimensions:
> time = UNLIMITED ; // (108 currently)
> lat = 62 ;
> lon = 58 ;
> variables:
> float pre(time, lat, lon) ;
> pre:time = 841512. ;
> pre:units = "mm/day" ;
> pre:long_name = "Average Monthly Rate of Precipitation" ;
> pre:_FillValue = -9999.f ;
> float lat2d(time, lat, lon) ;
> lat2d:_FillValue = 4.079829e-38f ;
> lat2d:time = 841512. ;
> lat2d:long_name = "2D Latitude" ;
> lat2d:units = "degrees north" ;
> lat2d:missing_value = 4.079829e-38f ;
> float lon2d(time, lat, lon) ;
> lon2d:_FillValue = 4.079829e-38f ;
> lon2d:time = 841512. ;
> lon2d:long_name = "2D Longitude" ;
> lon2d:units = "degrees east" ;
> lon2d:missing_value = 4.079829e-38f ;
> float time(time) ;
> time:long_name = "days since 199601" ;
> float lon(lon) ;
> lon:long_name = "Longitude" ;
> lon:units = "degrees_east" ;
> lon:actual_range = 23.55206f, 37.69477f ;
> float lat(lat) ;
> lat:long_name = "Latitude" ;
> lat:units = "degrees_north" ;
> lat:actual_range = 20.47314f, 33.98017f ;
> }
>
> so where is the problem?
> can you help me?
> I attached my interpolation script
> THANKS
> Mohamed OMAR
> Egyptian Meteorological Authority
> Scientific Research Department
> Climate Studies Centre
>
>
> --- On *Fri, 10/15/10, Dennis Shea /<shea@ucar.edu>/* wrote:
>
>
> From: Dennis Shea <shea@ucar.edu>
> Subject: Re: help
> To: "Mohamad Ibrahim OMAR" <egy_mohamad@yahoo.com>
> Cc: ncl-talk@ucar.edu
> Date: Friday, October 15, 2010, 9:31 PM
>
> The interpolation function requires that the latitude
> order be South-to-North. Your 'ncdump -h' seems to indicate that
> the data are ordered North-to-South
>
> 88.75f, -88.75f
>
>
> p = f->precip
> p = p(:,::-1,:) ; reverse the latitude order
>
> lat = p&lat
> print(lat)
>
>
>
> On 10/14/10 3:54 AM, Mohamad Ibrahim OMAR wrote:
> > I want to interpolate the precip.mon.mean.nc which has the ncdump as
> > following
> >
> > netcdf precip.mon.mean {
> > dimensions:
> > lat = 72 ;
> > lon = 144 ;
> > time = UNLIMITED ; // (369 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., 113103. ;
> > 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",
> > " GPCP Polar Satellite Precipitation Data Centre - Scattering (SSM/I
> > scattering estimates).\n",
> > " GPCP Geostationary Satellite Precipitation Data Centre (GPI and OPI
> > estimates and rain gauge analyses). NASA/GSFC Satellite Applications
> > Office (TOVS estimates).\n",
> > " GPCP Global Precipitation Climatology Centre (rain gauge
> analyses)\n",
> > " NASA ftp://precip.gsfc.nasa.gov/pub/gpcp-v2.1-v2/psg/" ;
> > :documentation =
> > "http://www.esrl.noaa.gov/psd/data/gridded/data.gpcp.html" ;
> > :version = "V2.1" ;
> > }
> >
> > to the regcm3 output file as I want to compare it with the other
> I did
> > this for CRU version 3 but when I did the sam with the GPCP I
> have the
> > following :
> >
> >
> > Segmentation fault (core dumped)
> > can anyone help me
> >
> > THANKS
> >
> > Mohamed OMAR
> > Egyptian Meteorological Authority
> > Scientific Research Department
> > Climate Studies Centre
> >
> >
> >
> >
> > _______________________________________________
> > 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 Mon Oct 25 06:45:15 2010