<p><b>laura@ucar.edu</b> 2013-04-18 12:54:41 -0600 (Thu, 18 Apr 2013)</p><p>limits the Richardson number to -100 in unstable layers in the free atmosphere, as indicated in Hong et al. (2006). Adding that limit avoids MPAS to crash when using lower resolution meshes with longer time-steps.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_atmos_physics/physics_wrf/module_bl_ysu.F
===================================================================
--- branches/atmos_physics/src/core_atmos_physics/physics_wrf/module_bl_ysu.F        2013-04-18 18:20:43 UTC (rev 2771)
+++ branches/atmos_physics/src/core_atmos_physics/physics_wrf/module_bl_ysu.F        2013-04-18 18:54:41 UTC (rev 2772)
@@ -380,6 +380,9 @@
 !               ==&gt; consider thermal z0 when differs from mechanical z0
 !              a bug fix in wscale computation in stable bl, sukanta basu, jun 2012
 !               ==&gt; wscale becomes small with height, and less mixing in stable bl
+!               ==&gt; ri = max(ri,rimin). limits the richardson number to -100 in
+!               unstable layers, following Hong et al. 2006.
+!               Laura D. Fowler (2013-04-18).
 !
 !     references:
 !
@@ -1002,6 +1005,7 @@
          rlamdz = min(dza(i,k+1),rlamdz)
          rl2 = (zk*rlamdz/(rlamdz+zk))**2
          dk = rl2*sqrt(ss)
+         ri = max(ri,rimin)
          if(ri.lt.0.)then
 ! unstable regime
            sri = sqrt(-ri)

</font>
</pre>