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

[Sri nandini]

Sloved.

Thank you. That worked perfectly.

Sri

On 05.11.21 23:13, Dennis Shea wrote:
> 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
>     <http://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
>
-- 
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20211108/77da9554/attachment.html>