diff --git a/dyn_em/module_advect_em.F b/dyn_em/module_advect_em.F
index ff03d1a695..2913ac5a9c 100644
--- a/dyn_em/module_advect_em.F
+++ b/dyn_em/module_advect_em.F
@@ -1,7 +1,7 @@
 
 !WRF:MODEL_LAYER:DYNAMICS
 !
-#if ( defined(ADVECT_KERNEL) )
+#ifdef ADVECT_KERNEL
 ! cpp -traditional-cpp -P -DADVECT_KERNEL module_advect_em.F > advection_kernel.f90
 ! gfortran -ffree-form -ffree-line-length-none advection_kernel.f90
 ! ./a.out
@@ -111,7 +111,7 @@ SUBROUTINE column (loop , data_list, its,ite)
    
 END SUBROUTINE column
 !----------------------------------------------------------------
-#elif ( ! defined(ADVECT_KERNEL) )
+#else
 
 MODULE module_advect_em
 
@@ -4357,7 +4357,7 @@ SUBROUTINE advect_scalar   ( field, field_old, tendency,    &
    ENDIF vert_order_test
 
 END SUBROUTINE advect_scalar
-#if ( ! defined(ADVECT_KERNEL) )
+#ifndef ADVECT_KERNEL
 
 !---------------------------------------------------------------------------------
 
@@ -10543,7 +10543,7 @@ END SUBROUTINE advect_scalar_mono
 
 !-----------------------------------------------------------
 
-#if ( defined(ADVECT_KERNEL) )
+#ifdef ADVECT_KERNEL
 
 END MODULE advection_kernel
 !================================================================
@@ -10851,7 +10851,7 @@ PROGRAM feeder
 
 END PROGRAM feeder
 #endif
-#if ( !defined(ADVECT_KERNEL) )
+#ifndef ADVECT_KERNEL
 
 !---------------------------------------------------------------------------------
 
diff --git a/external/atm_ocn/cmpcomm.F b/external/atm_ocn/cmpcomm.F
index a78e285337..89cd554e1c 100644
--- a/external/atm_ocn/cmpcomm.F
+++ b/external/atm_ocn/cmpcomm.F
@@ -1,4 +1,4 @@
-#if defined( DM_PARALLEL ) 
+#ifdef DM_PARALLEL
       MODULE CMP_COMM
 
       implicit none
diff --git a/external/io_adios2/wrf_io.F90 b/external/io_adios2/wrf_io.F90
index 3d5fdd6844..d53ad88481 100644
--- a/external/io_adios2/wrf_io.F90
+++ b/external/io_adios2/wrf_io.F90
@@ -702,9 +702,9 @@ subroutine Transpose(IO,MemoryOrder,di, Field,l1,l2,m1,m2,n1,n2 &
  
    call LowerCase(MemoryOrder,MemOrd)
    select case (MemOrd)
- 
- !#define XDEX(A,B,C) A-A ## 1+1+(A ## 2-A ## 1+1)*((B-B ## 1)+(C-C ## 1)*(B ## 2-B ## 1+1))
- ! define(`XDEX',($1-``$1''1+1+(``$1''2-``$1''1+1)*(($2-``$2''1)+($3-``$3''1)*(``$2''2-``$2''1+1))))
+! Cannot use following define due to gfortran cpp traditional mode concatenation limitations 
+!#define XDEX(A,B,C) A-A ## 1+1+(A ## 2-A ## 1+1)*((B-B ## 1)+(C-C ## 1)*(B ## 2-B ## 1+1))
+! define(`XDEX',($1-``$1''1+1+(``$1''2-``$1''1+1)*(($2-``$2''1)+($3-``$3''1)*(``$2''2-``$2''1+1))))
      case ('xzy')
 #undef  DFIELD
 #define DFIELD XField(1:di,XDEX(i,k,j))
diff --git a/external/io_netcdf/wrf_io.F90 b/external/io_netcdf/wrf_io.F90
index ec2162d2d5..8863e4e29f 100644
--- a/external/io_netcdf/wrf_io.F90
+++ b/external/io_netcdf/wrf_io.F90
@@ -754,7 +754,7 @@ subroutine Transpose(IO,MemoryOrder,di, Field,l1,l2,m1,m2,n1,n2 &
 
   call LowerCase(MemoryOrder,MemOrd)
   select case (MemOrd)
-
+! Cannot use following define due to gfortran cpp traditional mode concatenation limitations
 !#define XDEX(A,B,C) A-A ## 1+1+(A ## 2-A ## 1+1)*((B-B ## 1)+(C-C ## 1)*(B ## 2-B ## 1+1))
 ! define(`XDEX',($1-``$1''1+1+(``$1''2-``$1''1+1)*(($2-``$2''1)+($3-``$3''1)*(``$2''2-``$2''1+1))))
 
@@ -940,7 +940,7 @@ subroutine TransposeToR4(IO,MemoryOrder,di, Field,l1,l2,m1,m2,n1,n2 &
 
   call LowerCase(MemoryOrder,MemOrd)
   select case (MemOrd)
-
+! Cannot use following define due to gfortran cpp traditional mode concatenation limitations
 !#define XDEX(A,B,C) A-A ## 1+1+(A ## 2-A ## 1+1)*((B-B ## 1)+(C-C ## 1)*(B ## 2-B ## 1+1))
 ! define(`XDEX',($1-``$1''1+1+(``$1''2-``$1''1+1)*(($2-``$2''1)+($3-``$3''1)*(``$2''2-``$2''1+1))))
 
diff --git a/external/io_netcdfpar/wrf_io.F90 b/external/io_netcdfpar/wrf_io.F90
index a76ec5d82d..86e25dd2cb 100644
--- a/external/io_netcdfpar/wrf_io.F90
+++ b/external/io_netcdfpar/wrf_io.F90
@@ -767,7 +767,7 @@ subroutine Transpose(IO,MemoryOrder,di, Field,l1,l2,m1,m2,n1,n2 &
 
   call LowerCase(MemoryOrder,MemOrd)
   select case (MemOrd)
-
+! Cannot use following define due to gfortran cpp traditional mode concatenation limitations
 !#define XDEX(A,B,C) A-A ## 1+1+(A ## 2-A ## 1+1)*((B-B ## 1)+(C-C ## 1)*(B ## 2-B ## 1+1))
 ! define(`XDEX',($1-``$1''1+1+(``$1''2-``$1''1+1)*(($2-``$2''1)+($3-``$3''1)*(``$2''2-``$2''1+1))))
 
@@ -953,7 +953,7 @@ subroutine TransposeToR4a(IO,MemoryOrder,di, Field,l1,l2,m1,m2,n1,n2 &
 
   call LowerCase(MemoryOrder,MemOrd)
   select case (MemOrd)
-
+! Cannot use following define due to gfortran cpp traditional mode concatenation limitations
 !#define XDEX(A,B,C) A-A ## 1+1+(A ## 2-A ## 1+1)*((B-B ## 1)+(C-C ## 1)*(B ## 2-B ## 1+1))
 ! define(`XDEX',($1-``$1''1+1+(``$1''2-``$1''1+1)*(($2-``$2''1)+($3-``$3''1)*(``$2''2-``$2''1+1))))
 
diff --git a/external/io_pnetcdf/wrf_io.F90 b/external/io_pnetcdf/wrf_io.F90
index 18f6ac078a..9d9c3733b4 100644
--- a/external/io_pnetcdf/wrf_io.F90
+++ b/external/io_pnetcdf/wrf_io.F90
@@ -740,7 +740,7 @@ subroutine Transpose(IO,MemoryOrder,di, Field,l1,l2,m1,m2,n1,n2 &
 
   call LowerCase(MemoryOrder,MemOrd)
   select case (MemOrd)
-
+! Cannot use following define due to gfortran cpp traditional mode concatenation limitations
 !#define XDEX(A,B,C) A-A ## 1+1+(A ## 2-A ## 1+1)*((B-B ## 1)+(C-C ## 1)*(B ## 2-B ## 1+1))
 ! define(`XDEX',($1-``$1''1+1+(``$1''2-``$1''1+1)*(($2-``$2''1)+($3-``$3''1)*(``$2''2-``$2''1+1))))
 
diff --git a/frame/module_configure.F b/frame/module_configure.F
index 4e0ae808c3..8554a7d92a 100644
--- a/frame/module_configure.F
+++ b/frame/module_configure.F
@@ -15,7 +15,8 @@ SUBROUTINE init_module_scalar_tables
   END SUBROUTINE init_module_scalar_tables
 END MODULE module_scalar_tables
 
-#if( WRF_CHEM == 1 && WRF_KPP == 1 )
+#ifdef WRF_CHEM
+#ifdef WRF_KPP
 MODULE module_irr_diag
 
    INTEGER, parameter :: max_eqn = 1200
@@ -45,6 +46,7 @@ END SUBROUTINE init_module_irr_diag
 
 END MODULE module_irr_diag
 #endif
+#endif
 
 MODULE module_configure
 
diff --git a/phys/module_mp_SBM_polar_radar.F b/phys/module_mp_SBM_polar_radar.F
index 4f94129271..a6ba4e4cc1 100644
--- a/phys/module_mp_SBM_polar_radar.F
+++ b/phys/module_mp_SBM_polar_radar.F
@@ -6,6 +6,11 @@ SUBROUTINE SBM_polar_radar
          dummy = 1
       END SUBROUTINE SBM_polar_radar
       END MODULE module_mp_SBM_polar_radar
+
+! Stub module
+module scatt_tables
+end module scatt_tables
+
 #else
 !******************
 module scatt_tables
diff --git a/phys/module_mp_fast_sbm.F b/phys/module_mp_fast_sbm.F
index 184a9220ea..757f05db8c 100644
--- a/phys/module_mp_fast_sbm.F
+++ b/phys/module_mp_fast_sbm.F
@@ -6,6 +6,20 @@ SUBROUTINE SBM_fast
          dummy = 1
       END SUBROUTINE SBM_fast
       END MODULE module_mp_fast_sbm
+
+! Stub modules
+module module_mp_SBM_BreakUp
+end module module_mp_SBM_BreakUp
+
+module module_mp_SBM_Collision
+end module module_mp_SBM_Collision
+
+module module_mp_SBM_Auxiliary
+end module module_mp_SBM_Auxiliary
+
+module module_mp_SBM_Nucleation
+end module module_mp_SBM_Nucleation
+
 #else
 ! +-----------------------------------------------------------------------------+
 ! +-----------------------------------------------------------------------------+