segment number difference error while coupling ROMS and SWAN

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
chirantan
Posts: 1
Joined: Wed Dec 23, 2015 2:05 am
Location: IIT Kharagpur

segment number difference error while coupling ROMS and SWAN

#1 Unread post by chirantan »

Hi,

I am trying to couple ROMS and SWAN with refinement grids. I have followed the steps for INLET_TEST/Refined. The ROMS and SWAN grid are formed using roms2swan and interpolation weights are obtained by modifying the file scrip_coawst_inlet_test_refined.in. I am getting segment number difference error in MCT. While running the SWAN only with refinement grids having same input files it works properly. The details of my error is shown below:

chirantan@chirantan-System-Product-Name:~/COAWST$ mpirun -np 8 ./coawstM Projects/BoB/coupling_inlet_test_ref3.in

Coupled Input File name = Projects/BoB/coupling_inlet_test_ref3.in Coupled Input File name = Projects/BoB/coupling_inlet_test_ref3.in

Model Coupling:

Ocean Model MPI nodes: 000 - 003

Waves Model MPI nodes: 004 - 007

Model Input Parameters: ROMS/TOMS version 3.7

Thursday - March 16, 2017 - 11:29:45 AM

-----------------------------------------------------------------------------

SWAN grid 1 is preparing computation

Bay of Bengal sediment transport with nesting

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 : Unversioned directory

Local Root : /home/omsairam/COAWST Header Dir : /home/omsairam/COAWST/Projects/BoB Header file : bob_refined.h Analytical Dir: /home/omsairam/COAWST/Projects/BoB

Resolution, Grid 01: 0287x0248x032, Parallel Nodes: 4, Tiling: 002x002

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

16992 ntimes Number of timesteps for 3-D equations.

300.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.

288 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).

1728 nHIS Number of timesteps between the writing fields

into history file.

8766 ndefHIS Number of timesteps between creation of new

history files.

1.0000E-03 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

5.0000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)

for tracer 01: temp

5.0000E-06 Akt_bak(02) Background vertical mixing coefficient (m2/s)

for tracer 02: salt

5.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. 1.5000E-02 Zob Bottom roughness (m). 5.0000E-01 Zos Surface roughness (m).

2 Vtransform S-coordinate transformation equation. 4 Vstretching S-coordinate stretching function.

7.0000E+00 theta_s S-coordinate surface control parameter. 2.0000E+00 theta_b S-coordinate bottom control parameter.

250.000 Tcline S-coordinate surface/bottom layer width (m) used

in vertical coordinate stretching.

1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.

1.000 dstart Time-stamp assigned to model initialization (days).

20040101.00 time_ref Reference time for units attribute (yyyymmdd.dd)

2.0000E+01 Tnudg(01) Nudging/relaxation time scale (days)

for tracer 01: temp

2.0000E+01 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.

2.0000E+01 M3nudg Nudging/relaxation time scale (days)

for 3D momentum.

2.0000E+01 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.

1027.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(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 Hout(idUbms) Write out bottom U-momentum stress. T Hout(idVbms) Write out bottom V-momentum stress. T Hout(idUbrs) Write out bottom U-current stress. T Hout(idVbrs) Write out bottom V-current stress. T Hout(idUbws) Write out wind-induced, bottom U-wave stress. T Hout(idVbws) Write out wind-induced, bottom V-wave stress. T Hout(idUbcs) Write out max wind + current, bottom U-wave stress. T Hout(idVbcs) Write out max wind + current, bottom V-wave stress. T Hout(idUVwc) Write out max wind + current, bottom UV-wave stress. T Hout(idUbot) Write out bed wave orbital U-velocity. T Hout(idVbot) Write out bed wave orbital V-velocity. T Hout(idUbur) Write out bottom U-momentum above bed. T Hout(idVbvr) Write out bottom V-momentum above bed. T Hout(idU2rs) Write out total 2D u-radiation stress. T Hout(idV2rs) Write out total 2D v-radiation stress. T Hout(idU2Sd) Write out 2D u-momentum stokes velocity. T Hout(idV2Sd) Write out 2D v-momentum stokes velocity. T Hout(idU3rs) Write out total 3D u-radiation stress. T Hout(idV3rs) Write out total 3D v-radiation stress. T Hout(idU3Sd) Write out 3D u-momentum stokes velocity. T Hout(idV3Sd) Write out 3D v-momentum stokes velocity. T Hout(idW3Sd) Write out 3D omega-momentum stokes velocity. T Hout(idW3St) Write out 3D w-momentum stokes velocity. T Hout(idWztw) Write out WEC quasi-static sea level adjustment. T Hout(idWqsp) Write out WEC quasi-static pressure. T Hout(idWbeh) Write out WEC Bernoulli head. T Hout(idWamp) Write out wave height. T Hout(idWlen) Write out wavelength. T Hout(idWdir) Write out wave direction. T Hout(idWptp) Write out wave surface period. T Hout(idWpbt) Write out wave bottom period. T Hout(idWorb) Write out wave bottom orbital velocity. T Hout(idWdif) Write out wave dissipation due to bottom friction. T Hout(idWdib) Write out wave dissipation due to breaking. T Hout(idWdiw) Write out wave dissipation due to white capping. T Hout(idUwav) Wave-avg surface u-velocity. T Hout(idVwav) Wave-avg surface v-velocity. T Hout(idVvis) Write out vertical viscosity: AKv. T Hout(idMtke) Write out turbulent kinetic energy. T Hout(idMtls) Write out turbulent generic length-scale.

Output/Input Files:

Output Restart File: Projects/BoB/out/ocean_rst_ref1.nc

Prefix for History Files: Projects/BoB/out/ocean_inlet_his_ref1 Physical parameters File:

Input Grid File: Projects/BoB/in/roms_grd.nc

Nesting grid connectivity File: Projects/BoB/in/roms_out.nc

Input Boundary File 01: Projects/BoB/in/mercator_bry.nc

Resolution, Grid 02: 0320x0170x032, Parallel Nodes: 4, Tiling: 002x002

Physical Parameters, Grid: 02 =============================

101952 ntimes Number of timesteps for 3-D equations. 50.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).

1728 nHIS Number of timesteps between the writing fields

into history file.

1.0000E-03 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

5.0000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)

for tracer 01: temp

5.0000E-06 Akt_bak(02) Background vertical mixing coefficient (m2/s)

for tracer 02: salt

5.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. 1.5000E-02 Zob Bottom roughness (m). 5.0000E-01 Zos Surface roughness (m).

2 Vtransform S-coordinate transformation equation. 4 Vstretching S-coordinate stretching function.

7.0000E+00 theta_s S-coordinate surface control parameter. 2.0000E+00 theta_b S-coordinate bottom control parameter.

250.000 Tcline S-coordinate surface/bottom layer width (m) used

in vertical coordinate stretching.

1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.

1.000 dstart Time-stamp assigned to model initialization (days).

20040101.00 time_ref Reference time for units attribute (yyyymmdd.dd)

2.0000E+01 Tnudg(01) Nudging/relaxation time scale (days)

for tracer 01: temp

2.0000E+01 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.

2.0000E+01 M3nudg Nudging/relaxation time scale (days)

for 3D momentum.

2.0000E+01 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.

1027.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(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 Hout(idUbms) Write out bottom U-momentum stress. T Hout(idVbms) Write out bottom V-momentum stress. T Hout(idUbrs) Write out bottom U-current stress. T Hout(idVbrs) Write out bottom V-current stress. T Hout(idUbws) Write out wind-induced, bottom U-wave stress. T Hout(idVbws) Write out wind-induced, bottom V-wave stress. T Hout(idUbcs) Write out max wind + current, bottom U-wave stress. T Hout(idVbcs) Write out max wind + current, bottom V-wave stress. T Hout(idUVwc) Write out max wind + current, bottom UV-wave stress. T Hout(idUbot) Write out bed wave orbital U-velocity. T Hout(idVbot) Write out bed wave orbital V-velocity. T Hout(idUbur) Write out bottom U-momentum above bed. T Hout(idVbvr) Write out bottom V-momentum above bed. T Hout(idU2rs) Write out total 2D u-radiation stress. T Hout(idV2rs) Write out total 2D v-radiation stress. T Hout(idU2Sd) Write out 2D u-momentum stokes velocity. T Hout(idV2Sd) Write out 2D v-momentum stokes velocity. T Hout(idU3rs) Write out total 3D u-radiation stress. T Hout(idV3rs) Write out total 3D v-radiation stress. T Hout(idU3Sd) Write out 3D u-momentum stokes velocity. T Hout(idV3Sd) Write out 3D v-momentum stokes velocity. T Hout(idW3Sd) Write out 3D omega-momentum stokes velocity. T Hout(idW3St) Write out 3D w-momentum stokes velocity. T Hout(idWztw) Write out WEC quasi-static sea level adjustment. T Hout(idWqsp) Write out WEC quasi-static pressure. T Hout(idWbeh) Write out WEC Bernoulli head. T Hout(idWamp) Write out wave height. T Hout(idWlen) Write out wavelength. T Hout(idWdir) Write out wave direction. T Hout(idWptp) Write out wave surface period. T Hout(idWpbt) Write out wave bottom period. T Hout(idWorb) Write out wave bottom orbital velocity. T Hout(idWdif) Write out wave dissipation due to bottom friction. T Hout(idWdib) Write out wave dissipation due to breaking. T Hout(idWdiw) Write out wave dissipation due to white capping. T Hout(idUwav) Wave-avg surface u-velocity. T Hout(idVwav) Wave-avg surface v-velocity. T Hout(idVvis) Write out vertical viscosity: AKv. T Hout(idMtke) Write out turbulent kinetic energy. T Hout(idMtls) Write out turbulent generic length-scale.

Output/Input Files:

Output Restart File: Projects/BoB/out/ocean_rst_ref3.nc Output History File: Projects/BoB/out/ocean_inlet_his_ref3.nc

Physical parameters File:

Input Grid File: Projects/BoB/in/roms_nesting.nc

Nesting grid connectivity File: Projects/BoB/in/roms_out.nc

Tile partition information for Grid 01: 0287x0248x0032 tiling: 002x002

tile Istr Iend Jstr Jend Npts

Number of tracers: 3

0 1 144 1 124 571392 1 145 287 1 124 567424 2 1 144 125 248 571392 3 145 287 125 248 567424

Tile partition information for Grid 02: 0320x0170x0032 tiling: 002x002

tile Istr Iend Jstr Jend Npts

Number of tracers: 3

0 1 160 1 85 435200 1 161 320 1 85 435200 2 1 160 86 170 435200 3 161 320 86 170 435200

Tile minimum and maximum fractional coordinates for Grid 01:

(interior points only)

tile Xmin Xmax Ymin Ymax grid

0 0.50 144.50 0.50 124.50 RHO-points 1 144.50 287.50 0.50 124.50 RHO-points 2 0.50 144.50 124.50 248.50 RHO-points 3 144.50 287.50 124.50 248.50 RHO-points

0 1.00 144.50 0.50 124.50 U-points 1 144.50 287.00 0.50 124.50 U-points 2 1.00 144.50 124.50 248.50 U-points 3 144.50 287.00 124.50 248.50 U-points

0 0.50 144.50 1.00 124.50 V-points 1 144.50 287.50 1.00 124.50 V-points 2 0.50 144.50 124.50 248.00 V-points 3 144.50 287.50 124.50 248.00 V-points

Tile minimum and maximum fractional coordinates for Grid 02:

(interior points only)

tile Xmin Xmax Ymin Ymax grid

0 0.50 160.50 0.50 85.50 RHO-points 1 160.50 320.50 0.50 85.50 RHO-points 2 0.50 160.50 85.50 170.50 RHO-points 3 160.50 320.50 85.50 170.50 RHO-points

0 1.00 160.50 0.50 85.50 U-points 1 160.50 320.00 0.50 85.50 U-points 2 1.00 160.50 85.50 170.50 U-points 3 160.50 320.00 85.50 170.50 U-points

0 0.50 160.50 1.00 85.50 V-points 1 160.50 320.50 1.00 85.50 V-points 2 0.50 160.50 85.50 170.00 V-points 3 160.50 320.50 85.50 170.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 468 HaloSizeJ(1) = 411

TileSide?(1) = 150 TileSize?(1) = 19650

Maximum halo size in XI and ETA directions:

HaloSizeI(2) = 525 HaloSizeJ(2) = 300

TileSide?(2) = 169 TileSize?(2) = 15886

Sediment Parameters, Grid: 01 =============================

Size Sd50 Csed Srho Wsed Erate poros Class (mm) (kg/m3) (kg/m3) (mm/s) (kg/m2/s) (nondim)

1 1.0000E-01 0.0000E+00 2.6500E+03 1.0000E+01 5.0000E-03 0.0000E+00

tau_ce tau_cd nl_tnu2 nl_tnu4 Akt_bak Tnudg (N/m2) (N/m2) (m2/s) (m4/s) (m2/s) (day)

1 1.0000E-01 1.0000E-02 0.0000E+00 0.0000E+00 5.0000E-06 0.0000E+00

morph_fac (nondim)

1 1.0000E+01

New bed layer formed when deposition exceeds 0.10000E-01 (m). Two first layers are combined when 2nd layer smaller than 0.00000E+00 (m). Rate coefficient for bed load transport = 0.15000E+00

F LtracerSponge?(03Turning OFF sponge on tracer 03: mud_01 F LtracerSrc?(03) Turning OFF point sources/Sink on tracer 03: mud_01 F LtracerCLM(03) Turning OFF processing of climatology tracer 03: mud_01 F LnudgeTCLM(03) Turning OFF nudging of climatology tracer 03: mud_01 T Hout(idTvar) Write out sediment01: mud_01 T Hout(idfrac) Write out bed fraction, sediment 01: mudfrac_01 T Hout(idfrac) Write out mass, sediment 01: mudmass_01 T Hout(idSbed) Write out BED property 01: bed_thickness T Hout(idSbed) Write out BED property 02: bed_age T Hout(idSbed) Write out BED property 03: bed_porosity T Hout(idSbed) Write out BED property 04: bed_biodiff T Hout(idBott) Write out BOTTOM property01: grain_diameter T Hout(idBott) Write out BOTTOM property02: grain_density T Hout(idBott) Write out BOTTOM property03: settling_vel T Hout(idBott) Write out BOTTOM property04: erosion_stress T Hout(idBott) Write out BOTTOM property05: ripple_length T Hout(idBott) Write out BOTTOM property06: ripple_height T Hout(idBott) Write out BOTTOM property07: bed_wave_amp T Hout(idBott) Write out BOTTOM property08: Zo_def T Hout(idBott) Write out BOTTOM property09: Zo_app T Hout(idBott) Write out BOTTOM property15: active_layer_thickness T Hout(idBott) Write out BOTTOM property17: bed_inund_depth T Hout(idBott) Write out BOTTOM property18: dep_net

Sediment Parameters, Grid: 02 =============================

Size Sd50 Csed Srho Wsed Erate poros Class (mm) (kg/m3) (kg/m3) (mm/s) (kg/m2/s) (nondim)

1 1.0000E-01 0.0000E+00 2.6500E+03 1.0000E+01 5.0000E-03 0.0000E+00

tau_ce tau_cd nl_tnu2 nl_tnu4 Akt_bak Tnudg (N/m2) (N/m2) (m2/s) (m4/s) (m2/s) (day)

1 1.0000E-01 1.0000E-02 0.0000E+00 0.0000E+00 5.0000E-06 0.0000E+00

morph_fac (nondim)

1 1.0000E+01

New bed layer formed when deposition exceeds 0.10000E-01 (m). Two first layers are combined when 2nd layer smaller than 0.00000E+00 (m). Rate coefficient for bed load transport = 0.15000E+00

F LtracerSponge?(03Turning OFF sponge on tracer 03: mud_01 F LtracerSrc?(03) Turning OFF point sources/Sink on tracer 03: mud_01 F LtracerCLM(03) Turning OFF processing of climatology tracer 03: mud_01 F LnudgeTCLM(03) Turning OFF nudging of climatology tracer 03: mud_01 T Hout(idTvar) Write out sediment01: mud_01 T Hout(idfrac) Write out bed fraction, sediment 01: mudfrac_01 T Hout(idfrac) Write out mass, sediment 01: mudmass_01 T Hout(idSbed) Write out BED property 01: bed_thickness T Hout(idSbed) Write out BED property 02: bed_age T Hout(idSbed) Write out BED property 03: bed_porosity T Hout(idSbed) Write out BED property 04: bed_biodiff T Hout(idBott) Write out BOTTOM property01: grain_diameter T Hout(idBott) Write out BOTTOM property02: grain_density T Hout(idBott) Write out BOTTOM property03: settling_vel T Hout(idBott) Write out BOTTOM property04: erosion_stress T Hout(idBott) Write out BOTTOM property05: ripple_length T Hout(idBott) Write out BOTTOM property06: ripple_height T Hout(idBott) Write out BOTTOM property07: bed_wave_amp T Hout(idBott) Write out BOTTOM property08: Zo_def T Hout(idBott) Write out BOTTOM property09: Zo_app T Hout(idBott) Write out BOTTOM property15: active_layer_thickness T Hout(idBott) Write out BOTTOM property17: bed_inund_depth T Hout(idBott) Write out BOTTOM property18: dep_net

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

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

zeta 1 Chapman Imp Chapman Imp Chapman Imp Closed

2

@ �� �33333 3�?`

ubar 1 Flather Flather Flather Closed

2

@ �� �33333 3�?`

vbar 1 Flather Flather Flather Closed

2

@ �� �33333 3�?`

u 1 Rad + Nud Rad + Nud Rad + Nud Closed

2

@ �� �33333 3�?`

v 1 Rad + Nud Rad + Nud Rad + Nud Closed

2

@ �� �33333 3�?`

temp 1 Rad + Nud Rad + Nud Rad + Nud Closed

2

@ �� �33333 3�?`

salt 1 Rad + Nud Rad + Nud Rad + Nud Closed

2

@ �� �33333 3�?`

mud_01 1 Gradient Closed Gradient Gradient

2 Nested Nested Nested Nested

tke 1 Gradient Gradient Gradient Closed

2 Nested Nested Nested Nested

ubar_stokes 1 Gradient Closed Gradient Gradient

2 Nested Nested Nested Nested

vbar_stokes 1 Gradient Closed Gradient Gradient

2 Nested Nested Nested Nested

u_stokes 1 Gradient Closed Gradient Gradient

2 Nested Nested Nested Nested

v_stokes 1 Gradient Closed Gradient Gradient

2 Nested Nested Nested Nested

Activated C-preprocessing Options:

BOB_REFINED Bay of Bengal sediment transport with nesting ANA_BPFLUX Analytical bottom passive tracers fluxes. ANA_BSFLUX Analytical kinematic bottom salinity flux. ANA_BTFLUX Analytical kinematic bottom temperature flux. ANA_FSOBC Analytical free-surface boundary conditions. ANA_INITIAL Analytical initial conditions. ANA_M2OBC Analytical 2D momentum boundary conditions. ANA_SEDIMENT Analytical sediment initial conditions. ANA_SMFLUX Analytical kinematic surface momentum flux. ANA_SPFLUX Analytical surface passive tracer fluxes. ANA_SSFLUX Analytical kinematic surface salinity flux. ANA_STFLUX Analytical kinematic surface temperature flux. ASSUMED_SHAPE Using assumed-shape arrays. 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. MCT_LIB Using Model Coupling Toolkit library. MIX_S_UV Mixing of momentum along constant S-surfaces. MPI MPI distributed-memory configuration. WEC_VF Vortex Force wave current interaction. WDISS_WAVEMOD Wave energy dissipation acquired from coupled wave model. NESTING Nesting grids: Composite and Refinement. NONLINEAR Nonlinear Model. !NONLIN_EOS Linear Equation of State for seawater. N2S2_HORAVG Horizontal smoothing of buoyancy and shear. !ONE_WAY Two-way nesting in refinement grids. 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. SEDIMENT Cohesive and noncohesive sediments. SUSPLOAD Activate suspended sediment transport. NONCOHESIVE_BED1 Activate non-cohesive (original) sediment bed. SOLVE3D Solving 3D Primitive Equations. SPLINES_VDIFF Parabolic Spline Reconstruction for Vertical Diffusion. SPLINES_VVISC Parabolic Spline Reconstruction for Vertical Viscosity. SSW_BBL Styles and Glenn Bottom Boundary Layer - modified. SSW_CALC_ZNOT Internal computation of bottom roughness. SWAN_COUPLING SWAN model coupling. TS_MPDATA Recursive flux corrected MPDATA 3D advection 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_KIRBY Compute uwave and vwave Kirby avg velocities. UV_VIS2 Harmonic mixing of momentum. VAR_RHO_2D Variable density barotropic mode. WAVES_OCEAN Two-way wave-ocean models coupling.

Process Information:

Node # 1 (pid= 17599) is active. Node # 2 (pid= 17600) is active. Node # 3 (pid= 17602) is active. Node # 0 (pid= 17598) is active. Node # 0 (pid= 17598) is active. Node # 3 (pid= 17602) is active. Node # 1 (pid= 17599) is active. Node # 2 (pid= 17600) is active.

SWAN grid 2 is preparing computation

MCT::m_SparseMatrixDecompComputeSegments_: segment number difference error, stat =-32541 000.MCT(MPEU)::die.: from MCT::m_SparseMatrixDecompComputeSegments_()

MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD with errorcode 2.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them.

I shall be thankful for your suggestions.

Chirantan

Post Reply