Fixes for ustar and small to encourage heat exchange

dgwyther · Nov 22, 2013 · 290264e · 290264e
1 parent 4934e36
commit 290264e
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 2 deletions.
diff --git a/README b/README
@@ -4,3 +4,4 @@
 20/11/2013: Eva's edits for 'define' to 'defined' in 2 ana_ files
 21/11/2013: Eva's edits to iceshelf_vbc EWperiodic to EW_periodic
 21/11/2013: Bad code in iceshelf_vbc that executes with ICETEST
+22/11/2013: Fixes for ustar and small to encourage heat exchange
diff --git a/ROMS/Nonlinear/IceShelf/iceshelf_mod.h b/ROMS/Nonlinear/IceShelf/iceshelf_mod.h
@@ -32,7 +32,7 @@
       real(r8), parameter :: Cd = 5.0e-3_r8
       real(r8), parameter :: visc = 1.95e-6_r8
       real(r8), parameter :: L = 3.33e5_r8
-      real(r8), parameter :: small = 1.0e-10_r8
+      real(r8), parameter :: small = 1.0e-4_r8
       real(r8), parameter :: dt_i = 1.54e-6
       real(r8), parameter :: cp_w = 3947.0_r8
       real(r8), parameter :: rho_i = 920.0_r8

diff --git a/ROMS/Nonlinear/IceShelf/iceshelf_vbc.F b/ROMS/Nonlinear/IceShelf/iceshelf_vbc.F
@@ -304,7 +304,7 @@ SUBROUTINE iceshelf_vbc_tile (ng, tile,                           &
             gammaS(i,j) = ustar/(turb + Scadj - 9.0_r8)
 #    elif defined GAMMAS_Mc1987
 ! Uses simplified version of McPhee 1987
-            IF (ustar.gt.small.and.ABS(f(i,j)).gt.1.0E-8) THEN
+            IF (ustar.ge.small.and.ABS(f(i,j)).gt.1.0E-8) THEN
               turb = 2.5*LOG(5300.0_r8*ustar*ustar/ABS(f(i,j))) + 7.12  
             ELSE
               turb = 0.0_r8