function ep_write(E)
% EP_WRITE - Create files of eqpac data to be sent to nasa
%

% Brian Schlining
% 23 Feb 1998
ep_ini
for i = 1:length(E)
   PathS = ['\\tsunami.shore.mbari.org\oasis\eqpac\pete\realtimeprc\ftp_nasa\' E(i).id '.ftp'];
   eval(['x = ' E(i).id ';']);
   
   %======================
   % Deal with Y2K issues
   %======================
   d                = E(i).dat(:,3);
   d = d + 1900;
   %j                = find(d >= 95);
   %d(j)             = d(j) + 1900;
   %i                = find(d < 95);
   %d(j)             = d(j) + 2000;
   
   %=============
   % Format Time
   %=============
   Year = floor(d);
   
   % The decimal day (d) needs to rounded to the nearest day (Bob's program isn't perfect)
   leapYear  = ~rem(Year,4);        % Find the leap years
   Leap      = find(leapYear == 1);
   notLeap   = find(leapYear == 0); % The other years are not leap years
   JulianDay = ones(size(d))*NaN;
   
   if ~isempty(Leap)
      JulianDay(Leap) = round((d(Leap) - Year(Leap))*366) + 1;
   end
   if ~isempty(notLeap)
      JulianDay(notLeap) = round((d(notLeap) - Year(notLeap))*365) + 1;
   end
   
   JTime = E(i).dat(:,1);
   t     = find(JTime == 1230);
   JTime(t) = 1250;
   JTime = JTime/2400;
   JulianDay = JulianDay + JTime;
   sdn              = jul2date(JulianDay, Year);   % Convert to a serial date number
   [Year Month Day] = datevec(sdn);                % Convert back to year Month Day
   ymd              = Year*10000+Month*100+Day;
   
   
   
   %===================
   % Add Quality Flags
   %===================
   badFlagPrr620 = prod((E(i).dat(:,6:12) > -2)')';
   badFlagPrr600 = prod((E(i).dat(:,21:27) > -2)')';
   badFlagEd490  = E(i).dat(:,35) < -2;
   
   fr = find(strcmp(ColumnS,'Lu683/Lu555[0m]'));
   badFlagPrr620 = badFlagPrr620.*(E(i).dat(:,fr) <= 0.1) - 1; 
   
   fr = find(strcmp(ColumnS,'Lu683/Lu555[20m]'));
   badFlagPrr600 = badFlagPrr600.*(E(i).dat(:,fr) <= 0.1) - 1; 
   
   
   %===============
   % Begin Writing
   %===============
   fid = fopen(PathS,'wt');         % Open the file for writing
   
   %==================
   % Write the header
   %==================
   fprintf(fid,'/begin_header\n');
   fprintf(fid,'/affiliations=Monterey_Bay_Aquarium_Research_Institute\n');
   fprintf(fid,'/experiment=TAO_Moorings\n');
   fprintf(fid,'/cruise=TAO_98\n');
   fprintf(fid,'/investigators=Francisco_Chavez,Pete_Strutton,Brian_Schlining\n');
   fprintf(fid,'/contact=chfr@mbari.org,stpe@mbari.org,brian@mbari.org\n');
   fprintf(fid,'/parameters=Lu,Ed,Lw,Kd,CHL\n');
   fprintf(fid,'/data_type=mooring\n');
   fprintf(fid,'/original_file_name=%s\n',x.original_file_name);
   fprintf(fid,'/delimiter=space\n');
   %Old missing code: fprintf(fid,'/missing=-99.99)
   fprintf(fid,'/missing=-999\n');
   %Eaxample: fprintf(fid,'/north_latitude=-2.000[DEG]\n');
   fprintf(fid,'/north_latitude=%-5.3f[DEG]\n',x.north_latitude);
   fprintf(fid,'/south_latitude=%-5.3f[DEG]\n',x.south_latitude);
   fprintf(fid,'/east_longitude=%-5.3f[DEG]\n',x.east_longitude);
   fprintf(fid,'/west_longitude=%-5.3f[DEG]\n',x.west_longitude);
   fprintf(fid,'/start_date=%8i\n',min(ymd));
   fprintf(fid,'/end_date=%8i\n',max(ymd));
   %Example: fprintf(fid,'/start_time=12:00:00[GMT]\n');
   %fprintf(fid,'/start_time=12:%02i:00[Local]\n',(E(i).dat(1,1) - 1200));
   %fprintf(fid,'/end_time=12:%02i:00[Local]\n',(E(i).dat(end,1) - 1200));
   fprintf(fid,'/start_time=02:00:00[GMT]\n');
   fprintf(fid,'/end_time=02:30:00[GMT]\n');
   fprintf(fid,'/documents=MbariSeabassMooring.readme\n');
   fprintf(fid,'/calibration_files=%s\n',x.calibration_files);
   fprintf(fid,'! Comments: Last column is quality flag; u=unchecked, g=good\n');
   fprintf(fid,'! Comments: MorelChl[0-20m] is mean chl over the upper 20m calculated according to Morel (1988), OC2V2Chl[0m] is calculated from the surface Satlantic OCR-100 data, and Lu443/Lu555Chl at 0m and 20m calculated from an algorithm similar to Abbott and Letelier (1996), but with coefficients for the equatorial Pacific.  See documentation.\n');
   fprintf(fid,'! Comments: ELw is Lw calculated from the 1.5m Lu data, using an empirical relationship derived from Satlantic SPMR data for the equatorial Pacific, from June 1998.\n');
   fprintf(fid,'! Comments: LLw is Lw calculated from a linear extrapolation of the 20m and 1.5m data, with a factor of 0.544 applied for transmission across the water-air interface.\n');
   fprintf(fid,'! Comments: KLw is Lw calculated using the attenuation coefficient of downwelling irradiance at the relevant wavelength to back-propagate the 20m Lu data.\n');
   fprintf(fid,'! Comments: Quality Flags: -1=bad 0=unchecked, 1=good\n');
   % Write the fields...They are listed in ep_ini.m. Changes to field names or
   % the Field order need to be made there.
   fprintf(fid,'/fields=Date1,Time,Date2,');
   for k = 3:length(ColumnS)
      fprintf(fid,'%s,',char(ColumnS(k)));
   end
   fprintf(fid,'Quality[0m],Quality[20m]\n');
   
   % Write the Units
   fprintf(fid,'/units=yyyymmdd,hh:mm:ss,mm/dd/yyyy,');
   for k = 3:length(ColumnS)
      fprintf(fid,'%s,',char(UnitsS(k)));
   end
   fprintf(fid,'unitless,unitless\n');
   
   fprintf(fid,'/end_header@\n');
   
   %============
   % Write Body
   %============
   [r c] = size(E(i).dat);
   BAD = find(E(i).dat > 100000);
   E(i).dat(BAD) = NaN;
   BAD = find(isnan(E(i).dat) | isinf(E(i).dat));
   E(i).dat(BAD) = -999;
   Hour = 12;
   Minute = (E(i).dat(:,1)-1200);
   Second = 0;
   GmtDate = datenum(Year, Month, Day, Hour, Minute, Second) + gmtoffset(x.east_longitude);
   [Year Month Day Hour Minute Second] = datevec(GmtDate);
   for n = 1:r
      fprintf(fid,'%8i ',ymd(n));
      
      fprintf(fid,'%02i:%02i:00 ',Hour(n), Minute(n));         % Time in GMT
      fprintf(fid,'%02i/%02i/%4i ', Month(n), Day(n), Year(n));
      fprintf(fid,'%9.6f ', E(i).dat(n,3));                    % DecimalYear
      fprintf(fid,'%9.3f ', E(i).dat(n,4:5));                  % Temp Sal
      fprintf(fid,'%12.6f', E(i).dat(n,6:end));               % The data
      fprintf(fid,'%3i%3i\n',badFlagPrr620(n),badFlagPrr600(n)); % Quality flags
   end
   
   
   
   fclose(fid);
   pc2unix(PathS);
   
end