gpcureton / sunglint-polyspectra

Simulates the elevation and slope of the ocean surface with a particular power spectrum, and applies the specular condition to the slope to calculate the surface sunglint. Correlations are added between various elevation phase components.

Sunglint Cumulant Function Simulation

This program simulates multiple elevation and slope realisations with given spectra, introducing quadratic phase coupling at selected wavenumbers. Glint realisations are computed from the slope realisations. The average elevation, slope and glint moments and cumulants, and moment and cumulant functions are computed. The results are written to a HDF5 file.

Dependencies

• pytables
• h5py
• numpy
• scipy

Usage:

usage: cumulantFunctionSimulate.py [-h] [-n N] [-N NN] [-d DELTA_X]
                                   [-r NUMREALS] [-S SPECTRUMTYPE] [-l]
                                   [-o OUTFILE] [-v]

optional arguments:
  -h, --help            show this help message and exit
  -n N, --numdatapoints N
                        Number of points of the glint dataset. Must be an
                        integer power of 2. [default: 1024]
  -N NN, --num2Dpoints NN
                        Number of points of the bispectrum array side. Must be
                        an integer power of 2. [default: 64]
  -d DELTA_X, --deltax DELTA_X
                        The spatial increment in meters. [default: 0.02]
  -r NUMREALS, --num_realisations NUMREALS
                        Number of realisations. [default: 100]
  -S SPECTRUMTYPE, --spectrum_type SPECTRUMTYPE
                        Form of the elevation power spectrum. Possible values
                        are... 'phillips_3', 'phillips_4', 'gaussian'.
                        [default: 'phillips_3']
  -l, --nonlinear_modes
                        Switch on nonlinear modes. [default: False]
  -o OUTFILE, --output_file OUTFILE
                        The full path of the output HDF5 file. [default:
                        'outSimulatedGlint.h5']
  -v, --verbose         each occurrence increases verbosity 1 level from
                        ERROR: -v=WARNING -vv=INFO -vvv=DEBUG