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

Julie Deshayes jdeshayes at whoi.edu
Thu Mar 1 10:39:23 EST 2007


Dear MITgcm users,
I discovered my mistake : OLx,OLy=1. now it works perfectly fine with 3 !
than you for your help,
Julie

Matthew Mazloff a écrit :
> 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
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>




More information about the MITgcm-support mailing list