Input parameter to impose symmetry on ky=0 modes

src/model_mod.F90

@@ -28,9 +28,12 @@ MODULE model
   ! Auxiliary variable
   LOGICAL,  PUBLIC, PROTECTED ::      EM =  .false.   ! Electromagnetic effects flag
   LOGICAL,  PUBLIC, PROTECTED ::  MHD_PD =  .false.   ! MHD pressure drift
-
   ! Removes Landau damping in temperature and higher equation (Ivanov 2022)
   LOGICAL,  PUBLIC, PROTECTED :: RM_LD_T_EQ = .false.
+  ! Flag to force the reality condition symmetry for the kx at ky=0
+  LOGICAL,  PUBLIC, PROTECTED :: FORCE_SYMMETRY = .false.
+
+  ! Module's routines
   PUBLIC :: model_readinputs, model_outputinputs
 
 CONTAINS
@@ -42,7 +45,7 @@ CONTAINS
     USE prec_const
     IMPLICIT NONE
 
-    NAMELIST /MODEL_PAR/ KERN, LINEARITY, RM_LD_T_EQ, &
+    NAMELIST /MODEL_PAR/ KERN, LINEARITY, RM_LD_T_EQ, FORCE_SYMMETRY, &
                          mu_x, mu_y, N_HD, HDz_h, mu_z, mu_p, mu_j, HYP_V, Na,&
                          nu, k_gB, k_cB, lambdaD, MHD_PD, beta, ADIAB_E, ADIAB_I, tau_i
 

src/stepon.F90

@@ -110,13 +110,13 @@ CONTAINS
       USE utility,ONLY: is_nan, is_inf
       USE fields, ONLY: phi
       USE MPI
-      USE model,            ONLY: LINEARITY
+      USE model,            ONLY: LINEARITY, FORCE_SYMMETRY
       IMPLICIT NONE
       LOGICAL :: checkf_
       REAL    :: sum_
-      !! Maybe useless
-      IF(LINEARITY .NE. 'linear') CALL anti_aliasing   ! ensure 0 mode for 2/3 rule
-      ! IF(LINEARITY .NE. 'linear') CALL enforce_symmetry ! Enforcing symmetry on kx = 0
+      ! filtering
+      IF(LINEARITY .NE. 'linear') CALL anti_aliasing    ! ensure 0 mode for 2/3 rule
+      IF(FORCE_SYMMETRY)          CALL enforce_symmetry ! Enforcing symmetry on kx = 0 (looks useless)
 
       mlend=.FALSE.
       IF(.NOT.nlend) THEN