diff --git a/src/model_mod.F90 b/src/model_mod.F90 index b65026949373568956d5c786c0524f81eca6e606..1dc386eb733bef86e833cefb9e63aeb011e9a5c1 100644 --- a/src/model_mod.F90 +++ b/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 diff --git a/src/stepon.F90 b/src/stepon.F90 index af24f83107232d8346f49b63e8b5caa80e751d9a..91a7298b1281a5072d2c7c0559e5e1d542b7e008 100644 --- a/src/stepon.F90 +++ b/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