IMVL

The data package comprises Fortran code (Fortran code.zip), resulting NetCDF data files, and Matlab code (Matlab code.zip). Executing the Fortran code generates NetCDF data files containing model outputs such as temperature, salinity, currents, grid bathymetry, and sea surface height. The Matlab code is used to generate the figures presented in the paper. The Fortran code (Fortran code.zip) contains 10 subdirectories. Each directory is named according to the value of the sill steepness parameter A and the inclusion (or exclusion) of the Coriolis force. For example: · The directory A=0-0cf contains the Fortran code for simulations with the sill steepness parameter A set to 0 and the Coriolis force set to zero (representing a strait geometry without Coriolis). · The directory A=0.8 contains the Fortran code for simulations with A set to 0.8 and the Coriolis force set to 10^-4 s⁻¹ (representing a sill geometry with Coriolis). Each directory contains the following files: 1. The original POM model source code (pom2k.f). 2. The modified POM source code (pom2kba.f) configured with the specific sill steepness and Coriolis force parameters corresponding to the directory name. 3. The enhanced POM source code (pomdsba.f) incorporating the variable-limit integration method for solving temperature and salinity. Within each directory: 1. Executing the command make ba compiles and runs pom2kba.f, generating the corresponding NetCDF data file (pom2kba.nc) for the parameter settings (sill steepness A and Coriolis force) indicated by the directory name. 2. Executing the command make dsba compiles and runs pomdsba.f, generating the corresponding NetCDF data file (pomdsba.nc) for the parameter settings indicated by the directory name. The NetCDF output files are named with the prefix pom2kba (standard POM) or pomdsba (variable-limit integration enhanced POM), followed by identifiers specifying the sill steepness parameter A and Coriolis force setting. For instance: · The file pom2kba-ibm-4-0cf.nc contains results from the standard POM model (pom2kba.f) run with A=0.4 and Coriolis force set to zero (sill geometry, no Coriolis). · The file pomdsba-imvl-2.nc contains results from the enhanced POM model (pomdsba.f) incorporating variable-limit integration, run with A=0.2 and Coriolis force set to 10^-4 s⁻¹ (sill geometry, with Coriolis). The Matlab code (Matlab code.zip) contains all scripts necessary to reproduce the figures in the paper. Scripts are named with the prefix f followed by a figure number. Executing a specific script, for example F11.m, will generate the corresponding Figure 11 from the paper.