Re: How to regrid CMIP5 ocean wmo data to regular 1X1 grid

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Sun Mar 02 2014 - 12:50:14 MST

Did you not look at NCL's "Application" page?

http://www.ncl.ucar.edu/Applications/ESMF.shtml
Look at say, Example 6. Likely you will have to change
the script to match your input data grid.
Other Examples 7, 24

Generate the weight file once and apply thereafter.

On 3/2/14, 8:27 AM, Yang Li wrote:
> Hi,all
>
> I have downloaded CMIP5 ocean wmo
> data (wmo_Omon_CMCC-CM_historical_r1i1p1_187001-187912.nc), but find that
> the horizontal grid is not regular grid.
>
> Therefore I think I need to regrid this data from these irregular grid
> points to the regular 1X1 grid.
>
> I am asking if NCL can do the interpolation? If so, Could you point out
> to me the key sentences used in NCL scripts?
>
> Thanks!
>
> ;-----------------------------------------------
>
> *ncdump -h wmo_Omon_CMCC-CM_historical_r1i1p1_187001-187912.nc*
>
> dimensions:
>
> time = UNLIMITED ; // (120 currently)
>
> lev = 31 ;
>
> j = 149 ;
>
> i = 182 ;
>
> bnds = 2 ;
>
> vertices = 4 ;
>
> variables:
>
> double time(time) ;
>
> time:bounds = "time_bnds" ;
>
> time:units = "days since 1870-01-01" ;
>
> time:calendar = "standard" ;
>
> time:axis = "T" ;
>
> time:long_name = "time" ;
>
> time:standard_name = "time" ;
>
> double time_bnds(time, bnds) ;
>
> double lev(lev) ;
>
> lev:bounds = "lev_bnds" ;
>
> lev:units = "m" ;
>
> lev:axis = "Z" ;
>
> lev:positive = "down" ;
>
> lev:long_name = "ocean depth coordinate" ;
>
> lev:standard_name = "depth" ;
>
> double lev_bnds(lev, bnds) ;
>
> int j(j) ;
>
> j:units = "1" ;
>
> j:long_name = "cell index along second dimension" ;
>
> int i(i) ;
>
> i:units = "1" ;
>
> i:long_name = "cell index along first dimension" ;
>
> float lat(j, i) ;
>
> lat:standard_name = "latitude" ;
>
> lat:long_name = "latitude coordinate" ;
>
> lat:units = "degrees_north" ;
>
> lat:bounds = "lat_vertices" ;
>
> float lon(j, i) ;
>
> lon:standard_name = "longitude" ;
>
> lon:long_name = "longitude coordinate" ;
>
> lon:units = "degrees_east" ;
>
> lon:bounds = "lon_vertices" ;
>
> float lat_vertices(j, i, vertices) ;
>
> lat_vertices:units = "degrees_north" ;
>
> float lon_vertices(j, i, vertices) ;
>
> lon_vertices:units = "degrees_east" ;
>
> float wmo(time, lev, j, i) ;
>
> wmo:standard_name = "upward_ocean_mass_transport" ;
>
> wmo:long_name = "Upward Ocean Mass Transport" ;
>
> wmo:comment = "differs from CMIP3, which only had upward velocity." ;
>
> wmo:units = "kg s-1" ;
>
> wmo:original_name = "wmo" ;
>
> wmo:cell_methods = "time: mean (interval: 1 month)" ;
>
> wmo:cell_measures = "area: areacello volume: volcello" ;
>
> wmo:missing_value = 1.e+20f ;
>
> wmo:_FillValue = 1.e+20f ;
>
>
>
> _______________________________________________
> 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 Sun Mar 2 12:50:20 2014

This archive was generated by hypermail 2.1.8 : Mon Mar 03 2014 - 14:26:18 MST