search

knaughten / mitgcm_python

Python scripts designed for my Weddell Sea and Amundsen Sea configurations of MITgcm.
25 stars 17 forks source link

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

Open arrran opened 3 years ago

arrran commented 3 years ago

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 commented 3 years ago

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: https://github.com/knaughten/mitgcm_python/blob/d4c9b00161d9fd0ccd3dac9d98ad2c308b9e4987/ics_obcs.py#L335 Good luck! Kaitlin

arrran commented 3 years ago

Ah great exactly this, thanks very much!

I will adjust it to output SHELFICEloadAnomalyFile without running MITgcm