<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" style="display:none"><!--P{margin-top:0;margin-bottom:0;} --></style>

</head>

<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">

<div>

<div>

<div>

<div>

<div>

<div>Hi everyone,<br>

<br>

</div>

I use mnc/exf/seaice packages. I want to "control" my freshwater flux, sflux, so I prescribe the precipitation and evaporation fields such that my sflux (=evap-precip) integrates to zero (or at least integrates to an order of 10^(-11)) at the surface.

<br>

<br>

If I look at my results, my salt in my basin increases all the time during my simulation. I compared the sflux I generate with Matlab in binary file (that is evap-precip) with my EXFempmr diag and they differ to an order of 10^(-23) (my sflux generated with

 Matlab integrates to 10^(-11) while my sflux from EXF integrates to 10^(-9)). <br>

<br>

</div>

To generate my forcing files, I have tried:<br>

fid=fopen('precip.bin','w','<wbr>ieee-be.l64'); fwrite(fid,precip,'real*8'); fclose(fid);<br>

</div>

and<br>

fid=fopen('precip.bin','w','<wbr>ieee-be'); fwrite(fid,precip,'real*8'); fclose(fid);<br>

<br>

</div>

In data, PARM01, I have:<br>

readBinaryPrec=64,<br>

writeBinaryPrec=64, <br>

<br>

</div>

In data.exf, EXF_NML_01, I have:<br>

exf_iprec         = 64,<br>

<br>

</div>

I have to admit that my difference is smaller since I have used double precision (first I was using real*4 and precision 32, and my EXFempmr integrated to 10^(-3) which is very bad!) but I was wondering if I could get exactly in my diags what I generated with

 Matlab, and where does this difference between my sflux generated with Matlab and my diags comes from.<br>

<div><br>

</div>

<div>If you have ideas/suggestions it would be very cool! By the way, I compile with ifort v15 using build-option/linux_ia64_ifort.<br>

<br>

</div>

<div>I wish you a wonderful year 2018 :)<br>

<br>

</div>

Sandy

<p><br>

</p>

<p>ps: I get same difference without mnc package<br>

</p>

</body>

</html>