## Loisel and Stramski algorithm application to other sensors

Use this Forum to find information on, or ask a question about, NASA Earth Science data.
tihicha
Posts: 6
Joined: Wed Jun 05, 2013 4:47 pm America/New_York

### Loisel and Stramski algorithm application to other sensors

Hello,
Regarding the Loisel and Stramski IOP code as implemented in SeaADS:
https://oceancolor.gsfc.nasa.gov/docs/ocssw/las__iop_8c_source.html

I was wondering what your approach is to apply the algorithm to sensors other than SeaWiFS, as the original code coefficients are provided for SeaWiFS only. The Loisel code as posted here (https://www.ioccg.org/groups/Software_OCA/IOP_inversion.f) is specified at 410, 440 , 490, 510 and 550 nm. Is anything done, and does anything need to be done to apply it even to SeaWiFS, which has slightly different center wavelengths, notably 555 instead of 550 nm.

I know this is asking about code that is quite old by now. Thanks very much for your time looking into it!
Tiho

Tags:

OB.DAAC - SeanBailey
User Services
Posts: 1461
Joined: Wed Sep 18, 2019 6:15 pm America/New_York
Been thanked: 4 times

### Loisel and Stramski algorithm application to other sensors

Tiho,

Generally, if code requires a specific set of wavelengths, l2gen will generalize by picking the closest wavelength (to a point).
For example:
`  302         wave = l1rec->l1file->fwave;  303         ib443 = windex(443, wave, nbands);  304         ib555 = windex(555, wave, nbands);`

The `windex` function finds the sensor wavelength closest to the requested value.

In the case of the LAS code you reference, it  gives another clue as to how it specifically handles the table:

`  161         // if the input wavelength is not one of the table wavelengths, call the   162         // function recursively for bounding table wavelengths and interpolate`

Hope this helps.
Sean