Re: Convert data for another resolution

Hi Sara,

When using the addfile option "w", you cannot read a variable out of the
file with lat/lon dimensions 94x192 and write a variable back into the
file with the exact same name with different dimensions. You need to
write the variable to the file under a different variable name.

I also made some changes to your linint2_Wrap line:

begin
   in = addfile("prate_media2.nc", "w")
   data = in->prate
   new_lon = fspan(1.25, 358.75, 144)
   new_lat = fspan(88.75, -88.75, 72)
   new_data = linint2_Wrap(data&lon,data&lat,data,True,new_lon,new_lat,0)

   new_data!0 ="lat"
   new_data!1 = "lon"
   new_data&lat = new_lat
   new_data&lon = new_lon
   new_data&lat_at_units = "degrees_north" ; needed for plotting in NCL
   new_data&lon_at_units = "degrees_east"
   in->prate_regrid = new_data
end

Good luck,
Adam

Sara Reis wrote:
> Hi, NCL users!
>
> I need to convert a data type of an archive "prate_media2.nc" with not
> constant resolutions (specified below) for the resolution of 2,5x2,5
> (LATITUDE = 88.75, 86.25, 83.75,...,-83.75, -86.25, -88.75 and LONGITUDE
> = 1.25, 3.75, 6.25,...,353.75, 356.25, 358.75).
>
> At archive "prate_media2.nc", the lat has 94 values and the lon has 192
> values.
> I want the new lat with 72 values and the new lon with 144 values
> (because of the new resolution).
>
> My script...
>
> begin
>
> in = addfile("prate_media2.nc", "w")
> data = in->prate
>
> new_lon = fspan(1.25, 358.75, 144)
> new_lat = fspan(88.75, -88.75, 72)
>
> new_data =
> linint2_Wrap(data&lon,data&lat,inmatrix,True,new_data,new_data,0)
>
> new_data!0 ="lat"
> new_data!1 = "lon"
> new_data&lat = new_lat
> new_data&lon = new_lon
>
> in->prate = new_data ; THE PROBLEM IS HERE
>
> end
>
> THE LAT AND LON OF THE ARCHIVE ".nc":
> lat = 88.542, 86.6531, 84.7532, 82.8508, 80.9473, 79.0435, 77.1394,
> 75.2351,
> 73.3307, 71.4262, 69.5217, 67.6171, 65.7125, 63.8079, 61.9033, 59.9986,
> 58.0939, 56.1893, 54.2846, 52.3799, 50.4752, 48.5705, 46.6658, 44.7611,
> ...
> -67.6171, -69.5217, -71.4262, -73.3307, -75.2351, -77.1394, -79.0435,
> -80.9473, -82.8508, -84.7532, -86.6531, -88.542 ;
>
> lon = 0, 1.875, 3.75, 5.625, 7.5, 9.375, 11.25, 13.125, 15, 16.875, 18.75,
> 20.625, 22.5, 24.375, 26.25, 28.125, 30, 31.875, 33.75, 35.625, 37.5,
> 39.375, 41.25, 43.125, 45, 46.875, 48.75, 50.625, 52.5, 54.375, 56.25,
> ...
> 331.875, 333.75, 335.625, 337.5, 339.375, 341.25, 343.125, 345,
> 346.875,
> 348.75, 350.625, 352.5, 354.375, 356.25, 358.125 ;
>
> Can`t I do this?
>
> Thaks a lot...
>
> Sara
>
> __________________________________________________
> Fale com seus amigos de graça com o novo Yahoo! Messenger
> http://br.messenger.yahoo.com/
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
--------------------------------------------------------------
Adam Phillips			             asphilli_at_ucar.edu
National Center for Atmospheric Research   tel: (303) 497-1726
ESSL/CGD/CAS                               fax: (303) 497-1333
P.O. Box 3000				
Boulder, CO 80307-3000	  http://www.cgd.ucar.edu/cas/asphilli
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk