criterion

[redhat007]

hi
would you please tell me what are zeta, salt and temp criterion for blowing up ?
and usually what parameters should be checked for finding blow up reason?
cheers
redhat

[kate]

You can look in diag.F to see the blow-up criteria. One is a velocity larger than max_speed, the other is a density larger than max_rho. Both max_speed and max_rho are set in mod_scalars.F. The logic used to find the maximum velocity can be copied to zeta, T and S, but it has not been done in the trunk code.

If you are running my ice branch, there's a check in ice_frazil.F that will trip over odd values of T like -300. I have also made diag.F more verbose so that I get a report like (just this morning, in fact):

Code: Select all

 239325 33208 21:30:00  3.608264E-03  1.703917E+04  1.703917E+04  3.289577E+16
         (543,0266,50)  6.686609E-04  1.023531E-01  9.435409E+00  7.383824E+01
 DIAG speed trouble     73.83824443803645     

 Blowing-up: Saving latest model state into  RESTART file 

 DIAG speed ijk     88.75332729210966               404           71           25