Understing the vertical levels

[amhayes231]

• CliMT version: 0.17.12
• Python version: 3.10.13
• Operating System: EndeavourOS

Hello again! I am looking at using climt with a realistic orography. If I am to update the surface_geopotential variable, will this automatically adjust accordingly for the sigma pressure hybrid coordinate system? If not, how should I proceed with updating the model to properly accommodate for the new surface geopotential.

For context, please see the code snippet below, will this suffice, or will further work need to be done to accommodate this update?

dycore = GFSDynamicalCore(moist=True)
grid = climt.get_grid(nx=128, ny=62, nz=13, n_ice_interface_levels=None)
my_state = climt.get_default_state([dycore], grid_state=grid, n_ice_interface_levels=None)

# ...

my_state['surface_geopotential'].values = new_surf_geo.values
my_state.update()