[MITgcm-support] cg2d_init_res=NaN, cg2d_res=NaN

Julie Deshayes jdeshayes at whoi.edu
Tue Feb 27 14:48:17 EST 2007


Dear MITgcm users,

I am currently trying to make my first run with MITgcm but all the 
output data are NaN, except for the first timestep which conforms to the 
initial conditions. Note that my configuration passed all the 
verification tests (./testreport -ieee) except for 
global_ocean.cs32x15.thsice . The model runs without error (STDERR.0000 
file is empty) and here are a few lines of the output.txt file :

 (PID.TID 0000.0001) Start initial hydrostatic pressure computation
(PID.TID 0000.0001) Pressure is predetermined for buoyancyRelation OCEANIC
(PID.TID 0000.0001)
(PID.TID 0000.0001)  MDSREADFIELD: opening global file: Qnet.bin
(PID.TID 0000.0001) // 
=======================================================
(PID.TID 0000.0001) // Model current state
(PID.TID 0000.0001) // 
=======================================================
(PID.TID 0000.0001)
(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                =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_eta_max              =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_eta_min              =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_eta_mean             =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_eta_sd               =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_eta_del2             =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_uvel_max             =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_uvel_min             =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_uvel_mean            =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_uvel_sd              =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_uvel_del2            =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_vvel_max             =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_vvel_min             =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_vvel_mean            =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_vvel_sd              =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_vvel_del2            =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_wvel_max             =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_wvel_min             =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_wvel_mean            =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_wvel_sd              =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_wvel_del2            =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_theta_max            =   
5.3500000000000E+00
(PID.TID 0000.0001) %MON dynstat_theta_min            =   
4.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_theta_mean           =   
4.7300718480863E+00
(PID.TID 0000.0001) %MON dynstat_theta_sd             =   
4.2475300502820E-01
(PID.TID 0000.0001) %MON dynstat_theta_del2           =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON dynstat_salt_max             =   
1.0000000000000E+01
(PID.TID 0000.0001) %MON dynstat_salt_min             =   
1.0000000000000E+01
(PID.TID 0000.0001) %MON dynstat_salt_mean            =   
1.0000000000000E+01
(PID.TID 0000.0001) %MON dynstat_salt_sd              =   
1.4743761766996E-13
(PID.TID 0000.0001) %MON dynstat_salt_del2            =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON advcfl_uvel_max              =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON advcfl_vvel_max              =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON advcfl_wvel_max              =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON advcfl_W_hf_max              =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON pe_b_mean                    =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON ke_max                       =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON ke_mean                      =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON ke_vol                       =   
4.4227997438074E+14
(PID.TID 0000.0001) %MON vort_r_min                   =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON vort_r_max                   =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON vort_a_mean                  =   
1.0000000000000E-04
(PID.TID 0000.0001) %MON vort_a_sd                    =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON vort_p_mean                  =   
1.0429309436991E-04
(PID.TID 0000.0001) %MON vort_p_sd                    =   
4.8356310181416E-05
(PID.TID 0000.0001) %MON surfExpan_theta_mean         =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON surfExpan_salt_mean          =   
0.0000000000000E+00
(PID.TID 0000.0001) // 
=======================================================
(PID.TID 0000.0001) // End MONITOR dynamic field statistics
(PID.TID 0000.0001) // 
=======================================================
(PID.TID 0000.0001)                    cg2d_init_res 
=                     NaN
(PID.TID 0000.0001)                       cg2d_iters =  1000
(PID.TID 0000.0001)                         cg2d_res 
=                     NaN
(PID.TID 0000.0001) // 
=======================================================
(PID.TID 0000.0001) // Begin MONITOR dynamic field statistics
(PID.TID 0000.0001) // 
=======================================================
(PID.TID 0000.0001) %MON time_tsnumber                =                  
1152
(PID.TID 0000.0001) %MON time_secondsf                =   
1.7280000000000E+06
(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
(PID.TID 0000.0001) %MON dynstat_uvel_min             
=                   NaN
(PID.TID 0000.0001) %MON dynstat_uvel_mean            
=                   NaN
(PID.TID 0000.0001) %MON dynstat_uvel_sd              
=                   NaN
(PID.TID 0000.0001) %MON dynstat_uvel_del2            
=                   NaN
(PID.TID 0000.0001) %MON dynstat_vvel_max             
=                   NaN
(PID.TID 0000.0001) %MON dynstat_vvel_min             
=                   NaN
(PID.TID 0000.0001) %MON dynstat_vvel_mean            
=                   NaN
(PID.TID 0000.0001) %MON dynstat_vvel_sd              
=                   NaN
(PID.TID 0000.0001) %MON dynstat_vvel_del2            
=                   NaN
(PID.TID 0000.0001) %MON dynstat_wvel_max             
=                   NaN
(PID.TID 0000.0001) %MON dynstat_wvel_min             
=                   NaN
(PID.TID 0000.0001) %MON dynstat_wvel_mean            
=                   NaN
(PID.TID 0000.0001) %MON dynstat_wvel_sd              
=                   NaN
(PID.TID 0000.0001) %MON dynstat_wvel_del2            
=                   NaN
(PID.TID 0000.0001) %MON dynstat_theta_max            
=                   NaN
(PID.TID 0000.0001) %MON dynstat_theta_min            
=                   NaN
(PID.TID 0000.0001) %MON dynstat_theta_mean           
=                   NaN
(PID.TID 0000.0001) %MON dynstat_theta_sd             
=                   NaN
(PID.TID 0000.0001) %MON dynstat_theta_del2           
=                   NaN
(PID.TID 0000.0001) %MON dynstat_salt_max             =   
1.0000000000000E+01
(PID.TID 0000.0001) %MON dynstat_salt_min             =   
1.0000000000000E+01
(PID.TID 0000.0001) %MON dynstat_salt_mean            =   
1.0000000000000E+01
(PID.TID 0000.0001) %MON dynstat_salt_sd              =   
1.4743761766996E-13
(PID.TID 0000.0001) %MON dynstat_salt_del2            =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON advcfl_uvel_max              =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON advcfl_vvel_max              =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON advcfl_wvel_max              =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON advcfl_W_hf_max              =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON pe_b_mean                    
=                   NaN
(PID.TID 0000.0001) %MON ke_max                       =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON ke_mean                      
=                   NaN
(PID.TID 0000.0001) %MON ke_vol                       =   
4.4227997438074E+14
(PID.TID 0000.0001) %MON vort_r_min                   =   
1.0000000000000E+20
(PID.TID 0000.0001) %MON vort_r_max                   =  
-1.0000000000000E+20
(PID.TID 0000.0001) %MON vort_a_mean                  
=                   NaN
(PID.TID 0000.0001) %MON vort_a_sd                    
=                   NaN
(PID.TID 0000.0001) %MON vort_p_mean                  
=                   NaN
(PID.TID 0000.0001) %MON vort_p_sd                    =   
0.0000000000000E+00
(PID.TID 0000.0001) %MON surfExpan_theta_mean         
=                   NaN
(PID.TID 0000.0001) %MON surfExpan_salt_mean          
=                   NaN
(PID.TID 0000.0001) // 
=======================================================
(PID.TID 0000.0001) // End MONITOR dynamic field statistics
(PID.TID 0000.0001) // 
=======================================================


why does cg2d_init_res = NaN and cg2d_res = NaN ? after this time step 
these two variables remain equal to NaN and so do most of the "dynstat" 
variables, which may explain why I only have NaN in the output files.

ANY help is more than welcome !
and please let me know if you want me to send other scripts regarding 
the simulation I am running.

thank you in advance for your help,
Julie




More information about the MITgcm-support mailing list