<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:Helvetica;
        panose-1:0 0 0 0 0 0 0 0 0 0;}
@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:"Times New Roman \(Body CS\)";
        panose-1:2 2 6 3 5 4 5 2 3 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:11.0pt;
        font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
span.apple-converted-space
        {mso-style-name:apple-converted-space;}
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=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-size:12.0pt'>Thank you Matt! I have not noticed this before in the seaice pkg. <o:p></o:p></span></p><p class=MsoNormal><span style='font-size:12.0pt'>I will check </span>seaice_tracer_phys.F to see if I have further questions. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Best,<o:p></o:p></p><p class=MsoNormal>Yilang <span style='font-size:12.0pt'><o:p></o:p></span></p><p class=MsoNormal><span style='font-size:12.0pt'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:12.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='font-size:12.0pt;color:black'>From: </span></b><span style='font-size:12.0pt;color:black'>MITgcm-support <mitgcm-support-bounces@mitgcm.org> on behalf of Matthew Mazloff <mmazloff@ucsd.edu><br><b>Reply-To: </b><mitgcm-support@mitgcm.org><br><b>Date: </b>Thursday, April 2, 2020 at 23:47<br><b>To: </b><mitgcm-support@mitgcm.org><br><b>Subject: </b>Re: [MITgcm-support] change ptracers surface BC during runtime<o:p></o:p></span></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><p class=MsoNormal>Hi Yilang<o:p></o:p></p><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>The sea ice can take up tracer when formed and release when melted. This is turned on with<o:p></o:p></p></div><div><p class=MsoNormal>define ALLOW_SITRACER in SEAICE_OPTIONS.h<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>You will likely want to customize seaice_tracer_phys.F<o:p></o:p></p></div><div><p class=MsoNormal>to tell it which tracer to take up and what fraction should be taken up vs rejected. <o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>If this is what you are interested in doing I can give more details about customizing seaice_tracer_phys.F <o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Matt<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><div><p class=MsoNormal><br><br><o:p></o:p></p><blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><div><p class=MsoNormal>On Apr 2, 2020, at 8:36 PM, Yilang Xu <<a href="mailto:yxu@whoi.edu">yxu@whoi.edu</a>> wrote:<o:p></o:p></p></div><p class=MsoNormal><o:p> </o:p></p><div><div><p class=MsoNormal><span style='font-size:12.0pt'>Hi everyone,<span class=apple-converted-space> </span><o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt'>I have been using ptracers, but only with RBCS to restore the surface BC as a source. However, this only allows a fixed location for surface input.<span class=apple-converted-space> </span><o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt'>I am wondering if the following is possible, potentially with/without RBCS.<span class=apple-converted-space> </span><o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt'>1. Change the locations of surface source of ptracers while the code is running. The locations are based on other output.<span class=apple-converted-space> </span><o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt'>2. If the above is plausible, link the ptracers surface BC with the continuous changes of sea ice cover.<span class=apple-converted-space> </span><o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt'>These might need modifications in the code. My goal is to enforce ptracers=1 on the surface, when the sea ice cover percentage at that location is below a certain level. Is there any reference to it? I see some words mentioning this in gchem pkg, but I am unfamiliar with it.<span class=apple-converted-space> </span><o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt'>Thank you very much.<span class=apple-converted-space> </span><o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt'>Best,<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt'>Yilang<span class=apple-converted-space> </span><o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:12.0pt'> <o:p></o:p></span></p></div><p class=MsoNormal><span style='font-size:9.0pt;font-family:Helvetica'>_______________________________________________<br>MITgcm-support mailing list<br></span><a href="mailto:MITgcm-support@mitgcm.org"><span style='font-size:9.0pt;font-family:Helvetica'>MITgcm-support@mitgcm.org</span></a><span style='font-size:9.0pt;font-family:Helvetica'><br></span><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support"><span style='font-size:9.0pt;font-family:Helvetica'>http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</span></a><o:p></o:p></p></div></blockquote></div><p class=MsoNormal><o:p> </o:p></p></div><p class=MsoNormal>_______________________________________________ MITgcm-support mailing list MITgcm-support@mitgcm.org http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support <o:p></o:p></p></div></body></html>