From 3c014c4e33f78253ed939f6ccc1f3799d116fdb4 Mon Sep 17 00:00:00 2001 From: gmao-yzhu Date: Wed, 10 Mar 2021 13:55:10 -0500 Subject: [PATCH] Closes #70. Commit changes for PBL height --- src/Applications/GSI_App/GSI_GridComp.rc.tmpl | 4 ++-- src/Applications/GSI_App/GSIsa_Gaussian.F90 | 19 ++++++++++++------- 2 files changed, 14 insertions(+), 9 deletions(-) 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