[MITgcm-support] NaN output
liujinliang.627
liujinliang.627 at gmail.com
Sun Feb 12 20:29:39 EST 2012
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120213/b6b7aaf6/attachment-0001.htm>
More information about the MITgcm-support
mailing list