Datasets for "Europa's ocean translates interior tidal heating patterns to the ice-ocean boundary" by D. G. Lemasquerier, C. J. Bierson and K. M. Soderlund

<strong>heatflux-velocities.zip </strong> Heat flux, radial velocity, and azimuthal velocity for the simulations presented in the main text (Figure 3), i.e.<br> - 6 simulations with stress-free boundary conditions (folder "NeumannSF")<br> - 6 simulations with no slip boundary conditions (folder "NeumannNS")<br> Each folder contains the 6 simulations corresponding to 6 \(q^*\) values, \(q^*\)=[0.00,0.46,0.87,1.30,1.72,2.16]. For each simulation, we provide the data for the heat fluxes at the top and bottom boundaries, and the radial and azimuthal velocities in the ocean midplane. - "dsTdr_bottom.txt" and "dsTdr_top.txt" contain heat flux data at the bottom and top boundaries respectively. Each 2D map is an array of size Ntheta*Nphi where Ntheta is the number of grid points in latitude (-90° to 90°) and Nphi=2*Ntheta is the number of grid points in longitude (-180° to 180°). In the present files, this array is stored in column-major order. For instance, the first Ntheta lines correspond to latitudes from -90° to 90°, and a fixed longitude of -180°. The next Ntheta lines correspond to the next longitude.<br> The heat flux is given in units of the absolute value of the average bottom heat flux. - "vr_mid.txt" and "vp_mid.txt" contain radial and azimuthal velocity data at mid-height in the ocean. The ordering is the same as described above. Each 2D map is an array of size Ntheta*Nphi where Ntheta is the number of grid points in latitude (-90° to 90°) and Nphi=2*Ntheta is the number of grid points in longitude (-180° to 180°). In the present files, this array is stored in column-major order.<br> Velocities are provided in dimensionless units as a Reynolds number, i.e. they are normalized by \(\nu/D\) where \(\nu\) is the kinematic viscosity of the fluid and \(D\) the spherical shell thickness. <strong>icethickness.zip</strong> Matlab files for the ice thickness model input and outputs for four cases. The filenames refer to the simulation from which the relative heat flux variations were extracted:<br> - "NeumannSF_q0.mat": Figure 5(c), and Figure C2(c) (Neumann, q*=0.0, stress-free)<br> - "NeumannSF_q087.mat": Figure 5(f), and Figure C2(i) (Neumann, q*=0.87, stress-free)<br> - "NeumannNS_q0.mat": Figure C2(f) (Neumann, q*=0.0, no-slip)<br> - "NeumannNS_q087.mat": Figure C2(l) (Neumann, q*=0.87, no-slip) Each file contains:<br> 1) The ice thickness model inputs (Figure C1):<br> - Lambda: 100*50 array, latitude in radians.<br> - Phi: 100*50 array, longitude in radians.<br> - epsdot2avg: 100*50 array (longitude*latitude), quantity \(\overline{\dot{\epsilon}_{ij}^2}\) in \(s^{-2}\), where \(\epsilon_{ij}\) is the tidal strain within the ice. Computed from Ojakangas and Stevenson (1989), Appendix B.<br> - Fbmat: 100*50 array, oceanic heat flux in \(W/m^2\) .<br> - Ts: 100*50 array, surface temperature in K, computed from Ojakangas and Stevenson (1989), Appendix A.<br> <br> 2) The ice thickness model outputs:<br> - D: 100*50 array, ice thickness in meters.<br> - etafinal: 100*50*50 array (longitude*latitude*depth), ice viscosity in Pa s. The third dimension is the depth, which can be represented by a linearly spaced vector going from 0 to D(i,j).<br> - Hfinal: 100*50*50 array (longitude*latitude*depth), volumetric heating in the ice in W/m^3.<br> - TT: 100*50*50 array (longitude*latitude*depth), temperature in the ice in K.