<div dir="ltr">HI Barry,<div><br></div><div>I'm not entirely sure, but I'd guess that ice mixing ratio is needed to detect upper level cirrus, but anywhere that snow mixing ratio is being produced, the water vapor mixing ratio already has an optical depth of 1.0 above where that's occurring. The author probably just kept the algorithm simple. Keep in mind that this is just a diagnostic, which is simplified with things like "optical depth = 1.0". For some cloud types, an optical depth of 1.0 might be too thick (e.g. cirrus). A satellite observation is going to be more accurate. </div><div><br></div><div>Bill</div></div><br><div class="gmail_quote"><div dir="ltr">On Sat, Nov 3, 2018 at 9:46 AM Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com">barry.h.lynn@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Bill:<div><br></div><div>Thanks. Much appreciated.</div><div><br></div><div>Is there any particular reason why the scheme does not include snow in its calculations? Some schemes produce more snow than cloud ice....</div><div><br></div><div>Barry</div></div><br><div class="gmail_quote"><div dir="ltr">On Fri, Nov 2, 2018 at 7:31 PM Bill Ladwig <<a href="mailto:ladwig@ucar.edu" target="_blank">ladwig@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="ltr">Hi Barry,<div><br></div><div>In an upcoming release (or you can use wrf-python now), you can adjust the optical thickness to trigger when a cloud is detected, so that might help some. By default, that value is 1.0. The next version will also allow you to fill the no-cloud regions with fill values instead of the surface temperature, which is how it works now. Here is the source code for the algorithm:</div><div><br></div><div><a href="https://github.com/NCAR/wrf-python/blob/develop/fortran/wrf_fctt.f90" target="_blank">https://github.com/NCAR/wrf-python/blob/develop/fortran/wrf_fctt.f90</a><br></div><div><br></div><div>The algorithm itself is pretty basic. It looks from the top down, and integrates optical thickness until the threshold is met (1.0), then the cloud top temperature is set to the temperature of that level. Optical thickness at each level (if cloud ice mixing ratio is present in your files [note: converted to g/kg]) is:</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<p class="m_-8216277528516779165m_5278753536042180537gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Monaco;color:rgb(3,38,204)">opdepthd<span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(0,0,0)"> = </span>opdepthd<span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(0,0,0)"> + (</span>ABSCOEF<span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(0,0,0)">*</span>qcw<span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(0,0,0)">(</span>i<span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(0,0,0)">,</span>j<span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(0,0,0)">,</span>ripk<span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(0,0,0)">) + </span>ABSCOEFI<span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(0,0,0)">*</span>qci<span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(0,0,0)">(</span>i<span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(0,0,0)">,</span>j<span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(0,0,0)">,</span>ripk<span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(0,0,0)">))*</span>dp<span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(0,0,0)">/</span>G</p></blockquote><div><br></div><div>with the absorption coefficients:</div><div><br></div><div>
<blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote"><span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(147,26,104)">REAL</span><span class="m_-8216277528516779165m_5278753536042180537gmail-s2" style="color:rgb(0,0,0)">(</span><span class="m_-8216277528516779165m_5278753536042180537gmail-s3" style="color:rgb(117,42,204)">KIND</span><span class="m_-8216277528516779165m_5278753536042180537gmail-s2" style="color:rgb(0,0,0)">=8), </span><span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(147,26,104)">PARAMETER</span><span class="m_-8216277528516779165m_5278753536042180537gmail-s2" style="color:rgb(0,0,0)"> :: </span><span class="m_-8216277528516779165m_5278753536042180537gmail-s4" style="color:rgb(3,38,204)">ABSCOEFI</span><span class="m_-8216277528516779165m_5278753536042180537gmail-s2" style="color:rgb(0,0,0)"> = .272</span><span class="m_-8216277528516779165m_5278753536042180537gmail-s4" style="color:rgb(3,38,204)">D0</span><span class="m_-8216277528516779165m_5278753536042180537gmail-s2" style="color:rgb(0,0,0)"><span class="m_-8216277528516779165m_5278753536042180537gmail-Apple-converted-space"> </span></span>! cloud ice absorption coefficient in m^2/g<br><span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(147,26,104)">REAL</span><span class="m_-8216277528516779165m_5278753536042180537gmail-s2" style="color:rgb(0,0,0)">(</span><span class="m_-8216277528516779165m_5278753536042180537gmail-s3" style="color:rgb(117,42,204)">KIND</span><span class="m_-8216277528516779165m_5278753536042180537gmail-s2" style="color:rgb(0,0,0)">=8), </span><span class="m_-8216277528516779165m_5278753536042180537gmail-s1" style="color:rgb(147,26,104)">PARAMETER</span><span class="m_-8216277528516779165m_5278753536042180537gmail-s2" style="color:rgb(0,0,0)"> :: </span><span class="m_-8216277528516779165m_5278753536042180537gmail-s4" style="color:rgb(3,38,204)">ABSCOEF</span><span class="m_-8216277528516779165m_5278753536042180537gmail-s2" style="color:rgb(0,0,0)"> = .145</span><span class="m_-8216277528516779165m_5278753536042180537gmail-s4" style="color:rgb(3,38,204)">D0</span><span class="m_-8216277528516779165m_5278753536042180537gmail-s2" style="color:rgb(0,0,0)"> <span class="m_-8216277528516779165m_5278753536042180537gmail-Apple-converted-space"> </span></span>! cloud water absorption coefficient in m^2/g</blockquote>
</div><div><br></div><div>I'm not sure what the source is for those absorption constant values, but I imagine the author chose something safe, since it seems like something that could depend on cloud type, number concentration, etc. So, aside from changing the optical thickness threshold to trigger a cloud detection, the absorption coefficients could be changed from constants to variables, but that's really the only way to modify this algorithm in to something that uses more ice.</div><div><br></div><div>Hopefully this new behavior will be out in NCL in the next few months.</div><div><br></div><div>Hope this helps,</div><div><br></div><div>Bill</div><div><br></div><div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Fri, Nov 2, 2018 at 8:18 AM Barry Lynn <<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi:<div><br></div><div>I have used the "ctt" wrf get user calculation on two different case studies with WRF with 4 km grid spacing). In both, the cloud top height is much lower than the observations, much lower than the GEFS (and much lower than with WRF with a CU parameterization on a 12 km grid). </div><div><br></div><div>Hence, I think that there might be an error when calculating cloud top temperatures with (4 km) explicit microphysics. </div><div><br></div><div>Is there some easy way to modify the source code to be more sensitive to ice concentration? </div><div><br></div><div>Thanks,</div><div><br></div><div>Barry<br><br><div class="gmail_quote"><div dir="ltr">On Tue, Jul 24, 2018 at 11:02 PM Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div style="font-size:small">Barry,</div><div style="font-size:small"><br></div><div style="font-size:small">You can always get the latest version of the NCL source code from our github repo:</div><div style="font-size:small"><br></div><div><a href="https://github.com/ncar/ncl" target="_blank">https://github.com/ncar/ncl</a><br></div><div><br></div><div>If you click on the "Find File" button, then you can search for WRFUserARW.ncl individually.</div><div><br></div><div>From there you should be able to right click on the file name to "download linked file", or you can click on it to view the contents.</div><div><br></div><div>If you need the file from a particular version of NCL, like NCL V6.4.0, then from the <a href="http://github.com/ncar/ncl" target="_blank">github.com/ncar/ncl</a> page click on the downward triangle in the "Branch:master" button and select "Tags". You will then see different versions of NCL and can search the file in any one of those.</div><div><br></div><div>--Mary</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Jul 22, 2018 at 5:40 AM, Barry Lynn <span dir="ltr"><<a href="mailto:barry.h.lynn@gmail.com" target="_blank">barry.h.lynn@gmail.com</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">Hi:<div><br></div><div>I have been searching the web in vain to find the latest version of this file.</div><div><br></div><div>I must have an old version because it doesn't contain "ctt."</div><div><br></div><div>Could someone please point me in the direction of the download for this file.</div><div><br></div><div>(I am using version 6.4, but I must still be using an old file as it is not present.)</div><div><br></div><div>Thanks</div><div><div><br></div>-- <br><div dir="ltr" class="m_-8216277528516779165m_5278753536042180537m_-4936604952754196484m_-699335057869194960m_-1457772064199890876gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior Associate Scientist, Lecturer,</div><div><div><span style="color:rgb(136,136,136)">The Institute of the Earth Science, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">The Hebrew University of Jerusalem, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Givat Ram, Jerusalem 91904, Israel </span><br style="color:rgb(136,136,136)"></div><span style="color:rgb(136,136,136)">Tel: 972 547 231 170</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Fax: (972)-25662581</span></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>C.E.O, Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">http://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br><br></div></div></div></div></div></div></div>
<br>_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="m_-8216277528516779165m_5278753536042180537m_-4936604952754196484gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior Associate Scientist, Lecturer,</div><div><div><span style="color:rgb(136,136,136)">The Institute of the Earth Science, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">The Hebrew University of Jerusalem, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Givat Ram, Jerusalem 91904, Israel </span><br style="color:rgb(136,136,136)"></div><span style="color:rgb(136,136,136)">Tel: 972 547 231 170</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Fax: (972)-25662581</span></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>C.E.O, Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">http://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br><br></div></div></div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="m_-8216277528516779165gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr">Barry H. Lynn, Ph.D<div><div>Senior Associate Scientist, Lecturer,</div><div><div><span style="color:rgb(136,136,136)">The Institute of the Earth Science, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">The Hebrew University of Jerusalem, </span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Givat Ram, Jerusalem 91904, Israel </span><br style="color:rgb(136,136,136)"></div><span style="color:rgb(136,136,136)">Tel: 972 547 231 170</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Fax: (972)-25662581</span></div></div><div><span style="color:rgb(136,136,136)"><br></span></div><div>C.E.O, Weather It Is, LTD<br>Weather and Climate Focus<br><a href="http://weather-it-is.com" target="_blank">http://weather-it-is.com</a><br>Jerusalem, Israel<br>Local: 02 930 9525<br>Cell: 054 7 231 170<br>Int-IS: x972 2 930 9525<br><br></div></div></div></div></div>
</blockquote></div>