add the remove landau damping parameter

matlab/setup.m

@@ -23,6 +23,7 @@ GEOM.shift_y  = SHIFT_Y;
 GEOM.Npol  = NPOL;
 % Model parameters
 MODEL.LINEARITY = ['''',LINEARITY,''''];
+if RM_LD_T_EQ; MODEL.RM_LD_T_EQ = '.true.'; else; MODEL.RM_LD_T_EQ = '.false.'; end;
 MODEL.Na        = NA;
 if ADIAB_E; MODEL.ADIAB_E = '.true.'; else; MODEL.ADIAB_E = '.false.';end;
 MODEL.beta    = BETA;

matlab/write_fort90.m

@@ -53,7 +53,8 @@ fprintf(fid,['  write_temp  = ', OUTPUTS.write_temp,'\n']);
 fprintf(fid,'/\n');
 
 fprintf(fid,'&MODEL_PAR\n');
-fprintf(fid,['  LINEARITY = ', MODEL.LINEARITY,'\n']);
+fprintf(fid,['LINEARITY = ', MODEL.LINEARITY,'\n']);
+fprintf(fid,['RM_LD_T_EQ= ', MODEL.RM_LD_T_EQ,'\n']);
 fprintf(fid,['  Na      = ', num2str(MODEL.Na),'\n']);
 fprintf(fid,['  mu_x    = ', num2str(MODEL.mu_x),'\n']);
 fprintf(fid,['  mu_y    = ', num2str(MODEL.mu_y),'\n']);

wk/lin_Ivanov.m

@@ -31,6 +31,7 @@ SIGMA_E = 0.0233380;        % mass ratio sqrt(m_a/m_i) (correct = 0.0233380)
 NA = 1;                     % number of kinetic species
 ADIAB_E = (NA==1);          % adiabatic electron model
 BETA = 0.0;                 % electron plasma beta
+RM_LD_T_EQ = 1;             % to remove Landau damping in temperature and higher eq.
 %% Set up grid parameters
 P = 2;
 J = 1;%P/2;