diff --git a/matlab/setup.m b/matlab/setup.m
index 3236997da61f6536305bb378ab658ff2d20eb8d2..9e0490c8a1902dee8396db7962a8a94a9c9b88a7 100644
--- a/matlab/setup.m
+++ b/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;
diff --git a/matlab/write_fort90.m b/matlab/write_fort90.m
index eddf37789f23ce81d1fa0118d5a3d9b016750899..67520b5a0a8516b5e631ce9a0cf5a4623fb7f1ad 100644
--- a/matlab/write_fort90.m
+++ b/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']);
diff --git a/wk/lin_Ivanov.m b/wk/lin_Ivanov.m
index 505fe95e2246945b01035c824a2e5e584bdfd2ee..5ecf60b0002a993319ed9fe98ae466292fc7544e 100644
--- a/wk/lin_Ivanov.m
+++ b/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;