<div dir="ltr">Both worked.<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Aug 9, 2021 at 2:35 PM Martin Losch <<a href="mailto:Martin.Losch@awi.de">Martin.Losch@awi.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">which one?<br>
<br>
> On 9. Aug 2021, at 12:24, Mainak Mondal <<a href="mailto:mm10845@nyu.edu" target="_blank">mm10845@nyu.edu</a>> wrote:<br>
> <br>
> Dear Martin,<br>
> <br>
> The solution worked for me. Thanks a lot.<br>
> <br>
> <br>
> On Mon, Aug 9, 2021 at 12:28 PM Mainak Mondal <<a href="mailto:mm10845@nyu.edu" target="_blank">mm10845@nyu.edu</a>> wrote:<br>
> Great! thanks<br>
> <br>
> On Mon, Aug 9, 2021 at 12:20 PM Martin Losch <<a href="mailto:Martin.Losch@awi.de" target="_blank">Martin.Losch@awi.de</a>> wrote:<br>
> PS. I found it: <a href="https://github.com/MITgcm/MITgcm/issues/327" rel="noreferrer" target="_blank">https://github.com/MITgcm/MITgcm/issues/327</a><br>
> <a href="http://mailman.mitgcm.org/pipermail/mitgcm-support/2020-January/012328.html" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/pipermail/mitgcm-support/2020-January/012328.html</a><br>
> <br>
> > On 9. Aug 2021, at 10:16, Martin Losch <<a href="mailto:Martin.Losch@awi.de" target="_blank">Martin.Losch@awi.de</a>> wrote:<br>
> > <br>
> > Hi Mainak,<br>
> > <br>
> > I remember faintly, that in situations like this, the dynamics solver may have difficulties to get started (and there’s an MITgcm-support thread about this, but I could not find it, so I don’t if it got fixed). Please try this:<br>
> > 1. add LSR_mixIniGuess = 2, (or 4). This computes a free drift solution and initialises the LSR solver with this free drift solution and should get it going<br>
> > 2. just use free drift. To do so set<br>
> > SEAICEuseFREEDRIFT = .TRUE., (#define SEAICE_ALLOW_FREEDRIFT only compiles the corresponding code, but doesn’t turn it on).<br>
> > LSR_mixIniGuess = 0, (or comment out again)<br>
> > <br>
> > For your melange, a free drift solution is probably appropriate, but I would be particularly interested in what happens in the first option.<br>
> > <br>
> > Martin<br>
> > <br>
> > <br>
> >> On 9. Aug 2021, at 09:57, Mainak Mondal <<a href="mailto:mm10845@nyu.edu" target="_blank">mm10845@nyu.edu</a>> wrote:<br>
> >> <br>
> >> Hi Martin,<br>
> >> <br>
> >> There is a subglacial plume at the eastern boundary. It is coming up to the top and making a surface expression. There is the u component of velocity that should push the ice-melange. In the SEAICE_OPTIONS , I have already turned on the <br>
> >> #define SEAICE_ALLOW_FREEDRIFT<br>
> >> <br>
> >> In the data.seaice I used<br>
> >> # SEAICE parameters<br>
> >> &SEAICE_PARM01<br>
> >> SEAIcewriteState =.TRUE.,<br>
> >> # SEAICE_initialHEFF= 0.5,<br>
> >> SEAICE_EPS =1.E-8,<br>
> >> SEAICEuseDYNAMICS =.TRUE.,<br>
> >> HeffFile='Heff2.bin',<br>
> >> # AreaFile='Area1.bin',<br>
> >> SEAICEadvHeff=.TRUE.,<br>
> >> SEAICE_no_slip = .TRUE.,<br>
> >> &<br>
> >> <br>
> >> &SEAICE_PARM03<br>
> >> &<br>
> >> <br>
> >> Regards<br>
> >> <br>
> >> On Mon, Aug 9, 2021 at 11:37 AM Martin Losch <<a href="mailto:Martin.Losch@awi.de" target="_blank">Martin.Losch@awi.de</a>> wrote:<br>
> >> Hi Mainak,<br>
> >> <br>
> >> You need some kind of forcing, either by wind or be water. Can you confirm that the water in the grid cells covered by your ice-melange (not clear how you define that in the sea ice model) is moving in some way?<br>
> >> <br>
> >> Martin<br>
> >> <br>
> >>> On 8. Aug 2021, at 10:53, Mainak Mondal <<a href="mailto:mm10845@nyu.edu" target="_blank">mm10845@nyu.edu</a>> wrote:<br>
> >>> <br>
> >>> Hi,<br>
> >>> <br>
> >>> Recently I was trying to advect ice-melange in MITgcm. I used sea-ice package and created an ice-melange over a 5 by 5 grid cell with an Heff=1m while the rest Heff=0. I used the following options in data.seaice<br>
> >>> <br>
> >>> # SEAICE parameters<br>
> >>> &SEAICE_PARM01<br>
> >>> SEAIcewriteState =.TRUE.,<br>
> >>> SEAICE_EPS =1.E-8,<br>
> >>> SEAICEuseDYNAMICS =.TRUE.,<br>
> >>> HeffFile='Heff1.bin',<br>
> >>> SEAICEadvHeff=.TRUE.,<br>
> >>> &<br>
> >>> I am not putting any wind forcing but the melange is placed next to a subglacial plume, which was forced by the obcs package. The ice-melange is melting and sea-ice is growing around it.<br>
> >>> However, the ice-melange remains static at its location throughout. Any clue how I shall be able to advect it? I would really appreciate your help.<br>
> >>> <br>
> >>> Regards,<br>
> >>> <br>
> >>> <br>
> >>> -- <br>
> >>> Mainak Mondal<br>
> >>> Post-Doctoral Associate<br>
> >>> A2 104-C<br>
> >>> Center for Global Sea Level Change<br>
> >>> New York University Abu Dhabi<br>
> >>> ph: +971 26285339<br>
> >>> <br>
> >>> _______________________________________________<br>
> >>> MITgcm-support mailing list<br>
> >>> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
> >>> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
> >> <br>
> >> _______________________________________________<br>
> >> MITgcm-support mailing list<br>
> >> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
> >> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
> >> <br>
> >> <br>
> >> -- <br>
> >> Mainak Mondal<br>
> >> Post-Doctoral Associate<br>
> >> A2 104-C<br>
> >> Center for Global Sea Level Change<br>
> >> New York University Abu Dhabi<br>
> >> ph: +971 26285339<br>
> >> <br>
> >> _______________________________________________<br>
> >> MITgcm-support mailing list<br>
> >> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
> >> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
> > <br>
> > _______________________________________________<br>
> > MITgcm-support mailing list<br>
> > <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
> > <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
> <br>
> _______________________________________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
> <br>
> <br>
> -- <br>
> Mainak Mondal<br>
> Post-Doctoral Associate<br>
> A2 104-C<br>
> Center for Global Sea Level Change<br>
> New York University Abu Dhabi<br>
> ph: +971 26285339<br>
> <br>
> <br>
> <br>
> -- <br>
> Mainak Mondal<br>
> Post-Doctoral Associate<br>
> A2 104-C<br>
> Center for Global Sea Level Change<br>
> New York University Abu Dhabi<br>
> ph: +971 26285339<br>
> <br>
> _______________________________________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
<br>
_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><b>Mainak Mondal</b><div>Post-Doctoral Associate</div><div>A2 104-C</div><div>Center for Global Sea Level Change</div><div>New York University Abu Dhabi</div><div>ph: +971 26285339</div><div><br></div></div></div></div></div>