Merge pull request #10 from eclare108213/set_parameters

Initialize Icepack constants using values in ice_constants_colpkg.F90
E3SM-Ocean-Discussion · Apr 3, 2023 · 82fe7ad · 82fe7ad
2 parents d913a2b + 0c952ae
commit 82fe7ad
Show file tree

Hide file tree

Showing 2 changed files with 340 additions and 164 deletions.
diff --git a/components/mpas-seaice/src/shared/mpas_seaice_constants.F b/components/mpas-seaice/src/shared/mpas_seaice_constants.F
@@ -19,13 +19,16 @@ module seaice_constants
 
   ! fundamental constants
   real (kind=RKIND), parameter, public :: &
-       pii = 3.141592653589793_RKIND, &
+       pii = 3.141592653589793_RKIND, &          ! CESM uses SHR_CONST_PI
        seaiceDegreesToRadians = pii / 180.0_RKIND, &
        seaiceRadiansToDegrees = 180.0_RKIND / pii, &
        seaiceSecondsPerYear = 24.0_RKIND * 3600.0_RKIND * 365.0_RKIND, &
        seaiceSecondsPerDay  = 24.0_RKIND * 3600.0_RKIND, &
        seaiceDaysPerSecond = 1.0_RKIND/seaiceSecondsPerDay
 
+  real (kind=RKIND), public :: &
+       seaicePi                           ! pi
+
   ! Earth constants
   real (kind=RKIND), public :: &
        seaiceGravity                      ! gravitational acceleration (m/s^2)
@@ -50,16 +53,22 @@ module seaice_constants
        seaiceStefanBoltzmann, &         ! J m-2 K-4 s-1
        seaiceIceSnowEmissivity, &       ! emissivity of snow and ice
        seaiceFreshWaterFreezingPoint, & ! freezing temp of fresh ice (K)
+       seaiceFreshIceSpecificHeat, &    ! specific heat of fresh ice (J/kg/K)
        seaiceAirSpecificHeat, &         ! specific heat of air (J/kg/K)
+       seaiceWaterVaporSpecificHeat, &  ! specific heat of water vapor (J/kg/K)
+       seaiceZvir, &                    ! rh2o/rair - 1.0
        seaiceLatentHeatSublimation, &   ! latent heat, sublimation freshwater (J/kg)
        seaiceLatentHeatMelting, &       ! latent heat of melting of fresh ice (J/kg)
+       seaiceIceSurfaceMeltingTemperature, &  ! melting temp. ice top surface  (C)
+       seaiceSnowSurfaceMeltingTemperature, & ! melting temp. snow top surface  (C)
        seaiceOceanAlbedo, &             ! Ocean albedo
        seaiceVonKarmanConstant, &       ! Von Karman constant
        seaiceIceSurfaceRoughness, &     ! ice surface roughness (m)
        seaiceSeaWaterSpecificHeat, &    ! specific heat of ocn (J/kg/K)
        seaiceLatentHeatVaporization, &  ! latent heat, vaporization freshwater (J/kg)
        seaiceReferenceSalinity, &       ! ice reference salinity (ppt)
-       seaiceStabilityReferenceHeight   ! stability reference height (m)
+       seaiceStabilityReferenceHeight, &! stability reference height (m)
+       seaiceSnowPatchiness             ! snow patchiness parameter
 
   ! dynamics constants
   real(kind=RKIND), public :: &