[ncl-talk] converting generic POP grid to regular LatLon grid

[Dennis Shea]

The pop2lat_1.ncl i script approach was an older method to regrid.
I would suggest ESMF regridding:
*http://www.ncl.ucar.edu/Applications/ESMF.shtml*
<http://www.ncl.ucar.edu/Applications/ESMF.shtml>
Example 24

On Fri, Nov 5, 2021 at 7:33 AM Sri nandini via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hello dear ncl-users,
>
> Can someone advice me whether the following method is correct?
>
> I have an initial CESM POP file with lat and lon in generic grid which
> looks like the below:
>
> netcdf zos_Omon_CESM1-CAM5_historical_rcp85_r38i1p1_192001-210012 {
> dimensions:
>     time = UNLIMITED ; // (2172 currently)
>     i = 320 ;
>     j = 384 ;
> variables:
>     double time(time) ;
>         time:standard_name = "time" ;
>         time:long_name = "time" ;
>         time:units = "days since 0000-01-01 00:00:00" ;
>         time:calendar = "365_day" ;
>         time:axis = "T" ;
>     int i(i) ;
>         i:long_name = "cell index along first dimension" ;
>         i:units = "1" ;
>         i:axis = "X" ;
>     int j(j) ;
>         j:long_name = "cell index along second dimension" ;
>         j:units = "1" ;
>         j:axis = "Y" ;
>     float lat(j, i) ;
>         lat:standard_name = "latitude" ;
>         lat:long_name = "latitude coordinate" ;
>         lat:units = "degrees_north" ;
>     float lon(j, i) ;
>         lon:standard_name = "longitude" ;
>         lon:long_name = "longitude coordinate" ;
>         lon:units = "degrees_east" ;
>     float zos(time, j, i) ;
>         zos:standard_name = "sea_surface_height_above_geoid" ;
>         zos:long_name = "sea_surface_height_above_geoid" ;
>         zos:units = "m" ;
>         zos:_FillValue = 1.e+20f ;
>         zos:missing_value = 1.e+20f ;
>         zos:comment = "Created using NCL code" ;
>         zos:cell_methods = "time: mean (interval: 1 month)" ;
>         zos:history = "SSH/100" ;
>         zos:original_units = "cm" ;
>         zos:original_name = "SSH" ;
>
> I tried several methods to convert this into a regular LatLon grid,
> finally using the method from here:
> https://www.ncl.ucar.edu/Applications/pop2lat.shtml
>
> I used the pop2lat_1.ncl
> <https://www.ncl.ucar.edu/Applications/Scripts/pop2lat_1.ncl>: to
> interpolate with this POP weights file:
> map_gx1v6_to_1x1d_bilin_da_100716.nc.
>
> Attached here is the script i modified as well as the plot. I do not get
> any error, but i wish to know if this was the correct method? There are
> many POP weights available and i choose this one.
>
> Best
>
> Sri
>
> --
> Dr. Sri Nandini-Weiss
>
> Center for Earth System Research and Sustainability (CEN)
> Cluster of Excellence 'Climate, Climatic Change, and Society' (CLICCS)
>
> Universität Hamburg
> Institute of Oceanography
> Bundesstraße 53
> 20146 Hamburg
> Germany
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20211105/0702ce1f/attachment.html>