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 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221
| clc, clear
%% Dossier où sont stockés les fichiers
folder = 'D:\BARBEAU\BARBEAU-EDY\LI7200-HS50\CALC\_SMARTFLUX\'
%% Déclaration du fichier de sortie
fidout = fopen('D:\BARBEAU\BARBEAU-EDY\LI7200-HS50\CALC\_SMARTFLUX\Li7200Fluxes_EP-Express_Smartflux.dat', 'wt');
%% Création de l'entête Campbell Sci dans le fichier de sortie
%fprintf(fidout,'%s\n','"TOA5","EP-SmartFlux","CR3000","5488","CR3000.Std.24","CPU:Tower1_1.1_CR3","43371","EP-SmartFlux"');
%fprintf(fidout,'%s\n','"TIMESTAMP","RECORD","DATAH","filename","date","time","DOY","daytime","file_records","used_records","Tau","qc_Tau","rand_err_Tau","H","qc_H","rand_err_H","LE","qc_LE","rand_err_LE","co2_flux","qc_co2_flux","rand_err_co2_flux","h2o_flux","qc_h2o_flux","rand_err_h2o_flux","ch4_flux","qc_ch4_flux","rand_err_ch4_flux","none_flux","qc_none_flux","rand_err_none_flux","H_strg","LE_strg","co2_strg","h2o_strg","ch4_strg","none_strg","co2_v-adv","h2o_v-adv","ch4_v-adv","none_v-adv","co2_molar_density","co2_mole_fraction","co2_mixing_ratio","co2_time_lag","co2_def_timelag","h2o_molar_density","h2o_mole_fraction","h2o_mixing_ratio","h2o_time_lag","h2o_def_timelag","ch4_molar_density","ch4_mole_fraction","ch4_mixing_ratio","ch4_time_lag","ch4_def_timelag","none_molar_density","none_mole_fraction","none_mixing_ratio","none_time_lag","none_def_timelag","sonic_temperature","air_temperature","air_pressure","air_density","air_heat_capacity","air_molar_volume","ET","water_vapor_density","e","es","specific_humidity","RH","VPD","Tdew","u_unrot","v_unrot","w_unrot","u_rot","v_rot","w_rot","wind_speed","max_wind_speed","wind_dir","yaw","pitch","roll","u*","TKE","L","(z-d)/L","bowen_ratio","T*","model","x_peak","x_offset","x_10%","x_30%","x_50%","x_70%","x_90%","un_Tau","Tau_scf","un_H","H_scf","un_LE","LE_scf","un_co2_flux","co2_scf","un_h2o_flux","h2o_scf","un_ch4_flux","ch4_scf","un_none_flux","un_none_scf","spikes_hf","amplitude_resolution_hf","drop_out_hf","absolute_limits_hf","skewness_kurtosis_hf","skewness_kurtosis_sf","discontinuities_hf","discontinuities_sf","timelag_hf","timelag_sf","attack_angle_hf","non_steady_wind_hf","u_spikes","v_spikes","w_spikes","ts_spikes","co2_spikes","h2o_spikes","ch4_spikes","none_spikes","head_detect_LI-7200","t_out_LI-7200","t_in_LI-7200","aux_in_LI-7200","delta_p_LI-7200","chopper_LI-7200","detector_LI-7200","pll_LI-7200","sync_LI-7200","chopper_LI-7500","detector_LI-7500","pll_LI-7500","sync_LI-7500","not_ready_LI-7700","no_signal_LI-7700","re_unlocked_LI-7700","bad_temp_LI-7700","laser_temp_unregulated_LI-7700","block_temp_unregulated_LI-7700","motor_spinning_LI-7700","pump_on_LI-7700","top_heater_on_LI-7700","bottom_heater_on_LI-7700","calibrating_LI-7700","motor_failure_LI-7700","bad_aux_tc1_LI-7700","bad_aux_tc2_LI-7700","bad_aux_tc3_LI-7700","box_connected_LI-7700","mean_value_RSSI_LI-7200","mean_value_LI-7500","u_var","v_var","w_var","ts_var","co2_var","h2o_var","ch4_var","none_var","w/ts_cov","w/co2_cov","w/h2o_cov","w/ch4_cov","w/none_cov","co2_mean","h2o_mean","vin_sf_mean","dew_point_mean","avg_signal_strength_7200_mean","co2_signal_strength_7200_mean","h2o_signal_strength_7200_mean","delta_signal_strength_7200_mean","flowrate_mean","hit_power_mean","hit_vin_mean","CHK"');
%fprintf(fidout,'%s\n','"TS","RN","DATAU","","[yyyy-mm-dd]","[HH:MM]","[ddd.ddd]","[1=daytime]","[#]","[#]","[kg+1m-1s-2]","[#]","[kg+1m-1s-2]","[W+1m-2]","[#]","[W+1m-2]","[W+1m-2]","[#]","[W+1m-2]","[µmol+1s-1m-2]","[#]","[µmol+1s-1m-2]","[mmol+1s-1m-2]","[#]","[mmol+1s-1m-2]","[µmol+1s-1m-2]","[#]","[µmol+1s-1m-2]","[µmol+1s-1m-2]","[#]","[µmol+1s-1m-2]","[W+1m-2]","[W+1m-2]","[µmol+1s-1m-2]","[mmol+1s-1m-2]","[µmol+1s-1m-2]","[µmol+1s-1m-2]","[µmol+1s-1m-2]","[mmol+1s-1m-2]","[µmol+1s-1m-2]","[µmol+1s-1m-2]","[mmol+1m-3]","[µmol+1mol_a-1]","[µmol+1mol_d-1]","[s]","[1=default]","[mmol+1m-3]","[mmol+1mol_a-1]","[mmol+1mol_d-1]","[s]","[1=default]","[mmol+1m-3]","[µmol+1mol_a-1]","[µmol+1mol_d-1]","[s]","[1=default]","[mmol+1m-3]","[µmol+1mol_a-1]","[µmol+1mol_d-1]","[s]","[1=default]","[K]","[K]","[Pa]","[kg+1m-3]","[J+1kg-1K-1]","[m+3mol-1]","[mm]","[kg+1m-3]","[Pa]","[Pa]","[kg+1kg-1]","[%]","[Pa]","[K]","[m+1s-1]","[m+1s-1]","[m+1s-1]","[m+1s-1]","[m+1s-1]","[m+1s-1]","[m+1s-1]","[m+1s-1]","[deg_from_north]","[deg]","[deg]","[deg]","[m+1s-1]","[m+2s-2]","[m]","[#]","[#]","[K]","[0=KJ/1=KM/2=HS]","[m]","[m]","[m]","[m]","[m]","[m]","[m]","[kg+1m-1s-2]","[#]","[W+1m-2]","[#]","[W+1m-2]","[#]","[µmol+1s-1m-2]","[#]","[mmol+1s-1m-2]","[#]","[µmol+1s-1m-2]","[#]","[µmol+1s-1m-2]","[#]","8u/v/w/ts/co2/h2o/ch4/none","8u/v/w/ts/co2/h2o/ch4/none","8u/v/w/ts/co2/h2o/ch4/none","8u/v/w/ts/co2/h2o/ch4/none","8u/v/w/ts/co2/h2o/ch4/none","8u/v/w/ts/co2/h2o/ch4/none","8u/v/w/ts/co2/h2o/ch4/none","8u/v/w/ts/co2/h2o/ch4/none","8co2/h2o/ch4/none","8co2/h2o/ch4/none","8aa","8U","[#]","[#]","[#]","[#]","[#]","[#]","[#]","[#]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#_flagged_recs]","[#]","[#]","[m+2s-2]","[m+2s-2]","[m+2s-2]","[K+2]","--","--","--","--","[m+1s-1K+1]","--","--","--","--","--","--","--","--","--","--","--","--","--","--","--","026"');
fprintf(fidout,'%s\n','"TOA5","EP-SmartFlux","CR3000","5488","CR3000.Std.24","CPU:Tower1_1.1_CR3","43371","EP-SmartFlux"');
fprintf(fidout,'%s\n','"TIMESTAMP","RECORD","date" ,"time" ,"daytime" ,"Tau" ,"qc_Tau" ,"H" ,"qc_H" ,"LE" ,"qc_LE" ,"co2_flux" ,"qc_co2_flux","h2o_flux" ,"qc_h2o_flux"');
fprintf(fidout,'%s\n','"TS" ,"RN" ,"[yyyy-mm-dd]","[HH:MM]","[1=daytime]","[kg+1m-1s-2]","[#]" ,"[kg+1m-1s-2]","[#]" ,"[W+1m-2]","[#]" ,"[µmol+1s-1m-2]","[#]" ,"[mmol+1s-1m-2]","[#]" ');
%% Travail sur la liste des fichiers stockés dans le dossier
filePattern = fullfile(folder, '*.txt');
textFiles = dir(filePattern)
delimiter = '\t';
startRow = 3;
formatSpec = '%s%s%s%s%f%f%f%f%f%f%s%f%f%s%f%f%s%f%f%s%f%f%s%s%s%s%s%s%s%s%s%s%s%s%s%f%f%s%s%f%f%f%f%f%f%f%f%f%f%s%s%s%s%f%s%s%s%s%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%s%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%s%[^\n\r]';
for k = 1:1 %length(textFiles) % Boucle de lecture des fichiers dans leur ordre de création
%% Lecture nom du fichier et création de la table
baseFileName = textFiles(k).name;
fullFileName = fullfile(folder, textFiles(k).name);
fprintf('Now reading %s\n', fullFileName);
fileID = fopen(fullFileName, 'r');
dataArray = textscan(fileID, formatSpec, 'Delimiter', delimiter, 'EmptyValue' ,NaN,'HeaderLines' ,startRow-1, 'ReturnOnError', false);
%% Allocate imported array to column variable names
DATAH = dataArray{:, 1};
filename = dataArray{:, 2};
date = dataArray{:, 3};
time = dataArray{:, 4};
DOY = dataArray{:, 5};
daytime = dataArray{:, 6};
file_records = dataArray{:, 7};
used_records = dataArray{:, 8};
Tau = dataArray{:, 9};
qc_Tau = dataArray{:, 10};
rand_err_Tau = dataArray{:, 11};
H = dataArray{:, 12};
qc_H = dataArray{:, 13};
rand_err_H = dataArray{:, 14};
LE = dataArray{:, 15};
qc_LE = dataArray{:, 16};
rand_err_LE = dataArray{:, 17};
co2_flux = dataArray{:, 18};
qc_co2_flux = dataArray{:, 19};
rand_err_co2_flux = dataArray{:, 20};
h2o_flux = dataArray{:, 21};
qc_h2o_flux = dataArray{:, 22};
rand_err_h2o_flux = dataArray{:, 23};
ch4_flux = dataArray{:, 24};
qc_ch4_flux = dataArray{:, 25};
rand_err_ch4_flux = dataArray{:, 26};
none_flux = dataArray{:, 27};
qc_none_flux = dataArray{:, 28};
rand_err_none_flux = dataArray{:, 29};
H_strg = dataArray{:, 30};
LE_strg = dataArray{:, 31};
co2_strg = dataArray{:, 32};
h2o_strg = dataArray{:, 33};
ch4_strg = dataArray{:, 34};
none_strg = dataArray{:, 35};
co2_vadv = dataArray{:, 36};
h2o_vadv = dataArray{:, 37};
ch4_vadv = dataArray{:, 38};
none_vadv = dataArray{:, 39};
co2_molar_density = dataArray{:, 40};
co2_mole_fraction = dataArray{:, 41};
co2_mixing_ratio = dataArray{:, 42};
co2_time_lag = dataArray{:, 43};
co2_def_timelag = dataArray{:, 44};
h2o_molar_density = dataArray{:, 45};
h2o_mole_fraction = dataArray{:, 46};
h2o_mixing_ratio = dataArray{:, 47};
h2o_time_lag = dataArray{:, 48};
h2o_def_timelag = dataArray{:, 49};
ch4_molar_density = dataArray{:, 50};
ch4_mole_fraction = dataArray{:, 51};
ch4_mixing_ratio = dataArray{:, 52};
ch4_time_lag = dataArray{:, 53};
ch4_def_timelag = dataArray{:, 54};
none_molar_density = dataArray{:, 55};
none_mole_fraction = dataArray{:, 56};
none_mixing_ratio = dataArray{:, 57};
none_time_lag = dataArray{:, 58};
none_def_timelag = dataArray{:, 59};
sonic_temperature = dataArray{:, 60};
air_temperature = dataArray{:, 61};
air_pressure = dataArray{:, 62};
air_density = dataArray{:, 63};
air_heat_capacity = dataArray{:, 64};
air_molar_volume = dataArray{:, 65};
ET = dataArray{:, 66};
water_vapor_density = dataArray{:, 67};
e = dataArray{:, 68};
es = dataArray{:, 69};
specific_humidity = dataArray{:, 70};
RH = dataArray{:, 71};
VPD = dataArray{:, 72};
Tdew = dataArray{:, 73};
u_unrot = dataArray{:, 74};
v_unrot = dataArray{:, 75};
w_unrot = dataArray{:, 76};
u_rot = dataArray{:, 77};
v_rot = dataArray{:, 78};
w_rot = dataArray{:, 79};
wind_speed = dataArray{:, 80};
max_wind_speed = dataArray{:, 81};
wind_dir = dataArray{:, 82};
yaw = dataArray{:, 83};
pitch = dataArray{:, 84};
roll = dataArray{:, 85};
ustar = dataArray{:, 86};
TKE = dataArray{:, 87};
L = dataArray{:, 88};
VarName89 = dataArray{:, 89};
bowen_ratio = dataArray{:, 90};
Tstar = dataArray{:, 91};
model = dataArray{:, 92};
x_peak = dataArray{:, 93};
x_offset = dataArray{:, 94};
x_10 = dataArray{:, 95};
x_30 = dataArray{:, 96};
x_50 = dataArray{:, 97};
x_70 = dataArray{:, 98};
x_90 = dataArray{:, 99};
un_Tau = dataArray{:, 100};
Tau_scf = dataArray{:, 101};
un_H = dataArray{:, 102};
H_scf = dataArray{:, 103};
un_LE = dataArray{:, 104};
LE_scf = dataArray{:, 105};
un_co2_flux = dataArray{:, 106};
co2_scf = dataArray{:, 107};
un_h2o_flux = dataArray{:, 108};
h2o_scf = dataArray{:, 109};
un_ch4_flux = dataArray{:, 110};
ch4_scf = dataArray{:, 111};
un_none_flux = dataArray{:, 112};
un_none_scf = dataArray{:, 113};
spikes_hf = dataArray{:, 114};
amplitude_resolution_hf = dataArray{:, 115};
drop_out_hf = dataArray{:, 116};
absolute_limits_hf = dataArray{:, 117};
skewness_kurtosis_hf = dataArray{:, 118};
skewness_kurtosis_sf = dataArray{:, 119};
discontinuities_hf = dataArray{:, 120};
discontinuities_sf = dataArray{:, 121};
timelag_hf = dataArray{:, 122};
timelag_sf = dataArray{:, 123};
attack_angle_hf = dataArray{:, 124};
non_steady_wind_hf = dataArray{:, 125};
u_spikes = dataArray{:, 126};
v_spikes = dataArray{:, 127};
w_spikes = dataArray{:, 128};
ts_spikes = dataArray{:, 129};
co2_spikes = dataArray{:, 130};
h2o_spikes = dataArray{:, 131};
ch4_spikes = dataArray{:, 132};
none_spikes = dataArray{:, 133};
head_detect_LI7200 = dataArray{:, 134};
t_out_LI7200 = dataArray{:, 135};
t_in_LI7200 = dataArray{:, 136};
aux_in_LI7200 = dataArray{:, 137};
delta_p_LI7200 = dataArray{:, 138};
chopper_LI7200 = dataArray{:, 139};
detector_LI7200 = dataArray{:, 140};
pll_LI7200 = dataArray{:, 141};
sync_LI7200 = dataArray{:, 142};
chopper_LI7500 = dataArray{:, 143};
detector_LI7500 = dataArray{:, 144};
pll_LI7500 = dataArray{:, 145};
sync_LI7500 = dataArray{:, 146};
not_ready_LI7700 = dataArray{:, 147};
no_signal_LI7700 = dataArray{:, 148};
re_unlocked_LI7700 = dataArray{:, 149};
bad_temp_LI7700 = dataArray{:, 150};
laser_temp_unregulated_LI7700 = dataArray{:, 151};
block_temp_unregulated_LI7700 = dataArray{:, 152};
motor_spinning_LI7700 = dataArray{:, 153};
pump_on_LI7700 = dataArray{:, 154};
top_heater_on_LI7700 = dataArray{:, 155};
bottom_heater_on_LI7700 = dataArray{:, 156};
calibrating_LI7700 = dataArray{:, 157};
motor_failure_LI7700 = dataArray{:, 158};
bad_aux_tc1_LI7700 = dataArray{:, 159};
bad_aux_tc2_LI7700 = dataArray{:, 160};
bad_aux_tc3_LI7700 = dataArray{:, 161};
box_connected_LI7700 = dataArray{:, 162};
mean_value_RSSI_LI7200 = dataArray{:, 163};
mean_value_LI7500 = dataArray{:, 164};
u_var = dataArray{:, 165};
v_var = dataArray{:, 166};
w_var = dataArray{:, 167};
ts_var = dataArray{:, 168};
co2_var = dataArray{:, 169};
h2o_var = dataArray{:, 170};
ch4_var = dataArray{:, 171};
none_var = dataArray{:, 172};
wts_cov = dataArray{:, 173};
wco2_cov = dataArray{:, 174};
wh2o_cov = dataArray{:, 175};
wch4_cov = dataArray{:, 176};
wnone_cov = dataArray{:, 177};
co2_mean = dataArray{:, 178};
h2o_mean = dataArray{:, 179};
vin_sf_mean = dataArray{:, 180};
dew_point_mean = dataArray{:, 181};
avg_signal_strength_7200_mean = dataArray{:, 182};
co2_signal_strength_7200_mean = dataArray{:, 183};
h2o_signal_strength_7200_mean = dataArray{:, 184};
delta_signal_strength_7200_mean = dataArray{:, 185};
flowrate_mean = dataArray{:, 186};
hit_power_mean = dataArray{:, 187};
hit_vin_mean = dataArray{:, 188};
CHK = dataArray{:, 189};
%% Copie des variables à la suite du header
fprintf(fidout,formatSpec,DATAH,filename,date,time,DOY,daytime,file_records,used_records,Tau,qc_Tau,rand_err_Tau,H,qc_H,rand_err_H,LE,qc_LE,rand_err_LE,co2_flux,qc_co2_flux,rand_err_co2_flux,h2o_flux,qc_h2o_flux,rand_err_h2o_flux,ch4_flux,qc_ch4_flux,rand_err_ch4_flux,none_flux,qc_none_flux,rand_err_none_flux,H_strg,LE_strg,co2_strg,h2o_strg,ch4_strg,none_strg,co2_vadv,h2o_vadv,ch4_vadv,none_vadv,co2_molar_density,co2_mole_fraction,co2_mixing_ratio,co2_time_lag,co2_def_timelag,h2o_molar_density,h2o_mole_fraction,h2o_mixing_ratio,h2o_time_lag,h2o_def_timelag,ch4_molar_density,ch4_mole_fraction,ch4_mixing_ratio,ch4_time_lag,ch4_def_timelag,none_molar_density,none_mole_fraction,none_mixing_ratio,none_time_lag,none_def_timelag,sonic_temperature,air_temperature,air_pressure,air_density,air_heat_capacity,air_molar_volume,ET,water_vapor_density,e,es,specific_humidity,RH,VPD,Tdew,u_unrot,v_unrot,w_unrot,u_rot,v_rot,w_rot,wind_speed,max_wind_speed,wind_dir,yaw,pitch,roll,ustar,TKE,L,VarName89,bowen_ratio,Tstar,model,x_peak,x_offset,x_10,x_30,x_50,x_70,x_90,un_Tau,Tau_scf,un_H,H_scf,un_LE,LE_scf,un_co2_flux,co2_scf,un_h2o_flux,h2o_scf,un_ch4_flux,ch4_scf,un_none_flux,un_none_scf,spikes_hf,amplitude_resolution_hf,drop_out_hf,absolute_limits_hf,skewness_kurtosis_hf,skewness_kurtosis_sf,discontinuities_hf,discontinuities_sf,timelag_hf,timelag_sf,attack_angle_hf,non_steady_wind_hf,u_spikes,v_spikes,w_spikes,ts_spikes,co2_spikes,h2o_spikes,ch4_spikes,none_spikes,head_detect_LI7200,t_out_LI7200,t_in_LI7200,aux_in_LI7200,delta_p_LI7200,chopper_LI7200,detector_LI7200,pll_LI7200,sync_LI7200,chopper_LI7500,detector_LI7500,pll_LI7500,sync_LI7500,not_ready_LI7700,no_signal_LI7700,re_unlocked_LI7700,bad_temp_LI7700,laser_temp_unregulated_LI7700,block_temp_unregulated_LI7700,motor_spinning_LI7700,pump_on_LI7700,top_heater_on_LI7700,bottom_heater_on_LI7700,calibrating_LI7700,motor_failure_LI7700,bad_aux_tc1_LI7700,bad_aux_tc2_LI7700,bad_aux_tc3_LI7700,box_connected_LI7700,mean_value_RSSI_LI7200,mean_value_LI7500,u_var,v_var,w_var,ts_var,co2_var,h2o_var,ch4_var,none_var,wts_cov,wco2_cov,wh2o_cov,wch4_cov,wnone_cov,co2_mean,h2o_mean,vin_sf_mean,dew_point_mean,avg_signal_strength_7200_mean,co2_signal_strength_7200_mean,h2o_signal_strength_7200_mean,delta_signal_strength_7200_mean,flowrate_mean,hit_power_mean,hit_vin_mean,CHK)
end
fclose(fidout) |
Partager