[MITgcm-support] NaN output
Jean-Michel Campin
jmc at ocean.mit.edu
Sun Feb 12 22:11:48 EST 2012
Hi Liu,
I would suggest 2 things:
1) to check your input files, specially the heat flux,
since a global mean heat flux (forcing_qnet_mean) out of the ocean
of -4800 W/m2 seem unrealistic to me. The min & max (forcing_qnet_min
and forcing_qnet_max), respectively -18 kW/m2 and 32 kW/m2 are also suspicious.
2) in parameter file "data", to comment out most of the different time step:
> deltaTmom = 3600.,
> deltaTtracer= 3600.,
> deltaTClock = 3600.,
> deltaTfreesurf= 86400.,
and to keep only
> deltaT = 3600.,
Cheers,
Jean-Michel
On Mon, Feb 13, 2012 at 09:29:39AM +0800, liujinliang.627 wrote:
> Hi all,
> I'm trying to have a 10 years run for climatology with the example "\verification\tutorial_global_oce_latlon"
> (Global Ocean Simulation, 1 degree * 1 degree, Nx=360, Ny=160,Nr=15 ),
>
> but all the output data are NaN, except for the first timestep which conforms to the initial conditions.
>
> The model runs without error (STDERR.0000 file is empty) and here are a few lines of the output file(STDOUT.0000 file) :
>
> (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 = 3.0168800354004E+01
> (PID.TID 0000.0001) %MON dynstat_theta_min = -1.9000000000000E+00
> (PID.TID 0000.0001) %MON dynstat_theta_mean = 1.0071849421279E+01
> (PID.TID 0000.0001) %MON dynstat_theta_sd = 8.4344918791565E+00
> (PID.TID 0000.0001) %MON dynstat_theta_del2 = 7.3227199086626E-03
> (PID.TID 0000.0001) %MON dynstat_salt_max = 3.9133998870850E+01
> (PID.TID 0000.0001) %MON dynstat_salt_min = 5.5317001342773E+00
> (PID.TID 0000.0001) %MON dynstat_salt_mean = 3.4885608051622E+01
> (PID.TID 0000.0001) %MON dynstat_salt_sd = 6.2946557329227E-01
> (PID.TID 0000.0001) %MON dynstat_salt_del2 = 8.6282517248615E-04
> (PID.TID 0000.0001) %MON forcing_qnet_max = 3.2767000000000E+04
> (PID.TID 0000.0001) %MON forcing_qnet_min = -1.8822000000000E+04
> (PID.TID 0000.0001) %MON forcing_qnet_mean = -4.8783877681316E+03
> (PID.TID 0000.0001) %MON forcing_qnet_sd = 1.7581518323022E+04
> (PID.TID 0000.0001) %MON forcing_qnet_del2 = 7.0948941598268E+01
> (PID.TID 0000.0001) %MON forcing_qsw_max = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_qsw_min = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_qsw_mean = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_qsw_sd = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_qsw_del2 = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_empmr_max = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_empmr_min = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_empmr_mean = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_empmr_sd = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_empmr_del2 = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_fu_max = 2.7340000867844E-01
> (PID.TID 0000.0001) %MON forcing_fu_min = -2.0585998892784E-01
> (PID.TID 0000.0001) %MON forcing_fu_mean = 9.6384596145960E-03
> (PID.TID 0000.0001) %MON forcing_fu_sd = 3.6706517880918E-02
> (PID.TID 0000.0001) %MON forcing_fu_del2 = 4.1447118221703E-05
> (PID.TID 0000.0001) %MON forcing_fv_max = 1.3098999857903E-01
> (PID.TID 0000.0001) %MON forcing_fv_min = -1.2700000405312E-01
> (PID.TID 0000.0001) %MON forcing_fv_mean = -1.3944381280516E-02
> (PID.TID 0000.0001) %MON forcing_fv_sd = 2.2271252100874E-02
> (PID.TID 0000.0001) %MON forcing_fv_del2 = 3.6574661145645E-05
> (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 = 1.3146436338751E+18
> (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 = -2.1798011507750E-05
> (PID.TID 0000.0001) %MON vort_a_sd = 7.3751593981362E-05
> (PID.TID 0000.0001) %MON vort_p_mean = -2.4771337589814E-05
> (PID.TID 0000.0001) %MON vort_p_sd = 1.1397521446282E-04
> (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) // =======================================================
> (PID.TID 0000.0001) // Begin MONITOR ptracer field statistics
> (PID.TID 0000.0001) // =======================================================
> (PID.TID 0000.0001) %MON trcstat_ptracer01_max = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON trcstat_ptracer01_min = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON trcstat_ptracer01_mean = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON trcstat_ptracer01_sd = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON trcstat_ptracer01_del2 = 0.0000000000000E+00
> (PID.TID 0000.0001) // =======================================================
> (PID.TID 0000.0001) // End MONITOR ptracers field statistics
> (PID.TID 0000.0001) // =======================================================
> S/R EXTERNAL_FIELDS_LOAD: Reading new data: 12 1 0 0.000000000000E+00
> (PID.TID 0000.0001) MDS_READ_FIELD: opening global file: trenberth_taux.bin
> (PID.TID 0000.0001) MDS_READ_FIELD: opening global file: trenberth_taux.bin
> (PID.TID 0000.0001) MDS_READ_FIELD: opening global file: trenberth_tauy.bin
> (PID.TID 0000.0001) MDS_READ_FIELD: opening global file: trenberth_tauy.bin
> (PID.TID 0000.0001) MDS_READ_FIELD: opening global file: ncep_qnet.bin
> (PID.TID 0000.0001) MDS_READ_FIELD: opening global file: ncep_qnet.bin
> (PID.TID 0000.0001) MDS_READ_FIELD: opening global file: lev_sst.bin
> (PID.TID 0000.0001) MDS_READ_FIELD: opening global file: lev_sst.bin
> (PID.TID 0000.0001) MDS_READ_FIELD: opening global file: lev_sss.bin
> (PID.TID 0000.0001) MDS_READ_FIELD: opening global file: lev_sss.bin
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 0.0000E+00 1.2969E+01 3.4598E+01 1.4760E-02 -1.0825E-02 -1.4859E+04 0.0000E+00 12 1 5.0000E-01 5.0000E-01
> time,fu0,fu1,fu = 0.0000E+00 1.5620E-02 1.3900E-02 1.4760E-02 5.000000000000000E-01 5.000000000000000E-01
> (PID.TID 0000.0001) SOLVE_FOR_PRESSURE: putPmEinXvector = F
> cg2d: Sum(rhs),rhsMax = -8.88178419700125E-16 7.39588767266614E+01
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 3.6000E+03 1.2972E+01 3.4598E+01 1.4758E-02 -1.0828E-02 -1.4858E+04 0.0000E+00 12 1 5.0139E-01 4.9861E-01
> time,fu0,fu1,fu = 3.6000E+03 1.5620E-02 1.3900E-02 1.4758E-02 5.013888888888889E-01 4.986111111111111E-01
> cg2d: Sum(rhs),rhsMax = -9.76996261670138E-15 7.01807415543902E+01
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 7.2000E+03 1.2974E+01 3.4598E+01 1.4755E-02 -1.0831E-02 -1.4858E+04 0.0000E+00 12 1 5.0278E-01 4.9722E-01
> time,fu0,fu1,fu = 7.2000E+03 1.5620E-02 1.3900E-02 1.4755E-02 5.027777777777778E-01 4.972222222222222E-01
> cg2d: Sum(rhs),rhsMax = 8.07687250414801E-15 6.27097245563579E+01
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 1.0800E+04 1.2977E+01 3.4598E+01 1.4753E-02 -1.0834E-02 -1.4857E+04 0.0000E+00 12 1 5.0417E-01 4.9583E-01
> time,fu0,fu1,fu = 1.0800E+04 1.5620E-02 1.3900E-02 1.4753E-02 5.041666666666667E-01 4.958333333333333E-01
> cg2d: Sum(rhs),rhsMax = -5.41233724504764E-16 1.60822321321171E+03
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 1.4400E+04 1.2979E+01 3.4598E+01 1.4750E-02 -1.0836E-02 -1.4856E+04 0.0000E+00 12 1 5.0556E-01 4.9444E-01
> time,fu0,fu1,fu = 1.4400E+04 1.5620E-02 1.3900E-02 1.4750E-02 5.055555555555555E-01 4.944444444444445E-01
> cg2d: Sum(rhs),rhsMax = -2.87530416143156E-16 4.62166536906764E+07
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 1.8000E+04 1.2982E+01 3.4598E+01 1.4748E-02 -1.0839E-02 -1.4855E+04 0.0000E+00 12 1 5.0694E-01 4.9306E-01
> time,fu0,fu1,fu = 1.8000E+04 1.5620E-02 1.3900E-02 1.4748E-02 5.069444444444444E-01 4.930555555555556E-01
> cg2d: Sum(rhs),rhsMax = 3.52270594556211E-16 7.46011535631751E+26
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 2.1600E+04 1.2984E+01 3.4598E+01 1.4746E-02 -1.0842E-02 -1.4855E+04 0.0000E+00 12 1 5.0833E-01 4.9167E-01
> time,fu0,fu1,fu = 2.1600E+04 1.5620E-02 1.3900E-02 1.4746E-02 5.083333333333333E-01 4.916666666666667E-01
> cg2d: Sum(rhs),rhsMax = -1.40844458627650E-17 5.10453559128791E+58
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 2.5200E+04 1.2987E+01 3.4597E+01 1.4743E-02 -1.0845E-02 -1.4854E+04 0.0000E+00 12 1 5.0972E-01 4.9028E-01
> time,fu0,fu1,fu = 2.5200E+04 1.5620E-02 1.3900E-02 1.4743E-02 5.097222222222222E-01 4.902777777777778E-01
> cg2d: Sum(rhs),rhsMax = 3.53820825968914E-16 1.77043873723336+187
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 2.8800E+04 1.2989E+01 3.4597E+01 1.4741E-02 -1.0848E-02 -1.4853E+04 0.0000E+00 12 1 5.1111E-01 4.8889E-01
> time,fu0,fu1,fu = 2.8800E+04 1.5620E-02 1.3900E-02 1.4741E-02 5.111111111111111E-01 4.888888888888889E-01
> cg2d: Sum(rhs),rhsMax = NaN 5.27053278474446+265
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 3.2400E+04 1.2992E+01 3.4597E+01 1.4739E-02 -1.0851E-02 -1.4853E+04 0.0000E+00 12 1 5.1250E-01 4.8750E-01
> time,fu0,fu1,fu = 3.2400E+04 1.5620E-02 1.3900E-02 1.4739E-02 5.125000000000000E-01 4.875000000000000E-01
> cg2d: Sum(rhs),rhsMax = NaN 0.00000000000000E+00
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 3.6000E+04 1.2994E+01 3.4597E+01 1.4736E-02 -1.0853E-02 -1.4852E+04 0.0000E+00 12 1 5.1389E-01 4.8611E-01
> time,fu0,fu1,fu = 3.6000E+04 1.5620E-02 1.3900E-02 1.4736E-02 5.138888888888888E-01 4.861111111111112E-01
> cg2d: Sum(rhs),rhsMax = NaN 0.00000000000000E+00
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 3.9600E+04 1.2997E+01 3.4597E+01 1.4734E-02 -1.0856E-02 -1.4851E+04 0.0000E+00 12 1 5.1528E-01 4.8472E-01
> time,fu0,fu1,fu = 3.9600E+04 1.5620E-02 1.3900E-02 1.4734E-02 5.152777777777777E-01 4.847222222222223E-01
> cg2d: Sum(rhs),rhsMax = NaN 0.00000000000000E+00
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 4.3200E+04 1.2999E+01 3.4597E+01 1.4731E-02 -1.0859E-02 -1.4851E+04 0.0000E+00 12 1 5.1667E-01 4.8333E-01
> time,fu0,fu1,fu = 4.3200E+04 1.5620E-02 1.3900E-02 1.4731E-02 5.166666666666667E-01 4.833333333333333E-01
> cg2d: Sum(rhs),rhsMax = NaN 0.00000000000000E+00
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 4.6800E+04 1.3002E+01 3.4597E+01 1.4729E-02 -1.0862E-02 -1.4850E+04 0.0000E+00 12 1 5.1806E-01 4.8194E-01
> time,fu0,fu1,fu = 4.6800E+04 1.5620E-02 1.3900E-02 1.4729E-02 5.180555555555556E-01 4.819444444444444E-01
> cg2d: Sum(rhs),rhsMax = NaN 0.00000000000000E+00
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 5.0400E+04 1.3005E+01 3.4597E+01 1.4727E-02 -1.0865E-02 -1.4849E+04 0.0000E+00 12 1 5.1944E-01 4.8056E-01
> time,fu0,fu1,fu = 5.0400E+04 1.5620E-02 1.3900E-02 1.4727E-02 5.194444444444445E-01 4.805555555555555E-01
> cg2d: Sum(rhs),rhsMax = NaN 0.00000000000000E+00
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 5.4000E+04 1.3007E+01 3.4597E+01 1.4724E-02 -1.0868E-02 -1.4848E+04 0.0000E+00 12 1 5.2083E-01 4.7917E-01
> time,fu0,fu1,fu = 5.4000E+04 1.5620E-02 1.3900E-02 1.4724E-02 5.208333333333334E-01 4.791666666666666E-01
> cg2d: Sum(rhs),rhsMax = NaN 0.00000000000000E+00
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 5.7600E+04 1.3010E+01 3.4597E+01 1.4722E-02 -1.0871E-02 -1.4848E+04 0.0000E+00 12 1 5.2222E-01 4.7778E-01
> time,fu0,fu1,fu = 5.7600E+04 1.5620E-02 1.3900E-02 1.4722E-02 5.222222222222223E-01 4.777777777777777E-01
> cg2d: Sum(rhs),rhsMax = NaN 0.00000000000000E+00
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 6.1200E+04 1.3012E+01 3.4597E+01 1.4719E-02 -1.0873E-02 -1.4847E+04 0.0000E+00 12 1 5.2361E-01 4.7639E-01
> time,fu0,fu1,fu = 6.1200E+04 1.5620E-02 1.3900E-02 1.4719E-02 5.236111111111111E-01 4.763888888888889E-01
> cg2d: Sum(rhs),rhsMax = NaN 0.00000000000000E+00
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 6.4800E+04 1.3015E+01 3.4597E+01 1.4717E-02 -1.0876E-02 -1.4846E+04 0.0000E+00 12 1 5.2500E-01 4.7500E-01
> time,fu0,fu1,fu = 6.4800E+04 1.5620E-02 1.3900E-02 1.4717E-02 5.250000000000000E-01 4.750000000000000E-01
> cg2d: Sum(rhs),rhsMax = NaN 0.00000000000000E+00
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 6.8400E+04 1.3017E+01 3.4597E+01 1.4715E-02 -1.0879E-02 -1.4846E+04 0.0000E+00 12 1 5.2639E-01 4.7361E-01
> time,fu0,fu1,fu = 6.8400E+04 1.5620E-02 1.3900E-02 1.4715E-02 5.263888888888889E-01 4.736111111111111E-01
> cg2d: Sum(rhs),rhsMax = NaN 0.00000000000000E+00
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 7.2000E+04 1.3020E+01 3.4596E+01 1.4712E-02 -1.0882E-02 -1.4845E+04 0.0000E+00 12 1 5.2778E-01 4.7222E-01
> time,fu0,fu1,fu = 7.2000E+04 1.5620E-02 1.3900E-02 1.4712E-02 5.277777777777778E-01 4.722222222222222E-01
> cg2d: Sum(rhs),rhsMax = NaN 0.00000000000000E+00
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 7.5600E+04 1.3022E+01 3.4596E+01 1.4710E-02 -1.0885E-02 -1.4844E+04 0.0000E+00 12 1 5.2917E-01 4.7083E-01
> time,fu0,fu1,fu = 7.5600E+04 1.5620E-02 1.3900E-02 1.4710E-02 5.291666666666667E-01 4.708333333333333E-01
> cg2d: Sum(rhs),rhsMax = NaN 0.00000000000000E+00
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 7.9200E+04 1.3025E+01 3.4596E+01 1.4707E-02 -1.0888E-02 -1.4844E+04 0.0000E+00 12 1 5.3056E-01 4.6944E-01
> time,fu0,fu1,fu = 7.9200E+04 1.5620E-02 1.3900E-02 1.4707E-02 5.305555555555556E-01 4.694444444444444E-01
> cg2d: Sum(rhs),rhsMax = NaN 0.00000000000000E+00
> time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 8.2800E+04 1.3027E+01 3.4596E+01 1.4705E-02 -1.0890E-02 -1.4843E+04 0.0000E+00 12 1 5.3194E-01 4.6806E-01
> time,fu0,fu1,fu = 8.2800E+04 1.5620E-02 1.3900E-02 1.4705E-02 5.319444444444444E-01 4.680555555555556E-01
> cg2d: Sum(rhs),rhsMax = NaN 0.00000000000000E+00
> (PID.TID 0000.0001) cg2d_init_res = NaN
> (PID.TID 0000.0001) cg2d_iters = 500
> (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 = 24
> (PID.TID 0000.0001) %MON time_secondsf = 8.6400000000000E+04
> (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 = NaN
> (PID.TID 0000.0001) %MON dynstat_salt_min = NaN
> (PID.TID 0000.0001) %MON dynstat_salt_mean = NaN
> (PID.TID 0000.0001) %MON dynstat_salt_sd = NaN
> (PID.TID 0000.0001) %MON dynstat_salt_del2 = NaN
> (PID.TID 0000.0001) %MON forcing_qnet_max = 3.2767000000000E+04
> (PID.TID 0000.0001) %MON forcing_qnet_min = -1.9266184722222E+04
> (PID.TID 0000.0001) %MON forcing_qnet_mean = -4.9552740240484E+03
> (PID.TID 0000.0001) %MON forcing_qnet_sd = 1.7515051018205E+04
> (PID.TID 0000.0001) %MON forcing_qnet_del2 = 7.1535887255918E+01
> (PID.TID 0000.0001) %MON forcing_qsw_max = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_qsw_min = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_qsw_mean = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_qsw_sd = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_qsw_del2 = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_empmr_max = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_empmr_min = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_empmr_mean = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_empmr_sd = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_empmr_del2 = 0.0000000000000E+00
> (PID.TID 0000.0001) %MON forcing_fu_max = 2.6788163023690E-01
> (PID.TID 0000.0001) %MON forcing_fu_min = -1.9432710289127E-01
> (PID.TID 0000.0001) %MON forcing_fu_mean = 1.1124061701982E-02
> (PID.TID 0000.0001) %MON forcing_fu_sd = 3.7653689848615E-02
> (PID.TID 0000.0001) %MON forcing_fu_del2 = 3.6409954209718E-05
> (PID.TID 0000.0001) %MON forcing_fv_max = 1.2617838831825E-01
> (PID.TID 0000.0001) %MON forcing_fv_min = -1.1307534921086E-01
> (PID.TID 0000.0001) %MON forcing_fv_mean = -1.3345185202209E-02
> (PID.TID 0000.0001) %MON forcing_fv_sd = 2.1435827725580E-02
> (PID.TID 0000.0001) %MON forcing_fv_del2 = 3.2260748992731E-05
> (PID.TID 0000.0001) %MON advcfl_uvel_max = NaN
> (PID.TID 0000.0001) %MON advcfl_vvel_max = NaN
> (PID.TID 0000.0001) %MON advcfl_wvel_max = NaN
> (PID.TID 0000.0001) %MON advcfl_W_hf_max = NaN
> (PID.TID 0000.0001) %MON pe_b_mean = NaN
> (PID.TID 0000.0001) %MON ke_max = NaN
> (PID.TID 0000.0001) %MON ke_mean = NaN
> (PID.TID 0000.0001) %MON ke_vol = 1.3146436338751E+18
> (PID.TID 0000.0001) %MON vort_r_min = NaN
> (PID.TID 0000.0001) %MON vort_r_max = NaN
> (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) // =======================================================
> What's the meaning of variables begin with "dynstat_", why did they become NaN from the second step ?
>
> Here is my data file:
> # ====================
> # | Model parameters |
> # ====================
> #
> # Continuous equation parameters
> &PARM01
> tRef = 16.0, 15.2, 14.5, 13.9, 13.3,
> 12.4, 11.3, 9.9, 8.4, 6.7,
> 5.2, 3.8, 2.9, 2.3, 1.8,
> sRef = 34.65, 34.75, 34.82, 34.87, 34.90,
> 34.90, 34.86, 34.78, 34.69, 34.60,
> 34.58, 34.62, 34.68, 34.72, 34.73,
> viscAr=1.E-3,
> viscAh=5.E5,
> diffKhT=0,
> diffKrT=3.E-5,
> diffKhS=0,
> diffKrS=3.E-5,
> rhonil=1035.,
> rhoConstFresh=1000.,
> eosType = 'JMD95Z',
> ivdc_kappa=100.,
> implicitDiffusion=.TRUE.,
> allowFreezing=.TRUE.,
> exactConserv=.TRUE.,
> useRealFreshWaterFlux=.FALSE.,
> useCDscheme=.TRUE.,
> useSingleCpuIO=.TRUE.,
> # turn on looped cells
> hFacMin=.05,
> hFacMindr=50.,
> # set precision of data files
> readBinaryPrec=32,
> &
>
> # Elliptic solver parameters
> &PARM02
> cg2dMaxIters=500,
> cg2dTargetResidual=1.E-13,
> &
>
> # Time stepping parameters
> &PARM03
> # nIter0= 0,
> # nTimeSteps = 8640,
> # 100 years of integration will yield a reasonable flow field
> startTime = 0.,
> endTime = 31104000.,
> deltaT = 3600.,
> deltaTmom = 3600.,
> tauCD = 321428.,
> deltaTtracer= 3600.,
> deltaTClock = 3600.,
> deltaTfreesurf= 86400.,
> abEps = 0.1,
> pChkptFreq= 2592000.,
> dumpFreq= 2592000.,
> taveFreq= 2592000.,
> monitorFreq=86400.,
>
> tauThetaClimRelax= 0.,
> tauSaltClimRelax = 0.,
>
> periodicExternalForcing=.TRUE.,
> externForcingPeriod=2592000.,
> externForcingCycle=31104000.,
> &
>
> # Gridding parameters
> &PARM04
> usingSphericalPolarGrid=.TRUE.,
> delR= 50., 70., 100., 140., 190.,
> 240., 290., 340., 390., 440.,
> 490., 540., 590., 640., 690.,
> ygOrigin=-79.5,
> xgOrigin=0.5,
> dySpacing=1.,
> dxSpacing=1.,
> &
>
> # Input datasets
> &PARM05
>
>
> bathyFile= 'bathymetry.bin',
> hydrogThetaFile='lev_t.bin',
> hydrogSaltFile= 'lev_s.bin',
> zonalWindFile= 'trenberth_taux.bin',
> meridWindFile= 'trenberth_tauy.bin',
> thetaClimFile= 'lev_sst.bin',
> saltClimFile= 'lev_sss.bin',
> surfQFile= 'ncep_qnet.bin',
> the_run_name= 'global_oce_latlon',
> # fresh water flux is turned on, comment next line to it turn off
> # (maybe better with surface salinity restoring)
> # EmPmRFile= 'ncep_emp.bin',
> &
>
> Any help is more than welcome !
>
> And please let me know if you want me to send other scripts regarding.
>
> Thank you!
>
> Liu Jinliang
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list