Merge branch 'concepts/atmo-boundary-umin-parameter' into 'concepts/m…

…ain' (!10) icepack_parameters: make 'umin' a parameter
phil-blain · Feb 12, 2024 · b232242 · b232242
2 parents c28ab64 + 4066dcf
commit b232242
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 7 deletions.
diff --git a/columnphysics/icepack_atmo.F90 b/columnphysics/icepack_atmo.F90
@@ -65,7 +65,8 @@ subroutine atmo_boundary_layer (sfctype,            &
                                       uvel,     vvel,     &
                                       Uref,     zlvs      )
 
-      use icepack_parameters, only: highfreq, natmiter, atmiter_conv, zTrf
+      use icepack_parameters, only: highfreq, natmiter, atmiter_conv, zTrf, &
+                                    umin
 
       character (len=3), intent(in) :: &
          sfctype      ! ice or ocean
@@ -134,8 +135,7 @@ subroutine atmo_boundary_layer (sfctype,            &
          qqq   , & ! for qsat, dqsfcdt
          TTT   , & ! for qsat, dqsfcdt
          qsat  , & ! the saturation humidity of air (kg/m^3)
-         Lheat , & ! Lvap or Lsub, depending on surface type
-         umin      ! minimum wind speed (m/s)
+         Lheat     ! Lvap or Lsub, depending on surface type
 
       real (kind=dbl_kind) :: &
          ustar , & ! ustar (m/s)
@@ -172,8 +172,6 @@ subroutine atmo_boundary_layer (sfctype,            &
 
       if (highfreq) then
        umin  = p5 ! minumum allowable wind-ice speed difference of 0.5 m/s
-      else
-       umin  = c1 ! minumum allowable wind speed of 1m/s
       endif
 
       Tref = c0

diff --git a/columnphysics/icepack_parameters.F90 b/columnphysics/icepack_parameters.F90
@@ -253,6 +253,7 @@ module icepack_parameters
          zvir   = 0.606_dbl_kind     ,&! rh2o/rair - 1.0
          zref   = 10._dbl_kind       ,&! reference height for stability (m)
          zTrf   = c2                 ,&! reference height for {T,Q}ref (m)
+         umin   = c1                 ,&! min wind speed for turbulent fluxes (m/s)
          iceruf = 0.0005_dbl_kind    ,&! ice surface roughness (m)
          qqqice = 11637800._dbl_kind ,&! for qsat over ice
          TTTice = 5897.8_dbl_kind    ,&! for qsat over ice
@@ -444,7 +445,7 @@ subroutine icepack_init_parameters(   &
          stefan_boltzmann_in, ice_ref_salinity_in, &
          Tffresh_in, Lsub_in, Lvap_in, Timelt_in, Tsmelt_in, &
          iceruf_in, Cf_in, Pstar_in, Cstar_in, kappav_in, &
-         kice_in, ksno_in, zTrf_in, &
+         kice_in, ksno_in, zTrf_in, umin_in, &
          zref_in, hs_min_in, snowpatch_in, rhosi_in, sk_l_in, &
          saltmax_in, phi_init_in, min_salin_in, salt_loss_in, &
          min_bgc_in, dSin0_frazil_in, hi_ssl_in, hs_ssl_in, &
@@ -663,6 +664,7 @@ subroutine icepack_init_parameters(   &
          zvir_in,       & ! rh2o/rair - 1.0
          zref_in,       & ! reference height for stability (m)
          zTrf_in,       & ! reference height for {T,Q}ref (m)
+         umin_in,       & ! min wind speed for turbulent fluxes (m/s)
          qqqice_in,     & ! for qsat over ice
          TTTice_in,     & ! for qsat over ice
          qqqocn_in,     & ! for qsat over ocn
@@ -874,6 +876,7 @@ subroutine icepack_init_parameters(   &
       if (present(ksno_in)              ) ksno             = ksno_in
       if (present(zref_in)              ) zref             = zref_in
       if (present(zTrf_in)              ) zTrf             = zTrf_in
+      if (present(umin_in)              ) umin             = umin_in
       if (present(hs_min_in)            ) hs_min           = hs_min_in
       if (present(snowpatch_in)         ) snowpatch        = snowpatch_in
       if (present(rhosi_in)             ) rhosi            = rhosi_in
@@ -1155,7 +1158,7 @@ subroutine icepack_query_parameters(   &
          stefan_boltzmann_out, ice_ref_salinity_out, &
          Tffresh_out, Lsub_out, Lvap_out, Timelt_out, Tsmelt_out, &
          iceruf_out, Cf_out, Pstar_out, Cstar_out, kappav_out, &
-         kice_out, ksno_out, zTrf_out, &
+         kice_out, ksno_out, zTrf_out, umin_out, &
          zref_out, hs_min_out, snowpatch_out, rhosi_out, sk_l_out, &
          saltmax_out, phi_init_out, min_salin_out, salt_loss_out, &
          min_bgc_out, dSin0_frazil_out, hi_ssl_out, hs_ssl_out, &
@@ -1384,6 +1387,7 @@ subroutine icepack_query_parameters(   &
          zvir_out,       & ! rh2o/rair - 1.0
          zref_out,       & ! reference height for stability (m)
          zTrf_out,       & ! reference height for {T,Q}ref (m)
+         umin_out,       & ! min wind speed for turbulent fluxes (m/s)
          qqqice_out,     & ! for qsat over ice
          TTTice_out,     & ! for qsat over ice
          qqqocn_out,     & ! for qsat over ocn
@@ -1631,6 +1635,7 @@ subroutine icepack_query_parameters(   &
       if (present(ksno_out)              ) ksno_out         = ksno
       if (present(zref_out)              ) zref_out         = zref
       if (present(zTrf_out)              ) zTrf_out         = zTrf
+      if (present(umin_out)              ) umin_out         = umin
       if (present(hs_min_out)            ) hs_min_out       = hs_min
       if (present(snowpatch_out)         ) snowpatch_out    = snowpatch
       if (present(rhosi_out)             ) rhosi_out        = rhosi
@@ -1826,6 +1831,7 @@ subroutine icepack_write_parameters(iounit)
         write(iounit,*) "  ksno       = ",ksno
         write(iounit,*) "  zref       = ",zref
         write(iounit,*) "  zTrf       = ",zTrf
+        write(iounit,*) "  umin       = ",umin
         write(iounit,*) "  hs_min     = ",hs_min
         write(iounit,*) "  snowpatch  = ",snowpatch
         write(iounit,*) "  rhosi      = ",rhosi