I am seting up experiments looking at future changes in SST and their impact on the atmospheric circulation. The component set I generated is: RCP8_CAM5_CLM40%CN_DICE%PRES_DOCN%DOM_RTM_SGLC_SWAV. and I am conducting hybrid runs from output from the CESM-LE at 2046. My boundary conditions are created using the sea ice and sea surface temperatures from the CESM-LE (ensemble member 1). So my results should begin from a very similar state. When comparing my hybrid run the day after the run starts to the original CESM-LE we see anomalies in surface temperature immediately over the sea ice with colder anomalies where sea ice should be thin and warm anomalies where it should be thick. This seems to be caused by the prescribed mode assuming that the sea ice thickness is a constant at 2m. I found the location in the source code (ice_prescribed_mod.F90 line 459-463) where 2m is prescribed and I have done test changing the value and the results support this idea. These biases due to constant sea ice thickness, I believe, might be important sources for errors for my experiments so I am looking to use sea ice thickness instead.My question is whether this has ever been done before? Is it feasible to replace the contant thickness value with the sea ice thickness at each location by loading in a new variable or might there be other potential issues?Thanks so much!Melissa For reference sake I have uploaded images of the model run with 2m (Cnt_to_CESM-LETS_20460102.png), 1m (Cnt_hi1_to_CESM-LETS_20460102.png), and 4m (Cnt_hi4_to_CESM-LETS_20460102.png) ice thickness. The top figure is the CESM-LE the middle is my hybrid run and the bottom is the difference.