Tolosa-lct Outputs

Tolosa-lct supports several output formats. All files are written to the bin/ directory. Write intervals can be specified as formatted strings ("N seconds/minutes/hours/days") or as raw second values.

Tolosa-lct extends the Tolosa-sw output system with additional diagnostic fields specific to the non-hydrostatic solver (NH fields, wave breaking, atmospheric forcings).

Output Formats

Tecplot (w_tecplot)

7-digit flag controlling Tecplot .dat output:

Digit	Controls
D1	0=off, 1=on, 2=original format
D2	NH fields (w, p, s_x, s_y, surge tensors)
D3	Longitude/latitude fields
D4	Friction coefficients (cl, cb, z0 or Manning)
D5	Local acoustic speed a_loc
D6	Low-Mach coefficients coef_sw / coef_a (only with lm scheme)
D7	Wave breaking fields

w_tecplot   = 1100000    ! Tecplot on + NH fields
dt_tecplot  = 10 minutes

VTK (w_vtk)

8-digit flag. Note: D2 differs from Tecplot — for VTK, the surge fields (phi, psi, tensors) are written by D7, not D2.

Digit	Controls
D1	0=off, 1=on, 2=original format
D2	NH fields w, p, s_x, s_y (dispersive strain, only with DISP=yes)
D3	Longitude/latitude fields
D4	Friction coefficients (cl, cb, z0 or Manning)
D5	Local acoustic speed a_loc
D6	Low-Mach coefficients coef_sw / coef_a (only with lm scheme)
D7	Wave breaking fields — model 1: visc_break_B, visc_break_nu_t; model 2: phi, psi (or tensors for SURGE=full) + Br
D8	External forcings (p_atm, wind_u/v, tbf)

w_vtk      = 10000000
dt_vtk     = 30 minutes

Binary (w_bin)

Fortran unformatted binary output, readable by PyTolosa. 6-digit flag:

Digit	Controls
D1	1=separate file per snapshot, 2=concatenated (all timesteps in one file)
D2	NH fields w and p
D3	Longitude/latitude fields
D4	Friction coefficients
D5	Local acoustic speed a_loc
D6	Low-Mach coefficients coef_sw / coef_a

w_bin      = 100000
dt_binary  = 1 hour

Post-Treatment (w_post)

5-digit flag. When digit 1 is > 0, enables global diagnostics written at interval dt_stat / dtp:

Digit	Controls
D1	0=off, >0=on
D2	dt-weighted SSH statistics → VTK in res/statistics/ every dt_stat. Tide/mean sea level subtracted before computing variance. Outputs: ssh_mean (weighted mean) and H_m0 = 4σ (significant wave height, H_m0 = 4√m₀). D2 ≥ 2 also writes dt-weighted mean energy density maps: Epot_h_mean, Ekin_u_mean, Ekin_w_mean, Epot_p_mean.
D3	Energy time series: Ekin_h, Ekin_w, Epot_h, Epot_p, Etot (→ Tecplot)
D4	Volume, h_max, u_max time series (→ Tecplot)
D5	Hydrostatic/NH indicator field (→ Tecplot)

w_post     = 10100    ! enable post + energy + h/u time series
dtp        = 60.      ! post-treatment diagnostic interval (digits 3–5)
dt_stat    = 1 hour   ! SSH statistics export interval (digit 2)

Write Interval Parameters

Parameter	Default	Description
dt_tecplot	null	Tecplot write interval (formatted string)
dt_vtk	null	VTK write interval (formatted string)
dt_binary	null	Binary write interval (formatted string)
dt_restart	null	Restart file write interval
dt_stat	null	SSH statistics (w_post D2) VTK export interval — independent of dtp
t_min_stat	null	Simulation time at which statistics accumulation begins
t_max_stat	null	Simulation time at which statistics accumulation ends
dt_window	1e15	Statistics window duration (s); accumulator is reset after each export
dtw	-1.	Tecplot and VTK interval in seconds (fallback)
dtb	-1.	Binary interval in seconds (fallback)
dtrst	-1.	Restart interval in seconds (fallback)
dtp	-1.	Post-treatment diagnostic interval in seconds (fallback)

Restart Files

Parameter	Default	Description
restart	0	0=cold start, 1=restart from restart.bin
dt_restart	null	Write interval for restart snapshots

Verbosity

Parameter	Default	Description
verbose	1	Terminal output level: 0=quiet, higher=more detail

---

Specific Outputs — savefield.yaml

Tolosa-lct uses the same savefield.yaml mechanism as Tolosa-sw — see Tolosa-sw Outputs — savefield.yaml for the full reference (gauge, field, transect keys, location by x/y, lon/lat or cell, formats, t_min/t_max).

Tolosa-lct addition — variables: full

For gauges and fields, full is an additional variable set that includes the NH and wave fields on top of the standard shallow-water variables:

gauge:
-
  name: WG1
  x: 100.
  y: 500.
  variables: full   # ssh + u + v + bathy + w + p + phi/psi + s1/s2
  dt: 0.01

The exact columns written depend on the compilation flags:

Compiled with	Additional columns
SURGE=full	phi_xx, phi_yy, phi_xy, psi_xx, psi_yy, psi_xy
SURGE=scalarw or SURGE=scalar	phi (+ psi with scalar)
DISP=yes	s_x, s_y

full is always written in the column order: julian_cnes / ssh / u / v / bathy / w / p followed by the active wave/disp fields above.