NASA Ocean Color

ocssw V2022

 #include <stdio.h>
 #include <stdlib.h>
 #include "l2_struc.h"
 #include "l12_parms.h"
 #include "l12_proto.h"
  
 /* --------------------------------------------------------- */
 /* init_l2() - initialize a Level-2 record                   */
  
 /* --------------------------------------------------------- */
 void init_l2(l2str *l2rec, int32_t nbands) {
     int32_t ip, ib, ipb;
     int32_t npix = l2rec->l1rec->npix;
  
     for (ip = 0; ip < npix; ip++) {
  
         l2rec->num_iter[ip] = BAD_INT;
         l2rec->aermodmin[ip] = BAD_INT;
         l2rec->aermodmax[ip] = BAD_INT;
         l2rec->aermodmin2[ip] = BAD_INT;
         l2rec->aermodmax2[ip] = BAD_INT;
  
         l2rec->chl[ip] = BAD_FLT;
         l2rec->eps[ip] = BAD_FLT;
         l2rec->chi2[ip] = BAD_FLT;
         l2rec->aerratio[ip] = BAD_FLT;
         l2rec->aerratio2[ip] = BAD_FLT;
         l2rec->aerindex[ip] = BAD_FLT;
  
         for (ib = 0; ib < nbands; ib++) {
             ipb = ip * nbands + ib;
             l2rec->taua[ipb] = BAD_FLT;
             l2rec->La[ipb] = BAD_FLT;
             l2rec->Lw[ipb] = BAD_FLT;
             l2rec->nLw[ipb] = BAD_FLT;
             if(l2rec->l1rec->uncertainty){
                 l2rec->Rrs_unc[ipb] =BAD_FLT;
             }
  
             l2rec->brdf[ipb] = BAD_FLT;
             l2rec->Rrs[ipb] = BAD_FLT;
             l2rec->outband_correction[ipb] = 1.0;
             l2rec->a[ipb] = BAD_FLT;
             l2rec->bb[ipb] = BAD_FLT;
         }
         l2rec->chl_unc[ip]=BAD_FLT;
  
         if(input->proc_uncertainty==2){
             ipb=ip*nbands*nbands;
             for(ib=0;ib<nbands*nbands;ib++)
                 l2rec->covariance_matrix[ipb+ib]=0.;
         }
  
     }
 }
  