First step in adding Tf

columnphysics/icepack_itd.F90

@@ -27,7 +27,7 @@ module icepack_itd
 
       use icepack_kinds
       use icepack_parameters, only: c0, c1, c2, c3, c15, c25, c100, p1, p01, p001, p5, puny
-      use icepack_parameters, only: Lfresh, rhos, ice_ref_salinity, hs_min, cp_ice, Tocnfrz, rhoi
+      use icepack_parameters, only: Lfresh, rhos, ice_ref_salinity, hs_min, cp_ice, rhoi
       use icepack_parameters, only: rhosi, sk_l, hs_ssl, min_salin, rsnw_fall
       use icepack_tracers,    only: nt_Tsfc, nt_qice, nt_qsno, nt_aero, nt_isosno, nt_isoice
       use icepack_tracers,    only: nt_apnd, nt_hpnd, nt_fbri, tr_brine, nt_bgc_S, bio_index
@@ -108,7 +108,7 @@ subroutine rebin (ntrcr,    trcr_depend,     &
                         nt_strata,                 &
                         aicen,    trcrn,           &
                         vicen,    vsnon,           &
-                        ncat,     hin_max          )
+                        ncat,     hin_max, Tf      )
 
       integer (kind=int_kind), intent(in) :: &
          ntrcr , & ! number of tracers in use
@@ -135,6 +135,9 @@ subroutine rebin (ntrcr,    trcr_depend,     &
 
       real (kind=dbl_kind), dimension(0:ncat), intent(in) :: &
          hin_max   ! category limits (m)
+
+      real (kind=dbl_kind), intent(in) :: &
+         Tf                ! freezing temperature
 
       ! local variables
 
@@ -221,7 +224,7 @@ subroutine rebin (ntrcr,    trcr_depend,     &
                             aicen,    trcrn,      &
                             vicen,    vsnon,      &
                             hicen,    donor,      &
-                            daice,    dvice       )
+                            daice,    dvice, Tf   )
             if (icepack_warnings_aborted(subname)) return
 
       !-----------------------------------------------------------------
@@ -269,7 +272,7 @@ subroutine rebin (ntrcr,    trcr_depend,     &
                             aicen,    trcrn,      &
                             vicen,    vsnon,      &
                             hicen,    donor,      &
-                            daice,    dvice       )
+                            daice,    dvice, Tf   )
             if (icepack_warnings_aborted(subname)) return
 
       !-----------------------------------------------------------------
@@ -361,7 +364,7 @@ subroutine shift_ice (ntrcr,    ncat,        &
                             aicen,    trcrn,       &
                             vicen,    vsnon,       &
                             hicen,    donor,       &
-                            daice,    dvice        )
+                            daice,    dvice, Tf    )
 
       integer (kind=int_kind), intent(in) :: &
          ncat  , & ! number of thickness categories
@@ -396,6 +399,9 @@ subroutine shift_ice (ntrcr,    ncat,        &
          dvice         , & ! ice volume transferred across boundary
          hicen             ! ice thickness for each cat        (m)
 
+      real (kind=dbl_kind), intent(in) :: &
+         Tf                ! freezing temperature
+
       ! local variables
 
       integer (kind=int_kind) :: &
@@ -663,7 +669,7 @@ subroutine shift_ice (ntrcr,    ncat,        &
                                       atrcrn(:,n), aicen(n),    &
                                       vicen(n),    vsnon(n),    &
                                       trcr_base,   n_trcr_strata,  &
-                                      nt_strata,   trcrn(:,n))
+                                      nt_strata,   trcrn(:,n), Tf)
          if (icepack_warnings_aborted(subname)) return
 
       enddo                     ! ncat
@@ -773,8 +779,8 @@ subroutine cleanup_itd (dt,          ntrcr,      &
                               fpond,       fresh,      &
                               fsalt,       fhocn,      &
                               faero_ocn,   fiso_ocn,   &
-                              fzsal,                   &
-                              flux_bio,    limit_aice_in)
+                              fzsal,       &
+                              flux_bio,    Tf, limit_aice_in)
 
       integer (kind=int_kind), intent(in) :: &
          ncat  , & ! number of thickness categories
@@ -784,9 +790,12 @@ subroutine cleanup_itd (dt,          ntrcr,      &
          ntrcr , & ! number of tracers in use
          nbtrcr, & ! number of bio tracers in use
          n_aero    ! number of aerosol tracers
-
-      real (kind=dbl_kind), intent(in) :: &
-         dt        ! time step
+
+      real (kind=dbl_kind), intent(in) :: & 
+         dt        ! time step 
+
+      real (kind=dbl_kind), intent(in) :: & 
+         Tf        ! Freezing temperature
 
       real (kind=dbl_kind), dimension(0:ncat), intent(in) :: &
          hin_max   ! category boundaries (m)
@@ -934,7 +943,7 @@ subroutine cleanup_itd (dt,          ntrcr,      &
                      nt_strata,               &
                      aicen,      trcrn,       &
                      vicen,      vsnon,       &
-                     ncat,       hin_max      )
+                     ncat,       hin_max, Tf  )
          if (icepack_warnings_aborted(subname)) return
 
       endif ! aice > puny
@@ -959,7 +968,7 @@ subroutine cleanup_itd (dt,          ntrcr,      &
                                tr_aero,                     &
                                tr_pond_topo,                &
                                first_ice,    nbtrcr,        &
-                               dfzsal,       dflux_bio      )
+                               dfzsal,       dflux_bio, Tf )
 
          if (icepack_warnings_aborted(subname)) then
             write(warnstr,*) subname, 'aice:', aice
@@ -1061,7 +1070,7 @@ subroutine zap_small_areas (dt,        ntrcr,        &
                                   tr_aero,                 &
                                   tr_pond_topo, &
                                   first_ice, nbtrcr,       &
-                                  dfzsal,    dflux_bio     )
+                                  dfzsal,    dflux_bio, Tf)
 
       integer (kind=int_kind), intent(in) :: &
          ncat     , & ! number of thickness categories
@@ -1103,6 +1112,9 @@ subroutine zap_small_areas (dt,        ntrcr,        &
       real (kind=dbl_kind), dimension (:), intent(inout) :: &
          dflux_bio     ! zapped bio tracer flux from biology (mmol/m^2/s)
 
+      real (kind=dbl_kind), intent(in) :: &
+         Tf            ! Freezing temperature
+
       logical (kind=log_kind), intent(in) :: &
          tr_aero, &   ! aerosol flag
          tr_pond_topo ! pond flag
@@ -1220,7 +1232,7 @@ subroutine zap_small_areas (dt,        ntrcr,        &
             aice0 = aice0 + aicen(n)
             aicen(n) = c0
             vicen(n) = c0
-            trcrn(nt_Tsfc,n) = Tocnfrz
+            trcrn(nt_Tsfc,n) = Tf
 
       !-----------------------------------------------------------------
       ! Zap snow
@@ -2002,7 +2014,7 @@ subroutine icepack_aggregate (ncat,               &
                                    trcr_depend,        &
                                    trcr_base,          &
                                    n_trcr_strata,      &
-                                   nt_strata)
+                                   nt_strata, Tf)
 
       integer (kind=int_kind), intent(in) :: &
          ncat  , & ! number of thickness categories
@@ -2036,6 +2048,9 @@ subroutine icepack_aggregate (ncat,               &
       real (kind=dbl_kind), dimension (:), intent(out) :: &
          trcr      ! ice tracers
 
+      real (kind=dbl_kind), intent(in) :: &
+         Tf                ! freezing temperature
+
 !autodocument_end
 
       ! local variables
@@ -2097,7 +2112,7 @@ subroutine icepack_aggregate (ncat,               &
                                    atrcr,     aice,          &
                                    vice ,     vsno,          &
                                    trcr_base, n_trcr_strata, &
-                                   nt_strata, trcr)
+                                   nt_strata, trcr, Tf)   
       if (icepack_warnings_aborted(subname)) return
 
       deallocate (atrcr)

columnphysics/icepack_mechred.F90

@@ -114,7 +114,7 @@ subroutine ridge_ice (dt,          ndtd,       &
                             aparticn,    krdgn,      &
                             aredistn,    vredistn,   &
                             dardg1ndt,   dardg2ndt,  &
-                            dvirdgndt,               &
+                            dvirdgndt,   Tf,         &
                             araftn,      vraftn,     &
                             closing_flag,closing )
 
@@ -130,6 +130,9 @@ subroutine ridge_ice (dt,          ndtd,       &
          mu_rdg , & ! gives e-folding scale of ridged ice (m^.5)
          dt             ! time step
 
+      real (kind=dbl_kind), intent(in) :: &
+         Tf             ! freezing temperature
+
       real (kind=dbl_kind), dimension(0:ncat), intent(inout) :: &
          hin_max   ! category limits (m)
 
@@ -404,8 +407,8 @@ subroutine ridge_ice (dt,          ndtd,       &
                            nslyr,       n_aero,      &
                            msnow_mlt,   esnow_mlt,   &
                            maero,       miso,        &
-                           mpond,       &
-                           aredistn,    vredistn)
+                           mpond,       Tf,          &
+                           aredistn,    vredistn)    
          if (icepack_warnings_aborted(subname)) return
 
       !-----------------------------------------------------------------
@@ -1093,7 +1096,7 @@ subroutine ridge_shift (ntrcr,       dt,              &
                               nslyr,       n_aero,          &
                               msnow_mlt,   esnow_mlt,       &
                               maero,       miso,            &
-                              mpond,           &
+                              mpond,       Tf,              &
                               aredistn,    vredistn)
 
       integer (kind=int_kind), intent(in) :: &
@@ -1106,6 +1109,9 @@ subroutine ridge_shift (ntrcr,       dt,              &
       real (kind=dbl_kind), intent(in) :: &
          dt             ! time step (s)
 
+      real (kind=dbl_kind), intent(in) :: &
+         Tf             ! freezing temperature
+
       integer (kind=int_kind), dimension (:), intent(in) :: &
          trcr_depend, & ! = 0 for aicen tracers, 1 for vicen, 2 for vsnon
          n_trcr_strata  ! number of underlying tracer layers
@@ -1167,7 +1173,7 @@ subroutine ridge_shift (ntrcr,       dt,              &
 
       real (kind=dbl_kind), dimension(:), intent(inout) :: &
          miso           ! isotope mass added to ocean (kg m-2)
-
+ 
       real (kind=dbl_kind), dimension (:), intent(inout), optional :: &
          aredistn   , & ! redistribution function: fraction of new ridge area
          vredistn       ! redistribution function: fraction of new ridge volume
@@ -1580,7 +1586,7 @@ subroutine ridge_shift (ntrcr,       dt,              &
                                        atrcrn(:,n), aicen(n),      &
                                        vicen(n),    vsnon(n),      &
                                        trcr_base,   n_trcr_strata, &
-                                       nt_strata,   trcrn(:,n))
+                                       nt_strata,   trcrn(:,n), Tf)
          if (icepack_warnings_aborted(subname)) return
       enddo
 
@@ -1749,11 +1755,14 @@ subroutine icepack_step_ridge (dt,           ndtd,         &
                                     araftn,       vraftn,        &
                                     aice,         fsalt,         &
                                     first_ice,    fzsal,         &
-                                    flux_bio,     closing )
+                                    flux_bio,     closing, Tf )
 
       real (kind=dbl_kind), intent(in) :: &
          dt           ! time step
 
+      real (kind=dbl_kind), intent(in) :: &
+         Tf           ! freezing temperature
+
       integer (kind=int_kind), intent(in) :: &
          ncat  , & ! number of thickness categories
          ndtd  , & ! number of dynamics supercycles
@@ -1892,7 +1901,7 @@ subroutine icepack_step_ridge (dt,           ndtd,         &
                       aparticn,     krdgn,          &
                       aredistn,     vredistn,       &
                       dardg1ndt,    dardg2ndt,      &
-                      dvirdgndt,                    &
+                      dvirdgndt,    Tf,             &
                       araftn,       vraftn,         &
                       l_closing_flag,               &
                       l_closing )
@@ -1924,8 +1933,8 @@ subroutine icepack_step_ridge (dt,           ndtd,         &
                         fpond,                fresh,            &
                         fsalt,                fhocn,            &
                         faero_ocn,            l_fiso_ocn,       &
-                        fzsal,            &
-                        flux_bio)
+                        fzsal,                &
+                        flux_bio, Tf)
       if (icepack_warnings_aborted(subname)) return
 
       if (present(fiso_ocn)) fiso_ocn = l_fiso_ocn

columnphysics/icepack_parameters.F90

@@ -171,6 +171,7 @@ module icepack_parameters
       character(len=char_len), public :: &
          tfrz_option  = 'mushy'   ! form of ocean freezing temperature
                                   ! 'minus1p8' = -1.8 C
+                                  ! 'constant' = Tocnfrz from namelist
                                   ! 'linear_salt' = -depressT * sss
                                   ! 'mushy' conforms with ktherm=2
 
@@ -589,6 +590,7 @@ subroutine icepack_init_parameters(   &
         character(len=*), intent(in), optional :: &
              tfrz_option_in              ! form of ocean freezing temperature
                                          ! 'minus1p8' = -1.8 C
+                                         ! 'constant' = Tocnfrz from namelist
                                          ! 'linear_salt' = -depressT * sss
                                          ! 'mushy' conforms with ktherm=2
 
@@ -1313,6 +1315,7 @@ subroutine icepack_query_parameters(   &
       character(len=*), intent(out), optional :: &
          tfrz_option_out              ! form of ocean freezing temperature
                                       ! 'minus1p8' = -1.8 C
+                                      ! 'constant' = Tocnfrz from namelist
                                       ! 'linear_salt' = -depressT * sss
                                       ! 'mushy' conforms with ktherm=2
 

columnphysics/icepack_therm_itd.F90

@@ -101,11 +101,11 @@ subroutine linear_itd (ncat,        hin_max,     &
                              ntrcr,       trcr_depend, &
                              trcr_base,   n_trcr_strata,&
                              nt_strata,                &
-                             aicen_init,  vicen_init,  &
-                             aicen,       trcrn,       &
-                             vicen,       vsnon,       &
-                             aice,        aice0,       &
-                             fpond                     )
+                             aicen_init,  vicen_init,  & 
+                             aicen,       trcrn,       & 
+                             vicen,       vsnon,       & 
+                             aice,        aice0,       & 
+                             fpond,       Tf           )
 
       integer (kind=int_kind), intent(in) :: &
          ncat    , & ! number of thickness categories
@@ -127,6 +127,9 @@ subroutine linear_itd (ncat,        hin_max,     &
       integer (kind=int_kind), dimension (:,:), intent(in) :: &
          nt_strata      ! indices of underlying tracer layers
 
+      real (kind=dbl_kind), intent(in) :: &
+         Tf             ! freezing temperature
+
       real (kind=dbl_kind), dimension(:), intent(in) :: &
          aicen_init, & ! initial ice concentration (before vertical thermo)
          vicen_init    ! initial ice volume               (m)
@@ -595,7 +598,7 @@ subroutine linear_itd (ncat,        hin_max,     &
                          aicen,    trcrn,       &
                          vicen,    vsnon,       &
                          hicen,    donor,       &
-                         daice,    dvice        )
+                         daice,    dvice, Tf    )
          if (icepack_warnings_aborted(subname)) return
 
          ! maintain qsno negative definiteness
@@ -2177,7 +2180,7 @@ subroutine icepack_step_therm2 (dt, ncat, nltrcr,           &
                              vsnon,                 &
                              aice      ,         &
                              aice0     ,         &
-                             fpond       )
+                             fpond, Tf       )
             if (icepack_warnings_aborted(subname)) return
 
          endif ! aice > puny
@@ -2287,7 +2290,7 @@ subroutine icepack_step_therm2 (dt, ncat, nltrcr,           &
                         fpond,                fresh,            &
                         fsalt,                fhocn,            &
                         faero_ocn,            l_fiso_ocn,       &
-                        fzsal,                flux_bio)
+                        fzsal,                flux_bio, Tf)   
       if (icepack_warnings_aborted(subname)) return
 
       if (present(fiso_ocn)) then

columnphysics/icepack_therm_shared.F90

@@ -8,7 +8,7 @@ module icepack_therm_shared
 
       use icepack_kinds
 
-      use icepack_parameters, only: c0, c1, c2, c4, p5, pi, puny
+      use icepack_parameters, only: c0, c1, c2, c4, p5, pi, puny, Tocnfrz
       use icepack_parameters, only: cp_ocn, cp_ice, rhoi, rhos, Tffresh, TTTice, qqqice
       use icepack_parameters, only: stefan_boltzmann, emissivity, Lfresh, Tsmelt
       use icepack_parameters, only: saltmax, min_salin, depressT
@@ -409,7 +409,7 @@ function icepack_sea_freezing_temperature(sss) result(Tf)
 
         else
 
-           Tf = -1.8_dbl_kind
+           Tf = Tocnfrz
 
         endif