<div dir="ltr"><div><div><div><div><div><div><div>Thank you.<br><br></div>A JIRA ticket has been created: Issue <a class="" href="https://vets.development.ucar.edu/jira/browse/NCL-2077">NCL-2077 - omega_ccm_driver does not work with 3D array</a><br><br></div>This will be fixed in the 6.3.0 release<br><br>--<br></div>As you know, you can use an editor to extract the 'omega_ccm_driver' code from contributed.ncl.<br></div>Then modify it as you suggest (my_omega_ccm_driver.ncl) and load the new version, preceding it with <br>This will 'undefine' the original version and use the new version<br><br>So, my_omega_ccm_driver.ncl has<br><br>undef("omega_ccm_driver")<br>function omega_ccm_driver(p0,psfc,u,v,hyam[*],hybm[*],hyai[*],hybi[*])<br>;<br><br></div><div>Then <br><br></div><div>load "./my_omega_ccm_driver.ncl"<br><br><br></div><div><br></div><div>load<br></div><div><br>===<br></div>I will change the driver script but we must check the code to make sure it works.<br></div>I think Colin has done so.<br><br></div>D<br><div><div><div><br><br><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Sep 23, 2014 at 10:49 AM, Colin Zarzycki <span dir="ltr"><<a href="mailto:zarzycki@umich.edu" target="_blank">zarzycki@umich.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Omega_ccm_driver should accept rank 3 or 4 values of U and V (and 2 or 3 of PS). (with or without a time dimension)<div><br></div><div>However, when passing 3D wind (lev,lat,lon) instead of 4D wind into the function, it crashes because it searches for dimu(3) (4th index of u dimensions) which doesn't exist.</div><div><br></div><div>The simple fix is to replace:</div><div><br></div><div><div> ntim = dimu(0)</div><div> klev = dimu(1)</div><div> nlat = dimu(2)</div><div> mlon = dimu(3)</div></div><div><br></div><div>with...</div><div><br></div><div><div> if (ranku.eq.3) then</div><div> klev = dimu(0)</div><div> nlat = dimu(1)</div><div> mlon = dimu(2)</div><div> else</div><div> ntim = dimu(0)</div><div> klev = dimu(1)</div><div> nlat = dimu(2)</div><div> mlon = dimu(3)</div><div> end if</div></div><div><br></div><div>Cheers,</div><div><br></div><div>-Colin</div><span class="HOEnZb"><font color="#888888"><div><br>-- <br><div dir="ltr"><div dir="ltr" style="font-family:arial;font-size:small">Colin M. Zarzycki, ASP postdoctoral fellow</div><div dir="ltr" style="font-family:arial;font-size:small">Atmospheric Modeling and Predictability<br>National Center for Atmospheric Research<br><a href="http://www.colinzarzycki.com/" style="color:rgb(17,85,204)" target="_blank">http://www.colinzarzycki.com</a></div></div>
</div></font></span></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>