disp('prt_AMMA_flux_hr');
fclose('all');
clear jazz;
%This program takes data computed by the jas, nar, or kwaj da_red2.m progam and creates an ascii
%file appropriate for reading into a spreadsheet (e.g., lotus, quatropro, etc)
%To use it, you must first run the analysis program, or stick a call to thisin it

f=['C:\Data\AMMA\flux_data\' jdollars 'flux_hr.txt'];%jdollars is the 99 field prog name

flist=fopen(f,'w');
%flist=1;

npz=length(jda);
jazz(1,1:npz)=jda';						%date
jazz(2,1:npz)=yushp;						%uship
jazz(3,1:npz)=v1;							%u true
jazz(4,1:npz)=thet1;						%true dir
jazz(5,1:npz)=v1rel;						%u rel
jazz(6,1:npz)=thet1rel;					%rel dir
jazz(7,1:npz)=yhead;						%ship heading
jazz(8,1:npz)=yts';						%Ts
jazz(9,1:npz)=y1(:,19)';				%T
jazz(10,1:npz)=yqs';						%qs
jazz(11,1:npz)=y1(:,21)';				%qa	
jazz(12,1:npz)=yhsc;						%hs covariance
jazz(13,1:npz)=yhsib;					%hs ID
jazz(14,1:npz)=yhsb;						%hs bulk						
jazz(15,1:npz)=yhlc_bl;					%hl covariance
jazz(16,1:npz)=yhlib_lic;					%hl ID
jazz(17,1:npz)=yhlb;						%hl bulk
jazz(18,1:npz)=ytaucx;					%stress streamwise covariance
jazz(19,1:npz)=ytaucy;					%   "     xstream
jazz(20,1:npz)=ytauib;					%stress ID
jazz(21,1:npz)=ytaub;					%stress bulk
jazz(22,1:npz)=y1(:,22)';				%solar flux
jazz(23,1:npz)=y1(:,23)';				%IR flux
jazz(24,1:npz)=y1(:,24)';				%ORG rain rate
jazz(25,1:npz)=yj;						%ship plume contam index
%jazz(26,1:npz)=min(999,y1(:,133)');	%sigoph of good data
jazz(26,1:npz)=y1(:,44)';				%tilt of good data
jazz(27,1:npz)=yjm;						%maneuver index
jazz(28,1:npz)=y1(:,139)';				%ct
jazz(29,1:npz)=y1(:,141)';				%co
jazz(30,1:npz)=y1(:,135)';				%cu
jazz(31,1:npz)=y1(:,137)';				%cw
jazz(32,1:npz)=yhrain;					%Rain heat flux
jazz(33,1:npz)=hl_webar';				%Hl webb flux
jazz(34,1:npz)=y1(:,151)';				%lat
jazz(35,1:npz)=y1(:,150)';				%lon
jazz(36,1:npz)=y1(:,27)';				%zu
jazz(37,1:npz)=y1(:,28)';				%zt
jazz(38,1:npz)=y1(:,29)';				%zq
jazz(39,1:npz)=ysog; 					%uship  scs
jazz(40,1:npz)=v1_scs;					%u true	scs
jazz(41,1:npz)=thet1_scs;				%true dir scs
jazz(42,1:npz)=ycog;						%ship course scs 
jazz(43,1:npz)=y1(:,4)';				%Ts scs
jazz(44,1:npz)=y1(:,5)';				%T scs
yqs_scs=qsea(y1(:,4));
jazz(45,1:npz)=yqs_scs';				%qs scs
jazz(46,1:npz)=y1(:,7)';				%qa scs 
jazz(47,1:npz)=y1(:,8)';				%rs scs	 
jazz(48,1:npz)=y1(:,9)';				%rl scs
%jazz(50,1:npz)=hl_licbar;				%hl licor
jazz(49,1:npz)=wco2_licbar';				%wcCo2 licor
%jazz(52,1:npz)=hl_ophbar;             %Cov. Latent Heat Flux from Ophir
jazz(50,1:npz)=y3(:,2)';               %Specific Humidity from Licor (g/kg)
%jazz(54,1:npz)=q_ophbar;               %Specific Humidity from Ophir (g/kg)
%jazz(55,1:npz)=y1(:,134)';               %Median of Ophir Clear Channel (counts/s)
%jazz(56,1:npz)=y1(:,133)';               %Stv of Ophir Clear Channel (counts/s)
jazz(51,1:npz)=y3(:,15)';              %Stv of Specific Humidity from Licor (g/kg)
%jazz(58,1:npz)=y1(:,131)';              %Stv of Specific Humidity from Ophir (g/kg)
%jazz(59,1:npz)=yhlib_lic;					%hl ID
jazz(52,1:npz)=y1(:,124)';                    %variance of w' ((m/s)^2)
jazz(53,1:npz)=y1(:,121)';                % sonic w'T' covariance before the hunidity contribution
jazz(54,1:npz)=y5(:,2)';               % lagtime between w and O3 calculated by crosscovriance, 10 min lagtime
jazz(55,1:npz)=y5(:,5)';              % w'O3' calculated with crosscovariance, 10 min lagtime, no webb correction
jazz(56,1:npz)=y5(:,6)';              % w'O3' calculated with crosscovariance, flux at hourly mean lagtime, no webb correction
jazz(57,1:npz)=y5(:,7)';              % w'O3' calculated with crosscovariance, flux at dayly mean lagtime, no webb correction
jazz(58,1:npz)=y5(:,3)';           % mean O3 concentration cross (ppb)
jazz(59,1:npz)=y5(:,4)';            % stdev ozone cross (ppb)
jazz(60,1:npz)=y5(:,10)';                % O3'w' with constant lagtime of 5.1 s, and function cov, no webb correction
jazz(61,1:npz)=y5(:,8)';              % mean O3 concentration constant lagtime (ppb)
jazz(62,1:npz)=y5(:,9)';              % stdev ozone constant lagtime (ppb)
jazz(63,1:npz)=y5(:,5)'./y5(:,3)'*100;% deposition velocity 10 min lagtime (cm/s)
jazz(64,1:npz)=y5(:,6)'./y5(:,3)'*100;% deposition velocity hourly lagtime (cm/s)
jazz(65,1:npz)=y5(:,7)'./y5(:,3)'*100;% deposition velocity dayly lagtime (cm/s)
jazz(66,1:npz)=y5(:,10)'./y5(:,8)'*100;    % deposition velocity 5.1s lagtime (cm/s)
jazz(67,1:npz)=y5(:,13)';              % w'O3' calculated with crosscovariance, 10 min lagtime, webb correction
jazz(68,1:npz)=y5(:,14)';              % w'O3' calculated with crosscovariance, flux at hourly mean lagtime, webb correction
jazz(69,1:npz)=y5(:,15)';              % w'O3' calculated with crosscovariance, flux at dayly mean lagtime, webb correction
jazz(70,1:npz)=y5(:,16)';                 % O3'w' with constant lagtime of 5.1 s, and function cov, webb correction
jazz(71,1:npz)=y5(:,13)'./y5(:,3)'*100;% deposition velocity 10 min lagtime (cm/s)
jazz(72,1:npz)=y5(:,14)'./y5(:,3)'*100;% deposition velocity hourly lagtime (cm/s)
jazz(73,1:npz)=y5(:,14)'./y5(:,3)'*100;% deposition velocity dayly lagtime (cm/s)
jazz(74,1:npz)=y5(:,16)'./y5(:,8)'*100;     % deposition velocity 5.1s lagtime (cm/s)
jazz(75,1:npz)=y5(:,11)';
jazz(76,1:npz)=y5(:,12)';
jazz(77,1:npz)=y1(:,52)';               %density of dry air
jazz(78,1:npz)=y1(:,60)';               %ustr
%vectorized print, 76 columns
%                 1     2     3     4     5     6     7      8     9     10    11    12    13    14    15    16    17     18    19    20    21    22     23    24    25    26     27    28    29    30      31      32     33    34      35     36     37   38     39    40    41     42     43    44    45    46    47   48    49     50     51     52     53    54    55      56     57    58       59    60
fprintf(flist,'%11.5f %8.3f %8.3f %8.2f %8.3f %8.2f %8.2f  %8.3f %8.3f %8.3f %8.3f %8.2f %8.2f %8.2f %9.2f %9.2f %9.2f  %8.4f %8.4f %8.4f %8.4f %9.2f %9.2f %8.2f %8.1f %9.2f %8.3f %12.3e %12.3e %12.3e %12.3e  %9.2f %8.2f  %12.5f %12.5f  %9.2f %8.2f %9.2f %8.3f %8.3f  %9.2f  %9.2f %8.3f %8.3f %8.3f %8.3f %8.1f %8.1f %11.5f %11.5f %11.5f %11.5f %10.5f %10.5f %10.5f %10.5f %10.5f %10.5f  %10.5f %10.5f %10.5f %10.5f %10.5f %10.5f %10.5f %10.5f %10.5f %10.5f  %10.5f %10.5f %10.5f %10.5f %10.5f %10.5f %10.5f %10.5f %10.5f %10.5f \r\n',jazz);
%                jd    ushp utrue dirtr urel  dirrl  head  tsnk  tair   qs     qa    hsc   hsib  hsb   hlc   hlib  hlb   taucx taucy tauib  taub  psp   pir   org   J     tilt  jm    ct     cq      cu      cw     RF   webb    lat     lon   zu    zt    zq    sog utruim dirtrim  cog   tsg  taim  qsg   qaim  rsim  rlim  wclic   qlic   sgq_lic ww 
fclose('all');
disp(['File written as ' f]);	
