ROMS blows up after define MASKING

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
novice
Posts: 27
Joined: Fri Jan 18, 2013 6:02 pm
Location: Indian Institute of Technology

ROMS blows up after define MASKING

#1 Unread post by novice »

Hi Forum,

I am using ROMS model configured for the Arabian Sea in northern Indian Ocean. My model blows up soon after the start of simulation (climatological). To diagnose the problem, I run the model with different set of suitable cpp options and found that model runs fine without MASKING and blows up when #define MASKING given. I generated the masking files for u, v, rho, psi following the romstools guidelines i.e.

[Mp,Lp]=size(rfield);
M=Mp-1;
L=Lp-1;
%
vfield=0.5*(rfield(1:M,:)+rfield(2:Mp,:));
ufield=0.5*(rfield(:,1:L)+rfield(:,2:Lp));
pfield=0.5*(ufield(1:M,:)+ufield(2:Mp,:));

and my model grid is:--

netcdf roms_grd {
dimensions:
xi_u = 120 ;
eta_u = 85 ;
xi_v = 121 ;
eta_v = 84 ;
xi_rho = 121 ;
eta_rho = 85 ;
xi_psi = 120 ;
eta_psi = 84 ;
one = 1 ;
two = 2 ;
four = 4 ;
bath = 1 ;

Bathymetry (etopo2) smoothen by 0.1 r value to avoid any chances of bathymetry related issue.
Kindly suggest to solve the problem.
Following is my model configuration and simulation step....


Model Input Parameters: ROMS/TOMS version 3.7
Monday - February 3, 2014 - 4:38:14 PM
-----------------------------------------------------------------------------

Trial run for Indian Ocean

Operating system : Linux
CPU/hardware : x86_64
Compiler system : gfortran
Compiler command : /usr/bin/mpif90
Compiler flags : -frepack-arrays -O3 -ffast-math -ffree-form -ffree-line-length-none -ffree-form -ffree-line-length-none -ffree-form -ffree-line-length-none

Input Script :

SVN Root URL : https:://myroms.org/svn/src
SVN Revision : 669M

Local Root : /home/osf123/src0
Header Dir : /home/osf123/src0/test/test_as1
Header file : as_test1.h
Analytical Dir: /home/osf123/src0/test/test_as1

Resolution, Grid 01: 0119x0083x020, Parallel Nodes: 30, Tiling: 006x005


Physical Parameters, Grid: 01
=============================

2592000 ntimes Number of timesteps for 3-D equations.
60.000 dt Timestep size (s) for 3-D equations.
30 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.
864 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).
43200 nHIS Number of timesteps between the writing fields
into history file.
518400 ndefHIS Number of timesteps between creation of new
history files.
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
43200 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
518400 ndefAVG Number of timesteps between creation of new
time-averaged file.
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.
1.0000E-04 rdrg Linear bottom drag coefficient (m/s).
1.0000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
1 Vtransform S-coordinate transformation equation.
1 Vstretching S-coordinate stretching function.
6.0000E+00 theta_s S-coordinate surface control parameter.
0.0000E+00 theta_b S-coordinate bottom control parameter.
10.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1027.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 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.
-0.278 T0 Background potential temperature (C) constant.
34.649 S0 Background salinity (PSU) constant.
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
T LtracerCLM(01) Processing climatology on tracer 01: temp
T LtracerCLM(02) Processing climatology on tracer 02: salt
T Hout(idFsur) Write out free-surface.
T Hout(idUbar) Write out 2D U-momentum component.
T Hout(idVbar) Write out 2D V-momentum component.
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.
T Hout(idOvel) Write out omega vertical velocity.
T Hout(idTvar) Write out tracer 01: temp
T Hout(idTvar) Write out tracer 02: salt

T Aout(idFsur) Write out averaged free-surface.
T Aout(idUbar) Write out averaged 2D U-momentum component.
T Aout(idVbar) Write out averaged 2D V-momentum component.
T Aout(idUvel) Write out averaged 3D U-momentum component.
T Aout(idVvel) Write out averaged 3D V-momentum component.
T Aout(idWvel) Write out averaged W-momentum component.
T Aout(idOvel) Write out averaged omega vertical velocity.
T Aout(idTvar) Write out averaged tracer 01: temp
T Aout(idTvar) Write out averaged tracer 02: salt

Output/Input Files:

Output Restart File: ocean_rst.nc
Prefix for History Files: ocean_his
Prefix for Averages Files: ocean_avg
Input Grid File: roms_grd.nc
Input Nonlinear Initial File: roms_ini.nc
Input Forcing File 01: roms_frc.nc
Input Climatology File: roms_clm.nc

Tile partition information for Grid 01: 0119x0083x0020 tiling: 006x005

tile Istr Iend Jstr Jend Npts

0 1 20 1 16 6400
1 21 40 1 16 6400
2 41 60 1 16 6400
3 61 80 1 16 6400
4 81 100 1 16 6400
5 101 119 1 16 6080
6 1 20 17 33 6800
7 21 40 17 33 6800
8 41 60 17 33 6800
9 61 80 17 33 6800
10 81 100 17 33 6800
11 101 119 17 33 6460
12 1 20 34 50 6800
13 21 40 34 50 6800
14 41 60 34 50 6800
15 61 80 34 50 6800
16 81 100 34 50 6800
17 101 119 34 50 6460
18 1 20 51 67 6800
19 21 40 51 67 6800
20 41 60 51 67 6800
21 61 80 51 67 6800
22 81 100 51 67 6800
23 101 119 51 67 6460
24 1 20 68 83 6400
25 21 40 68 83 6400
26 41 60 68 83 6400
27 61 80 68 83 6400
28 81 100 68 83 6400
29 101 119 68 83 6080

Tile minimum and maximum fractional coordinates for Grid 01:
(interior points only)

tile Xmin Xmax Ymin Ymax grid

0 0.50 20.50 0.50 16.50 RHO-points
1 20.50 40.50 0.50 16.50 RHO-points
2 40.50 60.50 0.50 16.50 RHO-points
3 60.50 80.50 0.50 16.50 RHO-points
4 80.50 100.50 0.50 16.50 RHO-points
5 100.50 119.50 0.50 16.50 RHO-points
6 0.50 20.50 16.50 33.50 RHO-points
7 20.50 40.50 16.50 33.50 RHO-points
8 40.50 60.50 16.50 33.50 RHO-points
9 60.50 80.50 16.50 33.50 RHO-points
10 80.50 100.50 16.50 33.50 RHO-points
11 100.50 119.50 16.50 33.50 RHO-points
12 0.50 20.50 33.50 50.50 RHO-points
13 20.50 40.50 33.50 50.50 RHO-points
14 40.50 60.50 33.50 50.50 RHO-points
15 60.50 80.50 33.50 50.50 RHO-points
16 80.50 100.50 33.50 50.50 RHO-points
17 100.50 119.50 33.50 50.50 RHO-points
18 0.50 20.50 50.50 67.50 RHO-points
19 20.50 40.50 50.50 67.50 RHO-points
20 40.50 60.50 50.50 67.50 RHO-points
21 60.50 80.50 50.50 67.50 RHO-points
22 80.50 100.50 50.50 67.50 RHO-points
23 100.50 119.50 50.50 67.50 RHO-points
24 0.50 20.50 67.50 83.50 RHO-points
25 20.50 40.50 67.50 83.50 RHO-points
26 40.50 60.50 67.50 83.50 RHO-points
27 60.50 80.50 67.50 83.50 RHO-points
28 80.50 100.50 67.50 83.50 RHO-points
29 100.50 119.50 67.50 83.50 RHO-points

0 1.00 20.50 0.50 16.50 U-points
1 20.50 40.50 0.50 16.50 U-points
2 40.50 60.50 0.50 16.50 U-points
3 60.50 80.50 0.50 16.50 U-points
4 80.50 100.50 0.50 16.50 U-points
5 100.50 119.00 0.50 16.50 U-points
6 1.00 20.50 16.50 33.50 U-points
7 20.50 40.50 16.50 33.50 U-points
8 40.50 60.50 16.50 33.50 U-points
9 60.50 80.50 16.50 33.50 U-points
10 80.50 100.50 16.50 33.50 U-points
11 100.50 119.00 16.50 33.50 U-points
12 1.00 20.50 33.50 50.50 U-points
13 20.50 40.50 33.50 50.50 U-points
14 40.50 60.50 33.50 50.50 U-points
15 60.50 80.50 33.50 50.50 U-points
16 80.50 100.50 33.50 50.50 U-points
17 100.50 119.00 33.50 50.50 U-points
18 1.00 20.50 50.50 67.50 U-points
19 20.50 40.50 50.50 67.50 U-points
20 40.50 60.50 50.50 67.50 U-points
21 60.50 80.50 50.50 67.50 U-points
22 80.50 100.50 50.50 67.50 U-points
23 100.50 119.00 50.50 67.50 U-points
24 1.00 20.50 67.50 83.50 U-points
25 20.50 40.50 67.50 83.50 U-points
26 40.50 60.50 67.50 83.50 U-points
27 60.50 80.50 67.50 83.50 U-points
28 80.50 100.50 67.50 83.50 U-points
29 100.50 119.00 67.50 83.50 U-points

0 0.50 20.50 1.00 16.50 V-points
1 20.50 40.50 1.00 16.50 V-points
2 40.50 60.50 1.00 16.50 V-points
3 60.50 80.50 1.00 16.50 V-points
4 80.50 100.50 1.00 16.50 V-points
5 100.50 119.50 1.00 16.50 V-points
6 0.50 20.50 16.50 33.50 V-points
7 20.50 40.50 16.50 33.50 V-points
8 40.50 60.50 16.50 33.50 V-points
9 60.50 80.50 16.50 33.50 V-points
10 80.50 100.50 16.50 33.50 V-points
11 100.50 119.50 16.50 33.50 V-points
12 0.50 20.50 33.50 50.50 V-points
13 20.50 40.50 33.50 50.50 V-points
14 40.50 60.50 33.50 50.50 V-points
15 60.50 80.50 33.50 50.50 V-points
16 80.50 100.50 33.50 50.50 V-points
17 100.50 119.50 33.50 50.50 V-points
18 0.50 20.50 50.50 67.50 V-points
19 20.50 40.50 50.50 67.50 V-points
20 40.50 60.50 50.50 67.50 V-points
21 60.50 80.50 50.50 67.50 V-points
22 80.50 100.50 50.50 67.50 V-points
23 100.50 119.50 50.50 67.50 V-points
24 0.50 20.50 67.50 83.00 V-points
25 20.50 40.50 67.50 83.00 V-points
26 40.50 60.50 67.50 83.00 V-points
27 60.50 80.50 67.50 83.00 V-points
28 80.50 100.50 67.50 83.00 V-points
29 100.50 119.50 67.50 83.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 96
HaloSizeJ(1) = 87
TileSide(1) = 26
TileSize(1) = 598


Lateral Boundary Conditions: NLM
============================

Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------

zeta 1 Chapman Imp Periodic Chapman Imp Periodic

ubar 1 Flather Periodic Flather Periodic

vbar 1 Flather Periodic Flather Periodic

u 1 Radiation Periodic Radiation Periodic

v 1 Radiation Periodic Radiation Periodic

temp 1 Radiation Periodic Radiation Periodic

salt 1 Radiation Periodic Radiation Periodic

tke 1 Gradient Periodic Gradient Periodic

Activated C-preprocessing Options:

AS_TEST1 Trial run for Indian Ocean
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.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged nonlinear model fields.
CURVGRID Orthogonal curvilinear grid.
DEBUGGING Internal debugging switch activated.
DOUBLE_PRECISION Double precision arithmetic.
M3CLIMATOLOGY Processing 3D momentum climatology data.
MASKING Land/Sea masking.
MPI MPI distributed-memory configuration.
MY25_MIXING Mellor/Yamada Level-2.5 mixing closure.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
POWER_LAW Power-law shape time-averaging barotropic filter.
PRSGRD31 Standard density Jacobian formulation (Song, 1998).
PROFILE Time profiling activated .
QCORRECTION Surface net heat flux correction.
K_GSCHEME Third-order upstream advection of TKE fields.
RHO_SURF Include difference between rho0 and surface density.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLVE3D Solving 3D Primitive Equations.
TCLIMATOLOGY Processing tracer climatology data.
TS_C4HADVECTION Fourth-order centered horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection of tracers.
UV_ADV Advection of momentum.
UV_COR Coriolis term.
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum.
UV_C4VADVECTION Fourth-order centered vertical advection of momentum.
UV_LDRAG Linear bottom stress.
VAR_RHO_2D Variable density barotropic mode.
ZCLIMATOLOGY Processing sea surface height climatology data.

Process Information:

Node # 0 (pid= 13229) is active.
Node # 1 (pid= 13230) is active.
Node # 9 (pid= 13238) is active.
Node # 16 (pid= 13245) is active.
Node # 17 (pid= 13246) is active.
Node # 18 (pid= 13247) is active.
Node # 25 (pid= 13254) is active.
Node # 26 (pid= 13255) is active.
Node # 2 (pid= 13231) is active.
Node # 3 (pid= 13232) is active.
Node # 8 (pid= 13237) is active.
Node # 10 (pid= 13239) is active.
Node # 11 (pid= 13240) is active.
Node # 19 (pid= 13248) is active.
Node # 24 (pid= 13253) is active.
Node # 27 (pid= 13256) is active.
Node # 4 (pid= 13233) is active.
Node # 5 (pid= 13234) is active.
Node # 6 (pid= 13235) is active.
Node # 7 (pid= 13236) is active.
Node # 12 (pid= 13241) is active.
Node # 13 (pid= 13242) is active.
Node # 14 (pid= 13243) is active.
Node # 15 (pid= 13244) is active.
Node # 20 (pid= 13249) is active.
Node # 21 (pid= 13250) is active.
Node # 22 (pid= 13251) is active.
Node # 23 (pid= 13252) is active.
Node # 28 (pid= 13257) is active.
Node # 29 (pid= 13258) 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

20 0.0000000 0.0000000 0.000 0.000 0.000 0.000
19 -0.0500000 -0.0015097 -0.598 -0.500 -4.263 -7.928
18 -0.1000000 -0.0031562 -1.205 -1.000 -8.867 -16.529
17 -0.1500000 -0.0050890 -1.831 -1.500 -14.184 -26.538
16 -0.2000000 -0.0074832 -2.486 -2.000 -20.652 -38.817
15 -0.2500000 -0.0105560 -3.186 -2.500 -28.811 -54.435
14 -0.3000000 -0.0145859 -3.948 -3.000 -39.355 -74.762
13 -0.3500000 -0.0199385 -4.796 -3.500 -53.196 -101.597
12 -0.4000000 -0.0270990 -5.761 -4.000 -71.544 -137.326
11 -0.4500000 -0.0367168 -6.887 -4.500 -96.016 -185.145
10 -0.5000000 -0.0496640 -8.228 -5.000 -128.786 -249.344
9 -0.5500000 -0.0671145 -9.862 -5.500 -172.781 -335.700
8 -0.6000000 -0.0906508 -11.892 -6.000 -231.945 -451.998
7 -0.6500000 -0.1224070 -14.456 -6.500 -311.597 -608.737
6 -0.7000000 -0.1652627 -17.742 -7.000 -418.914 -820.085
5 -0.7500000 -0.2231040 -22.002 -7.500 -563.582 -1105.162
4 -0.8000000 -0.3011757 -27.576 -8.000 -758.674 -1489.771
3 -0.8500000 -0.4065570 -34.926 -8.500 -1021.834 -2008.742
2 -0.9000000 -0.5488038 -44.672 -9.000 -1376.881 -2709.090
1 -0.9500000 -0.7408145 -57.653 -9.500 -1855.963 -3654.274
0 -1.0000000 -1.0000000 -75.000 -10.000 -2502.478 -4929.955

Time Splitting Weights for Grid 01: ndtfast = 30 nfast = 42
==================================

Primary Secondary Accumulated to Current Step

1-0.0008094437383769 0.0333333333333333-0.0008094437383769 0.0333333333333333
2-0.0014053566728197 0.0333603147912792-0.0022148004111966 0.0666936481246126
3-0.0017877524645903 0.0334071600137066-0.0040025528757869 0.1001008081383191
4-0.0019566842408176 0.0334667517625262-0.0059592371166046 0.1335675599008453
5-0.0019122901320372 0.0335319745705535-0.0078715272486418 0.1670995344713988
6-0.0016548570247459 0.0335957175749547-0.0095263842733877 0.2006952520463536
7-0.0011849025289723 0.0336508794757796-0.0107112868023600 0.2343461315221331
8-0.0005032751608631 0.0336903762267453-0.0112145619632232 0.2680365077488784
9 0.0003887272597151 0.0337071520654408-0.0108258347035081 0.3017436598143192
10 0.0014892209965583 0.0336941944901169-0.0093366137069498 0.3354378543044362
11 0.0027955815694920 0.0336445537902317-0.0065410321374578 0.3690824080946679
12 0.0043042707117221 0.0335513677379153-0.0022367614257356 0.4026337758325831
13 0.0060106451121704 0.0334078920475245 0.0037738836864348 0.4360416678801076
14 0.0079087469427945 0.0332075372104522 0.0116826306292293 0.4692492050905598
15 0.0099910761708920 0.0329439123123590 0.0216737068001213 0.5021931174029188
16 0.0122483446563884 0.0326108764399960 0.0339220514565097 0.5348039938429148
17 0.0146692120341107 0.0322025982847830 0.0485912634906204 0.5670065921276978
18 0.0172400033810439 0.0317136245503127 0.0658312668716643 0.5987202166780105
19 0.0199444086685725 0.0311389577709445 0.0857756755402368 0.6298591744489550
20 0.0227631639997064 0.0304741441486588 0.1085388395399432 0.6603333185976138
21 0.0256737146312911 0.0297153720153352 0.1342125541712342 0.6900486906129490
22 0.0286498597812016 0.0288595815276255 0.1628624139524359 0.7189082721405746
23 0.0316613792205220 0.0279045862015855 0.1945237931729578 0.7468128583421600
24 0.0346736416507075 0.0268492068942347 0.2291974348236653 0.7736620652363948
25 0.0376471948657328 0.0256934188392112 0.2668446296893981 0.7993554840756060
26 0.0405373376992232 0.0244385123436867 0.3073819673886213 0.8237939964192927
27 0.0432936737565710 0.0230872677537126 0.3506756411451924 0.8468812641730054
28 0.0458596469320356 0.0216441452951603 0.3965352880772280 0.8685254094681656
29 0.0481720587108284 0.0201154903974257 0.4447073467880565 0.8886408998655914
30 0.0501605672561820 0.0185097551070648 0.4948679140442384 0.9071506549726561
31 0.0517471682814030 0.0168377361985254 0.5466150823256415 0.9239883911711815
32 0.0528456577069106 0.0151128305891453 0.5994607400325521 0.9391012217603267
33 0.0533610761022577 0.0133513086655816 0.6528218161348098 0.9524525304259084
34 0.0531891349131379 0.0115726061288397 0.7060109510479478 0.9640251365547481
35 0.0522156244733761 0.0097996349650684 0.7582265755213239 0.9738247715198165
36 0.0503158038019030 0.0080591141492892 0.8085423793232269 0.9818838856691057
37 0.0473537721847153 0.0063819206892258 0.8558961515079423 0.9882658063583315
38 0.0431818225418188 0.0048034616164019 0.8990779740497611 0.9930692679747334
39 0.0376397765791564 0.0033640675316746 0.9367177506289175 0.9964333355064080
40 0.0305543017255206 0.0021094083123694 0.9672720523544381 0.9985427438187774
41 0.0217382098544504 0.0010909315881854 0.9890102622088885 0.9996336754069628
42 0.0109897377911118 0.0003663245930371 1.0000000000000004 0.9999999999999999

ndtfast, nfast = 30 42 nfast/ndtfast = 1.40000

Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):

1.000000000000 1.047601458608 0.523800729304 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.18933

Metrics information for Grid 01:
===============================

Minimum X-grid spacing, DXmin = 2.51922285E+01 km
Maximum X-grid spacing, DXmax = 2.76666997E+01 km
Minimum Y-grid spacing, DYmin = 2.52153517E+01 km
Maximum Y-grid spacing, DYmax = 2.76720933E+01 km
Minimum Z-grid spacing, DZmin = 5.98128547E-01 m
Maximum Z-grid spacing, DZmax = 1.27568114E+03 m

Minimum barotropic Courant Number = 3.02208616E-03
Maximum barotropic Courant Number = 2.24802584E-02
Maximum Coriolis Courant Number = 3.67988396E-03


NLM: GET_STATE - Read state initial conditions, t = 0 00:00:00
(Grid 01, File: roms_ini.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.09958712E+00 Max = 2.89255675E+01)
- salinity
(Min = 3.23184533E+01 Max = 4.17184288E+01)
GET_2DFLD - surface u-momentum stress, t = 345 00:00:00
(Rec=0000012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -5.92220444E-05 Max = 2.50951080E-05)
GET_2DFLD - surface v-momentum stress, t = 345 00:00:00
(Rec=0000012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -8.02755791E-05 Max = -7.98899145E-06)
GET_2DFLD - surface net heat flux, t = 345 00:00:00
(Rec=0000012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -4.40694246E-05 Max = 1.33105372E-05)
GET_2DFLD - sea surface temperature climatology, t = 345 00:00:00
(Rec=0000012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.43573818E+01 Max = 2.83237804E+01)
GET_2DFLD - surface net heat flux sensitivity to SST, t = 345 00:00:00
(Rec=0000012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.23133173E-05 Max = -6.12059303E-06)
GET_2DFLD - surface net freswater flux, (E-P), t = 345 00:00:00
(Rec=0000012, Index=1, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -7.53600626E-08 Max = 9.57819544E-08)
GET_2DFLD - sea surface height climatology, t = 345 00:00:00
(Rec=0000012, Index=1, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.99058766E-01 Max = 2.64657134E-01)
GET_3DFLD - potential temperature climatology, t = 345 00:00:00
(Rec=0000012, Index=1, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.09958712E+00 Max = 2.89663960E+01)
GET_3DFLD - salinity climatology, t = 345 00:00:00
(Rec=0000012, Index=1, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.21366289E+01 Max = 4.12563504E+01)
GET_3DFLD - u-momentum component climatology, t = 345 00:00:00
(Rec=0000012, Index=1, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.61682382E-01 Max = 3.79378982E-01)
GET_3DFLD - v-momentum component climatology, t = 345 00:00:00
(Rec=0000012, Index=1, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -4.43708835E-01 Max = 2.21494382E-01)

Basin information for Grid 01:

Maximum grid stiffness ratios: rx0 = 9.429072E-01 (Beckmann and Haidvogel)
rx1 = 6.368397E+00 (Haney)

Initial basin volumes: TotVolume = 1.0928206657E+16 m3
MinVolume = 3.7995029515E+08 m3
MaxVolume = 9.7665573449E+11 m3
Max/Min = 2.5704828946E+03

NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 02592000)

GET_2DFLD - surface u-momentum stress, t = 15 00:00:00
(Rec=0000001, Index=2, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -6.29371451E-05 Max = 1.86384726E-05)
GET_2DFLD - surface v-momentum stress, t = 15 00:00:00
(Rec=0000001, Index=2, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -7.70138644E-05 Max = -4.20460912E-06)
GET_2DFLD - surface net heat flux, t = 15 00:00:00
(Rec=0000001, Index=2, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -3.04500745E-05 Max = 1.14627872E-05)
GET_2DFLD - sea surface temperature climatology, t = 15 00:00:00
(Rec=0000001, Index=2, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.19725246E+01 Max = 2.82567107E+01)
GET_2DFLD - surface net heat flux sensitivity to SST, t = 15 00:00:00
(Rec=0000001, Index=2, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.21069099E-05 Max = -5.86180488E-06)
GET_2DFLD - surface net freswater flux, (E-P), t = 15 00:00:00
(Rec=0000001, Index=2, File: roms_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -5.07285736E-08 Max = 9.24236958E-08)
GET_2DFLD - sea surface height climatology, t = 15 00:00:00
(Rec=0000001, Index=2, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -2.86289930E-01 Max = 2.95689316E-01)
GET_3DFLD - potential temperature climatology, t = 15 00:00:00
(Rec=0000001, Index=2, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 1.09958712E+00 Max = 2.89255675E+01)
GET_3DFLD - salinity climatology, t = 15 00:00:00
(Rec=0000001, Index=2, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.23184533E+01 Max = 4.17184288E+01)
GET_3DFLD - u-momentum component climatology, t = 15 00:00:00
(Rec=0000001, Index=2, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -4.17898139E-01 Max = 6.21943778E-01)
GET_3DFLD - v-momentum component climatology, t = 15 00:00:00
(Rec=0000001, Index=2, File: roms_clm.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -5.02392110E-01 Max = 3.89440053E-01)

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 1.559714E+04 1.559714E+04 1.127424E+16
(000,00,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_HIS - creating history file, Grid 01: ocean_his_0001.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
DEF_AVG - creating average file, Grid 01: ocean_avg_0001.nc
1 0 00:01:00 9.524038E-04 1.559708E+04 1.559708E+04 1.127421E+16
(021,83,01) 0.000000E+00 1.567634E-03 0.000000E+00 9.318216E-01
2 0 00:02:00 3.576918E-03 1.559692E+04 1.559693E+04 1.127412E+16
(021,83,20) 0.000000E+00 1.259230E-04 6.227519E+00 1.811713E+00
3 0 00:03:00 7.415215E-03 1.559668E+04 1.559669E+04 1.127398E+16
(021,83,20) 0.000000E+00 3.661471E-04 1.179486E+01 2.619368E+00
4 0 00:04:00 1.213967E-02 1.559637E+04 1.559638E+04 1.127380E+16
(021,83,20) 0.000000E+00 7.284765E-04 1.640104E+01 3.365025E+00
5 0 00:05:00 1.762348E-02 1.559601E+04 1.559603E+04 1.127359E+16
(021,83,20) 0.000000E+00 1.172137E-03 2.014017E+01 4.067640E+00
6 0 00:06:00 2.382833E-02 1.559563E+04 1.559565E+04 1.127335E+16
(021,83,20) 0.000000E+00 1.658619E-03 2.331828E+01 4.742009E+00
7 0 00:07:00 3.073024E-02 1.559522E+04 1.559525E+04 1.127310E+16
(021,83,20) 0.000000E+00 2.169071E-03 2.621270E+01 5.395613E+00
8 0 00:08:00 3.831151E-02 1.559480E+04 1.559484E+04 1.127283E+16
(021,83,20) 0.000000E+00 2.699173E-03 2.895292E+01 6.030531E+00
9 0 00:09:00 4.655473E-02 1.559437E+04 1.559441E+04 1.127255E+16
(021,83,20) 0.000000E+00 3.252172E-03 3.155294E+01 6.647103E+00
10 0 00:10:00 5.543465E-02 1.559392E+04 1.559398E+04 1.127226E+16
(021,83,20) 0.000000E+00 3.829872E-03 3.399906E+01 7.246263E+00
11 0 00:11:00 6.492676E-02 1.559346E+04 1.559353E+04 1.127196E+16
(021,83,20) 0.000000E+00 4.427987E-03 3.629865E+01 7.829489E+00
12 0 00:12:00 7.502018E-02 1.559299E+04 1.559307E+04 1.127164E+16
(021,83,20) 0.000000E+00 5.044217E-03 3.847580E+01 8.398327E+00
13 0 00:13:00 8.570729E-02 1.559251E+04 1.559259E+04 1.127132E+16
(021,83,20) 0.000000E+00 5.676931E-03 4.055026E+01 8.953969E+00
14 0 00:14:00 9.698193E-02 1.559201E+04 1.559211E+04 1.127100E+16
(021,83,20) 0.000000E+00 6.325238E-03 4.253392E+01 9.497639E+00
15 0 00:15:00 1.088335E-01 1.559150E+04 1.559161E+04 1.127066E+16
(021,83,20) 0.000000E+00 6.988407E-03 4.443437E+01 1.003026E+01
16 0 00:16:00 1.212453E-01 1.559099E+04 1.559111E+04 1.127032E+16
(021,83,20) 0.000000E+00 7.666602E-03 4.625522E+01 1.055217E+01
17 0 00:17:00 1.341968E-01 1.559046E+04 1.559059E+04 1.126997E+16
(021,83,20) 0.000000E+00 8.360840E-03 4.799163E+01 1.106306E+01
18 0 00:18:00 1.476655E-01 1.558992E+04 1.559007E+04 1.126962E+16
(021,83,20) 0.000000E+00 9.071715E-03 4.963209E+01 1.156252E+01
19 0 00:19:00 1.616293E-01 1.558938E+04 1.558954E+04 1.126926E+16
(021,83,20) 0.000000E+00 9.797689E-03 5.116890E+01 1.205082E+01
20 0 00:20:00 1.760664E-01 1.558883E+04 1.558901E+04 1.126889E+16
(021,83,20) 0.000000E+00 1.053555E-02 5.260682E+01 1.252878E+01
21 0 00:21:00 1.909552E-01 1.558829E+04 1.558848E+04 1.126853E+16
(021,83,20) 0.000000E+00 1.128005E-02 5.396726E+01 1.299817E+01
22 0 00:22:00 2.062711E-01 1.558774E+04 1.558794E+04 1.126815E+16
(021,83,20) 0.000000E+00 1.202647E-02 5.527468E+01 1.346040E+01
23 0 00:23:00 2.219875E-01 1.558719E+04 1.558742E+04 1.126778E+16
(021,83,20) 0.000000E+00 1.277101E-02 5.654734E+01 1.391618E+01
24 0 00:24:00 2.380750E-01 1.558665E+04 1.558689E+04 1.126740E+16
(021,83,20) 0.000000E+00 1.351097E-02 5.779234E+01 1.436528E+01
25 0 00:25:00 2.545029E-01 1.558611E+04 1.558637E+04 1.126701E+16
(021,83,20) 0.000000E+00 1.424439E-02 5.900692E+01 1.480686E+01
26 0 00:26:00 2.712379E-01 1.558558E+04 1.558585E+04 1.126662E+16
(021,83,20) 0.000000E+00 1.496933E-02 6.018373E+01 1.523999E+01
27 0 00:27:00 2.882460E-01 1.558505E+04 1.558534E+04 1.126623E+16
(021,83,20) 0.000000E+00 1.568342E-02 6.131664E+01 1.566411E+01
28 0 00:28:00 3.054883E-01 1.558453E+04 1.558483E+04 1.126584E+16
(021,83,20) 0.000000E+00 1.638385E-02 6.240348E+01 1.607901E+01
29 0 00:29:00 3.229282E-01 1.558401E+04 1.558433E+04 1.126544E+16
(021,83,20) 0.000000E+00 1.706808E-02 6.344426E+01 1.648455E+01
30 0 00:30:00 3.405275E-01 1.558349E+04 1.558383E+04 1.126504E+16
(021,83,20) 0.000000E+00 1.773450E-02 6.443721E+01 1.688043E+01
31 0 00:31:00 3.582535E-01 1.558299E+04 1.558334E+04 1.126463E+16
(021,83,20) 0.000000E+00 1.838279E-02 6.537702E+01 1.726625E+01
32 0 00:32:00 3.760680E-01 1.558249E+04 1.558286E+04 1.126422E+16
(021,83,20) 0.000000E+00 1.901352E-02 6.625729E+01 1.764180E+01
33 0 00:33:00 3.939393E-01 1.558199E+04 1.558238E+04 1.126381E+16
(021,83,20) 0.000000E+00 1.962758E-02 6.707430E+01 1.800721E+01
34 0 00:34:00 4.118357E-01 1.558150E+04 1.558191E+04 1.126339E+16
(021,83,20) 0.000000E+00 2.022577E-02 6.782855E+01 1.836277E+01
35 0 00:35:00 4.297313E-01 1.558102E+04 1.558145E+04 1.126298E+16
(021,83,20) 0.000000E+00 2.080873E-02 6.852329E+01 1.870887E+01
36 0 00:36:00 4.475903E-01 1.558055E+04 1.558100E+04 1.126256E+16
(021,83,20) 0.000000E+00 2.137681E-02 6.916250E+01 1.904615E+01
37 0 00:37:00 4.653914E-01 1.558009E+04 1.558055E+04 1.126213E+16
(021,83,20) 0.000000E+00 2.193008E-02 6.975025E+01 1.941748E+01
38 0 00:38:00 4.831067E-01 1.557964E+04 1.558012E+04 1.126171E+16
(021,83,20) 0.000000E+00 2.246836E-02 7.029068E+01 1.977916E+01
39 0 00:39:00 5.007131E-01 1.557919E+04 1.557969E+04 1.126128E+16
(021,83,20) 0.000000E+00 2.299134E-02 7.078714E+01 2.013169E+01

Blowing-up: Saving latest model state into RESTART file

WRT_RST - wrote re-start fields (Index=1,2) into time record = 0000001

Elapsed CPU time (seconds):

Node # 0 CPU: 1.172
Node # 1 CPU: 1.288
Node # 6 CPU: 1.288
Node # 7 CPU: 1.292
Node # 9 CPU: 1.288
Node # 12 CPU: 1.296
Node # 14 CPU: 1.296
Node # 16 CPU: 1.288
Node # 17 CPU: 1.292
Node # 18 CPU: 1.292
Node # 21 CPU: 1.292
Node # 22 CPU: 1.292
Node # 23 CPU: 1.292
Node # 25 CPU: 1.292
Node # 26 CPU: 1.288
Node # 29 CPU: 1.284
Node # 2 CPU: 1.248
Node # 3 CPU: 1.280
Node # 4 CPU: 1.272
Node # 5 CPU: 1.280
Node # 8 CPU: 1.256
Node # 10 CPU: 1.284
Node # 11 CPU: 1.236
Node # 13 CPU: 1.280
Node # 15 CPU: 1.280
Node # 19 CPU: 1.264
Node # 20 CPU: 1.268
Node # 24 CPU: 1.264
Node # 27 CPU: 1.256
Node # 28 CPU: 1.288

ROMS/TOMS - Output NetCDF summary for Grid 01:
number of time records written in HISTORY file = 00000001
number of time records written in RESTART file = 00000001

Analytical header files used:

ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_fsobc.h
ROMS/Functionals/ana_m2obc.h
ROMS/Functionals/ana_nudgcoef.h

ROMS/TOMS: DONE... Monday - February 3, 2014 - 4:38:15 PM
------------------------
Kindly help to resolve the problem.
Thanks in advance...

regards,
Sandy

ptimko
Posts: 35
Joined: Wed Mar 02, 2011 6:46 pm
Location: Environment and Climate Change Canada

Re: ROMS blows up after define MASKING

#2 Unread post by ptimko »

Question...

Why do you have periodic boundary conditions?

you say you're trying to model:
I am using ROMS model configured for the Arabian Sea in northern Indian Ocean.
but your output file says your boundary conditions are:
Lateral Boundary Conditions: NLM
============================

Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------

zeta 1 Chapman Imp Periodic Chapman Imp Periodic

ubar 1 Flather Periodic Flather Periodic

vbar 1 Flather Periodic Flather Periodic

u 1 Radiation Periodic Radiation Periodic

v 1 Radiation Periodic Radiation Periodic

temp 1 Radiation Periodic Radiation Periodic

salt 1 Radiation Periodic Radiation Periodic

tke 1 Gradient Periodic Gradient Periodic
I'm either misunderstanding what you are trying to model or I suspect from what you say that you are trying to model a system that you should not be using periodic boundary conditions.

Post Reply