I added code that initializes the snow, blue ice, and white ice thickness. This "hot start" is required for data assimilation applications. Some comments:
1) Initializing with non-zero values for snow thickness causes segmentation faults
2) This implementation requires values for the three thicknesses. Therefore the glm3.nml must have values for snow_thickness, white_ice_thickness, blue_ice_thickness. This would required an update to all the example glm3.nmls. We may want to make these values optional (and default to zero if not provided) but I am not sure how best to do that.
I added code that initializes the snow, blue ice, and white ice thickness. This "hot start" is required for data assimilation applications. Some comments:
1) Initializing with non-zero values for snow thickness causes segmentation faults 2) This implementation requires values for the three thicknesses. Therefore the glm3.nml must have values for snow_thickness, white_ice_thickness, blue_ice_thickness. This would required an update to all the example glm3.nmls. We may want to make these values optional (and default to zero if not provided) but I am not sure how best to do that.