disp('prt_OMECC_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\GOMECC\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 %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 \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]);