<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Wingdings;
panose-1:5 0 0 0 0 0 0 0 0 0;}
@font-face
{font-family:Wingdings;
panose-1:5 0 0 0 0 0 0 0 0 0;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
h3
{mso-style-priority:9;
mso-style-link:"Heading 3 Char";
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:13.5pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
code
{mso-style-priority:99;
font-family:"Courier New";}
p.MsoAcetate, li.MsoAcetate, div.MsoAcetate
{mso-style-priority:99;
mso-style-link:"Balloon Text Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:8.0pt;
font-family:"Tahoma","sans-serif";
mso-fareast-language:EN-US;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
span.Heading3Char
{mso-style-name:"Heading 3 Char";
mso-style-priority:9;
mso-style-link:"Heading 3";
font-family:"Times New Roman","serif";
mso-fareast-language:EN-GB;
font-weight:bold;}
p.author, li.author, div.author
{mso-style-name:author;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
span.BalloonTextChar
{mso-style-name:"Balloon Text Char";
mso-style-priority:99;
mso-style-link:"Balloon Text";
font-family:"Tahoma","sans-serif";}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
/* List Definitions */
@list l0
{mso-list-id:2039578349;
mso-list-template-ids:1439879786;}
@list l0:level1
{mso-level-number-format:bullet;
mso-level-text:\F0B7;
mso-level-tab-stop:36.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Symbol;}
@list l0:level2
{mso-level-number-format:bullet;
mso-level-text:o;
mso-level-tab-stop:72.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:"Courier New";
mso-bidi-font-family:"Times New Roman";}
@list l0:level3
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:108.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Wingdings;}
@list l0:level4
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:144.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Wingdings;}
@list l0:level5
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:180.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Wingdings;}
@list l0:level6
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:216.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Wingdings;}
@list l0:level7
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:252.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Wingdings;}
@list l0:level8
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:288.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Wingdings;}
@list l0:level9
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:324.0pt;
mso-level-number-position:left;
text-indent:-18.0pt;
mso-ansi-font-size:10.0pt;
font-family:Wingdings;}
ol
{margin-bottom:0cm;}
ul
{margin-bottom:0cm;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-GB" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal">Hi,<br>
<br>
I'm having an issue where WRF stops in upward_rad in phys/module_sf_bep.F as it's having an issue with the calculated roof temperature.<br>
<br>
I've run successfully before using BEP+NoahLSM for other time periods, but for some reason it now stops for this specific run. It looks like an element of the 'tr' (roof temperature) array is reaching less than 100 K, which is not allowed (according to the
code) and so the run stops. I'm not overly familiar with this routine so I'm not sure what the issue could be, or why tr can't be less than 100. Here’s the relevant snippet of code from the module it stops in.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:0cm;margin-right:0cm;margin-bottom:12.0pt;margin-left:36.0pt">
<span style="font-size:10.0pt;font-family:"Courier New";mso-fareast-language:EN-GB"> real tw(2*ndm,nz_um,nwr_u) ! Temperature in each layer of the wall [K]<br>
real tr(ndm,nz_um,nwr_u) ! Temperature in each layer of the roof [K]<br>
real tg(ndm,ng_u) ! Temperature in each layer of the ground [K]<br>
integer id ! street direction<br>
integer ndu ! number of street directions<br>
!OUTPUT/INPUT<br>
real rs_abs ! absrobed solar radiationfor this street direction<br>
real rl_up ! upward longwave radiation for this street direction<br>
real emiss ! mean emissivity<br>
real grdflx_urb ! ground heat flux<br>
!LOCAL<br>
integer iz,iw<br>
real rl_inc,rl_emit<br>
real gfl<br>
integer ix,iy,iwrong<br>
<br>
iwrong=1<br>
do iz=1,nzu+1<br>
do id=1,ndu<br>
do iw=1,nwr_u<br>
if(tr(id,iz,iw).lt.100.)then<br>
write(*,*)'in upward_rad ',iz,id,iw,tr(id,iz,iw)<br>
iwrong=0<br>
endif<br>
if(tw(2*id-1,iz,iw).lt.100.) then<br>
write(*,*)'in upward_rad ',iz,id,iw,tw(2*id-1,iz,iw)<br>
iwrong=0<br>
endif<br>
if(tw(2*id,iz,iw).lt.100.) then<br>
write(*,*)'in upward_rad ',iz,id,iw,tw(2*id,iz,iw)<br>
iwrong=0<br>
endif<br>
enddo<br>
enddo<br>
enddo<br>
do id=1,ndu<br>
do iw=1,ng_u<br>
if(tg(id,iw).lt.100.) then<br>
write(*,*)'in upward_rad ',id,iw,tg(id,iw)<br>
iwrong=0<br>
endif<br>
enddo<br>
enddo<br>
if(iwrong.eq.0)stop</span><span style="font-size:12.0pt;font-family:"Times New Roman","serif";mso-fareast-language:EN-GB"><o:p></o:p></span></p>
<p class="MsoNormal">I'm running for a period in the summer in the UK, so I wouldn't expect the temperature to be less than 100K(!), but the error file shows that this has happened in some of the roof layers:<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal" style="margin-left:36.0pt"><span style="font-size:10.0pt;font-family:"Courier New";mso-fareast-language:EN-GB"> in upward_rad 6 1 10 97.5907135<br>
in upward_rad 6 2 10 97.5907135<br>
in upward_rad 7 1 10 97.5907135<br>
in upward_rad 7 2 10 97.5907135<br>
in upward_rad 8 1 10 97.5907135<br>
in upward_rad 8 2 10 97.5907135</span><span style="font-size:12.0pt;font-family:"Times New Roman","serif";mso-fareast-language:EN-GB"><o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I have run this with no issues for other summer periods, so I'm very confused about what's going on. Any advice on how to track down the issue would be very helpful.<br>
<br>
Best,<br>
Helen<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b><span style="mso-fareast-language:EN-GB">Dr Helen Macintyre<o:p></o:p></span></b></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB">Senior Environmental Scientist (Air Quality)<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB">Public Health England<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB">Centre for Radiation, Chemical and Environmental Hazards<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB">Chilton, Didcot<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB">OX11 0RQ<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB"><a href="mailto:Helen.Macintyre@phe.gov.uk"><span style="color:blue">Helen.Macintyre@phe.gov.uk</span></a><o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB">+44 (0)1235 825388<o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal" style="line-height:12.0pt;mso-line-height-rule:exactly"><span lang="EN-US" style="font-size:10.0pt;font-family:"Arial","sans-serif""><a href="http://www.gov.uk/phe"><span style="color:blue">www.gov.uk/phe</span></a> Follow us on
<a href="https://twitter.com/PHE_uk"><span style="color:blue">Twitter</span><span style="font-size:12.0pt;color:blue"> @PHE_uk</span></a></span><span lang="EN-US" style="font-size:12.0pt;font-family:"Arial","sans-serif""><o:p></o:p></span></p>
<p class="MsoNormal"><b><span style="font-size:10.0pt;color:black;mso-fareast-language:EN-GB">Protecting and improving the nation’s health</span></b><span style="mso-fareast-language:EN-GB"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="mso-fareast-language:EN-GB"><o:p> </o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<br clear="both">
**************************************************************************<BR>
The information contained in the EMail and any attachments is confidential and intended solely and for the attention and use of the named addressee(s). It may not be disclosed to any other person without the express authority of Public Health England, or the intended recipient, or both. If you are not the intended recipient, you must not disclose, copy, distribute or retain this message or any part of it. This footnote also confirms that this EMail has been swept for computer viruses by Symantec.Cloud, but please re-sweep any attachments before opening or saving. http://www.gov.uk/PHE<BR>
**************************************************************************<BR>
</body>
</html>