[MITgcm-support] cg2d_init_res=NaN, cg2d_res=NaN
Matthew Mazloff
mmazloff at MIT.EDU
Tue Feb 27 17:11:19 EST 2007
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
More information about the MITgcm-support
mailing list