Merge pull request #855 from FESOM/bugfix/cavity_3eq_system

Update water_flux calculation in cavity_param.F90

Author: JanStreffing

setups/test_pi_cavity/setup.yml

@@ -59,12 +59,12 @@ namelist.io:
                 prec: 8
 
 fcheck:
-    a_ice: 0.8922816940404571
-    salt: 34.71518655581704
-    sst: 12.806226630290732
-    temp: 2.844705384878494
-    u: -0.0027710946118024962
-    v: -0.00025344388118151974
+    a_ice: 0.8922816930614722
+    salt: 34.71518654875407
+    sst: 12.806226998388576
+    temp: 2.8447054270089347
+    u: -0.002771103825917723
+    v: -0.00025341094981021055
 
 
 

src/cavity_param.F90

@@ -362,7 +362,7 @@ subroutine cavity_heat_water_fluxes_3eq(ice, dynamics, tracers, partit, mesh)
         !rt  s_surf_flux(i,j)=gas*(sf-(s(i,j,N,lrhs)+35.0))
 
         heat_flux(node)  = rhow*cpw*gat*(tin-tf)      ! [W/m2]  ! positive for upward
-        water_flux(node) =          gas*(sf-sal)/sf   ! [m/s]   !
+        water_flux(node) =          ep5*(sf-sal)/sf   ! [m/s]   !
 
         !      qo=-rhor*seta*oofw
         !      if(seta.le.0.) then