diff --git a/io/module_write_internal_state.F90 b/io/module_write_internal_state.F90
index 51c422227..6d0602a5f 100644
--- a/io/module_write_internal_state.F90
+++ b/io/module_write_internal_state.F90
@@ -98,6 +98,7 @@ module write_internal_state
       integer                  :: ncld !< Number of hydrometeors.
       integer                  :: nsoil !< Number of soil layers.
       integer                  :: imp_physics !< Choice of microphysics scheme.
+      integer                  :: landsfcmdl !< Choice of land surface model
       integer                  :: dtp !< Physics timestep.
       real,dimension(:),allocatable :: ak !< a parameter for sigma pressure level calculations.
       real,dimension(:),allocatable :: bk !< b parameter for sigma pressure level calculations.
diff --git a/io/post_fv3.F90 b/io/post_fv3.F90
index 97962bdd9..8296b7517 100644
--- a/io/post_fv3.F90
+++ b/io/post_fv3.F90
@@ -430,6 +430,7 @@ subroutine post_getattr_fv3(wrt_int_state,grid_id)
             if (trim(attName) == 'nsoil')  wrt_int_state%nsoil=varival
             if (trim(attName) == 'fhzero') wrt_int_state%fhzero=varival
             if (trim(attName) == 'imp_physics') wrt_int_state%imp_physics=varival
+            if (trim(attName) == 'landsfcmdl') wrt_int_state%landsfcmdl=varival
           endif
         else if (typekind==ESMF_TYPEKIND_R4) then
           if(n==1) then
@@ -624,7 +625,7 @@ subroutine set_postvars_fv3(wrt_int_state,grid_id,mype,mpicomp)
 !
       imp_physics = wrt_int_state%imp_physics       !set GFS mp physics to 99 for Zhao scheme
       dtp         = wrt_int_state%dtp
-      iSF_SURFACE_PHYSICS = 2
+      iSF_SURFACE_PHYSICS = wrt_int_state%landsfcmdl 
       spval = 9.99e20
 !
 ! nems gfs has zhour defined