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

Matthew Mazloff mmazloff at MIT.EDU
Tue Feb 27 17:27:16 EST 2007


Hi Julie,

I misspoke...it is unstable and blowing up somewhere.  Might be your  
cfl number...need to know more about your setup.

-Matt




On Feb 27, 2007, at 5:11 PM, Matthew Mazloff wrote:

> Hi Julie,
>
> The model is getting bad initial conditions somewhere/somehow.   
> Maybe you could send me your data file or explain your setup better.
>
> -Matt
>
>
> On Feb 27, 2007, at 2:48 PM, Julie Deshayes wrote:
>
>> 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
>>
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support




More information about the MITgcm-support mailing list