From 88a1401e47231895286efe0dd5475b79528e24d0 Mon Sep 17 00:00:00 2001
From: Ed Safford <62339196+EdwardSafford-NOAA@users.noreply.github.com>
Date: Mon, 18 Mar 2024 16:09:49 -0400
Subject: [PATCH] Ref #125 (#126)

Correct missing count value to 0 not rmiss (-999.0).
---
 .../image_gen/src/radmon_ig_time.fd/time.f90                | 6 +++++-
 .../nwprod/nam_radmon/fix/nam_radmon_satype.txt             | 2 +-
 .../nwprod/radmon_shared/sorc/verf_radbcoef.fd/bcoef.f90    | 3 +--
 .../nwprod/radmon_shared/sorc/verf_radbcor.fd/bcor.f90      | 2 +-
 .../nwprod/radmon_shared/sorc/verf_radtime.fd/time.f90      | 2 +-
 5 files changed, 9 insertions(+), 6 deletions(-)

diff --git a/src/Radiance_Monitor/image_gen/src/radmon_ig_time.fd/time.f90 b/src/Radiance_Monitor/image_gen/src/radmon_ig_time.fd/time.f90
index 5263ca6e..e9a1321f 100755
--- a/src/Radiance_Monitor/image_gen/src/radmon_ig_time.fd/time.f90
+++ b/src/Radiance_Monitor/image_gen/src/radmon_ig_time.fd/time.f90
@@ -190,7 +190,11 @@ program gatime
                 tot_cnt = tot_cnt + cnt(ftyp,cyc,j,k)
              end do
 
-             chi(ftyp,j,k) = (tot_pen/tot_cnt)
+             if( tot_cnt > 0.0 ) then
+                 chi(ftyp,j,k) = (tot_pen/tot_cnt)
+             else
+                 chi(ftyp,j,k) = rmiss
+             end if
          end do
       end do
    end do
diff --git a/src/Radiance_Monitor/nwprod/nam_radmon/fix/nam_radmon_satype.txt b/src/Radiance_Monitor/nwprod/nam_radmon/fix/nam_radmon_satype.txt
index cc4eb4ba..e56f10d9 100644
--- a/src/Radiance_Monitor/nwprod/nam_radmon/fix/nam_radmon_satype.txt
+++ b/src/Radiance_Monitor/nwprod/nam_radmon/fix/nam_radmon_satype.txt
@@ -1,2 +1,2 @@
-amsua_n18 amsua_n19 atms_npp hirs4_n19 amsua_metop-b iasi_metop-b mhs_metop-b mhs_n19
+amsua_n15 amsua_n18 amsua_n19 atms_npp hirs4_n19 amsua_metop-b iasi_metop-b mhs_metop-b mhs_n19
 
diff --git a/src/Radiance_Monitor/nwprod/radmon_shared/sorc/verf_radbcoef.fd/bcoef.f90 b/src/Radiance_Monitor/nwprod/radmon_shared/sorc/verf_radbcoef.fd/bcoef.f90
index 05dd18c8..cd6f6e82 100755
--- a/src/Radiance_Monitor/nwprod/radmon_shared/sorc/verf_radbcoef.fd/bcoef.f90
+++ b/src/Radiance_Monitor/nwprod/radmon_shared/sorc/verf_radbcoef.fd/bcoef.f90
@@ -218,7 +218,6 @@ program bcoef
 
 !       If observation was assimilated, accumulate sum
         if (data_chan(j)%errinv > 1.e-6) then
-!           pen        =  data_chan(j)%errinv*(data_chan(j)%omgbc)**2
            pen        =  (data_chan(j)%errinv*(data_chan(j)%omgbc))**2
            count(j)   = count(j) + 1.0 
            penalty(j) = penalty(j) + pen
@@ -238,7 +237,7 @@ program bcoef
      if (count(j)>0) then
         penalty(j)=penalty(j)/count(j)
      else
-        count(j)=rmiss
+        count(j)=0.0
         penalty(j)=rmiss
      endif
   end do
diff --git a/src/Radiance_Monitor/nwprod/radmon_shared/sorc/verf_radbcor.fd/bcor.f90 b/src/Radiance_Monitor/nwprod/radmon_shared/sorc/verf_radbcor.fd/bcor.f90
index 267109ea..96adb537 100755
--- a/src/Radiance_Monitor/nwprod/radmon_shared/sorc/verf_radbcor.fd/bcor.f90
+++ b/src/Radiance_Monitor/nwprod/radmon_shared/sorc/verf_radbcor.fd/bcor.f90
@@ -428,7 +428,7 @@ program bcor
         if (count(j,k)>0) then
            penalty(j,k)=penalty(j,k)/count(j,k)
         else
-           count(j,k)=rmiss
+           count(j,k) = 0.0
            penalty(j,k)=rmiss
         endif
      end do
diff --git a/src/Radiance_Monitor/nwprod/radmon_shared/sorc/verf_radtime.fd/time.f90 b/src/Radiance_Monitor/nwprod/radmon_shared/sorc/verf_radtime.fd/time.f90
index 3b8f6144..9ad63ca9 100755
--- a/src/Radiance_Monitor/nwprod/radmon_shared/sorc/verf_radtime.fd/time.f90
+++ b/src/Radiance_Monitor/nwprod/radmon_shared/sorc/verf_radtime.fd/time.f90
@@ -383,7 +383,7 @@ program time
               endif
 
            else
-              count(j,k)=rmiss
+              count(j,k)=0.0
               penalty(j,k)=rmiss
            endif