Dataset for "Boundary-integral solutions for freshwater lens volume in arbitrarily shaped islands"

Boundary-integral solutions for freshwater lens volume in arbitrarily shaped islands Summary This dataset provides Python code and example input files for implementing boundary-integral solutions to calculate freshwater lens volume in arbitrarily shaped islands. Description of files bis.pyMain Python script for calculating freshwater lens volume using the boundary-integral formulation.The boundary-point input file and the main model parameters can be specified in the script. For example: # CSV file containing boundary points boundary_file = 'gan.csv' # Baseline parameter values ---------------- K_base = 200 # Hydraulic conductivity (m/d) W_base = 1.86e-3 # Recharge rate (m/d) H0_base = 15 # Aquifer thickness (m) ep_base = 0.3 # Porosity (-) # Density ratio: delta = rho_f / (rho_s - rho_f) delta = 40 # Monte Carlo settings ---------------- # Total number of samples N_samples = 1000 # Standard deviations K_sigma = 0.1 * K_base W_sigma = 0.05 * W_base H0_sigma = 0.075 * H0_base ep_sigma = 0.1 * ep_base gan.csvExample boundary-point file for Gan Island with N=40 boundary points. velidhoo.csvCorresponding example file for Velidhoo Island.The boundary-point files should be provided in CSV format and contain the x and y coordinates of the island boundary points. The example files use the WGS 84 / UTM Zone 43N projected coordinate system. Other projected coordinate systems may also be used, provided that all coordinates are expressed consistently in a Cartesian coordinate system with length units such as meters.