Main menu


vertical thermo error

4 posts / 0 new
Last post
vertical thermo error

Dear all,

I am working with CESM1.2.2, E1850CN compset  (CAM4, CLM4, and a Slab-ocean).

I would like to remove the diurnal and seasonal cycle of solar forcing, so I modified files "zenith.F90", "radsw.F90", and "shr_orb_mod.F90" under case directory.  

Other settings are unchanged. I started with a "startup" run and the model crashed at year3, month 11 every time, with an error saying:


Starting thermo, Tsn < Tmin

 Tsn=  -100.007281731455     

 Tmin=  -100.000000000000     

 istep1, my_task, i, j:     1146098          83           3           6

 qsn  -179995997.739214     

 istep1, my_task, iblk =     1146098          83           2

 Global block:         986

 Global i and j:         252         365

 Lat, Lon:   72.6737451597182       -91.8653354459640     

(shr_sys_abort) ERROR: ice: Vertical thermo error

(shr_sys_abort) WARNING: calling shr_mpi_abort() and stopping

application called MPI_Abort(MPI_COMM_WORLD, 1001) - process 83


I found someone else also posted such error before, in which he suggested to increase ocean time step to get past this error. For slab-ocean, where could I modify the time step?  

And, i am thinking is there a code-based solution to solve this problem now?   





This message indicates that your surface temperature is dropping below -100C. So, there is something very wrong in the forcing and the modifications that you have done.



Dear dbailey,

Sorry to reply late.

I removed seasonal and diurnal cycle of solar insolation (no day and night, constant for each grid, only change with latitude) and fix SOLIN and zenith angle to their annual mean values. And, I replace all vegetations and urban type grid with bare ground.

Above is what I have done, and attached files are solar insolation and zenith angle I used.

Do you have any suggestions? Thanks! 


PS, I wonder why Tmin = -100 ?   Can I change this value lower?  or can I use prescribed sea ice mode to avoid this error?





Did you look at your sea ice? Sounds like you are cooling the Earth so much that the ice is running away. You can try an F compset with prescribed SST and ice fraction to see if the atmosphere and land can run with this forcing. You cannot run the slab ocean (SOM) with prescribed sea ice.

Log in or register to post comments

Who's new

  • yunqian.zhu@...
  • acostar@...
  • Daniele.visioni@...
  • federico
  • shreya.dhame@...