Re: help

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 Fri Oct 15 13:31:52 2010