Move rho_floatation from module to local variable

Author: cbegeman

components/mpas-ocean/src/shared/mpas_ocn_diagnostics.F

@@ -79,7 +79,6 @@ module ocn_diagnostics
    integer :: ke_cell_flag, ke_vertex_flag
    real (kind=RKIND) ::  fCoef
    real (kind=RKIND) ::  landIceTopDragCoeff
-   real (kind=RKIND) ::  rho_floatation
    real (kind=RKIND), pointer ::  coef_3rd_order
 
    ! Methods for computing thickness at edges for flux calculations
@@ -3679,7 +3678,8 @@ subroutine ocn_compute_land_ice_flux_input_fields(layerThickness, normalVelocity
 
       integer :: iCell, iEdge, cell1, cell2, iLevel, i
 
-      real (kind=RKIND) :: blThickness, dz, weightSum, h_nu, Gamma_turb, landIceEdgeFraction, velocityMagnitude
+      real (kind=RKIND) :: blThickness, dz, weightSum, h_nu, Gamma_turb, landIceEdgeFraction, &
+                           velocityMagnitude, rho_floatation
 
       ! Scratch Arrays
       real (kind=RKIND), dimension(:), allocatable ::  &
@@ -3705,6 +3705,8 @@ subroutine ocn_compute_land_ice_flux_input_fields(layerThickness, normalVelocity
 
       elseif ( trim(config_land_ice_draft_mode) == "pressure-dependent" ) then
 
+         rho_floatation = config_land_ice_rho_ocean
+
 #ifdef MPAS_OPENACC
          !$acc parallel loop present(landIceDraftForSsh, landIcePressure)
 #else
@@ -4693,8 +4695,6 @@ subroutine ocn_diagnostics_init(domain, err)!{{{
          landIceTopDragCoeff = config_land_ice_flux_explicit_topDragCoeff
       endif
 
-      rho_floatation = config_land_ice_rho_ocean
-
       if (trim(config_thickness_drag_type) == 'centered') then
          thickEdgeDragChoice = thickEdgeDragCenter
       elseif (trim(config_thickness_drag_type) == 'harmonic-mean') then