When I try to add a river input, I add to the header file: UV_PSOURCE, TS_PSOURCE ans ANA_PSOURCE. I've read previous problems.
In the ANA_PSOURCE I let the variables like this:
Nsrc=1
Dsrc=0
Isrc=1
Jsrc=6
Qshape= 1.0_r8/REAL(N(ng),r8)
Qbar=fac*1500.0_r8
Tsrc(is,k,itemp)=TO(ng)
Tsrc(is,k,isalt)=S0(ng)
As you can see, the values are the same as in the RIVERPLUME1 case.
With these values it compiles and it runs the model, but what I don't get it, is that Jsrc is supposed to be the location of j index, right? Well, I defined my grid analitically with Mm=4, so it shouldn't be more than 5 points, and when I put Jsrc=1 or 2 or 3 or 4 or 5, it compiles but blows up when I run the model. I don't understand why. Does somebody know?
Many thanks in advanced.
PROBLEM WITH ANA_PSOURCE
Re: PROBLEM WITH ANA_PSOURCE
Does ROMS give any kind of message on blowing up? That may or may not help us figure it out. Sometimes compiling with USE_DEBUG gives a little more information because the object files retain the line number information.
-
chaos4u2
Re: PROBLEM WITH ANA_PSOURCE
Thank you Kate, this is what I get with Jsrc(Nsrc)=4 for example, I've checked the USE_DEBUG and I have it on:
Operating system : Linux
CPU/hardware : i686
Compiler system : gfortran
Compiler command : /usr/bin/gfortran
Compiler flags : -frepack-arrays -g -fbounds-check -ffree-form -ffree-line-length-none -ffree-form -ffree-line-length-none -ffree-form -ffree-line-length-none
SVN Root URL : https://www.myroms.org/svn/src/trunk
SVN Revision : 575M
Local Root : /home/sergi/src
Header Dir : /home/sergi/Projects/CASPRACTIC
Header file : caspractic.h
Analytical Dir: /home/sergi/Projects/CASPRACTIC
Resolution, Grid 01: 0100x0004x010, Parallel Threads: 1, Tiling: 001x001
Physical Parameters, Grid: 01
=============================
4032 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).
288 nHIS Number of timesteps between the writing fields
into history file.
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.
1 ntsDIA Starting timestep for the accumulation of output
time-averaged diagnostics data.
576 nDIA Number of timesteps between the writing of
time-averaged data into diagnostics 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+00 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
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.
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).
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
3.0000E+00 theta_s S-coordinate surface control parameter.
0.0000E+00 theta_b S-coordinate bottom control parameter.
4.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)
5.0000E+01 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
5.0000E+01 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.
14.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).
T 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 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
T Dout(M2rate) Write out 2D momentun acceleration.
T Dout(M2pgrd) Write out 2D momentum pressure gradient.
T Dout(M2fcor) Write out 2D momentum Coriolis force.
T Dout(M2hadv) Write out 2D momentum horizontal advection.
T Dout(M2xadv) Write out 2D momentum horizontal X-advection.
T Dout(M2yadv) Write out 2D momentum horizontal Y-advection.
T Dout(M2hvis) Write out 2D momentum horizontal viscosity.
T Dout(M2xvis) Write out 2D momentum horizontal X-viscosity.
T Dout(M2yvis) Write out 2D momentum horizontal Y-viscosity.
T Dout(M2sstr) Write out 2D momentum surface stress.
T Dout(M2bstr) Write out 2D momentum bottom stress.
T Dout(M3rate) Write out 3D momentun acceleration.
T Dout(M3pgrd) Write out 3D momentum pressure gradient.
T Dout(M3fcor) Write out 3D momentum Coriolis force.
T Dout(M3hadv) Write out 3D momentum horizontal advection.
T Dout(M3xadv) Write out 3D momentum horizontal X-advection.
T Dout(M3yadv) Write out 3D momentum horizontal Y-advection.
T Dout(M3vadv) Write out 3D momentum vertical advection.
T Dout(M3hvis) Write out 3D momentum horizontal viscosity.
T Dout(M3xvis) Write out 3D momentum horizontal X-viscosity.
T Dout(M3yvis) Write out 3D momentum horizontal Y-viscosity.
T Dout(M3vvis) Write out 3D momentum vertical viscosity.
T Dout(iTrate) Write out rate of change of tracer 01: temp
T Dout(iTrate) Write out rate of change of tracer 02: salt
T Dout(iThadv) Write out horizontal advection, tracer 01: temp
T Dout(iThadv) Write out horizontal advection, tracer 02: salt
T Dout(iTxadv) Write out horizontal X-advection, tracer 01: temp
T Dout(iTxadv) Write out horizontal X-advection, tracer 02: salt
T Dout(iTyadv) Write out horizontal Y-advection, tracer 01: temp
T Dout(iTyadv) Write out horizontal Y-advection, tracer 02: salt
T Dout(iTvadv) Write out vertical advection, tracer 01: temp
T Dout(iTvadv) Write out vertical advection, tracer 02: salt
T Dout(iThdif) Write out horizontal diffusion, tracer 01: temp
T Dout(iThdif) Write out horizontal diffusion, tracer 02: salt
T Dout(iTxdif) Write out horizontal X-diffusion, tracer 01: temp
T Dout(iTxdif) Write out horizontal X-diffusion, tracer 02: salt
T Dout(iTydif) Write out horizontal Y-diffusion , tracer 01: temp
T Dout(iTydif) Write out horizontal Y-diffusion , tracer 02: salt
T Dout(iTvdif) Write out vertical diffusion, tracer 01: temp
T Dout(iTvdif) Write out vertical diffusion, tracer 02: salt
Output/Input Files:
Output Restart File: ocean_rst.nc
Output History File: ocean_his.nc
Output Averages File: ocean_avg.nc
Output Diagnostics File: ocean_dia.nc
Tile partition information for Grid 01: 0100x0004x0010 tiling: 001x001
tile Istr Iend Jstr Jend Npts
0 1 100 1 4 4000
Tile minimum and maximum fractional grid coordinates:
(interior points only)
tile Xmin Xmax Ymin Ymax grid
0 -1.50 102.50 0.50 5.50 RHO-points
0 -2.00 102.00 0.50 5.50 U-points
0 -1.50 102.50 0.00 5.00 V-points
Lateral Boundary Conditions: NLM
============================
Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------
zeta 1 Periodic Closed Periodic Closed
ubar 1 Periodic Closed Periodic Closed
vbar 1 Periodic Closed Periodic Closed
u 1 Periodic Closed Periodic Closed
v 1 Periodic Closed Periodic Closed
temp 1 Periodic Closed Periodic Closed
salt 1 Periodic Closed Periodic Closed
Activated C-preprocessing Options:
CASPRACTIC CASPRACTIC
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_GRID Analytical grid set-up.
ANA_INITIAL Analytical initial conditions.
ANA_PSOURCE Analytical point sources and sinks.
ANA_SMFLUX Analytical kinematic surface momentum flux.
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.
DIAGNOSTICS_TS Computing and writing tracer diagnostic terms.
DIAGNOSTICS_UV Computing and writing momentum diagnostic terms.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
LMD_MIXING Large/McWilliams/Doney interior mixing.
MIX_S_TS Mixing of tracers along constant S-surfaces.
MIX_S_UV Mixing of momentum along constant S-surfaces.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLVE3D Solving 3D Primitive Equations.
SPLINES Conservative parabolic spline reconstruction.
TS_U3HADVECTION Third-order upstream horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection of tracers.
TS_DIF2 Harmonic mixing of tracers.
TS_PSOURCE Tracers point sources and sinks.
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_QDRAG Quadratic bottom stress.
UV_PSOURCE Mass point sources and sinks.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
Process Information:
Thread # 0 (pid= 21489) is active.
INITIAL: Configuring and initializing forward nonlinear model ...
Vertical S-coordinate System:
level S-coord Cs-curve Z at hmin at hc half way at hmax
10 0.0000000 0.0000000 0.000 0.000 0.000 0.000
9 -0.1000000 -0.0050000 -0.321 -0.210 -0.321 -0.321
8 -0.2000000 -0.0204535 -0.718 -0.441 -0.718 -0.718
7 -0.3000000 -0.0477616 -1.198 -0.696 -1.198 -1.198
6 -0.4000000 -0.0894007 -1.781 -0.979 -1.781 -1.781
5 -0.5000000 -0.1491465 -2.494 -1.298 -2.494 -2.494
4 -0.6000000 -0.2324164 -3.374 -1.665 -3.374 -3.374
3 -0.7000000 -0.3467612 -4.477 -2.094 -4.477 -4.477
2 -0.8000000 -0.5025493 -5.875 -2.605 -5.875 -5.875
1 -0.9000000 -0.7139071 -7.671 -3.228 -7.671 -7.671
0 -1.0000000 -1.0000000 -10.000 -4.000 -10.000 -10.000
Time Splitting Weights: 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.0107112868023601 0.2343461315221332
8-0.0005032751608632 0.0336903762267453-0.0112145619632232 0.2680365077488785
9 0.0003887272597151 0.0337071520654408-0.0108258347035082 0.3017436598143193
10 0.0014892209965583 0.0336941944901169-0.0093366137069498 0.3354378543044362
11 0.0027955815694920 0.0336445537902317-0.0065410321374578 0.3690824080946679
12 0.0043042707117221 0.0335513677379153-0.0022367614257357 0.4026337758325832
13 0.0060106451121704 0.0334078920475245 0.0037738836864347 0.4360416678801077
14 0.0079087469427945 0.0332075372104522 0.0116826306292293 0.4692492050905599
15 0.0099910761708920 0.0329439123123590 0.0216737068001212 0.5021931174029189
16 0.0122483446563884 0.0326108764399960 0.0339220514565097 0.5348039938429149
17 0.0146692120341107 0.0322025982847830 0.0485912634906204 0.5670065921276979
18 0.0172400033810439 0.0317136245503127 0.0658312668716643 0.5987202166780106
19 0.0199444086685725 0.0311389577709445 0.0857756755402368 0.6298591744489551
20 0.0227631639997064 0.0304741441486588 0.1085388395399432 0.6603333185976139
21 0.0256737146312911 0.0297153720153352 0.1342125541712343 0.6900486906129492
22 0.0286498597812017 0.0288595815276255 0.1628624139524359 0.7189082721405747
23 0.0316613792205220 0.0279045862015855 0.1945237931729579 0.7468128583421602
24 0.0346736416507075 0.0268492068942347 0.2291974348236654 0.7736620652363949
25 0.0376471948657328 0.0256934188392112 0.2668446296893983 0.7993554840756061
26 0.0405373376992233 0.0244385123436867 0.3073819673886216 0.8237939964192929
27 0.0432936737565711 0.0230872677537126 0.3506756411451927 0.8468812641730055
28 0.0458596469320356 0.0216441452951603 0.3965352880772283 0.8685254094681657
29 0.0481720587108285 0.0201154903974257 0.4447073467880568 0.8886408998655915
30 0.0501605672561820 0.0185097551070648 0.4948679140442388 0.9071506549726562
31 0.0517471682814031 0.0168377361985254 0.5466150823256418 0.9239883911711816
32 0.0528456577069106 0.0151128305891453 0.5994607400325525 0.9391012217603268
33 0.0533610761022577 0.0133513086655816 0.6528218161348102 0.9524525304259084
34 0.0531891349131380 0.0115726061288397 0.7060109510479481 0.9640251365547481
35 0.0522156244733761 0.0097996349650684 0.7582265755213241 0.9738247715198165
36 0.0503158038019030 0.0080591141492892 0.8085423793232271 0.9818838856691057
37 0.0473537721847153 0.0063819206892258 0.8558961515079425 0.9882658063583314
38 0.0431818225418188 0.0048034616164019 0.8990779740497613 0.9930692679747333
39 0.0376397765791563 0.0033640675316746 0.9367177506289176 0.9964333355064079
40 0.0305543017255205 0.0021094083123694 0.9672720523544381 0.9985427438187773
41 0.0217382098544503 0.0010909315881854 0.9890102622088884 0.9996336754069627
42 0.0109897377911117 0.0003663245930371 1.0000000000000000 0.9999999999999998
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
Minimum X-grid spacing, DXmin = 1.00000000E-01 km
Maximum X-grid spacing, DXmax = 1.00000000E-01 km
Minimum Y-grid spacing, DYmin = 2.50000000E-01 km
Maximum Y-grid spacing, DYmax = 2.50000000E-01 km
Minimum Z-grid spacing, DZmin = 3.21428732E-01 m
Maximum Z-grid spacing, DZmax = 2.32923473E+00 m
Minimum barotropic Courant Number = 1.06675208E+00
Maximum barotropic Courant Number = 1.06675208E+00
Maximum Coriolis Courant Number = 3.00037500E-02
Maximum grid stiffness ratios: rx0 = 0.000000E+00 (Beckmann and Haidvogel)
rx1 = 0.000000E+00 (Haney)
Initial basin volumes: TotVolume = 1.0000000000E+08 m3
MinVolume = 8.0357183061E+03 m3
MaxVolume = 5.8230868312E+04 m3
Max/Min = 7.2465044311E+00
NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00004032)
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 8.750000E-01 1.236683E+02 1.245433E+02 2.000000E+08
(001,2,10) 3.000000E+00 1.200000E+00 0.000000E+00 1.414214E+00
DEF_HIS - creating history file: ocean_his.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
DEF_AVG - creating average file: ocean_avg.nc
DEF_DIAGS - creating diagnostics file: ocean_dia.nc
1 0 00:05:00 NaN NaN NaN NaN
(000,0,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
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):
Thread # 0 CPU: 1.372
Total: 1.372
Nonlinear model elapsed time profile:
Allocation and array initialization .............. 0.064 ( 4.6647 %)
Ocean state initialization ....................... 0.016 ( 1.1662 %)
Computation of vertical boundary conditions ...... 0.004 ( 0.2915 %)
Computation of global information integrals ...... 0.016 ( 1.1662 %)
Writing of output data ........................... 0.024 ( 1.7492 %)
Model 2D kernel .................................. 0.976 (71.1371 %)
2D/3D coupling, vertical metrics ................. 0.020 ( 1.4577 %)
Omega vertical velocity .......................... 0.004 ( 0.2915 %)
Equation of state for seawater ................... 0.040 ( 2.9155 %)
KPP vertical mixing parameterization ............. 0.004 ( 0.2916 %)
3D equations right-side terms .................... 0.008 ( 0.5831 %)
3D equations predictor step ...................... 0.012 ( 0.8747 %)
Pressure gradient ................................ 0.004 ( 0.2915 %)
Harmonic mixing of tracers, S-surfaces ........... 0.004 ( 0.2915 %)
Harmonic stress tensor, S-surfaces ............... 0.004 ( 0.2915 %)
Corrector time-step for 3D momentum .............. 0.096 ( 6.9971 %)
Corrector time-step for tracers .................. 0.076 ( 5.5394 %)
Total: 1.372 100.0001
All percentages are with respect to total time = 1.372
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_grid.h
ROMS/Functionals/ana_initial.h
ROMS/Functionals/ana_nudgcoef.h
ROMS/Functionals/ana_psource.h
ROMS/Functionals/ana_smflux.h
ROMS/Functionals/ana_stflux.h
Operating system : Linux
CPU/hardware : i686
Compiler system : gfortran
Compiler command : /usr/bin/gfortran
Compiler flags : -frepack-arrays -g -fbounds-check -ffree-form -ffree-line-length-none -ffree-form -ffree-line-length-none -ffree-form -ffree-line-length-none
SVN Root URL : https://www.myroms.org/svn/src/trunk
SVN Revision : 575M
Local Root : /home/sergi/src
Header Dir : /home/sergi/Projects/CASPRACTIC
Header file : caspractic.h
Analytical Dir: /home/sergi/Projects/CASPRACTIC
Resolution, Grid 01: 0100x0004x010, Parallel Threads: 1, Tiling: 001x001
Physical Parameters, Grid: 01
=============================
4032 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).
288 nHIS Number of timesteps between the writing fields
into history file.
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.
1 ntsDIA Starting timestep for the accumulation of output
time-averaged diagnostics data.
576 nDIA Number of timesteps between the writing of
time-averaged data into diagnostics 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+00 nl_visc2 NLM Horizontal, harmonic mixing coefficient
(m2/s) for momentum.
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.
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).
2 Vtransform S-coordinate transformation equation.
4 Vstretching S-coordinate stretching function.
3.0000E+00 theta_s S-coordinate surface control parameter.
0.0000E+00 theta_b S-coordinate bottom control parameter.
4.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)
5.0000E+01 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
5.0000E+01 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.
14.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).
T 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 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
T Dout(M2rate) Write out 2D momentun acceleration.
T Dout(M2pgrd) Write out 2D momentum pressure gradient.
T Dout(M2fcor) Write out 2D momentum Coriolis force.
T Dout(M2hadv) Write out 2D momentum horizontal advection.
T Dout(M2xadv) Write out 2D momentum horizontal X-advection.
T Dout(M2yadv) Write out 2D momentum horizontal Y-advection.
T Dout(M2hvis) Write out 2D momentum horizontal viscosity.
T Dout(M2xvis) Write out 2D momentum horizontal X-viscosity.
T Dout(M2yvis) Write out 2D momentum horizontal Y-viscosity.
T Dout(M2sstr) Write out 2D momentum surface stress.
T Dout(M2bstr) Write out 2D momentum bottom stress.
T Dout(M3rate) Write out 3D momentun acceleration.
T Dout(M3pgrd) Write out 3D momentum pressure gradient.
T Dout(M3fcor) Write out 3D momentum Coriolis force.
T Dout(M3hadv) Write out 3D momentum horizontal advection.
T Dout(M3xadv) Write out 3D momentum horizontal X-advection.
T Dout(M3yadv) Write out 3D momentum horizontal Y-advection.
T Dout(M3vadv) Write out 3D momentum vertical advection.
T Dout(M3hvis) Write out 3D momentum horizontal viscosity.
T Dout(M3xvis) Write out 3D momentum horizontal X-viscosity.
T Dout(M3yvis) Write out 3D momentum horizontal Y-viscosity.
T Dout(M3vvis) Write out 3D momentum vertical viscosity.
T Dout(iTrate) Write out rate of change of tracer 01: temp
T Dout(iTrate) Write out rate of change of tracer 02: salt
T Dout(iThadv) Write out horizontal advection, tracer 01: temp
T Dout(iThadv) Write out horizontal advection, tracer 02: salt
T Dout(iTxadv) Write out horizontal X-advection, tracer 01: temp
T Dout(iTxadv) Write out horizontal X-advection, tracer 02: salt
T Dout(iTyadv) Write out horizontal Y-advection, tracer 01: temp
T Dout(iTyadv) Write out horizontal Y-advection, tracer 02: salt
T Dout(iTvadv) Write out vertical advection, tracer 01: temp
T Dout(iTvadv) Write out vertical advection, tracer 02: salt
T Dout(iThdif) Write out horizontal diffusion, tracer 01: temp
T Dout(iThdif) Write out horizontal diffusion, tracer 02: salt
T Dout(iTxdif) Write out horizontal X-diffusion, tracer 01: temp
T Dout(iTxdif) Write out horizontal X-diffusion, tracer 02: salt
T Dout(iTydif) Write out horizontal Y-diffusion , tracer 01: temp
T Dout(iTydif) Write out horizontal Y-diffusion , tracer 02: salt
T Dout(iTvdif) Write out vertical diffusion, tracer 01: temp
T Dout(iTvdif) Write out vertical diffusion, tracer 02: salt
Output/Input Files:
Output Restart File: ocean_rst.nc
Output History File: ocean_his.nc
Output Averages File: ocean_avg.nc
Output Diagnostics File: ocean_dia.nc
Tile partition information for Grid 01: 0100x0004x0010 tiling: 001x001
tile Istr Iend Jstr Jend Npts
0 1 100 1 4 4000
Tile minimum and maximum fractional grid coordinates:
(interior points only)
tile Xmin Xmax Ymin Ymax grid
0 -1.50 102.50 0.50 5.50 RHO-points
0 -2.00 102.00 0.50 5.50 U-points
0 -1.50 102.50 0.00 5.00 V-points
Lateral Boundary Conditions: NLM
============================
Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------
zeta 1 Periodic Closed Periodic Closed
ubar 1 Periodic Closed Periodic Closed
vbar 1 Periodic Closed Periodic Closed
u 1 Periodic Closed Periodic Closed
v 1 Periodic Closed Periodic Closed
temp 1 Periodic Closed Periodic Closed
salt 1 Periodic Closed Periodic Closed
Activated C-preprocessing Options:
CASPRACTIC CASPRACTIC
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_GRID Analytical grid set-up.
ANA_INITIAL Analytical initial conditions.
ANA_PSOURCE Analytical point sources and sinks.
ANA_SMFLUX Analytical kinematic surface momentum flux.
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.
DIAGNOSTICS_TS Computing and writing tracer diagnostic terms.
DIAGNOSTICS_UV Computing and writing momentum diagnostic terms.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
LMD_MIXING Large/McWilliams/Doney interior mixing.
MIX_S_TS Mixing of tracers along constant S-surfaces.
MIX_S_UV Mixing of momentum along constant S-surfaces.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLVE3D Solving 3D Primitive Equations.
SPLINES Conservative parabolic spline reconstruction.
TS_U3HADVECTION Third-order upstream horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection of tracers.
TS_DIF2 Harmonic mixing of tracers.
TS_PSOURCE Tracers point sources and sinks.
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_QDRAG Quadratic bottom stress.
UV_PSOURCE Mass point sources and sinks.
UV_VIS2 Harmonic mixing of momentum.
VAR_RHO_2D Variable density barotropic mode.
Process Information:
Thread # 0 (pid= 21489) is active.
INITIAL: Configuring and initializing forward nonlinear model ...
Vertical S-coordinate System:
level S-coord Cs-curve Z at hmin at hc half way at hmax
10 0.0000000 0.0000000 0.000 0.000 0.000 0.000
9 -0.1000000 -0.0050000 -0.321 -0.210 -0.321 -0.321
8 -0.2000000 -0.0204535 -0.718 -0.441 -0.718 -0.718
7 -0.3000000 -0.0477616 -1.198 -0.696 -1.198 -1.198
6 -0.4000000 -0.0894007 -1.781 -0.979 -1.781 -1.781
5 -0.5000000 -0.1491465 -2.494 -1.298 -2.494 -2.494
4 -0.6000000 -0.2324164 -3.374 -1.665 -3.374 -3.374
3 -0.7000000 -0.3467612 -4.477 -2.094 -4.477 -4.477
2 -0.8000000 -0.5025493 -5.875 -2.605 -5.875 -5.875
1 -0.9000000 -0.7139071 -7.671 -3.228 -7.671 -7.671
0 -1.0000000 -1.0000000 -10.000 -4.000 -10.000 -10.000
Time Splitting Weights: 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.0107112868023601 0.2343461315221332
8-0.0005032751608632 0.0336903762267453-0.0112145619632232 0.2680365077488785
9 0.0003887272597151 0.0337071520654408-0.0108258347035082 0.3017436598143193
10 0.0014892209965583 0.0336941944901169-0.0093366137069498 0.3354378543044362
11 0.0027955815694920 0.0336445537902317-0.0065410321374578 0.3690824080946679
12 0.0043042707117221 0.0335513677379153-0.0022367614257357 0.4026337758325832
13 0.0060106451121704 0.0334078920475245 0.0037738836864347 0.4360416678801077
14 0.0079087469427945 0.0332075372104522 0.0116826306292293 0.4692492050905599
15 0.0099910761708920 0.0329439123123590 0.0216737068001212 0.5021931174029189
16 0.0122483446563884 0.0326108764399960 0.0339220514565097 0.5348039938429149
17 0.0146692120341107 0.0322025982847830 0.0485912634906204 0.5670065921276979
18 0.0172400033810439 0.0317136245503127 0.0658312668716643 0.5987202166780106
19 0.0199444086685725 0.0311389577709445 0.0857756755402368 0.6298591744489551
20 0.0227631639997064 0.0304741441486588 0.1085388395399432 0.6603333185976139
21 0.0256737146312911 0.0297153720153352 0.1342125541712343 0.6900486906129492
22 0.0286498597812017 0.0288595815276255 0.1628624139524359 0.7189082721405747
23 0.0316613792205220 0.0279045862015855 0.1945237931729579 0.7468128583421602
24 0.0346736416507075 0.0268492068942347 0.2291974348236654 0.7736620652363949
25 0.0376471948657328 0.0256934188392112 0.2668446296893983 0.7993554840756061
26 0.0405373376992233 0.0244385123436867 0.3073819673886216 0.8237939964192929
27 0.0432936737565711 0.0230872677537126 0.3506756411451927 0.8468812641730055
28 0.0458596469320356 0.0216441452951603 0.3965352880772283 0.8685254094681657
29 0.0481720587108285 0.0201154903974257 0.4447073467880568 0.8886408998655915
30 0.0501605672561820 0.0185097551070648 0.4948679140442388 0.9071506549726562
31 0.0517471682814031 0.0168377361985254 0.5466150823256418 0.9239883911711816
32 0.0528456577069106 0.0151128305891453 0.5994607400325525 0.9391012217603268
33 0.0533610761022577 0.0133513086655816 0.6528218161348102 0.9524525304259084
34 0.0531891349131380 0.0115726061288397 0.7060109510479481 0.9640251365547481
35 0.0522156244733761 0.0097996349650684 0.7582265755213241 0.9738247715198165
36 0.0503158038019030 0.0080591141492892 0.8085423793232271 0.9818838856691057
37 0.0473537721847153 0.0063819206892258 0.8558961515079425 0.9882658063583314
38 0.0431818225418188 0.0048034616164019 0.8990779740497613 0.9930692679747333
39 0.0376397765791563 0.0033640675316746 0.9367177506289176 0.9964333355064079
40 0.0305543017255205 0.0021094083123694 0.9672720523544381 0.9985427438187773
41 0.0217382098544503 0.0010909315881854 0.9890102622088884 0.9996336754069627
42 0.0109897377911117 0.0003663245930371 1.0000000000000000 0.9999999999999998
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
Minimum X-grid spacing, DXmin = 1.00000000E-01 km
Maximum X-grid spacing, DXmax = 1.00000000E-01 km
Minimum Y-grid spacing, DYmin = 2.50000000E-01 km
Maximum Y-grid spacing, DYmax = 2.50000000E-01 km
Minimum Z-grid spacing, DZmin = 3.21428732E-01 m
Maximum Z-grid spacing, DZmax = 2.32923473E+00 m
Minimum barotropic Courant Number = 1.06675208E+00
Maximum barotropic Courant Number = 1.06675208E+00
Maximum Coriolis Courant Number = 3.00037500E-02
Maximum grid stiffness ratios: rx0 = 0.000000E+00 (Beckmann and Haidvogel)
rx1 = 0.000000E+00 (Haney)
Initial basin volumes: TotVolume = 1.0000000000E+08 m3
MinVolume = 8.0357183061E+03 m3
MaxVolume = 5.8230868312E+04 m3
Max/Min = 7.2465044311E+00
NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00004032)
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 8.750000E-01 1.236683E+02 1.245433E+02 2.000000E+08
(001,2,10) 3.000000E+00 1.200000E+00 0.000000E+00 1.414214E+00
DEF_HIS - creating history file: ocean_his.nc
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000001
DEF_AVG - creating average file: ocean_avg.nc
DEF_DIAGS - creating diagnostics file: ocean_dia.nc
1 0 00:05:00 NaN NaN NaN NaN
(000,0,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
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):
Thread # 0 CPU: 1.372
Total: 1.372
Nonlinear model elapsed time profile:
Allocation and array initialization .............. 0.064 ( 4.6647 %)
Ocean state initialization ....................... 0.016 ( 1.1662 %)
Computation of vertical boundary conditions ...... 0.004 ( 0.2915 %)
Computation of global information integrals ...... 0.016 ( 1.1662 %)
Writing of output data ........................... 0.024 ( 1.7492 %)
Model 2D kernel .................................. 0.976 (71.1371 %)
2D/3D coupling, vertical metrics ................. 0.020 ( 1.4577 %)
Omega vertical velocity .......................... 0.004 ( 0.2915 %)
Equation of state for seawater ................... 0.040 ( 2.9155 %)
KPP vertical mixing parameterization ............. 0.004 ( 0.2916 %)
3D equations right-side terms .................... 0.008 ( 0.5831 %)
3D equations predictor step ...................... 0.012 ( 0.8747 %)
Pressure gradient ................................ 0.004 ( 0.2915 %)
Harmonic mixing of tracers, S-surfaces ........... 0.004 ( 0.2915 %)
Harmonic stress tensor, S-surfaces ............... 0.004 ( 0.2915 %)
Corrector time-step for 3D momentum .............. 0.096 ( 6.9971 %)
Corrector time-step for tracers .................. 0.076 ( 5.5394 %)
Total: 1.372 100.0001
All percentages are with respect to total time = 1.372
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_grid.h
ROMS/Functionals/ana_initial.h
ROMS/Functionals/ana_nudgcoef.h
ROMS/Functionals/ana_psource.h
ROMS/Functionals/ana_smflux.h
ROMS/Functionals/ana_stflux.h
-
arango
- Site Admin
- Posts: 1367
- Joined: Wed Feb 26, 2003 4:41 pm
- Location: DMCS, Rutgers University
- Contact:
Re: PROBLEM WITH ANA_PSOURCE
Do not use ANA_PSOURCE. This option is extremely difficult to code in parallel. It requires a lot of expertise. I don't think that it is a good idea for users learning about ROMS to start coding the ANA_PSOURCE option. This option is very deceiving. It may look trivial but it is actually very difficult even for me. Build a NetCDF file instead. I have mentioned this several times in this forum 
Re: PROBLEM WITH ANA_PSOURCE
The cause of your blowup could be here:
You want these numbers to all be less than one. Try a shorter timestep - like say half as long.
Code: Select all
Minimum barotropic Courant Number = 1.06675208E+00
Maximum barotropic Courant Number = 1.06675208E+00
Maximum Coriolis Courant Number = 3.00037500E-02