Looking for a python equivalent to MITgcm/verification/isomip/input/gendata.m

[arrran]

Hullo,

I use python to run MITgcm with shelf ice, but havent managed to find a python script to produce a SHELFICEloadAnomalyFile (aka input/phi0surf)
I have been running MITgcm/verification/isomip/input/gendata.m from python, which is slow.

Are you aware of a python script which produces a SHELFICEloadAnomalyFile from ice shelf topography and bathymetry?

Thanks!
arran

[knaughten]

I've previously written up this explanation of my domain generation tools which might help you: https://github.com/knaughten/mitgcm_python/wiki/Creating-a-new-MITgcm-domain
It includes the function calc_load_anomaly which I think is what you are looking for:

mitgcm_python/ics_obcs.py

Line 335 in d4c9b00

def calc_load_anomaly (grid, out_file, option='constant', ini_temp_file=None, ini_salt_file=None, ini_temp=None, ini_salt=None, constant_t=-1.9, constant_s=34.4, eosType='MDJWF', rhoConst=1035, tAlpha=None, sBeta=None, Tref=None, Sref=None, hfac=None, prec=64, check_grid=True):

Good luck!
Kaitlin

[arrran]

Ah great exactly this, thanks very much!

I will adjust it to output SHELFICEloadAnomalyFile without running MITgcm