Making RBD ustar filter switch more apparent

Author: Gregory E. Maurer

UNM_RemoveBadData.m

@@ -101,6 +101,17 @@
 
 draw_plots = args.Results.draw_plots;
 
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% USTAR filter switch
+% If this is set to false, the ustar filter for fluxes is turned off
+% (this is desireable for making ameriflux files)
+% See line 1272 to see where this takes affect (its pretty hacky)
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+ustar_filter_switch = false;
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% 
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
 
 data_for_analyses = 0; %1 to output file with data sorted for specific analyses
 ET_gap_filler = 0; %run ET gap-filler program
@@ -1258,8 +1269,10 @@
 if iteration > 1
 
     % Remove values with low U*
-    
-    ustar_lim = 0;
+    % Only happens if switch is set to true
+    if ~ustar_filter_switch
+        ustar_lim = 0;
+    end
     ustarflag = find(u_star < ustar_lim);
     removed_ustar = length(ustarflag);
     decimal_day_nan(ustarflag) = NaN;

scripts/script_make_ameriflux.m

@@ -1,20 +1,46 @@
 %close all;
 %clear all;
 %
-sitelist = {UNM_sites.MCon, UNM_sites.JSav, UNM_sites.PJ, UNM_sites.PJ_girdle, ...
-    UNM_sites.SLand, UNM_sites.GLand, UNM_sites.PPine, UNM_sites.New_GLand};
-%sitelist = {UNM_sites.JSav, ...
-%    UNM_sites.SLand, UNM_sites.GLand, UNM_sites.PPine, UNM_sites.MCon};
-sitelist={UNM_sites.MCon};
+
+ sitelist_1 = {UNM_sites.MCon, UNM_sites.JSav, UNM_sites.PJ, UNM_sites.PJ_girdle, ...
+     UNM_sites.SLand, UNM_sites.GLand, UNM_sites.PPine, UNM_sites.New_GLand};
+ yearlist_1 = 2010:2015;
+ 
+ sitelist_2 = {UNM_sites.JSav, ...
+     UNM_sites.SLand, UNM_sites.GLand, UNM_sites.PPine, UNM_sites.MCon};
+ yearlist_2 = 2007:2009;
+ 
+sitelist_3={UNM_sites.PJ};
 % Years to create files for
-yearlist = 2016;
+yearlist_3 = 2008:2009;
+
+sitelist_4 = {UNM_sites.PJ_girdle}
+yearlist = 2009 ;
 % Partitioned data source
 partmethod = 'eddyproc'; %'Reddyproc'
 % Make daily files? All AF files should be in $FLUXROOT$/Ameriflux_files
 make_daily = false;
-write_files = false;
+write_files = true;
 process_soil = false;
 
+for k = 1:4
+    % Automatically cycles through earlier site years to account for
+    % different startup dates.
+    switch k
+        case 1
+            sitelist = sitelist_1;
+            yearlist = yearlist_1;
+        case 2
+            sitelist = sitelist_2;
+            yearlist = yearlist_2;
+        case 3
+            sitelist = sitelist_3;
+            yearlist = yearlist_3;
+        case 4
+            sitelist = sitelist_4;
+            yearlist = yearlist_4;
+    end
+
 for i = 1:length(sitelist);
     %close all;
     for j = 1:length(yearlist);
@@ -47,4 +73,5 @@
     close all;
     clear sitecode;
 end
+end
 

scripts/script_make_qc_gf.m

@@ -4,7 +4,7 @@
 sitelist = {UNM_sites.GLand, UNM_sites.SLand, UNM_sites.JSav,...
     UNM_sites.PPine, UNM_sites.MCon, UNM_sites.PJ, UNM_sites.PJ_girdle, ...
     UNM_sites.New_GLand };
-sitelist = {UNM_sites.PPine};
+sitelist = {UNM_sites.New};
 yearlist = 2016;
 
 % True, overwrite files; False; do not overwrite