Understing the vertical levels
Opened this issue · 0 comments
amhayes231 commented
- 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()