<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN">
<html><body>
<p>Hi,</p>
<p><span><span> </span></span><span><span>I'm</span> <span>plotting</span> <span>sea </span><span>surface temperature</span> <span>from</span> <span>an output of the</span> <span>ROMS</span><span>, but a</span> <span>plot</span> <span>region</span> <span>is completely</span> <span>white</span><span>.</span> <span>I checked the</span> <span>file and it</span> <span>is correct</span><span>, I can</span> <span>usually</span> <span>plot</span> <span>in other</span> <span>languages (</span><span>matlab,</span> <span>for example), the problem occurs after calling the script Roms_utils.ncl (https://www.ncl.ucar.edu/Applications/Scripts/ROMS_utils.ncl).</span><br /></span></p>
<p><span><span>I had to make 2 modifications to be able to use this script:<br /></span></span></p>
<p><span><span>1 - comment some lines: 5-133 - <span id="result_box"><span>I do not know</span> <span>why, but</span> <span>the script</span> <span>does not run</span> <span>if these</span> <span>lines are</span> <span>uncommented</span><span>.</span></span><br /></span></span></p>
<p><span><span>2 - add another parameter to call function roms_3d_interp: in my case the variable h (Final bathymetry at RHO-points) and the variable that will be interpolated (temp: time-averaged potential temperature) don't be in the same file. h is in </span></span>my <strong>grid file</strong> and temp is in the <strong>output model file</strong>. So, <span id="result_box" class="short_text"><span>I added</span> <span>a new parameter</span> <span>when I calling</span> <span>the function - the file_grid</span></span>:</p>
<p><span style="color: #808080;">137</span> <span style="color: #ff0000;">function</span> roms_3d_interp( file_handle, <strong>file_grid</strong>, var:string, \<br /><span style="color: #808080;">138</span> rec:integer, z:numeric )</p>
<p>and I changed the lines 166-176 and 208: <span id="result_box" class="short_text"><span>to call</span> <span>the file</span> <span>that contains the variable h</span></span></p>
<p><span class="short_text"><span><span style="color: #808080;">166 </span>if(typeof(file_grid).eq."file") then<br /><span style="color: #808080;">167 </span> ISFILE = True<br /><span style="color: #808080;">168 </span> nc_grid = file_grid<br /><span class="short_text"><span><span style="color: #808080;">169 </span></span></span> else if(typeof(file_grid).eq."list") then<br /><span class="short_text"><span><span style="color: #808080;">170</span></span></span> ISFILE = False<br /><span class="short_text"><span><span style="color: #808080;">171 </span></span></span> nc_grid = file_grid[0]<br /><span class="short_text"><span><span style="color: #808080;">172</span></span></span> else<br /><span class="short_text"><span><span style="color: #808080;">173 </span></span></span> print("roms_interp_3d: error: the first argument must be a file or a list of files opened with addfile or addfiles")<br /><span class="short_text"><span><span style="color: #808080;">174</span></span></span> return<br /><span class="short_text"><span><span style="color: #808080;">175</span></span></span> end if<br /><span class="short_text"><span><span style="color: #808080;">176</span></span></span> end if</span></span></p>
<p><span class="short_text"><span>...<br /></span></span></p>
<p><span style="color: #808080;">208</span> if(isfilevar(<strong>nc_grid</strong>,"h"))</p>
<p> </p>
<p>After this function a region with constant temperature is interpolated with missing values. The figures and the scripts are attached. The files: nested_avg_19811216_LEILANE.nc and mesh_donor_grd_004nested_</p>
<div>LEILANE.nc are in the FTP server. My NCL version is 6.3.0.</div>
<div> </div>
<div><span id="result_box" class="short_text"><span>I'm new in NCL and I</span><span>'ll appreciate</span> <span>any help</span></span>.</div>
<div> </div>
<div> </div>
<div>Best regards.</div>
<div> </div>
<div>Leilane</div>
<p> </p>
<p><span><span> </span></span></p>
<p><span><span> </span></span></p>
<div> </div>
</body></html>