MorelMaritorena model
Instructions for use, includes example: text file
- The code computes above-water water reflectance Rw, defined as Rw=pi*Lw/Es where Lw is water-leaving radiance and Es downward solar irradiance just above the surface, as a function of wavelength (wav, nm), chlorophyll concentration (chl, mg/m^3), backscattering coefficient (bbs, m^-1), solar zenith angle (sza, degree), view zenith angle (vza, degree), relative azimuth angle (raa, degree), and wind speed (ws, m/s). bbs is introduced to account for backscattering by particles other than phytoplankton.
- Solar zenith angle and view zenith angle should be in the range (0,90). The f/Q table used in this code is built for solar zenith angle less than 75 degree, and view zenith angle less than 50 deg. Therefore, extension to higher solar / view zenith angles is based on nearest neighbor extrapolation.
- Relative azimuth angle should be in the range (0,180). raa=0: the sensor views along the retroreflecting direction; raa=180: sensor views along the mirror-reflecting direction.
- The model is valid for the range of wavelengths 350-2400 nm.
- Example python script to run the program: run_morel.py
Files needed to run the program:
- linux operating system, installation of python, and the python extensions NumPy and Cython.
- Software Files needed (Python,Cython,C):
- Data Files needed:
- Program needed to compile routine with cython: setup.py
Example:
- wav = [412., 443., 490., 510., 565., 670., 865., 1020.]
- chl for 0.1, 1.0, 10.0, bbs=0, wind speed = 5, sza=30, vza=20, raa=90
- Plot: Rw vs. wav
