Started migrating to formal indexing.

mwaxmonsky · Apr 12, 2024 · 2cd5f29 · 2cd5f29
1 parent 3f286f7
commit 2cd5f29
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 5 deletions.
diff --git a/tj2016/tj2016_sfc_pbl_hs.F90 b/tj2016/tj2016_sfc_pbl_hs.F90
@@ -16,7 +16,7 @@ module TJ2016_sfc_pbl_hs
 
     !> \section arg_table_tj2016_sfc_pbl_hs_run  Argument Table
     !! \htmlinclude tj2016_sfc_pbl_hs_run.html
-    subroutine tj2016_sfc_pbl_hs_run(ncol, pver, gravit, cappa, rairv,                        &
+    subroutine tj2016_sfc_pbl_hs_run(ncol, pver, index_surface, gravit, cappa, rairv,                 &
         cpairv, latvap, rh2o, epsilo, rhoh2o, zvirv, ps0, etamid, dtime, clat,                &
         PS, pmid, pint, lnpint, rpdel, T, U, dudt, V, dvdt, qv, shflx, lhflx, taux, tauy,     &
         evap, dqdt_vdiff, dtdt_vdiff, dtdt_heating, Km, Ke, Tsurf, tendency_of_air_enthalpy,  &
@@ -25,8 +25,9 @@ subroutine tj2016_sfc_pbl_hs_run(ncol, pver, gravit, cappa, rairv,
     !   Input / output parameters
     !------------------------------------------------
 
-    integer,  intent(in)    :: ncol                       ! number of columns
-    integer,  intent(in)    :: pver                       ! number of vertical levels
+    integer,  intent(in)           :: ncol                ! number of columns
+    integer,  intent(in)           :: pver                ! number of vertical levels
+    integer,  intent(in)           :: index_surface       ! index of surface
 
     real(kind_phys), intent(in)    :: gravit              ! g: gravitational acceleration (m/s2)
     real(kind_phys), intent(in)    :: cappa               ! Rd/cp
@@ -198,7 +199,7 @@ subroutine tj2016_sfc_pbl_hs_run(ncol, pver, gravit, cappa, rairv,
     ! Calculate hydrostatic height za of the lowermost model level
     !==========================================================================
     do i = 1, ncol
-        dlnpint = (lnpint(i,2) - lnpint(i,1))
+        dlnpint = (lnpint(i,index_surface+1) - lnpint(i,index_surface))
         za(i)   = rairv(i,pver)/gravit*T(i,pver)*(1._kind_phys+zvirv(i,pver)*qv(i,pver))*0.5_kind_phys*dlnpint
     end do
 
@@ -369,7 +370,7 @@ subroutine tj2016_sfc_pbl_hs_run(ncol, pver, gravit, cappa, rairv,
         !--------------------------------------------------------------------------
         kv  = kf*(etamid(pver) - sigmab)/onemsig                                 ! RF coefficient at the lowest level
         do i = 1, ncol
-            dlnpint = (lnpint(i,2) - lnpint(i,1))
+            dlnpint = (lnpint(i,index_surface+1) - lnpint(i,index_surface))
             za(i)   = rairv(i,pver)/gravit*T(i,pver)*(1._kind_phys+zvirv(i,pver)*qv(i,pver))*0.5_kind_phys*dlnpint ! height of lowest full model level
             rho(i)  = pmid(i,pver)/(rairv(i,pver) * T(i,pver) *(1._kind_phys+zvirv(i,pver)*qv(i,pver)))            ! air density at the lowest level rho = p/(Rd Tv)
             taux(i) = -kv * rho(i) * UCopy(i,pver) * za(i)                                                         ! U surface momentum flux in N/m2

diff --git a/tj2016/tj2016_sfc_pbl_hs.meta b/tj2016/tj2016_sfc_pbl_hs.meta
@@ -17,6 +17,12 @@
   type = integer
   dimensions = ()
   intent = in
+[ index_surface ]
+  standard_name = vertical_index_at_surface_adjacent_layer
+  units = index
+  type = integer
+  dimensions = ()
+  intent = in
 [ gravit ]
   standard_name = standard_gravitational_acceleration
   units = m s-2