Complete use_pbl usage change

jderber-NOAA · Dec 16, 2023 · 9f6103b · 9f6103b
1 parent 33d60d8
commit 9f6103b
Show file tree

Hide file tree

Showing 6 changed files with 6 additions and 12 deletions.
diff --git a/src/gsi/calctends.f90 b/src/gsi/calctends.f90
@@ -62,6 +62,7 @@ subroutine calctends(mype,teta,pri,guess,xderivative,yderivative,tendency)
   use gsi_bundlemod, only: gsi_bundlegetpointer
 
   use mpeu_util, only: die
+  use turblmod, only: use_pbl
   implicit none
 
 ! Declare passed variables
@@ -357,7 +358,7 @@ subroutine calctends(mype,teta,pri,guess,xderivative,yderivative,tendency)
       end do
     end do  !end do k
 
-    call turbl(u,v,pri,t,teta,z,u_t,v_t,t_t,jtstart(kk),jtstop(kk))
+    if(use_pbl)call turbl(u,v,pri,t,teta,z,u_t,v_t,t_t,jtstart(kk),jtstop(kk))
 
     if(.not.wrf_nmm_regional.and..not.nems_nmmb_regional)then
       do k=1,nsig

diff --git a/src/gsi/calctends_ad.f90 b/src/gsi/calctends_ad.f90
@@ -67,6 +67,7 @@ subroutine calctends_ad(fields,fields_dt,mype)
   use mpeu_util, only: die
   use derivsmod, only: gsi_xderivative_bundle
   use derivsmod, only: gsi_yderivative_bundle
+  use turblmod, only: use_pbl
   implicit none
 
 ! Declare passed variables
@@ -356,7 +357,7 @@ subroutine calctends_ad(fields,fields_dt,mype)
       end do                         
     end if
 
-    call turbl_ad(ges_prsi(1,1,1,it),ges_tv,ges_teta(1,1,1,it),&
+    if(use_pbl)call turbl_ad(ges_prsi(1,1,1,it),ges_tv,ges_teta(1,1,1,it),&
                 u,v,pri,t,u_t,v_t,t_t,jtstart(kk),jtstop(kk))
 
     do k=nsig,1,-1               

diff --git a/src/gsi/calctends_tl.f90 b/src/gsi/calctends_tl.f90
@@ -69,6 +69,7 @@ subroutine calctends_tl(fields,fields_dt,mype)
   use mpeu_util, only: die, getindex
   use derivsmod, only: gsi_xderivative_bundle
   use derivsmod, only: gsi_yderivative_bundle
+  use turblmod, only: use_pbl
   implicit none
 
 ! Declare passed variables
@@ -474,7 +475,7 @@ subroutine calctends_tl(fields,fields_dt,mype)
       end do   !end do j
     end do  !end do k
 
-    call turbl_tl(ges_prsi(1,1,1,it),ges_tv,ges_teta(1,1,1,it),&
+    if(use_pbl)call turbl_tl(ges_prsi(1,1,1,it),ges_tv,ges_teta(1,1,1,it),&
                u,v,pri,t,u_t,v_t,t_t,jtstart(kk),jtstop(kk))
 
     if(.not.wrf_nmm_regional.and..not.nems_nmmb_regional)then

diff --git a/src/gsi/turbl.f90 b/src/gsi/turbl.f90
@@ -41,7 +41,6 @@ subroutine turbl(uges,vges,pges,tges,oges,zges,termu,termv,termt,jstart,jstop)
   use kinds,only: r_kind,i_kind
   use constants,only: zero,one,two,half,rd_over_g,rd_over_cp,grav
   use gridmod,only: lat2,lon2,nsig,nsig_hlf
-  use turblmod, only: use_pbl
   use turblmod, only: dudz,dvdz,dodz,ri,rf,zi,km,kh,sm,sh
   use turblmod, only: lmix,dudtm,dvdtm,dtdtm,rdzi,rdzl
   use turblmod, only: kar0my20
@@ -70,8 +69,6 @@ subroutine turbl(uges,vges,pges,tges,oges,zges,termu,termv,termt,jstart,jstop)
   real(r_kind) px,rdzik,rdzlk,kmrdz,khrdz,ssq,aux,l0
   integer(i_kind) i,j,k
 
-  if(.not. use_pbl)return
-
   do k=1,nsig_hlf
      do j=jstart,jstop
         do i=1,lat2

diff --git a/src/gsi/turbl_ad.f90 b/src/gsi/turbl_ad.f90
@@ -40,7 +40,6 @@ subroutine turbl_ad(pges,tges,oges,u,v,prs,t,termu,termv,termt,jstart,jstop)
   use constants, only: rd_over_cp,two,rd_over_g,half,zero,one,three,grav
   use kinds, only: r_kind,i_kind 
   use gridmod, only: lat2,lon2,nsig,nsig_hlf
-  use turblmod, only: use_pbl
   use turblmod, only: dudz,dvdz,dodz,ri,rf,kar0my20,zi,km,kh,sm,sh
   use turblmod, only: lmix,dudtm,dvdtm,dtdtm,rdzi,rdzl
   use turblmod, only: a0my20,c0my20,d0my20, &
@@ -76,8 +75,6 @@ subroutine turbl_ad(pges,tges,oges,u,v,prs,t,termu,termv,termt,jstart,jstop)
   real(r_kind):: kmaz_bck,khaz_bck,kmaz_tl,khaz_tl
   integer(i_kind) i,j,k
 
-  if(.not. use_pbl)return
-
   do i=1,lat2
      do j=jstart,jstop
 

diff --git a/src/gsi/turbl_tl.f90 b/src/gsi/turbl_tl.f90
@@ -37,7 +37,6 @@ subroutine turbl_tl(pges,tges,oges,u,v,prs,t,termu,termv,termt,jstart,jstop)
   use constants,only: rd_over_cp,two,rd_over_g,half,zero,one,three,grav
   use kinds,only: r_kind,i_kind
   use gridmod, only: lat2,lon2,nsig,nsig_hlf
-  use turblmod, only: use_pbl
   use turblmod, only: dudz,dvdz,dodz,ri,rf,kar0my20,zi,km,kh,sm,sh
   use turblmod, only: lmix,dudtm,dvdtm,dtdtm,rdzi,rdzl
   use turblmod, only: a0my20,c0my20,d0my20,f7my20,f8my20,karmy20
@@ -72,8 +71,6 @@ subroutine turbl_tl(pges,tges,oges,u,v,prs,t,termu,termv,termt,jstart,jstop)
   integer(i_kind) i,j,k
   integer(i_kind),dimension(nsig):: lssq
 
-  if(.not. use_pbl)return
-
   do j=jstart,jstop
      do i=1,lat2