[MITgcm-devel] newest code giving nans

Matthew Mazloff mmazloff at MIT.EDU
Wed Jul 29 12:24:36 EDT 2009


Hello,

Just updated the newest code (older one was a couple months old) and  
ran it (with no sea-ice) and got nans at time-step 1.

One clue is that the time_tsnumber was not properly incrementing...for  
example, from timestep 0 to 1 STDOUT wrote:


(PID.TID 0000.0001) %MON vort_a_mean                  =    
7.7907401369469E-05
(PID.TID 0000.0001) %MON vort_a_sd                    =    
7.5333276615471E-06
(PID.TID 0000.0001) %MON vort_p_mean                  =    
8.1856709185886E-05
(PID.TID 0000.0001) %MON vort_p_sd                    =    
3.4115931535963E-05
(PID.TID 0000.0001) %MON surfExpan_theta_mean         =    
7.5944713547936E-06
(PID.TID 0000.0001) %MON surfExpan_salt_mean          =    
4.5528081632941E-06
(PID.TID 0000.0001) //  
=======================================================
(PID.TID 0000.0001) // End MONITOR dynamic field statistics
(PID.TID 0000.0001) //  
=======================================================
  cg2d: Sum(rhs),rhsMax =   1.45330722567394E-07  4.35515400006567E-01
(PID.TID 0000.0001) //  
=======================================================
(PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
(PID.TID 0000.0001) //  
=======================================================
(PID.TID 0000.0001) %MON time_tsnumber                 
=                     0
(PID.TID 0000.0001) %MON time_secondsf                 
=                   NaN
(PID.TID 0000.0001) %MON dynstat_eta_max               
=                   NaN
(PID.TID 0000.0001) %MON dynstat_eta_min               
=                   NaN
(PID.TID 0000.0001) %MON dynstat_eta_mean              
=                   NaN
(PID.TID 0000.0001) %MON dynstat_eta_sd                
=                   NaN
(PID.TID 0000.0001) %MON dynstat_eta_del2              
=                   NaN
(PID.TID 0000.0001) %MON dynstat_uvel_max              
=                   NaN

However the model ran the full simulation and properly incremented  
many things (e.g. dayrec) and the timer stats looked fine at the end  
of the run.
and notice initial (timestep 0) cg2d numbers are fine

So any idea what may have happened?  Could this possibly be a TAF  
thing...

Thanks
Matt



More information about the MITgcm-devel mailing list