Storing additional fields

[suhasdl]

• CliMT version: 0.16.9
• Python version: 3.7.3
• Operating System: Ubuntu

Description

Would it be possible to store additional fields of interest, say like geopotential or potential temperature during the run. The model uses surface geopotential in some of the modules but not the geopotential. So the model is able to save surface geopotential but not the 3D geopotential field.

Also is there a way to restart the model from certain time step. I was not able to find anything about restart files in the docs.

Thanks,
Suhas

[JoyMonteiro]

I have been thinking about this a bit, and the straightforward way would be to write a DiagnosticComponent which takes in the relevant fields and produces fields like geopotential and potential temperature. This way, the code will work for any model (1D, 3D) without having to replicate code.

Let me think a bit more how this can be done.

It should be quite easy to restart the model from the output. If you read in the data as a Dataset

fields = xarray.open_dataset('time_10.nc')
state = climt.get_default_state([...])
state.update(fields)

should work (not tested) since Dataset provides a dictionary like interface. Note that
this will only work if all fields required as inputs for all components of the model are present
in the nc file.