Re: help

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