Index '282' of dimension 2 of array 'finp' above upper bound of 281
- zhangtianhao
- Posts: 50
- Joined: Fri May 11, 2018 5:36 pm
- Location: Beijing Normal University
Index '282' of dimension 2 of array 'finp' above upper bound of 281
Hi there,
I met some problems when I ran a simple case:Index '282' of dimension 2 of array 'finp' above upper bound of 281. I checked interpolate.f90 and considered ‘finp’ may be related to grid file. Then I check my grid file, but the dimensions of the whole varibales are 281, no 282 appeared(my dimension is 179*281).
I can't figure out what happened to my grid file or other input files.
Maybe it's stupid to you but thanks.
Here is the result and grid file:
NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(Grid 01, File: palu_init.nc, Rec=0001, Index=1)
- free-surface
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- potential temperature
(Min = 1.00000000E+01 Max = 1.00000000E+01)
- salinity
(Min = 3.50000000E+01 Max = 3.50000000E+01)
At line 117 of file interpolate.f90
Fortran runtime error: Index '282' of dimension 2 of array 'finp' above upper bound of 281
Error termination. Backtrace:
#0 0x7fd5338322da in ???
#1 0x7fd533832ec5 in ???
#2 0x7fd533833297 in ???
#3 0x557ab0b8c343 in __interpolate_mod_MOD_linterp2d
at /home/ubuntu/app/COAWST_V3.2/Build/interpolate.f90:117
#4 0x557ab0b26fb6 in regrid_
at /home/ubuntu/app/COAWST_V3.2/Build/regrid.f90:180
#5 0x557ab0b1726b in __nf_fread2d_mod_MOD_nf_fread2d
at /home/ubuntu/app/COAWST_V3.2/Build/nf_fread2d.f90:310
#6 0x557ab0925c32 in get_2dfld_
at /home/ubuntu/app/COAWST_V3.2/Build/get_2dfld.f90:198
#7 0x557ab0533c01 in get_data_
at /home/ubuntu/app/COAWST_V3.2/Build/get_data.f90:104
#8 0x557ab0530791 in initial_
at /home/ubuntu/app/COAWST_V3.2/Build/initial.f90:296
#9 0x557ab034fbbd in __ocean_control_mod_MOD_roms_initialize
at /home/ubuntu/app/COAWST_V3.2/Build/ocean_control.f90:131
#10 0x557ab034ec15 in ocean
at /home/ubuntu/app/COAWST_V3.2/Build/master.f90:95
#11 0x557ab034ede1 in main
at /home/ubuntu/app/COAWST_V3.2/Build/master.f90:50
-------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code.. Per user-direction, the job has been aborted.
-------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:
Process name: [[4580,1],0]
Exit code: 2
I met some problems when I ran a simple case:Index '282' of dimension 2 of array 'finp' above upper bound of 281. I checked interpolate.f90 and considered ‘finp’ may be related to grid file. Then I check my grid file, but the dimensions of the whole varibales are 281, no 282 appeared(my dimension is 179*281).
I can't figure out what happened to my grid file or other input files.
Maybe it's stupid to you but thanks.
Here is the result and grid file:
NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(Grid 01, File: palu_init.nc, Rec=0001, Index=1)
- free-surface
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- potential temperature
(Min = 1.00000000E+01 Max = 1.00000000E+01)
- salinity
(Min = 3.50000000E+01 Max = 3.50000000E+01)
At line 117 of file interpolate.f90
Fortran runtime error: Index '282' of dimension 2 of array 'finp' above upper bound of 281
Error termination. Backtrace:
#0 0x7fd5338322da in ???
#1 0x7fd533832ec5 in ???
#2 0x7fd533833297 in ???
#3 0x557ab0b8c343 in __interpolate_mod_MOD_linterp2d
at /home/ubuntu/app/COAWST_V3.2/Build/interpolate.f90:117
#4 0x557ab0b26fb6 in regrid_
at /home/ubuntu/app/COAWST_V3.2/Build/regrid.f90:180
#5 0x557ab0b1726b in __nf_fread2d_mod_MOD_nf_fread2d
at /home/ubuntu/app/COAWST_V3.2/Build/nf_fread2d.f90:310
#6 0x557ab0925c32 in get_2dfld_
at /home/ubuntu/app/COAWST_V3.2/Build/get_2dfld.f90:198
#7 0x557ab0533c01 in get_data_
at /home/ubuntu/app/COAWST_V3.2/Build/get_data.f90:104
#8 0x557ab0530791 in initial_
at /home/ubuntu/app/COAWST_V3.2/Build/initial.f90:296
#9 0x557ab034fbbd in __ocean_control_mod_MOD_roms_initialize
at /home/ubuntu/app/COAWST_V3.2/Build/ocean_control.f90:131
#10 0x557ab034ec15 in ocean
at /home/ubuntu/app/COAWST_V3.2/Build/master.f90:95
#11 0x557ab034ede1 in main
at /home/ubuntu/app/COAWST_V3.2/Build/master.f90:50
-------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code.. Per user-direction, the job has been aborted.
-------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:
Process name: [[4580,1],0]
Exit code: 2
- Attachments
-
- simple_grid.nc
- (7.61 MiB) Downloaded 312 times
Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281
I suspect you are reading an atmospheric field here? The model wouldn't try to interpolate from your grid file, assuming all those fields are already on your grid. So, what atmospheric fields are you reading and are you expecting ROMS to do the interpolations? What grid are they on?at /home/ubuntu/app/COAWST_V3.2/Build/get_data.f90:104
- zhangtianhao
- Posts: 50
- Joined: Fri May 11, 2018 5:36 pm
- Location: Beijing Normal University
Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281
Thank you, I input uv-stress by stressTEST.nc file. But I have also checked this nc file, and the dimension of uv-stress is 178*281*time, not 282.kate wrote: ↑Wed Jan 08, 2020 11:40 amI suspect you are reading an atmospheric field here? The model wouldn't try to interpolate from your grid file, assuming all those fields are already on your grid. So, what atmospheric fields are you reading and are you expecting ROMS to do the interpolations? What grid are they on?at /home/ubuntu/app/COAWST_V3.2/Build/get_data.f90:104
Attached is stressTEST.nc file. I turn off ana_smflx and input it by nc file containing my group's data.
I am still confused of it.
![Crying or Very sad :cry:](./images/smilies/icon_cry.gif)
- Attachments
-
- stressTEST.zip
- (842.37 KiB) Downloaded 270 times
Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281
So you think you have already remapped the thing to your grid. You have to have the dimensions *exactly* right or else ROMS will decide to remap for you.
Some words on grid sizes. If you want powers of two or three in the grid for your decompositions without needing leftovers, you need Lm, Mm to have factors of two or three - something like 280 is good. Then the rho grid with halos would have say eta_rho=282. The v-stress is on the v-points of the grid, therefore dimensioned 281. Search the wiki for Arakawa C-grid to see the staggering used - not all variables are dimensioned the same.
Some words on grid sizes. If you want powers of two or three in the grid for your decompositions without needing leftovers, you need Lm, Mm to have factors of two or three - something like 280 is good. Then the rho grid with halos would have say eta_rho=282. The v-stress is on the v-points of the grid, therefore dimensioned 281. Search the wiki for Arakawa C-grid to see the staggering used - not all variables are dimensioned the same.
- zhangtianhao
- Posts: 50
- Joined: Fri May 11, 2018 5:36 pm
- Location: Beijing Normal University
Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281
Kate,kate wrote: ↑Wed Jan 08, 2020 12:53 pm So you think you have already remapped the thing to your grid. You have to have the dimensions *exactly* right or else ROMS will decide to remap for you.
Some words on grid sizes. If you want powers of two or three in the grid for your decompositions without needing leftovers, you need Lm, Mm to have factors of two or three - something like 280 is good. Then the rho grid with halos would have say eta_rho=282. The v-stress is on the v-points of the grid, therefore dimensioned 281. Search the wiki for Arakawa C-grid to see the staggering used - not all variables are dimensioned the same.
Thanks for your suggestion.
Maybe I am not clear about the exact meaning of my input forcing file. I set uv-stress on rho grid.(rho dimension is 178*281), but it seems that uv-stress should set on u-grid (dimension: 177*281) and v-grid (dimension: 178*280), I misunderstood it.
Thanks again.
- zhangtianhao
- Posts: 50
- Joined: Fri May 11, 2018 5:36 pm
- Location: Beijing Normal University
Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281
Sorry,kate wrote: ↑Wed Jan 08, 2020 12:53 pm So you think you have already remapped the thing to your grid. You have to have the dimensions *exactly* right or else ROMS will decide to remap for you.
Some words on grid sizes. If you want powers of two or three in the grid for your decompositions without needing leftovers, you need Lm, Mm to have factors of two or three - something like 280 is good. Then the rho grid with halos would have say eta_rho=282. The v-stress is on the v-points of the grid, therefore dimensioned 281. Search the wiki for Arakawa C-grid to see the staggering used - not all variables are dimensioned the same.
I have changed dimension of uv-stress, but I met another problem:Fortran runtime error: Index '1' of dimension 2 of array 'bry' above upper bound of 0.
It seems that there are some problems with my bry, but my forcing file only contains the information of sustrs, svstrs, lon_u, lon_v, lat_u, lat_v. I don't know it's relationship with bry.
NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(Grid 01, File: init.nc, Rec=0001, Index=1)
- free-surface
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- potential temperature
(Min = 1.00000000E+01 Max = 1.00000000E+01)
- salinity
(Min = 3.50000000E+01 Max = 3.50000000E+01)
GET_2DFLD - surface u-momentum stress, t = 0 00:00:00
(Rec=0000001, Index=1, File: stressTEST.nc)
(Tmin= 0.0000 Tmax= 0.0347)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
At line 139 of file get_data.f90
Fortran runtime error: Index '1' of dimension 2 of array 'bry' above upper bound of 0
Error termination. Backtrace:
GET_2DFLD - surface v-momentum stress, t = 0 00:00:00
(Rec=0000001, Index=1, File: stressTEST.nc)
(Tmin= 0.0000 Tmax= 0.0347)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
#0 0x7f162ca952da in ???
#1 0x7f162ca95ec5 in ???
#2 0x7f162ca96297 in ???
#3 0x55d9a7b3a288 in get_data_
at /home/ubuntu/app/COAWST_V3.2/Build/get_data.f90:139
#4 0x55d9a7b35791 in initial_
at /home/ubuntu/app/COAWST_V3.2/Build/initial.f90:296
#5 0x55d9a7954bbd in __ocean_control_mod_MOD_roms_initialize
at /home/ubuntu/app/COAWST_V3.2/Build/ocean_control.f90:131
#6 0x55d9a7953c15 in ocean
at /home/ubuntu/app/COAWST_V3.2/Build/master.f90:95
#7 0x55d9a7953de1 in main
at /home/ubuntu/app/COAWST_V3.2/Build/master.f90:50
-------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code.. Per user-direction, the job has been aborted.
-------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:
Process name: [[29868,1],0]
Exit code: 2
- Attachments
-
- stressTEST.zip
- (1.03 MiB) Downloaded 333 times
Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281
Since you are using COAWST, you need to set the number of boundary and climatology files, not just the number of forcing files. See the coawst roms.in for an example. Same for the river inputs.
- zhangtianhao
- Posts: 50
- Joined: Fri May 11, 2018 5:36 pm
- Location: Beijing Normal University
Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281
Attached is my bry condition
- Attachments
-
- bry.png (30.67 KiB) Viewed 6097 times
Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281
The Flather condition requires values from outside the domain for both ubar/vbar and zeta. You can provide them via ANA_xxx options or you can specify a file:
Code: Select all
NBCFILES == 1 ! number of boundary files
BRYNAME == /center1/AKWATERS/kshedstrom/Arctic/Files/Arctic4_bdry_1980_SODA3.3.1.nc |
/center1/AKWATERS/kshedstrom/Arctic/Files/Arctic4_bdry_1981_SODA3.3.1.nc |
/center1/AKWATERS/kshedstrom/Arctic/Files/Arctic4_bdry_1982_SODA3.3.1.nc
- zhangtianhao
- Posts: 50
- Joined: Fri May 11, 2018 5:36 pm
- Location: Beijing Normal University
Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281
Thanks for your reply, I use our group's server and I only use the part of ROMS. I will go to learn the input of COAWST.kate wrote: ↑Thu Jan 09, 2020 3:04 am The Flather condition requires values from outside the domain for both ubar/vbar and zeta. You can provide them via ANA_xxx options or you can specify a file:Code: Select all
NBCFILES == 1 ! number of boundary files BRYNAME == /center1/AKWATERS/kshedstrom/Arctic/Files/Arctic4_bdry_1980_SODA3.3.1.nc | /center1/AKWATERS/kshedstrom/Arctic/Files/Arctic4_bdry_1981_SODA3.3.1.nc | /center1/AKWATERS/kshedstrom/Arctic/Files/Arctic4_bdry_1982_SODA3.3.1.nc
Thanks again
- zhangtianhao
- Posts: 50
- Joined: Fri May 11, 2018 5:36 pm
- Location: Beijing Normal University
Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281
Hi all,
Now I can run this case but it blows up after several steps.
And it's weird that: usually results become large and model blows up, but this time results suddenly become NaN.
So I checked my forcing file, but there are no NaN in it.
One thing I noticed:
I turn on WET_DRY, and I set depth of land area < 0.
I find that before model blows up, zeta in land area changed bigger(changed up tp 20m). (Why land-mask has zeta value??)
I wonder if this matters? (!!!!I CHECKED U V W BEFORE BLOW-UP, they are fine! all=zero )
Thanks!!
Basin information for Grid 01:
Maximum grid stiffness ratios: rx0 = 1.500000E-01 (Beckmann and Haidvogel)
rx1 = 5.728206E+00 (Haney)
Initial basin volumes: TotVolume = 6.3638374556E+10 m3
MinVolume = 2.5571097553E+02 m3
MaxVolume = 1.4785505052E+06 m3
Max/Min = 5.7821159305E+03
NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00029000)
GET_2DFLD - surface u-momentum stress, t = 0 02:24:00
(Rec=0000002, Index=2, File: stressTEST.nc)
(Tmin= 0.0000 Tmax= 3001.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface v-momentum stress, t = 0 02:24:00
(Rec=0000002, Index=2, File: stressTEST.nc)
(Tmin= 0.0000 Tmax= 3001.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME
C => (i,j,k) Cu Cv Cw Max Speed
0 0 00:00:00 0.000000E+00 9.300949E+03 9.300949E+03 6.367420E+10
(000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_HIS - creating history file, Grid 01: simplecase_his.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
DEF_AVG - creating average file, Grid 01: simplecase_avg.nc
1 0 00:00:00 7.403984E-35 9.301714E+03 9.301714E+03 6.366897E+10
(073,204,01) 2.838148E-20 3.182801E-20 0.000000E+00 7.471913E-17
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000002
2 0 00:00:00 2.892330E-34 9.301714E+03 9.301714E+03 6.366897E+10
(058,226,25) 6.956945E-23 4.667068E-22 1.367662E-16 1.488681E-16
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000003
3 0 00:00:00 6.180541E-34 9.301714E+03 9.301714E+03 6.366897E+10
(058,226,25) 2.184770E-22 1.522074E-21 2.705075E-16 2.451863E-16
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000004
4 0 00:00:00 1.027599E-33 9.301714E+03 9.301714E+03 6.366897E+10
(048,263,25) 3.982611E-22 6.445907E-21 3.350156E-16 4.986654E-16
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000005
5 0 00:00:00 1.477361E-33 9.301714E+03 9.301714E+03 6.366897E+10
(057,240,25) 2.241456E-21 7.977925E-21 4.357218E-16 7.637114E-16
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000006
6 0 00:00:00 1.952023E-33 9.301714E+03 9.301714E+03 6.366897E+10
(086,192,25) 3.322778E-21 3.886765E-21 4.632618E-16 1.028624E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000007
7 0 00:00:00 2.444697E-33 9.301714E+03 9.301714E+03 6.366897E+10
(086,192,25) 4.586721E-21 2.804571E-21 5.329079E-16 1.293538E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000008
8 0 00:00:00 3.007035E-33 9.301714E+03 9.301714E+03 6.366897E+10
(047,269,25) 5.364246E-21 2.059773E-20 5.409963E-16 1.558241E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000009
9 0 00:00:00 3.690863E-33 9.301714E+03 9.301714E+03 6.366897E+10
(059,225,25) 1.215610E-20 3.419549E-21 6.021904E-16 1.822926E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000010
10 0 00:00:00 4.620259E-33 9.301714E+03 9.301714E+03 6.366897E+10
(058,238,25) 3.370067E-20 1.212422E-20 8.214504E-16 2.087500E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000011
11 0 00:00:01 5.830281E-33 9.301714E+03 9.301714E+03 6.366897E+10
(052,249,25) 2.251131E-21 7.650399E-21 1.012366E-15 2.351966E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000012
12 0 00:00:01 7.741074E-33 9.301714E+03 9.301714E+03 6.366897E+10
(057,270,25) 2.112356E-21 3.158259E-20 1.175885E-15 2.616385E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000013
13 0 00:00:01 1.033973E-32 9.301714E+03 9.301714E+03 6.366897E+10
(057,270,25) 2.272178E-21 3.745048E-20 1.355493E-15 2.880627E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000014
14 0 00:00:01 1.377428E-32 9.301714E+03 9.301714E+03 6.366897E+10
(055,247,25) 2.595464E-20 5.314083E-20 1.593740E-15 3.144802E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000015
15 0 00:00:01 1.782040E-32 9.301714E+03 9.301714E+03 6.366897E+10
(055,247,25) 2.649002E-20 4.711058E-20 1.979339E-15 3.408829E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000016
16 0 00:00:01 2.255626E-32 9.301714E+03 9.301714E+03 6.366897E+10
(055,247,25) 3.381972E-20 6.569763E-20 2.318928E-15 3.672661E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000017
17 0 00:00:01 2.762860E-32 9.301714E+03 9.301714E+03 6.366897E+10
(063,217,25) 1.388260E-19 1.283977E-19 2.858346E-15 3.936374E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000018
18 0 00:00:01 3.298347E-32 9.301714E+03 9.301714E+03 6.366897E+10
(063,217,25) 1.484557E-19 1.637055E-19 3.103611E-15 4.199973E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000019
19 0 00:00:01 NaN NaN NaN 6.366897E+10
(063,217,25) 1.554516E-19 1.820026E-19 3.170670E-15 1.839072E-15
Blowing-up: Saving latest model state into RESTART file
WRT_RST - wrote re-start fields (Index=2,2) into time record = 0000001
Elapsed CPU time (seconds):
Node # 0 CPU: 21.440
Total: 21.440
Nonlinear ocean model elapsed time profile, Grid: 01
Allocation and array initialization .............. 0.119 ( 0.5533 %)
Ocean state initialization ....................... 0.015 ( 0.0717 %)
Reading of input data ............................ 0.006 ( 0.0274 %)
Processing of input data ......................... 0.006 ( 0.0264 %)
Processing of output time averaged data .......... 0.011 ( 0.0502 %)
Computation of vertical boundary conditions ...... 0.031 ( 0.1456 %)
Computation of global information integrals ...... 0.135 ( 0.6315 %)
Writing of output data ........................... 2.427 (11.3187 %)
Model 2D kernel .................................. 3.937 (18.3644 %)
2D/3D coupling, vertical metrics ................. 0.145 ( 0.6742 %)
Omega vertical velocity .......................... 0.113 ( 0.5256 %)
Equation of state for seawater ................... 0.125 ( 0.5815 %)
GLS vertical mixing parameterization ............. 9.713 (45.3013 %)
3D equations right-side terms .................... 0.700 ( 3.2662 %)
3D equations predictor step ...................... 0.729 ( 3.4019 %)
Pressure gradient ................................ 0.358 ( 1.6705 %)
Harmonic mixing of tracers, geopotentials ........ 0.812 ( 3.7854 %)
Harmonic stress tensor, S-surfaces ............... 0.304 ( 1.4182 %)
Corrector time-step for 3D momentum .............. 0.555 ( 2.5868 %)
Corrector time-step for tracers .................. 0.628 ( 2.9313 %)
Total: 20.868 97.3319
Nonlinear model message Passage profile, Grid: 01
Message Passage: 2D halo exchanges ............... 0.011 ( 0.0507 %)
Message Passage: 3D halo exchanges ............... 0.001 ( 0.0032 %)
Message Passage: 4D halo exchanges ............... 0.000 ( 0.0007 %)
Message Passage: data broadcast .................. 0.001 ( 0.0060 %)
Message Passage: data reduction .................. 0.000 ( 0.0008 %)
Message Passage: data gathering .................. 0.310 ( 1.4437 %)
Message Passage: data scattering.................. 0.020 ( 0.0942 %)
Message Passage: boundary data gathering ......... 0.000 ( 0.0006 %)
Total: 0.343 1.5999
All percentages are with respect to total time = 21.440
ROMS/TOMS - Output NetCDF summary for Grid 01:
number of time records written in HISTORY file = 00000019
number of time records written in RESTART file = 00000001
Analytical header files used:
ROMS/Functionals/ana_btflux.h
/home/ubuntu/app/COAWST_V3.2/Projects/simplecase/ana_fsobc.h
/home/ubuntu/app/COAWST_V3.2/Projects/simplecase/ana_m2obc.h
ROMS/Functionals/ana_stflux.h
ROMS/TOMS: DONE... Thursday - January 9, 2020 - 6:57:05 PM
Now I can run this case but it blows up after several steps.
And it's weird that: usually results become large and model blows up, but this time results suddenly become NaN.
So I checked my forcing file, but there are no NaN in it.
One thing I noticed:
I turn on WET_DRY, and I set depth of land area < 0.
I find that before model blows up, zeta in land area changed bigger(changed up tp 20m). (Why land-mask has zeta value??)
I wonder if this matters? (!!!!I CHECKED U V W BEFORE BLOW-UP, they are fine! all=zero )
Thanks!!
Basin information for Grid 01:
Maximum grid stiffness ratios: rx0 = 1.500000E-01 (Beckmann and Haidvogel)
rx1 = 5.728206E+00 (Haney)
Initial basin volumes: TotVolume = 6.3638374556E+10 m3
MinVolume = 2.5571097553E+02 m3
MaxVolume = 1.4785505052E+06 m3
Max/Min = 5.7821159305E+03
NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00029000)
GET_2DFLD - surface u-momentum stress, t = 0 02:24:00
(Rec=0000002, Index=2, File: stressTEST.nc)
(Tmin= 0.0000 Tmax= 3001.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface v-momentum stress, t = 0 02:24:00
(Rec=0000002, Index=2, File: stressTEST.nc)
(Tmin= 0.0000 Tmax= 3001.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME
C => (i,j,k) Cu Cv Cw Max Speed
0 0 00:00:00 0.000000E+00 9.300949E+03 9.300949E+03 6.367420E+10
(000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_HIS - creating history file, Grid 01: simplecase_his.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
DEF_AVG - creating average file, Grid 01: simplecase_avg.nc
1 0 00:00:00 7.403984E-35 9.301714E+03 9.301714E+03 6.366897E+10
(073,204,01) 2.838148E-20 3.182801E-20 0.000000E+00 7.471913E-17
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000002
2 0 00:00:00 2.892330E-34 9.301714E+03 9.301714E+03 6.366897E+10
(058,226,25) 6.956945E-23 4.667068E-22 1.367662E-16 1.488681E-16
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000003
3 0 00:00:00 6.180541E-34 9.301714E+03 9.301714E+03 6.366897E+10
(058,226,25) 2.184770E-22 1.522074E-21 2.705075E-16 2.451863E-16
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000004
4 0 00:00:00 1.027599E-33 9.301714E+03 9.301714E+03 6.366897E+10
(048,263,25) 3.982611E-22 6.445907E-21 3.350156E-16 4.986654E-16
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000005
5 0 00:00:00 1.477361E-33 9.301714E+03 9.301714E+03 6.366897E+10
(057,240,25) 2.241456E-21 7.977925E-21 4.357218E-16 7.637114E-16
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000006
6 0 00:00:00 1.952023E-33 9.301714E+03 9.301714E+03 6.366897E+10
(086,192,25) 3.322778E-21 3.886765E-21 4.632618E-16 1.028624E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000007
7 0 00:00:00 2.444697E-33 9.301714E+03 9.301714E+03 6.366897E+10
(086,192,25) 4.586721E-21 2.804571E-21 5.329079E-16 1.293538E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000008
8 0 00:00:00 3.007035E-33 9.301714E+03 9.301714E+03 6.366897E+10
(047,269,25) 5.364246E-21 2.059773E-20 5.409963E-16 1.558241E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000009
9 0 00:00:00 3.690863E-33 9.301714E+03 9.301714E+03 6.366897E+10
(059,225,25) 1.215610E-20 3.419549E-21 6.021904E-16 1.822926E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000010
10 0 00:00:00 4.620259E-33 9.301714E+03 9.301714E+03 6.366897E+10
(058,238,25) 3.370067E-20 1.212422E-20 8.214504E-16 2.087500E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000011
11 0 00:00:01 5.830281E-33 9.301714E+03 9.301714E+03 6.366897E+10
(052,249,25) 2.251131E-21 7.650399E-21 1.012366E-15 2.351966E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000012
12 0 00:00:01 7.741074E-33 9.301714E+03 9.301714E+03 6.366897E+10
(057,270,25) 2.112356E-21 3.158259E-20 1.175885E-15 2.616385E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000013
13 0 00:00:01 1.033973E-32 9.301714E+03 9.301714E+03 6.366897E+10
(057,270,25) 2.272178E-21 3.745048E-20 1.355493E-15 2.880627E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000014
14 0 00:00:01 1.377428E-32 9.301714E+03 9.301714E+03 6.366897E+10
(055,247,25) 2.595464E-20 5.314083E-20 1.593740E-15 3.144802E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000015
15 0 00:00:01 1.782040E-32 9.301714E+03 9.301714E+03 6.366897E+10
(055,247,25) 2.649002E-20 4.711058E-20 1.979339E-15 3.408829E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000016
16 0 00:00:01 2.255626E-32 9.301714E+03 9.301714E+03 6.366897E+10
(055,247,25) 3.381972E-20 6.569763E-20 2.318928E-15 3.672661E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000017
17 0 00:00:01 2.762860E-32 9.301714E+03 9.301714E+03 6.366897E+10
(063,217,25) 1.388260E-19 1.283977E-19 2.858346E-15 3.936374E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000018
18 0 00:00:01 3.298347E-32 9.301714E+03 9.301714E+03 6.366897E+10
(063,217,25) 1.484557E-19 1.637055E-19 3.103611E-15 4.199973E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000019
19 0 00:00:01 NaN NaN NaN 6.366897E+10
(063,217,25) 1.554516E-19 1.820026E-19 3.170670E-15 1.839072E-15
Blowing-up: Saving latest model state into RESTART file
WRT_RST - wrote re-start fields (Index=2,2) into time record = 0000001
Elapsed CPU time (seconds):
Node # 0 CPU: 21.440
Total: 21.440
Nonlinear ocean model elapsed time profile, Grid: 01
Allocation and array initialization .............. 0.119 ( 0.5533 %)
Ocean state initialization ....................... 0.015 ( 0.0717 %)
Reading of input data ............................ 0.006 ( 0.0274 %)
Processing of input data ......................... 0.006 ( 0.0264 %)
Processing of output time averaged data .......... 0.011 ( 0.0502 %)
Computation of vertical boundary conditions ...... 0.031 ( 0.1456 %)
Computation of global information integrals ...... 0.135 ( 0.6315 %)
Writing of output data ........................... 2.427 (11.3187 %)
Model 2D kernel .................................. 3.937 (18.3644 %)
2D/3D coupling, vertical metrics ................. 0.145 ( 0.6742 %)
Omega vertical velocity .......................... 0.113 ( 0.5256 %)
Equation of state for seawater ................... 0.125 ( 0.5815 %)
GLS vertical mixing parameterization ............. 9.713 (45.3013 %)
3D equations right-side terms .................... 0.700 ( 3.2662 %)
3D equations predictor step ...................... 0.729 ( 3.4019 %)
Pressure gradient ................................ 0.358 ( 1.6705 %)
Harmonic mixing of tracers, geopotentials ........ 0.812 ( 3.7854 %)
Harmonic stress tensor, S-surfaces ............... 0.304 ( 1.4182 %)
Corrector time-step for 3D momentum .............. 0.555 ( 2.5868 %)
Corrector time-step for tracers .................. 0.628 ( 2.9313 %)
Total: 20.868 97.3319
Nonlinear model message Passage profile, Grid: 01
Message Passage: 2D halo exchanges ............... 0.011 ( 0.0507 %)
Message Passage: 3D halo exchanges ............... 0.001 ( 0.0032 %)
Message Passage: 4D halo exchanges ............... 0.000 ( 0.0007 %)
Message Passage: data broadcast .................. 0.001 ( 0.0060 %)
Message Passage: data reduction .................. 0.000 ( 0.0008 %)
Message Passage: data gathering .................. 0.310 ( 1.4437 %)
Message Passage: data scattering.................. 0.020 ( 0.0942 %)
Message Passage: boundary data gathering ......... 0.000 ( 0.0006 %)
Total: 0.343 1.5999
All percentages are with respect to total time = 21.440
ROMS/TOMS - Output NetCDF summary for Grid 01:
number of time records written in HISTORY file = 00000019
number of time records written in RESTART file = 00000001
Analytical header files used:
ROMS/Functionals/ana_btflux.h
/home/ubuntu/app/COAWST_V3.2/Projects/simplecase/ana_fsobc.h
/home/ubuntu/app/COAWST_V3.2/Projects/simplecase/ana_m2obc.h
ROMS/Functionals/ana_stflux.h
ROMS/TOMS: DONE... Thursday - January 9, 2020 - 6:57:05 PM
Last edited by zhangtianhao on Thu Jan 09, 2020 12:19 pm, edited 1 time in total.
Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281
With WET_DRY, zeta in the land area is always some DCRIT above the land, nothing to worry about. It should not be changing in time.
Have you tried fussing with your timestep?
Have you tried fussing with your timestep?
- zhangtianhao
- Posts: 50
- Joined: Fri May 11, 2018 5:36 pm
- Location: Beijing Normal University
Re: Index '282' of dimension 2 of array 'finp' above upper bound of 281
I changed dt from 0.1 to 0.05, and it delays the blow-up time from 1.9s to 44s. I am keeping trying.
Here is the complete result.out.
And btw , my zeta in the land area seems too big, it nearly equals to height of land ~ sevral hundreds?(My forcing file cover the area of land, I doubt if it matters)
Operating system : Linux
CPU/hardware : x86_64
Compiler system : gfortran
Compiler command : /usr/bin/mpif90
Compiler flags : -frepack-arrays -O3 -ffast-math -ftree-vectorize -ftree-loop-linear -funroll-loops -w -ffree-form -ffree-line-length-none -frecord-marker=4 -fconvert=big-endian
Input Script :
SVN Root URL : https:://myroms.org/svn/src
SVN Revision : 65M
Local Root : /home/ubuntu/app/COAWST_V3.2
Header Dir : /home/ubuntu/app/COAWST_V3.2/Projects/simplecase
Header file : palu.h
Analytical Dir: /home/ubuntu/app/COAWST_V3.2/Projects/simplecase
Resolution, Grid 01: 0176x0279x025, Parallel Nodes: 1, Tiling: 001x001
Physical Parameters, Grid: 01
=============================
29000 ntimes Number of timesteps for 3-D equations.
0.100 dt Timestep size (s) for 3-D equations.
20 ndtfast Number of timesteps for 2-D equations between
each 3D timestep.
1 ERstr Starting ensemble/perturbation run number.
1 ERend Ending ensemble/perturbation run number.
0 nrrec Number of restart records to read from disk.
T LcycleRST Switch to recycle time-records in restart file.
120 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
T ldefout Switch to create a new output NetCDF file(s).
1 nHIS Number of timesteps between the writing fields
into history file.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
360 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
2.0000E-01 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
2.0000E-01 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
1.0000E-01 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
F LuvSponge Turning OFF sponge on horizontal momentum.
F LtracerSponge(01) Turning OFF sponge on tracer 01: temp
F LtracerSponge(02) Turning OFF sponge on tracer 02: salt
1.0000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
1.0000E-06 Akt_bak(02) Background vertical mixing coefficient (m2/s)
for tracer 02: salt
1.0000E-05 Akv_bak Background vertical mixing coefficient (m2/s)
for momentum.
5.0000E-06 Akk_bak Background vertical mixing coefficient (m2/s)
for turbulent energy.
5.0000E-06 Akp_bak Background vertical mixing coefficient (m2/s)
for turbulent generic statistical field.
3.000 gls_p GLS stability exponent.
1.500 gls_m GLS turbulent kinetic energy exponent.
-1.000 gls_n GLS turbulent length scale exponent.
7.6000E-06 gls_Kmin GLS minimum value of turbulent kinetic energy.
1.0000E-12 gls_Pmin GLS minimum value of dissipation.
5.4770E-01 gls_cmu0 GLS stability coefficient.
1.4400E+00 gls_c1 GLS shear production coefficient.
1.9200E+00 gls_c2 GLS dissipation coefficient.
-4.0000E-01 gls_c3m GLS stable buoyancy production coefficient.
1.0000E+00 gls_c3p GLS unstable buoyancy production coefficient.
1.0000E+00 gls_sigk GLS constant Schmidt number for TKE.
1.3000E+00 gls_sigp GLS constant Schmidt number for PSI.
1400.000 charnok_alpha Charnok factor for Zos calculation.
0.500 zos_hsig_alpha Factor for Zos calculation using Hsig(Awave).
0.250 sz_alpha Factor for Wave dissipation surface tke flux .
100.000 crgban_cw Factor for Craig/Banner surface tke flux.
0.000 wec_alpha WEC factor for roller/breaking energy distribution.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
2.5000E-02 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
2.0000E-02 Zos Surface roughness (m).
1.0000E-01 Dcrit Minimum depth for wetting and drying (m).
1 levsfrc Deepest level to apply surface stress as a
bodyforce.
1 levbfrc Shallowest level to apply bottom stress as a
bodyforce.
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
6.0000E+00 theta_s S-coordinate surface control parameter.
4.0000E-01 theta_b S-coordinate bottom control parameter.
10.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
0.000 dstart Time-stamp assigned to model initialization (days).
0.00 time_ref Reference time for units attribute (yyyymmdd.dd)
0.0000E+00 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
0.0000E+00 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
0.0000E+00 Tnudg_SSS Nudging/relaxation time scale (days)
for sea surface salinity.
0.0000E+00 Znudg Nudging/relaxation time scale (days)
for free-surface.
0.0000E+00 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
0.0000E+00 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
0.0000E+00 obcfac Factor between passive and active
open boundary conditions.
F VolCons(1) NLM western edge boundary volume conservation.
F VolCons(2) NLM southern edge boundary volume conservation.
F VolCons(3) NLM eastern edge boundary volume conservation.
F VolCons(4) NLM northern edge boundary volume conservation.
10.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
1025.000 R0 Background density (kg/m3) used in linear Equation
of State.
1.7000E-04 Tcoef Thermal expansion coefficient (1/Celsius).
7.6000E-04 Scoef Saline contraction coefficient (1/PSU).
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
F LuvSrc Turning OFF momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
F LtracerSrc(01) Turning OFF point Sources/Sinks on tracer 01: temp
F LtracerSrc(02) Turning OFF point Sources/Sinks on tracer 02: salt
F LsshCLM Turning OFF processing of SSH climatology.
F Lm2CLM Turning OFF processing of 2D momentum climatology.
F Lm3CLM Turning OFF processing of 3D momentum climatology.
F LtracerCLM(01) Turning OFF processing of climatology tracer 01: temp
F LtracerCLM(02) Turning OFF processing of climatology tracer 02: salt
F LnudgeM2CLM Turning OFF nudging of 2D momentum climatology.
F LnudgeM3CLM Turning OFF nudging of 3D momentum climatology.
F LnudgeTCLM(01) Turning OFF nudging of climatology tracer 01: temp
F LnudgeTCLM(02) Turning OFF nudging of climatology tracer 02: salt
T Hout(idFsur) Write out free-surface.
T Hout(idUvel) Write out 3D U-momentum component.
T Hout(idVvel) Write out 3D V-momentum component.
T Hout(idWvel) Write out W-momentum component.
Output/Input Files:
Output Restart File: simplecase_rst.nc
Output History File: simplecase_his.nc
Output Averages File: simplecase_avg.nc
Input Grid File: /home/ubuntu/app/COAWST_V3.2/Projects/simplecase/grid.nc
Input Nonlinear Initial File: /home/ubuntu/app/COAWST_V3.2/Projects/simplecase/init.nc
Input Forcing File 01: /home/ubuntu/app/COAWST_V3.2/Projects/simplecase/stressTEST.nc
Tile partition information for Grid 01: 0176x0279x0025 tiling: 001x001
tile Istr Iend Jstr Jend Npts
Number of tracers: 2
0 1 176 1 279 1227600
Tile minimum and maximum fractional coordinates for Grid 01:
(interior points only)
tile Xmin Xmax Ymin Ymax grid
0 0.50 176.50 0.50 279.50 RHO-points
0 1.00 176.00 0.50 279.50 U-points
0 0.50 176.50 1.00 279.00 V-points
Maximum halo size in XI and ETA directions:
HaloSizeI(1) = 374
HaloSizeJ(1) = 578
TileSide(1) = 283
TileSize(1) = 51223
Lateral Boundary Conditions: NLM
============================
Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------
zeta 1 Chapman Imp Chapman Imp Chapman Imp Chapman Imp
ubar 1 Closed Flather Closed Flather
vbar 1 Closed Flather Closed Flather
u 1 Closed Gradient Closed Gradient
v 1 Closed Gradient Closed Gradient
temp 1 Closed Closed Closed Closed
salt 1 Closed Closed Closed Closed
tke 1 Closed Gradient Closed Gradient
Activated C-preprocessing Options:
PALU palu
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_FSOBC Analytical free-surface boundary conditions.
ANA_M2OBC Analytical 2D momentum boundary conditions.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ANA_STFLUX Analytical kinematic surface temperature flux.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged nonlinear model fields.
BODYFORCE Momentum stresses as body-forces.
CURVGRID Orthogonal curvilinear grid.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
GLS_MIXING Generic Length-Scale turbulence closure.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
MASKING Land/Sea masking.
MIX_GEO_TS Mixing of tracers along geopotential surfaces.
MIX_S_UV Mixing of momentum along constant S-surfaces.
MPI MPI distributed-memory configuration.
NONLINEAR Nonlinear Model.
!NONLIN_EOS Linear Equation of State for seawater.
N2S2_HORAVG Horizontal smoothing of buoyancy and shear.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
K_GSCHEME Third-order upstream advection of TKE fields.
RI_SPLINES Parabolic Spline Reconstruction for Richardson Number.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLVE3D Solving 3D Primitive Equations.
TS_U3HADVECTION Third-order upstream horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection of tracers.
TS_DIF2 Harmonic mixing of tracers.
UV_ADV Advection of momentum.
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum.
UV_C4VADVECTION Fourth-order centered vertical advection of momentum.
UV_LOGDRAG Logarithmic bottom stress.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
WET_DRY Wetting and drying activated.
Process Information:
Node # 0 (pid= 27497) is active.
INITIAL: Configuring and initializing forward nonlinear model ...
*******
Vertical S-coordinate System, Grid 01:
level S-coord Cs-curve Z at hmin at hc half way at hmax
25 0.0000000 0.0000000 0.000 0.000 0.000 0.000
24 -0.0400000 -0.0001749 -0.060 -0.201 -0.308 -0.538
23 -0.0800000 -0.0007098 0.582 -0.404 -0.410 -1.368
22 -0.1200000 -0.0016353 1.988 -0.608 -0.290 -2.516
21 -0.1600000 -0.0030046 4.260 -0.815 0.084 -4.025
20 -0.2000000 -0.0048963 7.553 -1.024 0.755 -5.960
19 -0.2400000 -0.0074189 12.079 -1.237 1.785 -8.407
18 -0.2800000 -0.0107165 18.118 -1.454 3.258 -11.485
17 -0.3200000 -0.0149770 26.037 -1.675 5.279 -15.347
16 -0.3600000 -0.0204424 36.311 -1.902 7.988 -20.189
15 -0.4000000 -0.0274214 49.540 -2.137 11.559 -26.262
14 -0.4400000 -0.0363057 66.492 -2.382 16.216 -33.886
13 -0.4800000 -0.0475899 88.131 -2.638 22.241 -43.462
12 -0.5200000 -0.0618955 115.671 -2.909 29.988 -55.497
11 -0.5600000 -0.0799996 150.631 -3.200 39.901 -70.622
10 -0.6000000 -0.1028685 194.898 -3.514 52.530 -89.623
9 -0.6400000 -0.1316956 250.803 -3.858 68.555 -113.473
8 -0.6800000 -0.1679410 321.197 -4.240 88.809 -143.358
7 -0.7200000 -0.2133714 409.533 -4.667 114.299 -180.716
6 -0.7600000 -0.2700910 519.919 -5.150 146.224 -227.259
5 -0.8000000 -0.3405555 657.152 -5.703 185.984 -284.985
4 -0.8400000 -0.4275481 826.669 -6.338 235.166 -356.158
3 -0.8800000 -0.5340940 1034.379 -7.070 295.494 -443.241
2 -0.9200000 -0.6632759 1286.303 -7.916 368.726 -548.740
1 -0.9600000 -0.8179050 1587.932 -8.890 456.463 -674.944
0 -1.0000000 -1.0000000 1943.209 -10.000 559.857 -823.494
Time Splitting Weights for Grid 01: ndtfast = 20 nfast = 29
==================================
Primary Secondary Accumulated to Current Step
1-0.0009651193358779 0.0500000000000000-0.0009651193358779 0.0500000000000000
2-0.0013488780126037 0.0500482559667939-0.0023139973484816 0.1000482559667939
3-0.0011514592651644 0.0501156998674241-0.0034654566136460 0.1501639558342180
4-0.0003735756740661 0.0501732728306823-0.0038390322877122 0.2003372286649003
5 0.0009829200513762 0.0501919516143856-0.0028561122363360 0.2505291802792859
6 0.0029141799764308 0.0501428056118168 0.0000580677400949 0.3006719858911027
7 0.0054132615310267 0.0499970966129953 0.0054713292711215 0.3506690825040980
8 0.0084687837865133 0.0497264335364439 0.0139401130576348 0.4003955160405419
9 0.0120633394191050 0.0493029943471183 0.0260034524767398 0.4496985103876601
10 0.0161716623600090 0.0486998273761630 0.0421751148367488 0.4983983377638231
11 0.0207585511322367 0.0478912442581626 0.0629336659689855 0.5462895820219857
12 0.0257765478740990 0.0468533167015507 0.0887102138430846 0.5931428987235364
13 0.0311633730493853 0.0455644893078458 0.1198735868924699 0.6387073880313822
14 0.0368391158442262 0.0440063206553765 0.1567127027366961 0.6827137086867586
15 0.0427031802506397 0.0421643648631652 0.1994158829873358 0.7248780735499238
16 0.0486309868367616 0.0400292058506332 0.2480468698240974 0.7649072794005570
17 0.0544704302037592 0.0375976565087951 0.3025173000278565 0.8025049359093521
18 0.0600380921294286 0.0348741349986072 0.3625553921572851 0.8373790709079593
19 0.0651152103984763 0.0318722303921357 0.4276706025557614 0.8692513013000950
20 0.0694434033194839 0.0286164698722119 0.4971140058752453 0.8978677711723069
21 0.0727201499285570 0.0251442997062377 0.5698341558038023 0.9230120708785446
22 0.0745940258796570 0.0215082922098099 0.6444281816834592 0.9445203630883545
23 0.0746596950216179 0.0177785909158270 0.7190878767050771 0.9622989540041815
24 0.0724526566618460 0.0140456061647461 0.7915405333669231 0.9763445601689277
25 0.0674437485167025 0.0104229733316538 0.8589842818836255 0.9867675335005816
26 0.0590334053485720 0.0070507859058187 0.9180176872321975 0.9938183194064003
27 0.0465456732896125 0.0040991156383901 0.9645633605218099 0.9979174350447905
28 0.0292219798521905 0.0017718319739095 0.9937853403740003 0.9996892670187000
29 0.0062146596259993 0.0003107329813000 0.9999999999999997 1.0000000000000000
ndtfast, nfast = 20 29 nfast/ndtfast = 1.45000
Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):
1.000000000000 1.060707743385 0.530353871693 1.000000000000 1.000000000000
Power filter parameters, Fgamma, gamma = 0.28400 0.14200
Metrics information for Grid 01:
===============================
Minimum X-grid spacing, DXmin = 1.00256452E-01 km
Maximum X-grid spacing, DXmax = 1.00262565E-01 km
Minimum Y-grid spacing, DYmin = 9.95247380E-02 km
Maximum Y-grid spacing, DYmax = 9.95366015E-02 km
Minimum Z-grid spacing, DZmin = -7.95450909E+03 m
Maximum Z-grid spacing, DZmax = 3.92962302E+03 m
Minimum barotropic Courant Number = 1.83413880E-04
Maximum barotropic Courant Number = 6.36202595E-03
Maximum Coriolis Courant Number = 2.26677842E-07
NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(Grid 01, File: init.nc, Rec=0001, Index=1)
- free-surface
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- vertically integrated v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- u-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- v-momentum component
(Min = 0.00000000E+00 Max = 0.00000000E+00)
- potential temperature
(Min = 1.00000000E+01 Max = 1.00000000E+01)
- salinity
(Min = 3.50000000E+01 Max = 3.50000000E+01)
GET_2DFLD - surface u-momentum stress, t = 0 00:00:00
(Rec=0000001, Index=1, File: stressTEST.nc)
(Tmin= 0.0000 Tmax= 3001.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface v-momentum stress, t = 0 00:00:00
(Rec=0000001, Index=1, File: stressTEST.nc)
(Tmin= 0.0000 Tmax= 3001.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
Basin information for Grid 01:
Maximum grid stiffness ratios: rx0 = 1.500000E-01 (Beckmann and Haidvogel)
rx1 = 5.728206E+00 (Haney)
Initial basin volumes: TotVolume = 6.3638374556E+10 m3
MinVolume = 2.5571097553E+02 m3
MaxVolume = 1.4785505052E+06 m3
Max/Min = 5.7821159305E+03
NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00029000)
GET_2DFLD - surface u-momentum stress, t = 0 02:24:00
(Rec=0000002, Index=2, File: stressTEST.nc)
(Tmin= 0.0000 Tmax= 3001.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface v-momentum stress, t = 0 02:24:00
(Rec=0000002, Index=2, File: stressTEST.nc)
(Tmin= 0.0000 Tmax= 3001.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME
C => (i,j,k) Cu Cv Cw Max Speed
0 0 00:00:00 0.000000E+00 9.300949E+03 9.300949E+03 6.367420E+10
(000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_HIS - creating history file, Grid 01: palu_his.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
DEF_AVG - creating average file, Grid 01: palu_avg.nc
1 0 00:00:00 7.403984E-35 9.301714E+03 9.301714E+03 6.366897E+10
(073,204,01) 2.838148E-20 3.182801E-20 0.000000E+00 7.471913E-17
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000002
2 0 00:00:00 2.892330E-34 9.301714E+03 9.301714E+03 6.366897E+10
(058,226,25) 6.956945E-23 4.667068E-22 1.367662E-16 1.488681E-16
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000003
3 0 00:00:00 6.180541E-34 9.301714E+03 9.301714E+03 6.366897E+10
(058,226,25) 2.184770E-22 1.522074E-21 2.705075E-16 2.451863E-16
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000004
4 0 00:00:00 1.027599E-33 9.301714E+03 9.301714E+03 6.366897E+10
(048,263,25) 3.982611E-22 6.445907E-21 3.350156E-16 4.986654E-16
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000005
5 0 00:00:00 1.477361E-33 9.301714E+03 9.301714E+03 6.366897E+10
(057,240,25) 2.241456E-21 7.977925E-21 4.357218E-16 7.637114E-16
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000006
6 0 00:00:00 1.952023E-33 9.301714E+03 9.301714E+03 6.366897E+10
(086,192,25) 3.322778E-21 3.886765E-21 4.632618E-16 1.028624E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000007
7 0 00:00:00 2.444697E-33 9.301714E+03 9.301714E+03 6.366897E+10
(086,192,25) 4.586721E-21 2.804571E-21 5.329079E-16 1.293538E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000008
8 0 00:00:00 3.007035E-33 9.301714E+03 9.301714E+03 6.366897E+10
(047,269,25) 5.364246E-21 2.059773E-20 5.409963E-16 1.558241E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000009
9 0 00:00:00 3.690863E-33 9.301714E+03 9.301714E+03 6.366897E+10
(059,225,25) 1.215610E-20 3.419549E-21 6.021904E-16 1.822926E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000010
10 0 00:00:00 4.620259E-33 9.301714E+03 9.301714E+03 6.366897E+10
(058,238,25) 3.370067E-20 1.212422E-20 8.214504E-16 2.087500E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000011
11 0 00:00:01 5.830281E-33 9.301714E+03 9.301714E+03 6.366897E+10
(052,249,25) 2.251131E-21 7.650399E-21 1.012366E-15 2.351966E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000012
12 0 00:00:01 7.741074E-33 9.301714E+03 9.301714E+03 6.366897E+10
(057,270,25) 2.112356E-21 3.158259E-20 1.175885E-15 2.616385E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000013
13 0 00:00:01 1.033973E-32 9.301714E+03 9.301714E+03 6.366897E+10
(057,270,25) 2.272178E-21 3.745048E-20 1.355493E-15 2.880627E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000014
14 0 00:00:01 1.377428E-32 9.301714E+03 9.301714E+03 6.366897E+10
(055,247,25) 2.595464E-20 5.314083E-20 1.593740E-15 3.144802E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000015
15 0 00:00:01 1.782040E-32 9.301714E+03 9.301714E+03 6.366897E+10
(055,247,25) 2.649002E-20 4.711058E-20 1.979339E-15 3.408829E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000016
16 0 00:00:01 2.255626E-32 9.301714E+03 9.301714E+03 6.366897E+10
(055,247,25) 3.381972E-20 6.569763E-20 2.318928E-15 3.672661E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000017
17 0 00:00:01 2.762860E-32 9.301714E+03 9.301714E+03 6.366897E+10
(063,217,25) 1.388260E-19 1.283977E-19 2.858346E-15 3.936374E-15
WRT_HIS - wrote history fields (Index=2,2) into time record = 0000018
18 0 00:00:01 3.298347E-32 9.301714E+03 9.301714E+03 6.366897E+10
(063,217,25) 1.484557E-19 1.637055E-19 3.103611E-15 4.199973E-15
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000019
19 0 00:00:01 NaN NaN NaN 6.366897E+10
(063,217,25) 1.554516E-19 1.820026E-19 3.170670E-15 1.839072E-15
Blowing-up: Saving latest model state into RESTART file
WRT_RST - wrote re-start fields (Index=2,2) into time record = 0000001
Elapsed CPU time (seconds):
Node # 0 CPU: 21.440
Total: 21.440
Nonlinear ocean model elapsed time profile, Grid: 01
Allocation and array initialization .............. 0.119 ( 0.5533 %)
Ocean state initialization ....................... 0.015 ( 0.0717 %)
Reading of input data ............................ 0.006 ( 0.0274 %)
Processing of input data ......................... 0.006 ( 0.0264 %)
Processing of output time averaged data .......... 0.011 ( 0.0502 %)
Computation of vertical boundary conditions ...... 0.031 ( 0.1456 %)
Computation of global information integrals ...... 0.135 ( 0.6315 %)
Writing of output data ........................... 2.427 (11.3187 %)
Model 2D kernel .................................. 3.937 (18.3644 %)
2D/3D coupling, vertical metrics ................. 0.145 ( 0.6742 %)
Omega vertical velocity .......................... 0.113 ( 0.5256 %)
Equation of state for seawater ................... 0.125 ( 0.5815 %)
GLS vertical mixing parameterization ............. 9.713 (45.3013 %)
3D equations right-side terms .................... 0.700 ( 3.2662 %)
3D equations predictor step ...................... 0.729 ( 3.4019 %)
Pressure gradient ................................ 0.358 ( 1.6705 %)
Harmonic mixing of tracers, geopotentials ........ 0.812 ( 3.7854 %)
Harmonic stress tensor, S-surfaces ............... 0.304 ( 1.4182 %)
Corrector time-step for 3D momentum .............. 0.555 ( 2.5868 %)
Corrector time-step for tracers .................. 0.628 ( 2.9313 %)
Total: 20.868 97.3319
Nonlinear model message Passage profile, Grid: 01
Message Passage: 2D halo exchanges ............... 0.011 ( 0.0507 %)
Message Passage: 3D halo exchanges ............... 0.001 ( 0.0032 %)
Message Passage: 4D halo exchanges ............... 0.000 ( 0.0007 %)
Message Passage: data broadcast .................. 0.001 ( 0.0060 %)
Message Passage: data reduction .................. 0.000 ( 0.0008 %)
Message Passage: data gathering .................. 0.310 ( 1.4437 %)
Message Passage: data scattering.................. 0.020 ( 0.0942 %)
Message Passage: boundary data gathering ......... 0.000 ( 0.0006 %)
Total: 0.343 1.5999
All percentages are with respect to total time = 21.440
ROMS/TOMS - Output NetCDF summary for Grid 01:
number of time records written in HISTORY file = 00000019
number of time records written in RESTART file = 00000001
Analytical header files used:
ROMS/Functionals/ana_btflux.h
/home/ubuntu/app/COAWST_V3.2/Projects/simplecase/ana_fsobc.h
/home/ubuntu/app/COAWST_V3.2/Projects/simplecase/ana_m2obc.h
ROMS/Functionals/ana_stflux.h
ROMS/TOMS: DONE... Thursday - January 9, 2020 - 6:57:05 PM