diff --git a/src/Applications/GSI_App/GSI_GridComp.rc.tmpl b/src/Applications/GSI_App/GSI_GridComp.rc.tmpl
index 8226e895..a4d686bf 100644
--- a/src/Applications/GSI_App/GSI_GridComp.rc.tmpl
+++ b/src/Applications/GSI_App/GSI_GridComp.rc.tmpl
@@ -74,8 +74,8 @@ post-analysis_obs_filename: %s.ana.obs.%y4%m2%d2_%h2z.ods
 ensemble_upaprgb_filename: %s.prog.eta.%y4%m2%d2_%h2%n2z.nc4
 forecast_errorB_filename: fsens.eta.nc4
 
-vars_anaeta_file: phis,ps,ts,frland,frlandice,frlake,frocean,frseaice,delp,u,v,tv,sphu,ozone,qitot,qltot,qrtot,qstot
-vars_bkgeta_file: phis,ps,ts,frland,frlandice,frlake,frocean,frseaice,u,v,tv,sphu,ozone,qitot,qltot,qrtot,qstot
+vars_anaeta_file: phis,ps,ts,PBLH,frland,frlandice,frlake,frocean,frseaice,delp,u,v,tv,sphu,ozone,qitot,qltot,qrtot,qstot
+vars_bkgeta_file: phis,ps,ts,PBLH,frland,frlandice,frlake,frocean,frseaice,u,v,tv,sphu,ozone,qitot,qltot,qrtot,qstot
 vars_bkgsfc_file: U10M,V10M,SNOWDP,GWETTOP,TSOIL1,Z0M,DCOOL,DWARM,TDROP,TS_FOUND,TDEL
 
 # used in read_guess
diff --git a/src/Applications/GSI_App/GSIsa_Gaussian.F90 b/src/Applications/GSI_App/GSIsa_Gaussian.F90
index c64202b3..7a67be4a 100644
--- a/src/Applications/GSI_App/GSIsa_Gaussian.F90
+++ b/src/Applications/GSI_App/GSIsa_Gaussian.F90
@@ -587,7 +587,7 @@ subroutine ReadFirstGuess (dynfile, sfcfile, impSt)
 ! Open file
 !----------
 
-!!  only_vars='phis,ps,uwnd,vwnd,theta,sphu,ozone,qitot,qltot')
+!!  only_vars='phis,ps,PBLH,uwnd,vwnd,theta,sphu,ozone,qitot,qltot')
 
 !I/O
 
@@ -604,6 +604,11 @@ subroutine ReadFirstGuess (dynfile, sfcfile, impSt)
    call FieldRead (unit, Efield, rc=status)
    VERIFY_(STATUS)
 
+   call ESMF_StateGetField(impSt, 'PBLH', Efield, rc=status)
+   VERIFY_(STATUS)
+   call FieldRead (unit, Efield, rc=status) !PBLH
+   VERIFY_(STATUS)
+
    call ESMF_StateGetField(impSt, 'theta', Efield, rc=status)
    VERIFY_(STATUS)
    call FieldRead (unit, Efield, rc=status) !PS3
@@ -958,15 +963,15 @@ subroutine FillBundle( BUNDLE, NVARS, stType, RC)
     real, pointer                :: lons_in(:), lats_in(:)
     real, pointer                ::  PTR2(:,:),  PTR3(:,:,:)
     character(len=ESMF_MAXSTR)   :: NAME
-    character(len=ESMF_MAXSTR) :: difnam(16), defnam(7)
-    integer                    :: difdim(16), defdim(7)
+    character(len=ESMF_MAXSTR) :: difnam(17), defnam(8)
+    integer                    :: difdim(17), defdim(8)
 
-   data difnam /'phis','ps','uwnd','vwnd','theta','sphu','ozone',&
+   data difnam /'phis','ps','PBLH','uwnd','vwnd','theta','sphu','ozone',&
               'qitot','qltot','TSKIN','U10M','V10M','SNOWDP', &
               'GWETTOP','TSOIL1','ORO'/
-   data difdim /2,2,3,3,3,3,3,3,3,2,2,2,2,2,2,2/
-   data defnam /'dps','du','dv','dt','dq','doz','dql'/
-   data defdim /2,3,3,3,3,3,3/
+   data difdim /2,2,2,3,3,3,3,3,3,3,2,2,2,2,2,2,2/
+   data defnam /'dps','dpblh','du','dv','dt','dq','doz','dql'/
+   data defdim /2,2,3,3,3,3,3,3/
 
 ! Start