From 14b3e8e7b63cd20eee59775cec4da2ded3c872ea Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Fri, 26 May 2023 15:37:14 +0200
Subject: [PATCH] add the remove landau damping parameter

---
 matlab/setup.m        | 1 +
 matlab/write_fort90.m | 3 ++-
 wk/lin_Ivanov.m       | 1 +
 3 files changed, 4 insertions(+), 1 deletion(-)

diff --git a/matlab/setup.m b/matlab/setup.m
index 3236997d..9e0490c8 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 eddf3778..67520b5a 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 505fe95e..5ecf60b0 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;
-- 
GitLab