Add checking for prior calculation of virtual temperture before running 'prevents'

src/Applications/NCEP_Paqc/GMAOprev/CMakeLists.txt

@@ -3,6 +3,11 @@ ecbuild_add_executable (
   SOURCES gblevents_gmao.f prevents.f
   LIBS NCEP_bufr_r8i4 NCEP_w3_r8i4 GMAO_hermes)
 
+ecbuild_add_executable (
+  TARGET check_virtmp.x
+  SOURCES check_virtmp.f
+  LIBS NCEP_bufr_r4i4)
+
 string(REPLACE " " ";" tmp ${FREAL8})
 foreach (flag ${tmp})
    target_compile_options (gmao_prevents.x PRIVATE $<$<COMPILE_LANGUAGE:Fortran>:${flag}>)

src/Applications/NCEP_Paqc/GMAOprev/check_virtmp.f

@@ -0,0 +1,58 @@
+        implicit none
+        real*8 hdr(1), tpc(1,20)
+        real*8 bmiss, getbmiss
+        CHARACTER(len=180) filn
+        CHARACTER(len=8)  SUBSET
+        integer  iret, idate, ilev, j, tvflg, lubfi
+        integer nargs, iargc
+        real vtcd
+
+        data lubfi /10/
+
+        nargs = iargc()
+        if (nargs .lt. 1) then
+           print *, 'Usage: check_virtmp.x prepbufr-file'
+           stop
+        end if
+        call getarg(1,filn)
+        open(unit=lubfi,file=filn,form='unformatted')
+
+        CALL OPENBF(LUBFI,'IN',LUBFI)
+        call ufbqcd(lubfi,'VIRTMP',vtcd)
+        tvflg = 1
+        subset = ''
+        iret = 0
+        do while (iret == 0 .and. subset /= 'SFCSHP')
+           CALL READMG(LUBFI,SUBSET,IDATE,IRET)
+        end do
+
+
+C  LOOP THROUGH THE INPUT MESSAGES - READ THE NEXT SUBSET
+C  ------------------------------------------------------
+
+        do while(iret .eq. 0) 
+          CALL READSB(LUBFI,IRET)
+          IF(IRET .eq. 0) then
+            call ufbint(lubfi,hdr,1,1,ilev,'TYP')
+            call ufbevn(lubfi,tpc,1,1,20,ilev,'TPC')
+            do j = 1,20
+               if(tpc(1,j) == vtcd) tvflg = 0
+               if(tpc(1,j) >= bmiss) exit
+            end do
+            if (tvflg == 0) exit
+          else
+              subset = ''
+              iret = 0
+              do while (iret == 0 .and. subset /= 'SFCSHP')
+                 CALL READMG(LUBFI,SUBSET,IDATE,IRET)
+              end do
+          endif
+        end do
+        if (tvflg == 0) then
+           print *,'VIRTUAL'
+        else
+           print *,'DRY'
+        end if
+
+        stop
+        end

src/Applications/NCEP_Paqc/oiqc/gmao_prepqc

@@ -26,6 +26,9 @@
 #                       can be found by DAS, few other tweaks
 #  18Oct2019  Meta      Removed reference to 'Shell'; cleaned out 
 #                        code for old ACQC, ACARSQC
+#  16Nov2021  Meta      Add routine to check for previous virtual temperature
+#                        calculation for SFCSHP in input prepbufr and skip
+#                        recalculation of VIRTMP if found
 #------------------------------------------------------------------
 
 # make env vars readily available
@@ -302,13 +305,30 @@ sub gmaoprevents {
    Assign("prepdv.$nymd.$hh",       51);
    Assign("prevents.out.$nymd.$hh", 52);
 
+   $VTEMP = `check_virtmp.x $pref`;
+   open(PRVPARM, '>prepobs_prevents.parm');
+
+   if ( $VTEMP  =~ /DRY/) {     # input has sensible temperature   
+       print "$0: PREVENTS calculate VIRTMP\n";
+       print PRVPARM " &PREVDATA DOVTMP= TRUE, DOFCST= TRUE DOBERR= TRUE / \n";
+   } else  {                   # already has virtual temperature
+       print "$0: PREVENTS without VIRTMP\n";
+       print PRVPARM " &PREVDATA DOVTMP= FALSE, RECALC_Q= FALSE, DOFCST= TRUE DOBERR= TRUE  / \n";
+   }
+
+   close(PRVPARM);
+
 
    $FORT_CONVERT12 = "BIG_ENDIAN";
    $FORT_CONVERT13 = "BIG_ENDIAN";
 
+# Run check for virtual temperature
+# ---------------------------------
+
+
 # Run prevents
 # -------------
-  $cmd = "gmao_prevents.x $dynf < $rcdir/prepobs_prevents.merra.parm";
+  $cmd = "gmao_prevents.x $dynf < prepobs_prevents.parm";
   print "$0: $cmd\n" unless ( $opt_q );
 #  $rc = System ( $cmd, "/dev/null" ) unless ( $opt_n ) ;    
  $rc = system ( $cmd ) unless ( $opt_n ) ;