-
Notifications
You must be signed in to change notification settings - Fork 54
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Temp/salinity and thickness bug when initializing with an ice shelf #379
Comments
Verifying that this issue was resolved by 88d9eb7 and 25989ad With the current dev/gfdl branch, the following parameters work well for me for ALE ice shelf initialisation:
(may need to increase this tolerance depending on the geometry if model crashes (#609), but want it to be as small as possible) Additionally, I use:
Folder with some plots and configs showing the initialisation verification here if useful to someone in the future! |
Hello, I've come across strange features in the initialization of temperature/salinity and layer thickness/interface levels when an ice shelf is present. I initialized the model with a simple triangular ice shelf and
![image](https://private-user-images.githubusercontent.com/61528379/245993927-59895c49-778f-4bb5-a059-bea1a65220bd.png?jwt=eyJhbGciOiJIUzI1NiIsInR5cCI6IkpXVCJ9.eyJpc3MiOiJnaXRodWIuY29tIiwiYXVkIjoicmF3LmdpdGh1YnVzZXJjb250ZW50LmNvbSIsImtleSI6ImtleTUiLCJleHAiOjE3MjU3NDU2NDEsIm5iZiI6MTcyNTc0NTM0MSwicGF0aCI6Ii82MTUyODM3OS8yNDU5OTM5MjctNTk4OTVjNDktNzc4Zi00YmI1LWEwNTktYmVhMWE2NTIyMGJkLnBuZz9YLUFtei1BbGdvcml0aG09QVdTNC1ITUFDLVNIQTI1NiZYLUFtei1DcmVkZW50aWFsPUFLSUFWQ09EWUxTQTUzUFFLNFpBJTJGMjAyNDA5MDclMkZ1cy1lYXN0LTElMkZzMyUyRmF3czRfcmVxdWVzdCZYLUFtei1EYXRlPTIwMjQwOTA3VDIxNDIyMVomWC1BbXotRXhwaXJlcz0zMDAmWC1BbXotU2lnbmF0dXJlPTdmNzNjNWFhYTY2NzY0YzVkOWU4ZGQ2ZDJjMjZmNDZkMjM2NDI3NTc0ZjA3ZWU3YWJhYjVjMzE4ZDgzMDkzYTgmWC1BbXotU2lnbmVkSGVhZGVycz1ob3N0JmFjdG9yX2lkPTAma2V5X2lkPTAmcmVwb19pZD0wIn0.-Gc4ztlqajUWHY-YEOSdEHDwFgOV2Rp9qx1qhZXTaHY)
TS_CONFIG = "ISOMIP"
andTHICKNESS_CONFIG = "ISOMIP"
(though the same happens with other config options, e.g. TS=linear or Thickness=uniform) which is meant to produce a linear stratification according to given temperature and salinity bounds. However, although it looks mostly linear, there are strange periodic deviations from linearity in the initial conditions file which are highlighted below in this figure, showing the anomaly of the salt initial conditions from the theoretical linear profile.Sea level/layer thicknesses also have very slight perturbations from their expected smooth distributions (for sigma_shelf_zstar coordinates) when
TRIM_IC_FOR_P_SURF
is not used (turning it on seems to fix those discontinuities but not the T/S ones). When I run the model, horizontal velocities emerge very quickly with the same spatial frequency as the salt anomalies suggesting that the anomalous horizontal density gradients drive this unwanted flow.A temporary fix: I can get around this issue by using netcdf files as input for the thicknesses (with ice shelf accounted for), and by asking the model to skip the
calc_sfc_displacement
step inMOM_state_initialization.F90
. This somehow results in a linear T/S profile and very low (order 10^(-8) m/s) velocities for the simple triangular ice shelf, however requires running the model twice. If it is helpful, I have linked a github repo that includes the files/MOM_input/override used in my test cases.The text was updated successfully, but these errors were encountered: