<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hello all,<br class=""><br class="">I am trying to run a channel model configuration where I have a cold, freshwater input over half the domain that is constant in time. To do this I have turned on EmPmR with <br class="">EmPmRFile = 'M_dT5e+00C.bin’,<div class="">The file is a 2-D array of flux values in m/s where values<=0 to indicate the addition of freshwater to the system.<br class=""><br class="">I also want my surface fresh water flux to have a temperature of 0 deg C so I include this line in my data file:<br class="">temp_EvPrRn=0.,<div class=""><br class=""></div><div class="">Over the first 5 days, salinity in the surface layer looks ok, with values of 0 psu, however salinity in the top grid cells become increasingly negative after this point. At the same time, the temperature of the surface layer in the region of forcing drops below 0 deg C after a day or so becoming unrealistically cold (-15 deg C) by day 5. </div><div class=""><br class=""></div><div class="">If any one can offer some guidance for where I’ve gone wrong, I would greatly appreciate it. The full set up for the simulation can be found here:</div><div class=""><a href="https://github.com/marionalberty/MITgcm/tree/front_relax_SODA2/verification/front_surfFlux_highres" class="">https://github.com/marionalberty/MITgcm/tree/front_relax_SODA2/verification/front_surfFlux_highres</a></div><div class=""><br class=""></div><div class="">Thanks in advance for any assistance,</div><div class=""><br class=""><div class="">Marion Alberty (she/her), Ph.D.<br class="">Postdoctoral Research Associate<br class="">Princeton University<br class="">NOAA Geophysical Fluid Dynamics Laboratory<br class=""></div><br class=""></div></div></body></html>