Fix virtual temperature

ADCollard · Oct 10, 2024 · 0fe60bb · 0fe60bb
1 parent d9e7c54
commit 0fe60bb
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/gsi/setupt.f90 b/src/gsi/setupt.f90
@@ -610,6 +610,8 @@ subroutine setupt(obsLL,odiagLL,lunin,mype,bwork,awork,nele,nobs,is,conv_diagsav
      prsltmp2 = exp(prsltmp)
      call tintrp2a1(ges_prsi,prsitmp,dlat,dlon,dtime,hrdifsig,&
           nsig+1,mype,nfldsig)
+     call tintrp2a1(ges_tv,tvtmp,dlat,dlon,dtime,hrdifsig,&
+          nsig,mype,nfldsig)
      call tintrp2a1(ges_tsen,ttmp,dlat,dlon,dtime,hrdifsig,&
           nsig,mype,nfldsig)
      call tintrp2a1(ges_q,qtmp,dlat,dlon,dtime,hrdifsig,&
@@ -741,6 +743,8 @@ subroutine setupt(obsLL,odiagLL,lunin,mype,bwork,awork,nele,nobs,is,conv_diagsav
           nsig+1,mype,nfldsig)
      call tintrp2a1(ges_tsen,ttmp,dlat,dlon,dtime,hrdifsig,&
           nsig,mype,nfldsig)
+     call tintrp2a1(ges_tv,tvtmp,dlat,dlon,dtime,hrdifsig,&
+          nsig,mype,nfldsig)
      call tintrp2a1(ges_q,qtmp,dlat,dlon,dtime,hrdifsig,&
           nsig,mype,nfldsig)
      call tintrp2a1(ges_u,utmp,dlat,dlon,dtime,hrdifsig,&