[MITgcm-support] advection of ice-melange

[Mainak Mondal]

Both worked.

On Mon, Aug 9, 2021 at 2:35 PM Martin Losch <Martin.Losch at awi.de> wrote:

> which one?
>
> > On 9. Aug 2021, at 12:24, Mainak Mondal <mm10845 at nyu.edu> wrote:
> >
> > Dear Martin,
> >
> > The solution worked for me. Thanks a lot.
> >
> >
> > On Mon, Aug 9, 2021 at 12:28 PM Mainak Mondal <mm10845 at nyu.edu> wrote:
> > Great! thanks
> >
> > On Mon, Aug 9, 2021 at 12:20 PM Martin Losch <Martin.Losch at awi.de>
> wrote:
> > PS. I found it: https://github.com/MITgcm/MITgcm/issues/327
> >
> http://mailman.mitgcm.org/pipermail/mitgcm-support/2020-January/012328.html
> >
> > > On 9. Aug 2021, at 10:16, Martin Losch <Martin.Losch at awi.de> wrote:
> > >
> > > Hi Mainak,
> > >
> > > 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:
> > > 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
> > > 2. just use free drift. To do so set
> > > SEAICEuseFREEDRIFT = .TRUE., (#define SEAICE_ALLOW_FREEDRIFT only
> compiles the corresponding code, but doesn’t turn it on).
> > > LSR_mixIniGuess = 0, (or comment out again)
> > >
> > > For your melange, a free drift solution is probably appropriate, but I
> would be particularly interested in what happens in the first option.
> > >
> > > Martin
> > >
> > >
> > >> On 9. Aug 2021, at 09:57, Mainak Mondal <mm10845 at nyu.edu> wrote:
> > >>
> > >> Hi Martin,
> > >>
> > >> 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
> > >> #define SEAICE_ALLOW_FREEDRIFT
> > >>
> > >> In the data.seaice I used
> > >> # SEAICE parameters
> > >> &SEAICE_PARM01
> > >> SEAIcewriteState =.TRUE.,
> > >> # SEAICE_initialHEFF= 0.5,
> > >> SEAICE_EPS =1.E-8,
> > >> SEAICEuseDYNAMICS =.TRUE.,
> > >> HeffFile='Heff2.bin',
> > >> # AreaFile='Area1.bin',
> > >> SEAICEadvHeff=.TRUE.,
> > >> SEAICE_no_slip = .TRUE.,
> > >> &
> > >>
> > >> &SEAICE_PARM03
> > >> &
> > >>
> > >> Regards
> > >>
> > >> On Mon, Aug 9, 2021 at 11:37 AM Martin Losch <Martin.Losch at awi.de>
> wrote:
> > >> Hi Mainak,
> > >>
> > >> 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?
> > >>
> > >> Martin
> > >>
> > >>> On 8. Aug 2021, at 10:53, Mainak Mondal <mm10845 at nyu.edu> wrote:
> > >>>
> > >>> Hi,
> > >>>
> > >>> 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
> > >>>
> > >>> # SEAICE parameters
> > >>> &SEAICE_PARM01
> > >>> SEAIcewriteState =.TRUE.,
> > >>> SEAICE_EPS =1.E-8,
> > >>> SEAICEuseDYNAMICS =.TRUE.,
> > >>> HeffFile='Heff1.bin',
> > >>> SEAICEadvHeff=.TRUE.,
> > >>> &
> > >>> 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.
> > >>> 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.
> > >>>
> > >>> Regards,
> > >>>
> > >>>
> > >>> --
> > >>> Mainak Mondal
> > >>> Post-Doctoral Associate
> > >>> A2 104-C
> > >>> Center for Global Sea Level Change
> > >>> New York University Abu Dhabi
> > >>> ph: +971 26285339
> > >>>
> > >>> _______________________________________________
> > >>> MITgcm-support mailing list
> > >>> MITgcm-support at mitgcm.org
> > >>> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> > >>
> > >> _______________________________________________
> > >> MITgcm-support mailing list
> > >> MITgcm-support at mitgcm.org
> > >> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> > >>
> > >>
> > >> --
> > >> Mainak Mondal
> > >> Post-Doctoral Associate
> > >> A2 104-C
> > >> Center for Global Sea Level Change
> > >> New York University Abu Dhabi
> > >> ph: +971 26285339
> > >>
> > >> _______________________________________________
> > >> MITgcm-support mailing list
> > >> MITgcm-support at mitgcm.org
> > >> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> > >
> > > _______________________________________________
> > > MITgcm-support mailing list
> > > MITgcm-support at mitgcm.org
> > > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> >
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
> >
> >
> > --
> > Mainak Mondal
> > Post-Doctoral Associate
> > A2 104-C
> > Center for Global Sea Level Change
> > New York University Abu Dhabi
> > ph: +971 26285339
> >
> >
> >
> > --
> > Mainak Mondal
> > Post-Doctoral Associate
> > A2 104-C
> > Center for Global Sea Level Change
> > New York University Abu Dhabi
> > ph: +971 26285339
> >
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>


-- 
*Mainak Mondal*
Post-Doctoral Associate
A2 104-C
Center for Global Sea Level Change
New York University Abu Dhabi
ph: +971 26285339
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20210809/371dc4a0/attachment.html>