Move inv_PI3 outside of k loop

NOAA-GFDL · Jan 24, 2019 · 581b30d · 581b30d
1 parent 462dfd8
commit 581b30d
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/parameterizations/lateral/MOM_lateral_mixing_coeffs.F90 b/src/parameterizations/lateral/MOM_lateral_mixing_coeffs.F90
@@ -777,9 +777,10 @@ subroutine calc_QG_Leith_viscosity(CS, G, GV, h, k, div_xx_dx, div_xx_dy, vort_x
   is  = G%isc  ; ie  = G%iec  ; js  = G%jsc  ; je  = G%jec
   Isq = G%IscB ; Ieq = G%IecB ; Jsq = G%JscB ; Jeq = G%JecB
 
+  inv_PI3 = 1.0/((4.0*atan(1.0))**3)
+
   if (k > 1) then
 
-    inv_PI3 = 1.0/((4.0*atan(1.0))**3)
   ! Add in stretching term for the QG Leith vsicosity
 !  if (CS%use_QG_Leith) then
 !    do j=js-1,Jeq+1 ; do I=is-2,Ieq+1