Grid and Computational TimeΒΆ

SPECIFICATION OF DIMENSION

  • Mglob: global dimension in x direction.
  • Nglob: global dimension in y direction.

PECIFICATION OF GRID SIZE

  • DX: grid size(m) in x direction.
  • DY: grid size(m) in y direction.

SPECIFICATION OF WATER DEPTH

  • DEPTH_TYPE: depth input type.

DEPTH_TYPE=DATA: from a depth file.

The program includes several simple bathymetry configurations such as

DEPTH_TYPE=FLAT: flat bottom, need DEPTH_FLAT

DEPTH_TYPE=SLOPE: plane beach along x direction. It needs three parameters: slope,SLP, slope starting point, Xslp and flat part of depth, DEPTH_FLAT

  • DEPTH_FILE: bathymetry file if DEPTH_TYPE=DATA, file dimension should be Mglob x Nglob with the first point as the south-west corner. The read format in the code is shown below.

    DO J=1,Nglob

    READ(1,*)(Depth(I,J),I=1,Mglob)

    ENDDO

  • DEPTH_FLAT: water depth of flat bottom if DEPTH_TYPE=FLAT or DEPTH_TYPE=SLOPE (flat part of a plane beach).

  • SLP: slope if DEPTH_TYPE=SLOPE

  • Xslp: starting x (m) of a slope, if DEPTH_TYPE=SLOPE

  • WaterLevel: Specify a water level which will be added to the input bathymetry and wavemaker depth such as DEP_WK (internal wave generator) and DepthWaveMaker (left boundary generator).

Note

IF you add surge or tide level using ‘WaterLevel’, please keep DEP_WK the same as the original depth in the depth file because the water level will be automatically added to the model bathymetry.

SPECIFICATION OF TIME

  • TOTAL_TIME: simulation time in seconds
  • PLOT_INTV: output interval in seconds (Note, output time is not exact because adaptive dt is used.)
  • SCREEN_INTV: time interval (s) of screen print.
  • PLOT_INTV_STATION: time interval (s) of gauge output
  • DT_fixed: time step (s) if use fixed DT. But DT_fixed will be checked by CFL. IF not satisfy CLF, DT/2, DT/4 ... will be checked until it satisfy CFL. Default is using variable DT based on CFL.