Fix the error for the gnucompiler

NOAA-GFDL · Oct 11, 2021 · 64b4a04 · 64b4a04
1 parent 81024dd
commit 64b4a04
Showing 1 changed file with 7 additions and 1 deletion.
diff --git a/model/fv_grid_utils.F90 b/model/fv_grid_utils.F90
@@ -201,6 +201,7 @@ subroutine grid_utils_init(Atm, npx, npy, npz, non_ortho, grid_type, c2l_order)
       sina   => Atm%gridstruct%sina_64
       cosa   => Atm%gridstruct%cosa_64
 ! MOLECULAR_DIFFUSION
+   if ( Atm%flagstruct%molecular_diffusion ) then
       area_u  => Atm%gridstruct%area_u_64
       area_v  => Atm%gridstruct%area_v_64
       dx6    => Atm%gridstruct%dx6_64
@@ -210,6 +211,7 @@ subroutine grid_utils_init(Atm, npx, npy, npz, non_ortho, grid_type, c2l_order)
       delv_6 => Atm%gridstruct%delv_6
       delu_5 => Atm%gridstruct%delu_5
       delv_5 => Atm%gridstruct%delv_5
+   endif
 
       divg_u => Atm%gridstruct%divg_u
       divg_v => Atm%gridstruct%divg_v
@@ -530,7 +532,9 @@ subroutine grid_utils_init(Atm, npx, npy, npz, non_ortho, grid_type, c2l_order)
            cosa = big_number
            sina = big_number
 ! MOLECULAR_DIFFUSION
+      if ( Atm%flagstruct%molecular_diffusion ) then
            sina_6 = big_number
+      endif
 
         do j=js,je+1
            do i=is,ie+1
@@ -726,7 +730,9 @@ subroutine grid_utils_init(Atm, npx, npy, npz, non_ortho, grid_type, c2l_order)
            cosa_s = 0.
            rsin_u = 1.
            rsin_v = 1.
-           sina_6 = 1.
+           if ( Atm%flagstruct%molecular_diffusion ) then
+               sina_6 = 1.
+           endif
    endif
 
    if ( grid_type < 3 ) then