Inialize uvel, vvel at every substep (#43)

apcraig · Nov 29, 2021 · c2c0c4a · c2c0c4a
1 parent 9e6e8f3
commit c2c0c4a
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/cicecore/cicedynB/dynamics/ice_dyn_evp.F90 b/cicecore/cicedynB/dynamics/ice_dyn_evp.F90
@@ -646,6 +646,13 @@ subroutine evp (dt)
 
          do ksub = 1,ndte        ! subcycling
 
+         ! shift velocity components from CD grid locations (N, E) to B grid location (U) for stress_U
+
+            if (grid_system == 'CD') then
+                call grid_average_X2Y('E2US',uvelE,uvel)
+                call grid_average_X2Y('N2US',vvelN,vvel)
+            endif
+
          !-----------------------------------------------------------------
          ! stress tensor equation, total surface stress
          !-----------------------------------------------------------------