function [Z_select, f_T1_select,f_T2_select] = RB_linear_init(Z, criteria)
%%% with deterministic initial data but the pde depends on z

%% set up
Nz = 50;
%initial data for all Z
%%%%%%%%%%%%only need to run once%%%%%%%%%%%%
% Z = Z_select_init(Z,Nz);
% save('initial_data_eq7','Z');
%%%%%%%%%%%%only need to run once%%%%%%%%%%%%
load('initial_data_eq7');
[Z_select, f_T1_select,f_T2_select] = RB_linear(Z, criteria);




