[MITgcm-support] advection of ice-melange

[Martin Losch]

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

-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 1665 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20210809/2a930aa2/attachment-0001.p7s>