From 711a613ec8ff1c70b8cd75f75cec1ac990a29636 Mon Sep 17 00:00:00 2001 From: Gustavo Marques Date: Fri, 23 Jun 2017 13:48:50 -0600 Subject: [PATCH] Removes frazil's contribution to evaporation Frazil formation was being converted to surface layer evaporation under ice sheves and this was double counting the effect of frazil on the system. --- src/ice_shelf/MOM_ice_shelf.F90 | 11 ----------- 1 file changed, 11 deletions(-) diff --git a/src/ice_shelf/MOM_ice_shelf.F90 b/src/ice_shelf/MOM_ice_shelf.F90 index f7de07a41c..1294ed065b 100644 --- a/src/ice_shelf/MOM_ice_shelf.F90 +++ b/src/ice_shelf/MOM_ice_shelf.F90 @@ -1053,17 +1053,6 @@ subroutine add_shelf_flux(G, CS, state, fluxes) endif - ! Add frazil formation diagnosed by the ocean model (J m-2) in the - ! form of surface layer evaporation (kg m-2 s-1). Update lprec in the - ! control structure for diagnostic purposes. - - if (associated(state%frazil)) then - fraz = state%frazil(i,j) / CS%time_step / CS%Lat_fusion - if (associated(fluxes%evap)) fluxes%evap(i,j) = fluxes%evap(i,j) - fraz - CS%lprec(i,j)=CS%lprec(i,j) - fraz - state%frazil(i,j) = 0.0 - endif - if (associated(fluxes%sens)) fluxes%sens(i,j) = -frac_area*CS%t_flux(i,j)*CS%flux_factor if (associated(fluxes%salt_flux)) fluxes%salt_flux(i,j) = frac_area * CS%salt_flux(i,j)*CS%flux_factor if (associated(fluxes%p_surf)) fluxes%p_surf(i,j) = frac_area * CS%g_Earth * CS%mass_shelf(i,j)