[MITgcm-support] advection of ice-melange

[Mainak Mondal]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20210809/91894282/attachment.html>