[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