<div dir="ltr"><font face="arial, sans-serif">Hi guys,</font><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">can you please clarify for me: </font></div><div><font face="arial, sans-serif">1) if I want to calculate the potential density relative to the surface using ECCO output (THETA and SALT fields....) and &quot;density.m&quot; function from gcmfaces, where &quot;<b>p</b>&quot; is the depth of tracer points, do I specify <b>pref=0, </b></font><font face="arial, sans-serif">and the potential density is given by <b>rhop, </b>as<b>:</b></font><span style="font-family:arial,sans-serif"> </span></div><div></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">[rhop,rhpis,rhor] = density(t(:),s(:),p(:),pref(:));<br></font></div><div><font face="arial, sans-serif"><br></font></div><div><br></div><div><font face="arial, sans-serif">2) If following Iudicone et al. (2008), I want to estimate <i>the equation of evolution of the locally referenced
potential density defined at a reference pressure chosen as the local pressure pr</i>, I get that by multiplying the temperature tendency equation by alpha, salt</font> tendency equation by beta, and combine them. In doing so, is this a correct way to obtain alpha?</div><div><font face="arial, sans-serif"><br></font></div><div> alpha=sw_alpha(s(:),t(:),p(:),&#39;ptmp&#39;),<font face="arial, sans-serif"><br></font></div><div>where again <b>s</b> is obtained from the SALT fields, <b>t</b> from THETA and <b>p</b> is the depth of the tracer points.</div><div><br></div><div><font face="arial, sans-serif">I would very much appreciate your help,</font></div><div><font face="arial, sans-serif">thanks,</font></div><div><font face="arial, sans-serif">Ivana </font></div><div><span style="font-size:9pt;font-family:TimesTen"><br></span></div>
        
        
                </div>