Hello,

I'm trying to run CESM1.2.2 with F1850W compset. And I want to increase the SST by 30 degrees centigrade by changing the SST_cpl and ice_cov data in "sst_HadOIBl_bc_1.9x2.5_clim_pi_c101028.nc".

But when I want to run this model, it crashes with the following error(see log_file attached for more imformation):

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

DEPTH exceeding table value

findsp not converging at i = 11

t, q, p 236.970523240216 3.735770877134584E-002

6152.05000000000

tsp, qsp 263.787360155026 2.971450329655875E-002

(shr_sys_abort) ERROR: wv_saturation::FINDSP -- not converging

(shr_sys_abort) WARNING: calling shr_mpi_abort() and stopping

application called MPI_Abort(MPI_COMM_WORLD, 1001) - process 37

yhrun: error: cn7068: task 37: Exited with exit code 233

yhrun: First task exited 60s ago

yhrun: tasks 0-36,38-119: running

yhrun: task 37: exited abnormally

yhrun: Terminating job step 4588727.0

yhrun: Job step aborted: Waiting up to 2 seconds for job step to finish.

slurmd[cn7067]: *** STEP 4588727.0 KILLED AT 2018-01-17T12:27:46 WITH SIGNAL 9 ***

slurmd[cn7067]: *** STEP 4588727.0 KILLED AT 2018-01-17T12:27:46 WITH SIGNAL 9 ***

yhrun: error: cn7067: task 22: Killed

Can anyone tell me why the model crashed with the error message "DEPTH exceeding table value", and I also run a model with SST increased by 10 degrees centigrade and it runs fine.Any help is highly appreciated.

Shaozhi Lin

Mike MillsWACCM LiaisonAtmospheric Chemistry Division

NCAR Foothills Lab

Boulder, Colorado USA

Hi Shoazhi,

I haven't looked in your exact code branch, but have some suggestions for you.

1) I am not sure about the "DEPTH" error massage. I don't see where it is generated in the code (at least looking at the latest branch), but I suspect it is coming from the estbl lookup table used to calculate water saturation. I suspect you are going past the ends of this table because of your extreme conditions.

2) If it wasn't for your extreme conditions and what I think are table lookup errors, you could try increasing the number of iterations to alow it to try to converge. However, once you hit the edge of the table, I don't think you will be able to converge.

3) The fundamental problem seems like that your model state has gotten very wet. It looks like it failed at 61 hPa and a temperature of 263K with a specific humidity of 3.7e-2 kg/kg. This is a very large value for the stratosphere and for this temperature. Do you see the model state getting wetter and wetter?

How abruptly did you change the surface temperature? If you did it quickly, you might try introducing the change more slowly to allow the model state to adjust. If the stratosphere is getting so much moisture that it is above saturation, you may need to increase the altitude at which cloud microphysics is applied. This is controlled by a namelist variable that I think is set to 1 hPa in your version. You might try setting it to a lower value (beyond where you see the relative humidity exceeding 0.8), to allow clouds and precipitation to remove some of the excess humidity. The varaible you should be looking for is trop_cloud_top_press. and it is set in Pa. So you might try adding the following to your user_nl_cam:

trop_cloud_top_press = 1.d-4