[MITgcm-support] ISOMIP melt rate settings

[Martin Losch]

Hi Kaylie,

the if’s decide between “ISOMIP” thermodynamics and the 3-equation thermodynamics. Unfortunately there is no code to prescribe the meltrate.

In data.shelfice you can select useISOMIPTD=.TRUE. (default is false).

I would useISOMIPTD=.TRUE., because that code is simpler, and then set shelfIceFreshWaterFlux(I,J,bi,bj) in l568 to your favourite value.

setting SHELFICEheatTransCoeff=0. turns off the heat flux,  and SHELFICEconserve=.FALSE. turns off the heat flux associated with the fresh water flux.

Martin



On 18. Nov 2021, at 19:12, Kaylie Cohanim <kcohanim at princeton.edu> wrote:
> 
> Hi all,
> 
> I'm a second-year graduate student, currently working on an ISOMIP model configuration. The configuration is running well, but I would now like to apply a constant melt rate, rather than use the 3-equations melt rate formula that is currently in place.
> 
> I found a few lines in mitgcm/pkg/shelfice/shelfice_thermodynamics.F where I can change the heat and freshwater flux,
> 
> shelfIceHeatFlux(I,J,bi,bj) = 0. _d 
> 0
> shelfIceFreshWaterFlux(I,J,bi,bj) 
> = 0. _d 
> 0
> 
> shelficeForcingT(I,J,bi,bj) = 0. _d 
> 0
> shelficeForcingS(I,J,bi,bj) 
> = 0. _d 0 .
> 
> 
> 
> However, these lines show up multiple times in different "if" statements and I don't quite understand which block should be changed. I'm also wondering if there's a simpler way to do this, like a line I can add to the data script, rather than manually setting it?
> 
> Thank you!
> 
> Kaylie
> 
>  
> _______________________________________________
> 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/20211119/c4ce3407/attachment.p7s>