<div dir="ltr"><div>The pop2lat_1.ncl i script approach was an older method to regrid.</div><div>I would suggest ESMF regridding:<br></div><div><a href="http://www.ncl.ucar.edu/Applications/ESMF.shtml"><b>http://www.ncl.ucar.edu/Applications/ESMF.shtml</b></a></div><div>Example 24<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Nov 5, 2021 at 7:33 AM Sri nandini via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<p>Hello dear ncl-users,</p>
<p>Can someone advice me whether the following method is correct?</p>
<p>I have an initial CESM POP file with lat and lon in generic grid
which looks like the below:</p>
<p>netcdf zos_Omon_CESM1-CAM5_historical_rcp85_r38i1p1_192001-210012
{<br>
dimensions:<br>
time = UNLIMITED ; // (2172 currently)<br>
i = 320 ;<br>
j = 384 ;<br>
variables:<br>
double time(time) ;<br>
time:standard_name = "time" ;<br>
time:long_name = "time" ;<br>
time:units = "days since 0000-01-01 00:00:00" ;<br>
time:calendar = "365_day" ;<br>
time:axis = "T" ;<br>
int i(i) ;<br>
i:long_name = "cell index along first dimension" ;<br>
i:units = "1" ;<br>
i:axis = "X" ;<br>
int j(j) ;<br>
j:long_name = "cell index along second dimension" ;<br>
j:units = "1" ;<br>
j:axis = "Y" ;<br>
float lat(j, i) ;<br>
lat:standard_name = "latitude" ;<br>
lat:long_name = "latitude coordinate" ;<br>
lat:units = "degrees_north" ;<br>
float lon(j, i) ;<br>
lon:standard_name = "longitude" ;<br>
lon:long_name = "longitude coordinate" ;<br>
lon:units = "degrees_east" ;<br>
float zos(time, j, i) ;<br>
zos:standard_name = "sea_surface_height_above_geoid" ;<br>
zos:long_name = "sea_surface_height_above_geoid" ;<br>
zos:units = "m" ;<br>
zos:_FillValue = 1.e+20f ;<br>
zos:missing_value = 1.e+20f ;<br>
zos:comment = "Created using NCL code" ;<br>
zos:cell_methods = "time: mean (interval: 1 month)" ;<br>
zos:history = "SSH/100" ;<br>
zos:original_units = "cm" ;<br>
zos:original_name = "SSH" ;</p>
<p>I tried several methods to convert this into a regular LatLon
grid, finally using the method from here:
<a href="https://www.ncl.ucar.edu/Applications/pop2lat.shtml" target="_blank">https://www.ncl.ucar.edu/Applications/pop2lat.shtml</a></p>
<p>I used the <a href="https://www.ncl.ucar.edu/Applications/Scripts/pop2lat_1.ncl" target="_blank">pop2lat_1.ncl</a>:
to interpolate with this POP weights file:
<a href="http://map_gx1v6_to_1x1d_bilin_da_100716.nc" target="_blank">map_gx1v6_to_1x1d_bilin_da_100716.nc</a>.</p>
<p>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. <br>
</p>
<p>Best</p>
<p>Sri<br>
</p>
<pre cols="72">--
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</pre>
</div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>