In preparation for a future PR with diagnostic plots, modify the code that extracts and saves the raw GLM output to extract ice thickness, evaporation, and the number of layers for each timestep in addition to the temperature predictions and ice flags.
Downstream, when the GLM output is compiled for each lake-gcm combo (across all 3 time periods), only write the depth-specific temperature predictions and ice flag variables to the final feather files.
In preparation for a future PR with diagnostic plots, modify the code that extracts and saves the raw GLM output to extract ice thickness, evaporation, and the number of layers for each timestep in addition to the temperature predictions and ice flags.
Downstream, when the GLM output is compiled for each lake-gcm combo (across all 3 time periods), only write the depth-specific temperature predictions and ice flag variables to the final feather files.