<html 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=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:DengXian;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"\@DengXian";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style></head><body lang=EN-US link="#0563C1" vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-size:11.0pt'>Hi all, <o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt'>I think the problem has been solved. I should have a reasonable value of input for shortwave and longwave downward radiation. <o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt'>Thanks,<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt'>Yilang<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Times New Roman",serif;color:black'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt'><o:p> </o:p></span></p><div style='border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in'><p class=MsoNormal><b><span style='color:black'>From: </span></b><span style='color:black'>MITgcm-support <mitgcm-support-bounces@mitgcm.org> on behalf of Yilang Xu <yxu@whoi.edu><br><b>Reply-To: </b><mitgcm-support@mitgcm.org><br><b>Date: </b>Friday, November 15, 2019 at 06:06<br><b>To: </b>"mitgcm-support@mitgcm.org" <mitgcm-support@mitgcm.org><br><b>Subject: </b>[MITgcm-support] Weird cooling with simple EXF setup<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:11.0pt'><o:p> </o:p></span></p></div><p class=MsoNormal><span style='font-size:11.0pt'>Hi all, </span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'> </span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>I am testing EXF package, and am using the default EXF_OPTIONS.h file. In my data.exf, I have the following input, </span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>Air temperature 10 degC all over the domain, with “exf_offset_atemp = 273.15” added;</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>Humidity 0.0001 all over the domain;</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>Vwind ramps from 0 to 20m/s within certain area; </span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>For precipfile, lwdownfile, swdownfile, runoffFile, I have input files with value 0 all over the domain. </span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'> </span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>The ocean is initialized with 5 degC in data. During the runs, the ocean keeps cooling, with strong positive hflux and lwflux in STDOUT.0000. I checked exf_atemp in STDOUT, it is 283.15. So I am not sure why the ocean is getting colder even the air is warmer. </span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'> </span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>My data.exf and STDOUT files are uploaded<b> <a href="https://drive.google.com/open?id=1H-Q-YZwYXdxay6jUXdfysjF1YHpL9mQ8">here</a></b> through Google drive. </span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>Appreciate it if anyone would help with this problem. </span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'> </span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>Thanks,</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>Yilang</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'> </span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>File input within data.exf: </span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'> hfluxfile = ' ',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>sfluxfile = ' ',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>ustressfile = ' ',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>vstressfile = ' ',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>atempfile = 'airtempfield.bin',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>aqhfile = 'humidfield.bin',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>uwindfile = 'empty2dfield.bin',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>vwindfile = 'vwind.bin',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'># evapfile = 'empty2dfield.bin',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>precipfile = 'empty2dfield.bin',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>lwfluxfile = ' ',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>swfluxfile = ' ',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>lwdownfile = 'empty2dfield.bin',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>swdownfile = 'empty2dfield.bin',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>runoffFile = 'empty2dfield.bin',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>climsstfile = ' ',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>climsssfile = ' ',</span><o:p></o:p></p><p class=MsoNormal><span style='font-size:11.0pt'>_______________________________________________ MITgcm-support mailing list MITgcm-support@mitgcm.org http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support <o:p></o:p></span></p></div></body></html>