diff --git a/matlab/setup.m b/matlab/setup.m
index 562293dd39382b95e3c3a8fb516fc80023a0d4c3..2cdb3cc20077c9614bf5316cf2a133a4dfcb501c 100644
--- a/matlab/setup.m
+++ b/matlab/setup.m
@@ -71,7 +71,6 @@ end
 MODEL.K_mB   = K_mB;      % artificial mirror force   tuner
 MODEL.K_tB   = K_tB;      % artificial trapping term  tuner
 MODEL.K_ldB  = K_ldB;     % artificial Landau damping tuner
-MODEL.lambdaD = LAMBDAD;
 % CLOSURE parameters
 CLOSURE.hierarchy_closure = ['''',HRCY_CLOS,''''];
 CLOSURE.nonlinear_closure = ['''',NLIN_CLOS,''''];
diff --git a/matlab/write_fort90.m b/matlab/write_fort90.m
index 01e41e080b09397b3f7b9b8795eb36fce078eef6..398628c6f1f8f9dc4c7e2f879551ee11d8bc6b21 100644
--- a/matlab/write_fort90.m
+++ b/matlab/write_fort90.m
@@ -74,7 +74,6 @@ fprintf(fid,['  k_cB    = ', num2str(MODEL.K_cB),'\n']);
 fprintf(fid,['  k_mB    = ', num2str(MODEL.K_mB),'\n']);
 fprintf(fid,['  k_tB    = ', num2str(MODEL.K_tB),'\n']);
 fprintf(fid,['  k_ldB   = ', num2str(MODEL.K_ldB),'\n']);
-fprintf(fid,['  lambdaD = ', num2str(MODEL.lambdaD),'\n']);
 fprintf(fid,['  beta    = ', num2str(MODEL.beta),'\n']);
 fprintf(fid,['  ExBrate = ', num2str(MODEL.ExBrate),'\n']);
 fprintf(fid,['  ADIAB_E = ', MODEL.ADIAB_E,'\n']);