pml.c

Go to the documentation of this file.

/*
 *  pml.c
 *
 *  Compute IOPs from rrs using PML IOP algorithm
 *
 *  Plymouth Marine Laboratory
 *  UK
 */
 
/*
 *  This code implements the IOP code of Smyth et al. "" Applied Optics
 */
 
#include <stdlib.h>
#include <math.h>
#include "pml.h"
#include "l12_proto.h"
 
static int iterate = 1;
static int idx410 = -1;
static int idx440 = -1;
static int idx490 = -1;
static int idx510 = -1;
static int idx555 = -1;
static int idx670 = -1;
static int initialized = 0;
 
static float *aw;
static float *bbw;
 
void pml_assert(int expr) {
    if (expr) return;
    else exit(1);
}
 
int pml_is_initialized(void) {
    return initialized;
}
 
int pml_init(int iter8, int i410, int i440, int i490, int i510, int i555, int i670, float *awptr, float *bbwptr) {
    iterate = iter8;
    idx410 = i410;
    idx440 = i440;
    idx490 = i490;
    idx510 = i510;
    idx555 = i555;
    idx670 = i670;
    aw = awptr;
    bbw = bbwptr;
    initialized = 1;
    return 0;
}