Hello. I have been running CESM2.1.3 with B1850 compset but reset to year 2000 on a large computing system of Korea Meteorological Administration. I have encountered the same problem mentioned in a previous message (ERROR: ice: Vertical thermo error) twice. At the first time the model run stopped after simulating 50 years, I resolved the problem by changing a clubb_gamma coefficient a bit as suggested on the board. But the model stopped again after simulating 190 years. This time, I am looking into the source code and tracking down the energy terms.

Among the terms, I found out that the advective heat flux to ocean, fadvocn, is initialized to zero and then not updated with POP or CICE source codes, but that the subroutine conservation_check_vthermo in the ice_therm_vertical.F90 module printed an error message with a value of -570.xxx for fadvocn when the model stopped recently. I wonder where the actual value of fadvocn is calculated/updated after the initialization. I attached the cesm.log file to this message. Any help will be greatly appreciated. Thanks in advance!

Among the terms, I found out that the advective heat flux to ocean, fadvocn, is initialized to zero and then not updated with POP or CICE source codes, but that the subroutine conservation_check_vthermo in the ice_therm_vertical.F90 module printed an error message with a value of -570.xxx for fadvocn when the model stopped recently. I wonder where the actual value of fadvocn is calculated/updated after the initialization. I attached the cesm.log file to this message. Any help will be greatly appreciated. Thanks in advance!