<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Dear Todd,<br>
<br>
My name is Breogan Gomez and I work in the same department as Eduardo
Penabad, who has started this discussion. We have tried your
suggestion, modifying the Registry and recompiling the code as wrfhelp
have suggested, and the model still crashes. This are the parameters we
have chosen for the adaptive time_step for our 36Km, 12Km, 4Km domain.
As you can see this are fairly conservative. <br>
<br>
<i>&nbsp;time_step&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 225,<br>
[...]<br>
&nbsp;use_adaptive_time_step&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = .true.,<br>
&nbsp;step_to_output_time&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = .true.,<br>
&nbsp;target_cfl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.1,<br>
&nbsp;max_step_increase_pct&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 5,51,51<br>
&nbsp;starting_time_step&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -1,<br>
&nbsp;max_time_step&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 475,<br>
&nbsp;min_time_step&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 180,</i><br>
<br>
<br>
As you have suggested we had a look to the last lines before crashing.
This are the values:<br>
<br>
<i>Timing for main: time 2008-06-24_06:23:45 on domain&nbsp;&nbsp; 1:&nbsp;&nbsp;&nbsp; 3.39700
elapsed seconds.<br>
Timing for main: time 2008-06-24_06:26:23 on domain&nbsp;&nbsp; 3:&nbsp;&nbsp;&nbsp; 0.23600
elapsed seconds.<br>
Timing for main: time 2008-06-24_06:26:23 on domain&nbsp;&nbsp; 2:&nbsp;&nbsp;&nbsp; 0.74010
elapsed seconds.<br>
Timing for main: time 2008-06-24_06:29:01 on domain&nbsp;&nbsp; 3:&nbsp;&nbsp;&nbsp; 0.29760
elapsed seconds.<br>
Timing for main: time 2008-06-24_06:29:01 on domain&nbsp;&nbsp; 2:&nbsp;&nbsp;&nbsp; 0.75270
elapsed seconds.<br>
Timing for main: time 2008-06-24_06:31:40 on domain&nbsp;&nbsp; 3:&nbsp;&nbsp;&nbsp; 0.27950
elapsed seconds.<br>
Timing for main: time 2008-06-24_06:31:40 on domain&nbsp;&nbsp; 2:&nbsp;&nbsp;&nbsp; 0.75810
elapsed seconds.<br>
Timing for main: time 2008-06-24_06:31:40 on domain&nbsp;&nbsp; 1:&nbsp;&nbsp;&nbsp; 2.72560
elapsed seconds.<br>
Timing for main: time 2008-06-24_06:34:18 on domain&nbsp;&nbsp; 3:&nbsp;&nbsp;&nbsp; 0.32230
elapsed seconds.<br>
Timing for main: time 2008-06-24_06:34:18 on domain&nbsp;&nbsp; 2:&nbsp;&nbsp;&nbsp; 0.80300
elapsed seconds.<br>
&nbsp; WOULD GO OFF TOP: KF_ETA_PARA I,J,DPTHMX,DPMIN&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 11&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
7 NaN</i><br>
<br>
<br>
As you can see model domain 3 (4Km) has reaches aproximated time step
of 150 seconds! It could be possible that model will be blowing because
the inner domains are reaching too high time-step values? Which yields
to another question.&nbsp; <br>
<br>
Are the time_steps from the inner domains limited by the
"max_time_step" and "min_time_step" parameters? If so, are this values
divided by the resolution ratio? If not, Could I modify the Regitry, in
the same way I did for the "max_step_increase_pct", to select different
"max_time_step" and "min_time_step" for each domain?<br>
<br>
Finally, do you have any other suggestion we could try?<br>
<br>
All the best<br>
<br>
<br>
Hutchinson, Todd escribi&oacute;:
<blockquote
 cite="mid:8536F69C1FCC294B859D07B179F069440E48F9F9@EXCHANGE.ad.wsicorp.com"
 type="cite">
  <meta http-equiv="Content-Type" content="text/html; ">
  <meta content="MSHTML 6.00.6000.16674" name="GENERATOR">
  <div dir="ltr" align="left"><span class="440004713-20062008"><font
 color="#0000ff" face="Arial" size="2">Eduardo,</font></span></div>
  <div dir="ltr" align="left"><span class="440004713-20062008"></span>&nbsp;</div>
  <div dir="ltr" align="left"><span class="440004713-20062008"><font
 color="#0000ff" face="Arial" size="2">Glad to see that you are working
with the adaptive time-step.&nbsp; As I'm sure you know, it is new in
WRF3.0.&nbsp; We developed and have been using the adaptive time-step in
operational forecasts&nbsp;with a modified version of WRF2.2 for quite some
time.&nbsp; Our operational simulations consist of an outer 12km nest and an
inner 4km nest.&nbsp;&nbsp; This configuration has been very stable for us.&nbsp;
However, we have only done limited testing with three levels of nesting.</font></span></div>
  <div dir="ltr" align="left"><span class="440004713-20062008"></span>&nbsp;</div>
  <div dir="ltr" align="left"><span class="440004713-20062008"><font
 color="#0000ff" face="Arial" size="2">I'm interested to see the length
of the time-steps just before the model crash.&nbsp; Could you post the last
hundred or so "Timing for main" lines from the rsl.out.0000 file?</font></span></div>
  <div dir="ltr" align="left"><span class="440004713-20062008"></span>&nbsp;</div>
  <div dir="ltr" align="left"><span class="440004713-20062008"><font
 color="#0000ff" face="Arial" size="2">Have you set any parameters in
the namelist file other than the use_adaptive_time_step setting?</font></span></div>
  <div dir="ltr" align="left"><span class="440004713-20062008"></span>&nbsp;</div>
  <div dir="ltr" align="left"><span class="440004713-20062008"><font
 color="#0000ff" face="Arial" size="2">I'd suggest that you add the
following to your namelist file:</font></span></div>
  <div dir="ltr" align="left"><span class="440004713-20062008"><font
 color="#0000ff" face="Arial" size="2">max_step_increase_pct = 5, 51,
51&nbsp;&nbsp; (the default is 5,5,5)</font></span></div>
  <div dir="ltr" align="left"><span class="440004713-20062008"></span>&nbsp;</div>
  <div dir="ltr" align="left"><span class="440004713-20062008"><font
 color="#0000ff" face="Arial" size="2">This will limit the time-step of
the outer nest from changing by more than 5%, but allow the time-step
on the inner two nests to change by up to 51%.&nbsp; Since the inner nest(s)
time steps must be evenly divisible into the outer nest time-step, with
the default of 5,5,5, the time step on the inner nests will
not&nbsp;increase unless the time-step on the outer nest changes.&nbsp; </font></span></div>
  <div dir="ltr" align="left"><span class="440004713-20062008"></span>&nbsp;</div>
  <div dir="ltr" align="left"><span class="440004713-20062008"><font
 color="#0000ff" face="Arial" size="2">Just FYI:</font></span></div>
  <div dir="ltr" align="left"><span class="440004713-20062008"><font
 color="#0000ff" face="Arial" size="2">I believe the 'WOULD GO OFF TOP"
lines that you see are a symptom of the model being unstable.&nbsp; I've
seen these types of errors in several situations:</font></span></div>
  <div dir="ltr" align="left"><span class="440004713-20062008"><font
 color="#0000ff" face="Arial" size="2">1) a couple of minutes into a
simulation -- often this is caused by bad input data</font></span></div>
  <div dir="ltr" align="left"><span class="440004713-20062008"><font
 color="#0000ff" face="Arial" size="2">2) later on during a simulation
-- this is most often caused by a time-step problem, or bad boundary
data.</font></span></div>
  <br>
  <div class="OutlookMessageHeader" dir="ltr" align="left" lang="en-us">
  <hr tabindex="-1"><font face="Tahoma" size="2"><b>From:</b>
<a class="moz-txt-link-abbreviated" href="mailto:wrf-users-bounces@ucar.edu">wrf-users-bounces@ucar.edu</a> [<a class="moz-txt-link-freetext" href="mailto:wrf-users-bounces@ucar.edu">mailto:wrf-users-bounces@ucar.edu</a>] <b>On
Behalf Of </b>Eduardo Penabad Ramos<br>
  <b>Sent:</b> Tuesday, June 17, 2008 11:37 AM<br>
  <b>To:</b> <a class="moz-txt-link-abbreviated" href="mailto:wrf-users@ucar.edu">wrf-users@ucar.edu</a><br>
  <b>Cc:</b> 'wrfhelp'<br>
  <b>Subject:</b> [Wrf-users] Tests with adaptive time step in WRF-ARW
(v3)<br>
  </font><br>
  </div>
  <div><span class="050511815-17062008"><font face="Arial" size="2">Hello!</font></span></div>
  <div><span class="050511815-17062008"></span>&nbsp;</div>
  <div><span class="050511815-17062008"><font face="Arial" size="2">we
are running a setup of 3 nested grids with WRF-ARW (v3) over Galicia
(NW Iberian Peninsula), with dx=36,12 and 4 km.</font></span></div>
  <div><span class="050511815-17062008"></span>&nbsp;</div>
  <div><span class="050511815-17062008"><font face="Arial" size="2">We
are testing the adaptive time step with this grid configuration, and
with the following set of physics:</font></span></div>
  <div><span class="050511815-17062008"></span>&nbsp;</div>
  <div><span class="050511815-17062008"><font face="Arial" size="2">Radiation:
Longwave RRTM, Shortwave MM5-Dudhia</font></span></div>
  <div><span class="050511815-17062008"><font face="Arial" size="2">PBL:
YSU (with sfclay: Monin-Obukhov from MM5 MRF)</font></span></div>
  <div><span class="050511815-17062008"><font face="Arial" size="2">Surface:
5 layer MM5 LSM</font></span></div>
  <div><span class="050511815-17062008"><font face="Arial" size="2">Cumulus:
Kain-Fritsch</font></span></div>
  <div><span class="050511815-17062008"><font face="Arial" size="2">Microphysics:
WSM6</font></span></div>
  <div><span class="050511815-17062008"></span>&nbsp;</div>
  <div><span class="050511815-17062008"><font face="Arial" size="2">With
this setup we have experienced some problems related with Kain-Fritsch
(we guess). After about 40-50 hours of simulation, error messages
appear in the WRF output textfiles:</font></span></div>
  <div><span class="050511815-17062008">
  <pre><i>WOULD GO OFF TOP: KF_ETA_PARA I,J,DPTHMX,DPMIN</i></pre>
  <pre><font face="Arial" size="2">where NaNs appeared for DPTHMAX</font></pre>
  <pre><font face="Arial" size="2">We have tested with Betts-Miller-Janjic cumulus parameterization, and although no error messages appeared, wrfout files also have some NaNs inside.</font></pre>
  <pre>&nbsp;</pre>
  <pre><font face="Arial" size="2">Do you have any suggestions? Thanks in advance!</font></pre>
  <pre>&nbsp;</pre>
  <pre><font face="Arial" size="2">Eduardo Penabad
MeteoGalicia</font></pre>
  <pre>&nbsp;</pre>
  <pre><em></em>&nbsp;</pre>
  </span></div>
  <div><span class="050511815-17062008"></span>&nbsp;</div>
  <pre wrap="">
<hr size="4" width="90%">
_______________________________________________
Wrf-users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Wrf-users@ucar.edu">Wrf-users@ucar.edu</a>
<a class="moz-txt-link-freetext" href="http://mailman.ucar.edu/mailman/listinfo/wrf-users">http://mailman.ucar.edu/mailman/listinfo/wrf-users</a>
  </pre>
</blockquote>
<br>
<div class="moz-signature">-- <br>
<meta http-equiv="Content-Type" content="text/html; ">
<meta name="ProgId" content="Word.Document">
<meta name="Generator" content="Microsoft Word 11">
<meta name="Originator" content="Microsoft Word 11">
<link rel="File-List" href="BreoganGomez_archivos/filelist.xml">
<title>Num&eacute;rico Firma</title>
<!--[if gte mso 9]><xml>
 <o:DocumentProperties>
  <o:Author>Eduardo Penabad</o:Author>
  <o:LastAuthor>breogan.gomez</o:LastAuthor>
  <o:Revision>5</o:Revision>
  <o:TotalTime>7</o:TotalTime>
  <o:Created>2006-12-20T08:35:00Z</o:Created>
  <o:LastSaved>2008-05-15T08:10:00Z</o:LastSaved>
  <o:Pages>1</o:Pages>
  <o:Words>65</o:Words>
  <o:Characters>359</o:Characters>
  <o:Lines>2</o:Lines>
  <o:Paragraphs>1</o:Paragraphs>
  <o:CharactersWithSpaces>423</o:CharactersWithSpaces>
  <o:Version>11.5606</o:Version>
 </o:DocumentProperties>
</xml><![endif]--><!--[if gte mso 9]><xml>
 <w:WordDocument>
  <w:Zoom>BestFit</w:Zoom>
  <w:SpellingState>Clean</w:SpellingState>
  <w:GrammarState>Clean</w:GrammarState>
  <w:HyphenationZone>21</w:HyphenationZone>
  <w:ValidateAgainstSchemas/>
  <w:SaveIfXMLInvalid>false</w:SaveIfXMLInvalid>
  <w:IgnoreMixedContent>false</w:IgnoreMixedContent>
  <w:AlwaysShowPlaceholderText>false</w:AlwaysShowPlaceholderText>
  <w:BrowserLevel>MicrosoftInterne
tExplorer4</w:BrowserLevel>
 </w:WordDocument>
</xml><![endif]--><!--[if gte mso 9]><xml>
 <w:LatentStyles DefLockedState="false" LatentStyleCount="156">
 </w:LatentStyles>
</xml><![endif]-->
<style>
<!--
 /* Font Definitions */
 @font-face
        {font-family:Verdana;
        panose-1:2 11 6 4 3 5 4 4 2 4;
        mso-font-charset:0;
        mso-generic-font-family:swiss;
        mso-font-pitch:variable;
        mso-font-signature:536871559 0 0 0 415 0;}
 /* Style Definitions */
 p.MsoNormal, li.MsoNormal, div.MsoNormal
        {mso-style-parent:"";
        margin:0cm;
        margin-bottom:.0001pt;
        mso-pagination:widow-orphan;
        font-size:12.0pt;
        font-family:"Times New Roman";
        mso-fareast-font-family:"Times New Roman";}
a:link, span.MsoHyperlink
        {color:blue;
        text-decoration:underline;
        text-underline:single;}
a:visited, span.MsoHyperlinkFollowed
        {color:purple;
        text-decoration:underline;
        text-underline:single;}
span.SpellE
        {mso-style-name:"";
        mso-spl-e:yes;}
span.GramE
        {mso-style-name:"";
        mso-gram-e:yes;}
@page Section1
        {size:595.3pt 841.9pt;
        margin:70.85pt 3.0cm 70.85pt 3.0cm;
        mso-header-margin:35.4pt;
        mso-footer-margin:35.4pt;
        mso-paper-source:0;}
div.Section1
        {page:Section1;}
-->
</style><!--[if gte mso 10]>
<style>
 /* Style Definitions */
 table.MsoNormalTable
        {mso-style-name:"Tabla normal";
        mso-tstyle-rowband-size:0;
        mso-tstyle-colband-size:0;
        mso-style-noshow:yes;
        mso-style-parent:"";
        mso-padding-alt:0cm 5.4pt 0cm 5.4pt;
        mso-para-margin:0cm;
        mso-para-margin-bottom:.0001pt;
        mso-pagination:widow-orphan;
        font-size:10.0pt;
        font-family:"Times New Roman";
        mso-ansi-language:#0400;
        mso-fareast-language:#0400;
        mso-bidi-language:#0400;}
</style>
<![endif]-->
<div class="Section1">
<p class="MsoNormal"><b><span style="font-family: Verdana;">Breog&aacute;n
G&oacute;mez Hombre</span></b></p>
<p class="MsoNormal">&nbsp;</p>
<p class="MsoNormal"><strong><span
 style="font-family: Verdana; color: blue;">Meteo</span></strong><strong><span
 style="font-family: Verdana; color: rgb(153, 153, 153);">Galicia &#8211; </span></strong><b
 style=""><span style="font-size: 10pt; font-family: Verdana;">Dpto.
Investigaci&oacute;n e <span class="SpellE">Predici&oacute;n</span> Num&eacute;rica</span></b><span
 style="font-size: 8pt;"><o:p></o:p></span></p>
<p class="MsoNormal"><span class="SpellE"><b style=""><i style=""><span
 style="font-size: 8pt; font-family: Verdana;">Conseller&iacute;a</span></i></b></span><b
 style=""><i style=""><span
 style="font-size: 8pt; font-family: Verdana;"> de Medio Ambiente e <span
 class="SpellE">Desenvolvemento</span> <span class="SpellE">Sostible</span><o:p></o:p></span></i></b></p>
<p class="MsoNormal"><span class="SpellE"><b style=""><i style=""><span
 style="font-size: 8pt; font-family: Verdana;">Xunta</span></i></b></span><b
 style=""><i style=""><span
 style="font-size: 8pt; font-family: Verdana;"> de Galicia<o:p></o:p></span></i></b></p>
<p class="MsoNormal"><b style=""><span
 style="font-size: 8pt; font-family: Verdana;"><o:p>&nbsp;</o:p></span></b></p>
<p class="MsoNormal"><span style="font-size: 8pt; font-family: Verdana;"
 lang="PT-BR">&Aacute;rea Central, Local 31-<span class="GramE">C</span><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size: 8pt; font-family: Verdana;"
 lang="PT-BR">Pol&iacute;gono de <span class="SpellE">Fonti&ntilde;as</span>. <span
 class="GramE">15703 Santiago</span> de Compostela<o:p></o:p></span></p>
<p class="MsoNormal"><span class="SpellE"><span
 style="font-size: 8pt; font-family: Verdana;" lang="PT-BR">Tfno</span></span><span
 style="font-size: 8pt; font-family: Verdana;" lang="PT-BR">:
(+34) 981 95 74 62<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size: 8pt; font-family: Verdana;"
 lang="PT-BR">FAX: (+34) 981 95 74 66<o:p></o:p></span></p>
<p class="MsoNormal"><b style=""><span
 style="font-size: 8pt; font-family: Verdana;" lang="PT-BR"><a
 href="http://www.meteogalicia.es/">www.meteogalicia.es</a> &#8211; <a
 href="mailto:numerico@meteogalicia.es">breogangomez@meteogalicia.es</a>
<o:p></o:p></span></b></p>
<p class="MsoNormal"><b style=""><span
 style="font-size: 8pt; font-family: Verdana;" lang="PT-BR"><o:p>&nbsp;</o:p></span></b></p>
</div>
</div>
</body>
</html>