remove restart flag from noahmpdrv_init

HelinWei-NOAA · Jan 10, 2024 · de58149 · de58149
1 parent cf2be29
commit de58149
Show file tree

Hide file tree

Showing 2 changed files with 15 additions and 18 deletions.
diff --git a/physics/module_soil_init.f90 b/physics/module_soil_init.f90
@@ -54,7 +54,7 @@ subroutine noahmp_soil_init(im                      , & ! in
       real (kind=kind_phys), dimension(     lsoil_lsm)                :: level_wilting_point 
       real (kind=kind_phys), dimension(     lsoil_lsm)                :: level_water 
 
-      integer :: iloc, ilev, lhave, lwant
+      integer :: iloc, ilev, lhave, lwant,styp
       real (kind=kind_phys) :: porosity, bexp, psisat, soil_matric_potential
       real (kind=kind_phys) :: supercool_water, field_capacity, wilting_point
 
@@ -95,6 +95,7 @@ subroutine noahmp_soil_init(im                      , & ! in
                  ( soil_temperature_interp(iloc,lsoil_input-1) - soil_temperature_interp(iloc,lsoil_input) ) /  &
                     ( interp_levels(lsoil_input-1) - interp_levels(lsoil_input) )
 
+
         soil_moisture_interp(iloc,0) = soil_moisture_interp(iloc,1) + &
                  ( interp_levels(0) - interp_levels(1) ) * &
                  ( soil_moisture_interp(iloc,1) - soil_moisture_interp(iloc,2) ) /  &
@@ -114,7 +115,7 @@ subroutine noahmp_soil_init(im                      , & ! in
             soil_temperature_output(iloc,lwant) = soil_temperature_interp(iloc,lsoil_input+1)
             soil_moisture_output   (iloc,lwant) = soil_moisture_interp(iloc,lsoil_input+1)
           end do
-          exit level_want1
+!         exit level_want1
         end if
 
         level_have1 : do lhave = 0 , lsoil_input
@@ -135,6 +136,7 @@ subroutine noahmp_soil_init(im                      , & ! in
         end do level_have1
       end do level_want1
 
+
 ! Some arbitrary limits to temperature
 
       where(soil_temperature_output < 200.0) soil_temperature_output = 200.0
@@ -148,10 +150,11 @@ subroutine noahmp_soil_init(im                      , & ! in
       end do
 
       do iloc = 1 , im
-
-        wilting_point  = smcwlt_table(soil_type(iloc))
-        field_capacity = smcref_table(soil_type(iloc))
-        porosity       = smcmax_table(soil_type(iloc))
+        styp = soil_type(iloc)
+        if(styp .le. 0) styp = 16
+        wilting_point  = smcwlt_table(styp)
+        field_capacity = smcref_table(styp)
+        porosity       = smcmax_table(styp)
         level_wilting_point  = level_thickness_output * wilting_point   ! meters of water
         level_field_capacity = level_thickness_output * field_capacity
         level_porosity       = level_thickness_output * porosity
@@ -182,9 +185,11 @@ subroutine noahmp_soil_init(im                      , & ! in
 ! Initialize liquid soil moisture from total soil moisture and soil temperature using Niu and Yang (2006)
 
       do iloc = 1 , im
-        porosity = smcmax_table(soil_type(iloc))
-        bexp     = bexp_table(soil_type(iloc))
-        psisat   = psisat_table(soil_type(iloc))
+        styp = soil_type(iloc)
+        if(styp .le. 0) styp = 16
+        porosity = smcmax_table(styp)
+        bexp     = bexp_table(styp)
+        psisat   = psisat_table(styp)
       do ilev = 1 , lsoil_lsm
         if(soil_temperature_output(iloc,ilev) >= temperature_freezing) then
           soil_liquid_output(iloc,ilev) = soil_moisture_output(iloc,ilev)

diff --git a/physics/noahmpdrv.F90 b/physics/noahmpdrv.F90
@@ -33,7 +33,7 @@ module noahmpdrv
 !!
       subroutine noahmpdrv_init(im,lsm, lsm_noahmp, me, isot, ivegsrc, &
                                 nlunit, pores, resid,               &
-                                lsm_cold_start,lsoil,zsi,lsoil_lsm,     &
+                                lsoil,zsi,lsoil_lsm,                &
                                 zs,soiltyp,vegtype,                 &
                                 tskin,tg3,                          &
                                 smc,slc,stc,                        &
@@ -53,7 +53,6 @@ subroutine noahmpdrv_init(im,lsm, lsm_noahmp, me, isot, ivegsrc, &
 
         real (kind=kind_phys), dimension(:), intent(out) :: pores, resid
 
-        logical,              intent(in) :: lsm_cold_start
         integer,              intent(in) :: lsoil, lsoil_lsm
         real (kind=kind_phys), dimension(:), intent(in) :: zsi,zs
 
@@ -165,13 +164,6 @@ subroutine noahmpdrv_init(im,lsm, lsm_noahmp, me, isot, ivegsrc, &
 
         ! -- call to init n-layer Noah MP soil layers from coldstart
 
-!       if ( lsoil /= lsoil_lsm) then
-
-!       call   noahmpsoilinit (lsm_cold_start, im, lsoil_lsm,lsoil,&
-!                      zsin,zsout,dzsout,tskin,tg3,smc,slc,stc,   &
-!                      sh2o,tslb,smois,soiltyp,vegtype,            &
-!                               errmsg, errflg)
-!       endif
         if ( lsoil /= lsoil_lsm) then
         call noahmp_soil_init ( im                      , & ! in
                        lsoil_lsm               , & ! in