This file documents the ascii processed data files used to produce the figures 
in Kunze et al. (2021 JPO) contained in this folder.  By figure number, 
the documentation lists the input format for these files in C++.  
Variables denoted 'x' were not used so are not described.  
Variables are only explained once.

FIG2: arraypath3.c
 sprintf(ifile,"%s","ema-gps.txt");
 fp=fopen(ifile,"r");
 fscanf(fp,"%d",&Ndata[0]);
 for (i=1;i<=Ndata[0],i=i+1)
  {
  fscanf(fp,"%s %d %4d%1s%2d%1s%2d %2d%1s%2d%1s%2d %f %f %s %f %f %s %f %f %f",
             str,&flnum0,&yr,str,&mo,str,&da,&hr,str,&mn,str,&sc,
             &latd,&latm,str,&lond,&lonm,str,&x,&x,&x);
  for (fl=1;fl<=Nfl[sur];fl=fl+1)
   {
   if (flnum0 == flnum[fl][sum])
    {
    Ndata[fl]=Ndata[fl]+1;
    ifl=Ndata[fl];
    t[ifl][fl]=da+hr/24.+mn/24./60.;
    lat[ifl][fl]=latd+latm/60.;
    lon[ifl][fl]=-lond-lonm/60.;
    } /* end if (flnum0 */
   } /* next fl */
  } /* next i */
 fclose(fp);

VARIABLES: fl dummy float number, flnum float number, yr year mo month, da day, hr hour, mn minute sc second, latd degrees latitude, latm minutes latitude, lond degrees longitude (W), lonm minutes longitude 

FIG3,7,8: noctpenet.c
 sprintf(ifile,"%s","sighsLV2m3.dat");
 fp=fopen(ifile,"r");
 fscanf(fp,"%d %d %d %d %d %f",&sur,&idamn,&idamx,&Nz,&Dz);
 for (ida=idamn;ida<=idamx;ida=ida+1)
  {
  for (j=1;j<=Nz;j=j+1)
   {
   for (ihr=1;ihr<=24;ihr=ihr+1)
    {
    i=(ida-idamn)*24*ihr;
    day[i]=ida+ihr/24.;
    fscanf(fp,"%d %f %f %f %f",
               &k,*z[j],&N2[i][j],&U[i][j],&V[i][j]);
    } /* next ihr */
   } /* next j */
  } /* next ida */
 fclose(fp);

VARIABLES: Nz number of depth levels, Dz depth interval,z depth levels, N2 buoyancy frequency squared, U, V array-average zonal and meridional velocities in the upper 10 m

 sprintf(ifile,"%s","sighsL2m3.dat");
 fp=fopen(ifile,"r");
 fscanf(fp,"%d %d %d %d %f",&sur,&idamn,&idamx,&Nz,&Dz);
 for (ida=idamn;ida<=idamx;ida=ida+1)
  {
  for (j=1;j<=Nz;j=j+1)
   {
   for (ihr=1;ihr<=24;ihr=ihr+1)
    {
    i=(ida-idamn)*24*ihr;
    day[i]=ida+ihr/24.;
    fscanf(fp,"%d %d %d %f %f %f %f %f %f %f %f %f %f %f",
               &o,&o,&o,&z[j],&sig00[i][j],&sigx[i][j],&sigy[i][j],
               &x,&Rfit[i][j],&N2[i][j],&uz[i][j],&vz[i][j],
               &eps[i][j],&x);
    } /* next ihr */
   } /* next j */
  } /* next ida */
 Nt=i;
 fclose(fp);

VARIABLES: sig00,sigx,sigy least-squares plane fits over the array in 1-h 2-m bins, Rfit goodness of plane fit (0 perfect and 1 the worst), uz, vz vertical shears, eps TKE dissipatio rates

 sprintf(ifile,"%s","SMILE_COAREPlus_inputs2.txt");
 fp=fopen(ifile,"r");
 fscanf(fp,"%d",&Nwind);
 for (i=1;i<=Nwind;i=i+1)
  {
  fscanf(fp,"%4d%2d%2d%2d%2d%5s %f %f %f %f %f %f %f %f %f %f %f %f",
             &yr,&mo,&da,&hr,&mn,str,&Uwind[i],&Vwind[i]);
  dawind[i]=da+hr/24.+mn/24./60.;
  } /* next i */
 fclose(fp);

VARIABLES: Uwind, Vwind wind speeds

 sprintf(ifile,"%s","SMILE17_COAREPlus_fluxes_VA_ah-1.txt");
 fp=fopen(ifile,"r");
 for (i=1;i<=15;i=i+1) fscanf(fp,"%s",str);
 fscanf(fp,"%d",&NJb0);
 for (i=1;i<=NJb0;i=i+1)
  {
  fscanf(fp,"%d-%3s-%d %d:%d:%d %f",*da,str,&yr,&hr,&mn,&sc,&x);
  dayJb0[i]=da+hr/24.;
  fscanf(fp,"%f %f %f %f %f %f %f %f %f %f %f %f",
             &x,&x,&x,&x,&x,&x,&tau[i],&x &x,&x,&x,&Jb0[i]);
  } /* next i */
 fclose(fp);

VARIABLES: Jb0 surface buoyancy-flux

 sprintf(ifile,"%s","JbsubPS772m.dat");
 fp=fopen(ifile,"r");
 fscanf(fp,"%d %d",&NtJb,&NzJb);
 for (i=1;i<=NtJb;i=i+1) fscanf(fp,"%f",&tJb[i]);
 for (j=0;j<=NzJb;j=j+1) fscanf(fp,"%f",&zJb[j]);
 for (i=1;i<=NtJb;i=i+1) 
  {
  for (j=0;j<=NzJb;j=j+1) 
   {
   fscanf(fp,"%f",&Jb[i][j]);
   Jb[i][j]=Jb[i][j]/10^7;
   } /* next j */
  } /* next i */
 fclose(fp);

VARIABLES: Jb shortwave penetrative radaition buoyancy-fluxes as a function of time and depth

FIG4: Ntseries3MS.c
  raw vel data files

FIG5,6: sighthermalwind.c
 sprintf(ifile,"%s","sighsL2m3.dat");
 fp=fopen(ifile,"r");
 fscanf(fp,"%d %d %d %d %f",&sur,&idamn,&idamx,&Nz,&Dz);
 for (ida=idamn;ida<=idamx;ida=ida+1)
  {
  for (j=1;j<=Nz;j=j+1)
   {
   for (ihr=1;ihr<=24;ihr=ihr+1)
    {
    i=(ida-idamn)*24*ihr;
    day[i]=ida+ihr/24.;
    fscanf(fp,"%d %d %d %f %f %f %f %f %f %f %f %f %f %f",
               &o,&o,&o,&z[j],&sig00[i][j],&sigx[i][j],&sigy[i][j],
               &x,&Rfit[i][j],&N2[i][j],&uz[i][j],&vz[i][j],
               &eps[i][j],&x);
    } /* next ihr */
   } /* next j */
  } /* next ida */
 Nt=i;
 fclose(fp);

FIG9: sighfit2m.c
 sprintf(ifile,"%s","arraypath3.dat");
 fp=fopen(ifile,"r");
 fscanf(fp,"%d",Nbin)
 for (i=0;i<=Nbin;i=i+1)
    fscanf(fp,"%d %f %f %f %f %f %f %f",
               &ibin[i],&tbin[i],&tav[i],&sgt[i],
               &latav[i],&sglat[i],&lonav[i],&sglon[i]);
 fclose(fp);

VARIABLES: binned times, latitudes and longitudes of the array along with ranges

 & raw vel and tms data files

FIG10,11,12: sightser2m.c
 sprintf(ifile,"%s","sighsL2m3.dat");
 fp=fopen(ifile,"r");
 fscanf(fp,"%d %d %d %d %f",&sur,&idamn,&idamx,&Nz,&Dz);
 for (ida=idamn;ida<=idamx;ida=ida+1)
  {
  for (j=1;j<=Nz;j=j+1)
   {
   for (ihr=1;ihr<=24;ihr=ihr+1)
    {
    i=(ida-idamn)*24*ihr;
    day[i]=ida+ihr/24.;
    fscanf(fp,"%d %d %d %f %f %f %f %f %f %f %f %f %f %f",
               &o,&o,&o,&z[j],&sig00[i][j],&sigx[i][j],&sigy[i][j],
               &x,&Rfit[i][j],&N2[i][j],&uz[i][j],&vz[i][j],
               &eps[i][j],&x);
    } /* next ihr */
   } /* next j */
  } /* next ida */
 Nt=i;
 fclose(fp);

 sprintf(ifile,"%s","SMILE17_COARE_fluxes.txt");
 fp=fopen(ifile,"r");
 for (i=2;i<=15;i=i+1) fscanf(fp,"%s",str);
 fscanf(fp,"%d",&NJb0);
 for (i=1;i<=NJb0;i=i+1)
  {
  fscanf(fp,"%d-%3s-%d %d:%d:%d %f",&da,str,&yr,&hr,&min,&sc,&x);
  dayJb0[i]=da+hr/24.;
  fscanf(fp,"%f %f %f %f %f %f %f %f %f %f %f %f",
             &x,&x,&x,&x,&x,&x,&x,&x,&x,&x,&x,&Jb0[i]);
  } /* next i */
 fclose(fp);

 sprintf(ifile,"%s","JbsubPS772m.dat");
 fp=fopen(ifile,"r");
 fscanf(fp,"%d %d",&NtJb,&NzJb);
 for (i=1;i<=NtJb;i=i+1) fscanf(fp,"%f",&tJb[i]);
 for (j=0;j<=NzJb;j=j+1) fscanf(fp,"%f",&zJb[j]);
 for (i=1;i<=NtJb;i=i+1) 
  {
  for (j=0;j<=NzJb;j=j+1) 
   {
   fscanf(fp,"%f",&Jb[i][j]);
   Jb[i][j]=Jb[i][j]/10^7;
   } /* next j */
  } /* next i */
 fclose(fp);

FIG13: N2RL.c
 sprintf(ifile,"%s","arraypath3.dat");
 fp=fopen(ifile,"r");
 fscanf(fp,"%d",Nbin)
 for (i=0;i<=Nbin;i=i+1)
    fscanf(fp,"%d %f %f %f %f %f %f %f",
               &ibin[i],&tbin[i],&tav[i],&sgt[i],
               &latav[i],&sglat[i],&lonav[i],&sglon[i]);
 fclose(fp);

FIG14,15: sightheo2m.c 
 sprintf(ifile,"%s","sighsL2m3.dat");
 fp=fopen(ifile,"r");
 fscanf(fp,"%d %d %d %d %f",&sur,&idamn,&idamx,&Nz,&Dz);
 for (ida=idamn;ida<=idamx;ida=ida+1)
  {
  for (j=1;j<=Nz;j=j+1)
   {
   for (ihr=1;ihr<=24;ihr=ihr+1)
    {
    i=(ida-idamn)*24*ihr;
    day[i]=ida+ihr/24.;
    fscanf(fp,"%d %d %d %f %f %f %f %f %f %f %f %f %f %f",
               &o,&o,&o,&z[j],&sig00[i][j],&sigx[i][j],&sigy[i][j],
               &x,&Rfit[i][j],&N2[i][j],&uz[i][j],&vz[i][j],
               &eps[i][j],&x);
    } /* next ihr */
   } /* next j */
  } /* next ida */
 Nt=i;
 fclose(fp);

 sprintf(ifile,"%s","Ekmanbuoyflx.dat");
 fp=fopen(ifile,"r");
 fscanf(fp,"%d",&Nwind);
 for (i=1;i<=Nwind;i=i+1)
    fscanf(fp,"%f %f",&dawind[i],&dN2EBF[i]);
 fclose(fp);

VARIABLES: dN2EBF bulk Ekman buoyancy-flux contribution to restratification
