Dear All,
i am running model for inter-annual for upto 2015. my model running upto 2013 is fine, but then 2014 intially its blow up. i can't fined the error. Please any body help..me.
herewith i am sending model details also.
Model Input Parameters: ROMS/TOMS version 3.6
Sunday - November 13, 2016 - 3:25:18 PM
-----------------------------------------------------------------------------
Indian Ocean ROMS #3.6, 0.1 Resolution
Operating system : Linux
CPU/hardware : x86_64
Compiler system : ifort
Compiler command : /apps/mvapich2-gcc/bin/mpif90
Compiler flags : -heap-arrays -fp-model precise -ip -O2 -msse4.2 -free -free -free -free -free
Input Script : /home/sreeja/romsp1/SMAG_run/interannual_2013-2014/indian_ocean.in
SVN Root URL : https://www.myroms.org/svn/src/trunk
SVN Revision : 583M
Local Root : /home/sreeja/ROMS/models/roms3.6_new
Header Dir : /home/sreeja/ROMS/models/roms3.6_new/ROMS/Include
Header file : indianocean.h
Analytical Dir: /home/sreeja/ROMS/models/roms3.6_new/ROMS/Functionals
Resolution, Grid 01: 0899x0629x040, Parallel Nodes: 288, Tiling: 024x012
Physical Parameters, Grid: 01
=============================
5469984 ntimes Number of timesteps for 3-D equations.
150.000 dt Timestep size (s) for 3-D equations.
15 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.
10519200 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).
17532 nHIS Number of timesteps between the writing fields
into history file.
17532 ndefHIS Number of timesteps between creation of new
1 ntsAVG Starting timestep for the accumulation of output
time-averaged data.
576 nAVG Number of timesteps between the writing of
time-averaged data into averages file.
576 ndefAVG Number of timesteps between creation of new
time-averaged file.
5.0000E+00 nl_tnu2(01) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 01: temp
5.0000E+00 nl_tnu2(02) NLM Horizontal, harmonic mixing coefficient
(m2/s) for tracer 02: salt
5.0000E-03 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
1.0000E+02 nl_visc4 NLM Horizontal, biharmonic mixing coefficient
(m4/s) for momentum.
5.0000E-07 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
5.0000E-07 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.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
3.0000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
1.0000E+01 blk_ZQ Height (m) of surface air humidity measurement.
1.0000E+01 blk_ZT Height (m) of surface air temperature measurement.
1.0000E+01 blk_ZW Height (m) of surface winds measurement.
1 lmd_Jwt Jerlov water type.
1 Vtransform S-coordinate transformation equation.
1 Vstretching S-coordinate stretching function.
7.0000E+00 theta_s S-coordinate surface control parameter.
2.0000E-01 theta_b S-coordinate bottom control parameter.
50.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).
19910101.00 time_ref Reference time for units attribute (yyyymmdd.dd)
7.0000E+00 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
7.0000E+00 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02: salt
7.0000E+00 Znudg Nudging/relaxation time scale (days)
for free-surface.
7.0000E+00 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
7.0000E+00 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
1.2000E+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.
T 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.
-1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
F LtracerSrc(01) Processing point sources/Sink on tracer 01: temp
T LtracerSrc(02) Processing point sources/Sink 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 Hout(idDano) Write out density anomaly.
T Hout(idVvis) Write out vertical viscosity: AKv.
T Hout(idTdif) Write out vertical diffusion: AKt(itemp).
T Hout(idSdif) Write out vertical diffusion: AKt(isalt).
T Hout(idHsbl) Write out depth of surface boundary layer.
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: output/ocean_rst_s.nc
Prefix for History Files: output/ocean_his_ia
Prefix for Averages Files: output/ocean_avg_ia
Input Grid File: /home/sreeja/romsp1/romsfiles_p1_newgrid_sep8_2012/roms_grd.nc
Input Nonlinear Initial File: ./output/ocean_his_ia_0300.nc
Input Forcing File 01: /home/sreeja/romsp1/interannual_inputs/14-16/airtemp_13_16.nc
Input Forcing File 02: /home/sreeja/romsp1/interannual_inputs/14-16/lwrad_13_16.nc
Input Forcing File 03: /home/sreeja/romsp1/interannual_inputs/14-16/precip_13_16.nc
Input Forcing File 04: /home/sreeja/romsp1/interannual_inputs/14-16/sh_13_16.nc
Input Forcing File 05: /home/sreeja/romsp1/interannual_inputs/14-16/slp_13_16.nc
Input Forcing File 06: /home/sreeja/romsp1/interannual_inputs/14-16/swrad_13_16.nc
Input Forcing File 07: /home/sreeja/romsp1/interannual_inputs/14-16/uwind_13_16.nc
Input Forcing File 08: /home/sreeja/romsp1/interannual_inputs/14-16/vwind_13_16.nc
Input Forcing File 09: /home/sreeja/romsp1/romsfiles_p1_newgrid_sep8_2012/roms_frc.nc
Input Forcing File 10: /home/sreeja/romsp1/river_input/large_gb_ir/roms_river_ir_gb_KG_MH_p1.nc
Input Boundary File: /home/sreeja/romsp1/romsfiles_p1_newgrid_sep8_2012/roms_bry.nc
Tile partition information for Grid 01: 0899x0629x0040 tiling: 024x012
tile Istr Iend Jstr Jend Npts
0 1 32 1 50 64000
1 33 70 1 50 76000
2 71 108 1 50 76000
3 109 146 1 50 76000
4 147 184 1 50 76000
5 185 222 1 50 76000
6 223 260 1 50 76000
7 261 298 1 50 76000
8 299 336 1 50 76000
9 337 374 1 50 76000
10 375 412 1 50 76000
11 413 450 1 50 76000
12 451 488 1 50 76000
Maximum halo size in XI and ETA directions:
HaloSizeI(1) = 156
HaloSizeJ(1) = 201
TileSide(1) = 61
TileSize(1) = 2806
Lateral Boundary Conditions: NLM
============================
Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------
zeta 1 Closed Chapman Chapman Closed
ubar 1 Closed Flather Flather Closed
vbar 1 Closed Flather Flather Closed
u 1 Closed Rad + Nud Rad + Nud Closed
v 1 Closed Rad + Nud Rad + Nud Closed
temp 1 Closed Rad + Nud Rad + Nud Closed
salt 1 Closed Rad + Nud Rad + Nud Closed
Activated C-preprocessing Options:
Minimum horizontal viscosity coefficient = 8.02177632E+01 m2/s
Maximum horizontal viscosity coefficient = 9.26106741E+01 m2/s
NLM: GET_STATE - Read state initial conditions, t = 9100 19:30:00
(File: ocean_his_ia_0300.nc, Rec=0001, Index=1)
- free-surface
(Min = -3.90019208E-01 Max = 1.18797386E+00)
- vertically integrated u-momentum component
(Min = -1.00571752E+00 Max = 1.38676095E+00)
- vertically integrated v-momentum component
(Min = -7.06835926E-01 Max = 1.10835421E+00)
- u-momentum component
(Min = -2.95794272E+00 Max = 4.63440323E+00)
- v-momentum component
(Min = -2.31057048E+00 Max = 2.22049809E+00)
- potential temperature
(Min = 5.52539825E-01 Max = 3.35481720E+01)
- salinity
(Min = 2.15598464E-01 Max = 6.07375221E+01)
- vertical viscosity coefficient
(Min = 0.00000000E+00 Max = 1.79698721E-01)
- temperature vertical diffusion coefficient
(Min = 0.00000000E+00 Max = 3.97602588E-01)
- salinity vertical diffusion coefficient
(Min = 0.00000000E+00 Max = 3.97602588E-01)
GET_NGFLD - river runoff XI-positions at RHO-points
(Min = 5.10000000E+02 Max = 6.53000000E+02)
GET_NGFLD - river runoff ETA-positions at RHO-points
(Min = 4.74000000E+02 Max = 5.48000000E+02)
GET_NGFLD - river runoff direction
(Min = 1.00000000E+00 Max = 1.00000000E+00)
GET_NGFLD - river runoff mass transport vertical profile
(Min = 0.00000000E+00 Max = 2.50000000E-02)
GET_NGFLD - river runoff mass transport, t = 75 00:00:00
(Rec=0003, Index=2, File: roms_river_ir_gb_KG_MH_p1.nc)
(Tmin= 15.0000 Tmax= 345.0000)
GET_NGFLD - river runoff salinity, t = 75 00:00:00
(Rec=0003, Index=2, File: roms_river_ir_gb_KG_MH_p1.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 1.00000000E+01)
GET_2DFLD - surface u-wind component, t = 9100 00:00:00
(Rec=0334, Index=1, File: uwind_13_16.nc)
(Tmin= 8767.0000 Tmax= 10135.0000)
(Min = -1.19670162E+01 Max = 1.04587536E+01)
GET_2DFLD - surface v-wind component, t = 9100 00:00:00
(Rec=0334, Index=1, File: vwind_13_16.nc)
(Tmin= 8767.0000 Tmax= 9861.0000)
(Min = -1.42567110E+01 Max = 1.20455198E+01)
GET_2DFLD - surface air pressure, t = 9100 00:00:00
(Rec=0334, Index=1, File: slp_13_16.nc)
(Tmin= 8767.0000 Tmax= 9861.0000)
(Min = 1.00501855E+03 Max = 1.03692053E+03)
GET_2DFLD - solar shortwave radiation flux, t = 9100 00:00:00
(Rec=0334, Index=1, File: swrad_13_16.nc)
(Tmin= 8767.0000 Tmax= 10135.0000)
(Min = 2.21972257E-05 Max = 8.04871280E-05)
GET_2DFLD - net longwave radiation flux, t = 9100 00:00:00
(Rec=0334, Index=1, File: lwrad_13_16.nc)
(Tmin= 8767.0000 Tmax= 10135.0000)
(Min = -3.62977720E-05 Max = -5.91071844E-06)
GET_2DFLD - surface air temperature, t = 9100 00:00:00
(Rec=0334, Index=1, File: airtemp_13_16.nc)
(Tmin= 8767.0000 Tmax= 10163.0000)
(Min = -1.46754036E+01 Max = 3.10997448E+01)
GET_2DFLD - surface air relative humidity, t = 9100 00:00:00
(Rec=0334, Index=1, File: sh_13_16.nc)
(Tmin= 8767.0000 Tmax= 10163.0000)
(Min = 7.69979262E-01 Max = 2.10090205E+01)
GET_2DFLD - rain fall rate, t = 9100 00:00:00
(Rec=0334, Index=1, File: precip_13_16.nc)
(Tmin= 8767.0000 Tmax= 10116.0000)
(Min = 0.00000000E+00 Max = 7.17369363E-02)
GET_NGFLD - free-surface eastern boundary condition, t = 315 00:00:00
(Rec=0011, Index=2, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 2.94756093E-01 Max = 3.34678628E+00)
Maximum grid stiffness ratios: rx0 = 2.168438E-01 (Beckmann and Haidvogel)
rx1 = 5.784310E+00 (Haney)
Initial basin volumes: TotVolume = 1.3803256035E+17 m3
MinVolume = 1.1592918981E+08 m3
MaxVolume = 7.8719566710E+10 m3
Max/Min = 6.7903145738E+02
NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 05242069 - 05469984)
GET_NGFLD - river runoff mass transport, t = 105 00:00:00
(Rec=0004, Index=1, File: roms_river_ir_gb_KG_MH_p1.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.67757143E+03 Max = 0.00000000E+00)
GET_NGFLD - river runoff salinity, t = 105 00:00:00
(Rec=0004, Index=1, File: roms_river_ir_gb_KG_MH_p1.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 1.00000000E+01)
GET_2DFLD - surface u-wind component, t = 9101 00:00:00
(Rec=0335, Index=2, File: uwind_13_16.nc)
(Min = 2.11732128E-05 Max = 7.93587780E-05)
GET_2DFLD - net longwave radiation flux, t = 9101 00:00:00
(Rec=0335, Index=2, File: lwrad_13_16.nc)
(Tmin= 8767.0000 Tmax= 10135.0000)
(Min = -3.77141866E-05 Max = -6.27035041E-06)
GET_2DFLD - surface air temperature, t = 9101 00:00:00
(Rec=0335, Index=2, File: airtemp_13_16.nc)
(Tmin= 8767.0000 Tmax= 10163.0000)
(Min = -1.70001659E+01 Max = 3.14050999E+01)
GET_2DFLD - surface air relative humidity, t = 9101 00:00:00
(Rec=0335, Index=2, File: sh_13_16.nc)
(Tmin= 8767.0000 Tmax= 10163.0000)
(Min = 8.45432281E-01 Max = 2.17764881E+01)
GET_2DFLD - rain fall rate, t = 9101 00:00:00
(Rec=0335, Index=2, File: precip_13_16.nc)
(Tmin= 8767.0000 Tmax= 10116.0000)
(Min = 0.00000000E+00 Max = 6.48545953E-02)
GET_NGFLD - free-surface eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 3.29263715E-01 Max = 3.32140055E+00)
GET_NGFLD - free-surface southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -1.94705785E-01 Max = 3.32140055E+00)
GET_NGFLD - 2D u-momentum eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -7.76264698E-02 Max = 1.65414579E-02)
GET_NGFLD - 2D v-momentum eastern boundary condition, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 2.42844598E-02)
GET_NGFLD - 2D u-momentum southern boundary condition, t = 345 00:00:00
(Rec=0012, Index=1, File: roms_bry.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = -5.21521725E-02 Max = 1.47307941E-03)
STEP Day HH:MM:SS KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME
C => (i,j,k) Cu Cv Cw Max Speed
5242068 9100 19:30:00 3.581106E-03 2.044979E+04 2.044980E+04 1.402491E+17
(648,306,40) 6.217847E-02 7.010756E-03 0.000000E+00 4.654828E+00
DEF_HIS - creating history file: output/ocean_his_ia_0300.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
5242069 9100 19:32:30 3.581186E-03 2.044979E+04 2.044980E+04 1.402491E+17
(683,311,30) 1.618281E-02 1.994954E-03 3.296907E-01 4.653886E+00
5242070 9100 19:35:00 3.581256E-03 2.044979E+04 2.044980E+04 1.402491E+17
(683,311,30) 1.618812E-02 2.021088E-03 3.296978E-01 4.652932E+00
5242071 9100 19:37:30 3.581327E-03 2.044979E+04 2.044980E+04 1.402491E+17
(683,311,30) 1.619315E-02 2.046804E-03 3.296941E-01 4.651968E+00
5242072 9100 19:40:00 3.581398E-03 2.044979E+04 2.044980E+04 1.402491E+17
(683,311,30) 1.619791E-02 2.072099E-03 3.296829E-01 4.650996E+00
5260024 9131 23:40:00 3.583920E-03 2.044923E+04 2.044923E+04 1.402475E+17
(432,335,33) 2.274356E-02 5.484634E-03 6.402293E-02 2.689458E+00
5260025 9131 23:42:30 3.583625E-03 2.044923E+04 2.044923E+04 1.402475E+17
(432,335,33) 2.273146E-02 5.488724E-03 6.405025E-02 2.688659E+00
5260026 9131 23:45:00 3.583322E-03 2.044923E+04 2.044923E+04 1.402475E+17
(432,335,33) 2.271922E-02 5.492786E-03 6.407983E-02 2.688070E+00
5260027 9131 23:47:30 3.583001E-03 2.044923E+04 2.044923E+04 1.402475E+17
(432,335,33) 2.270847E-02 5.496920E-03 6.411384E-02 2.686619E+00
5260028 9131 23:50:00 3.582654E-03 2.044923E+04 2.044923E+04 1.402475E+17
(432,335,33) 2.269777E-02 5.501373E-03 6.414321E-02 2.685188E+00
5260029 9131 23:52:30 3.582242E-03 2.044923E+04 2.044923E+04 1.402475E+17
(432,335,33) 2.272749E-02 5.569635E-03 6.417803E-02 2.683646E+00
5260030 9131 23:55:00 3.736959E-03 2.142642E+04 2.142642E+04 1.402475E+17
(131,077,40) 9.263582E-01 1.145389E+00 3.560067E-05 1.003550E+02
Elapsed CPU time (seconds):
Node # 0 CPU: 3956.663
Node # 2 CPU: 4295.943
Node # 3 CPU: 4281.912
Node # 4 CPU: 4292.011
Node #276 CPU: 4296.153
Node #204 CPU: 4296.482
Node # 12 CPU: 4296.407
Node # 84 CPU: 4296.447
Node #180 CPU: 4296.449
Node #120 CPU: 4296.624
Node #192 CPU: 4296.667
Node # 96 CPU: 4296.658
Node #240 CPU: 4296.354
Node # 60 CPU: 4296.464
Node #167 CPU: 4296.958
Total: 1236933.528
Nonlinear model elapsed time profile:
Allocation and array initialization .............. 14.062 ( 0.0011 %)
Ocean state initialization ....................... 4.144 ( 0.0003 %)
Reading of input data ............................ 8148.930 ( 0.6588 %)
Processing of input data ......................... 20431.687 ( 1.6518 %)
Processing of output time averaged data .......... 20016.105 ( 1.6182 %)
Computation of vertical boundary conditions ...... 987.961 ( 0.0799 %)
Computation of global information integrals ...... 22406.928 ( 1.8115 %)
Writing of output data ........................... 136156.966 (11.0076 %)
Model 2D kernel .................................. 213780.327 (17.2831 %)
2D/3D coupling, vertical metrics ................. 38901.712 ( 3.1450 %)
Omega vertical velocity .......................... 28658.835 ( 2.3169 %)
Equation of state for seawater ................... 57297.323 ( 4.6322 %)
Atmosphere-Ocean bulk flux parameterization ...... 31039.049 ( 2.5094 %)
KPP vertical mixing parameterization ............. 235358.994 (19.0276 %)
3D equations right-side terms .................... 36176.144 ( 2.9247 %)
3D equations predictor step ...................... 115755.457 ( 9.3583 %)
Pressure gradient ................................ 36970.898 ( 2.9889 %)
Harmonic mixing of tracers, S-surfaces ........... 22090.072 ( 1.7859 %)
Biharmonic mixing of tracers, geopotentials ...... 19596.064 ( 1.5842 %)
Harmonic stress tensor, S-surfaces ............... 29481.635 ( 2.3834 %)
Biharmonic stress tensor, S-surfaces ............. 24038.354 ( 1.9434 %)
Corrector time-step for 3D momentum .............. 65137.180 ( 5.2660 %)
Corrector time-step for tracers .................. 50523.094 ( 4.0845 %)
Total: 1212971.920 98.0628
Nonlinear model message Passage profile:
Message Passage: 2D halo exchanges ............... 103689.973 ( 8.3828 %)
Message Passage: 3D halo exchanges ............... 92052.591 ( 7.4420 %)
Message Passage: 4D halo exchanges ............... 39345.462 ( 3.1809 %)
Message Passage: data broadcast .................. 131819.603 (10.6570 %)
Message Passage: data reduction .................. 45003.488 ( 3.6383 %)
Message Passage: data gathering .................. 9571.233 ( 0.7738 %)
Message Passage: data scattering.................. 2831.388 ( 0.2289 %)
Total: 424313.737 34.3037
All percentages are with respect to total time = 1236933.528
ROMS/TOMS - Output NetCDF summary for Grid 01:
number of time records written in HISTORY file = 00000001
number of time records written in AVERAGE file = 00000001
Analytical header files used:
ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_hmixcoef.h
ROMS/Functionals/ana_nudgcoef.h
ROMS/TOMS - Blows up ................ exit_flag: 1
MAIN: Abnormal termination: BLOWUP.
Anybody help me..
MAIN: Abnormal termination: BLOWUP.
Re: MAIN: Abnormal termination: BLOWUP.
It looks like your v-velocity went nuts from one timestep to the next, perhaps at the point 131,77,40. Have you looked at the blowup record in the restart file to see what's what there? It's at the surface - what do the surface fluxes look like there? Where is this point in relation to the land mask or river sources or??? We don't know.5260029 9131 23:52:30 3.582242E-03 2.044923E+04 2.044923E+04 1.402475E+17
(432,335,33) 2.272749E-02 5.569635E-03 6.417803E-02 2.683646E+00
5260030 9131 23:55:00 3.736959E-03 2.142642E+04 2.142642E+04 1.402475E+17
(131,077,40) 9.263582E-01 1.145389E+00 3.560067E-05 1.003550E+02
Actually, you do have rivers and you provide transport and salinity. What about the river temperature? Better to provide some value to avoid the downwind advection instability.