Add FY3-E MWSFY3E obtype (#248)

* Adding MWSFY3E support to opsinputs. Merged from OPS branch 17656_mwsfy3e

* add cx writer fields

* Update Ops_SubTypeNameToNum.inc

Accidentally overwrote MWSSG!

* Update Ops_SubTypeNameToNum.inc

preserve alphabetical order

* Add test copied from mwsfy3

* add unit test netcdf

---------

Co-authored-by: Matt Shin <[email protected]>
Co-authored-by: Chris Thomas <[email protected]>

Author: 3 people

deps/ops/stubs/OpsMod_CXGenerate/Ops_GetDefaultCxFields.inc

@@ -211,7 +211,7 @@ SELECT CASE (ObsGroup)
                        StashItem_qcl,StashItem_p,StashItem_p_surface,          &
                        StashCode_u10,StashCode_v10,StashCode_t2,StashCode_rh2, &
                        StashCode_pmsl/)
-  CASE (ObsGroupMWSFY3)
+  CASE (ObsGroupMWSFY3, ObsGroupMWSFY3E)
       CxFields(1:14) = (/StashItem_theta,StashItem_q,StashItem_qcf,StashItem_SST,  &
                          StashItem_SeaIce,StashItem_orog,StashItem_qcl,            &
                          StashItem_p,StashItem_p_Surface,StashCode_u10,            &

deps/ops/stubs/OpsMod_Control/Ops_ReadObsGroupSwitches.inc

@@ -93,6 +93,7 @@ Using(ObsGroupSurfaceCloud) = ops_env_is_true ("OPS_SURFACECLOUD")
 Using(ObsGroupMWSFY3B) = ops_env_is_true ("OPS_MWSFY3B")
 Using(ObsGroupRainAccum) = ops_env_is_true ("OPS_RAINACCUM")
 Using(ObsGroupMWSFY3) = ops_env_is_true ("OPS_MWSFY3")
+Using(ObsGroupMWSFY3E) = ops_env_is_true ("OPS_MWSFY3E")
 Using(ObsGroupMWRI) = ops_env_is_true ("OPS_MWRI")
 Using(ObsGroupSAPHIR) = ops_env_is_true ("OPS_SAPHIR")
 Using(ObsGroupIN3DIClr) = ops_env_is_true ("OPS_IN3DICLR")

deps/ops/stubs/OpsMod_ODB/OpsFn_ObsToReportType.inc

@@ -719,6 +719,24 @@ SELECT CASE (Obs % Header % ObsGroup)
         END IF
       END DO
     END IF
+  CASE (ObsGroupMWSFY3E)
+    IF (Obs % Header % Satid % Present) THEN
+      DO i = ob_start, ob_end
+        IF (Obs % Satid(i) == SatId_FY3E) THEN
+          report_types(i - ob_start + 1) = rep_type_mwts3_fy3e
+        ELSE
+          found_report_type = .FALSE.
+          DO j = 1, SIZE (map, DIM = 1)
+            IF (map(j,1) == ObsGroupMWSFY3E .AND. map(j,2) == Obs % Satid(i)) THEN
+              report_types(i - ob_start + 1) = map(j,3)
+              found_report_type = .TRUE.
+              EXIT
+            END IF
+          END DO
+          CALL update_unrecognized_satids (Obs % Satid(i))
+        END IF
+      END DO
+    END IF
   CASE (ObsGroupOcean)
     IF (Obs % Header % ObsType % Present) THEN
       DO i = ob_start, ob_end

deps/ops/stubs/OpsMod_ODB/OpsODB_InitODBElemDesp_Body.inc

@@ -42,6 +42,7 @@ USE OpsMod_ObsGroupInfo, ONLY: &
     ObsGroupSAPHIR,            &
     ObsGroupMTSATImClr,        &
     ObsGroupMWSFY3,            &
+    ObsGroupMWSFY3E,           &
     ObsGroupMWRI,              &
     ObsGroupOcean,             &
     ObsGroupSatSST,            &
@@ -317,7 +318,7 @@ IF (ObsGroup == ObsGroupTOVS) THEN
   ElemDesp2ODBDesp(IndexMSU_SSU_BRGTS_TMPR) % Varno              = ops_odb_getval (h, "$rawbt_amsu")  ! brightness temperature
 END IF
 
-IF (ObsGroup == ObsGroupMWSFY3) THEN
+IF (ObsGroup == ObsGroupMWSFY3 .OR. ObsGroup == ObsGroupMWSFY3E) THEN
   ElemDesp2ODBDesp(IndexBRGTS_TMPR1) % TableName              = 'body'
   ElemDesp2ODBDesp(IndexBRGTS_TMPR1) % ColName                = 'initial_obsvalue'
   ElemDesp2ODBDesp(IndexBRGTS_TMPR1) % Varno                  = ops_odb_getval (h, "$rawbt_mwts")  ! brightness temperature

deps/ops/stubs/OpsMod_Varobs/Ops_GetDefaultVarfields.inc

@@ -219,7 +219,7 @@ SELECT CASE (ObsGroup)
                         Varfield_solzenith,Varfield_tskin,Varfield_tcozone, &
                         Varfield_surface,Varfield_numchans,Varfield_channum, &
                         Varfield_stratt,Varfield_iremiss,Varfield_qcinfo/)
-  CASE (ObsGroupMWSFY3)
+  CASE (ObsGroupMWSFY3, ObsGroupMWSFY3E)
     Varfields(1:13) = (/VarField_britemp,VarField_tskin,VarField_mwemiss, &
                         VarField_satzenith,VarField_surface, Varfield_stratt, &
                         Varfield_numchans,Varfield_channum,VarField_satid, &

deps/ops/stubs/OpsMod_Varobs/Ops_SetupVarobsLevDepC.inc

@@ -70,6 +70,7 @@ USE OpsMod_ObsGroupInfo, ONLY: &
     ObsGroupRainAccum,         &
     ObsGroupMWRI,              &
     ObsGroupMWSFY3,            &
+    ObsGroupMWSFY3E,           &
     ObsGroupGMIlow,            &
     ObsGroupGMIhigh,           &
     ObsGroupHLOSWIND,          &
@@ -311,7 +312,7 @@ SELECT CASE (Observations % header % ObsGroup)
     ObsLevelType = 0.0
     NumLevelLevs = 1
 
-  CASE (ObsGroupMWSFY3)
+  CASE (ObsGroupMWSFY3,ObsGroupMWSFY3E )
 
     ObsLevelType = 0.0
     NumLevelLevs = 1

deps/ops/stubs/Ops_Constants/OpsFn_IsSatRadGroup.inc

@@ -38,6 +38,7 @@ IF (obsgroup == ObsGroupABIClr .OR. &
     obsgroup == ObsGroupSAPHIR .OR. &
     obsgroup == ObsGroupMVIRIClr .OR. &
     obsgroup == ObsGroupMWSFY3 .OR. &
+    obsgroup == ObsGroupMWSFY3E .OR. &
     obsgroup == ObsGroupMWRI .OR. &
     obsgroup == ObsGroupSEVIRIClr .OR. &
     obsgroup == ObsGroupSEVIRIASR .OR. &

deps/ops/stubs/Ops_Constants/OpsFn_ObsGroupNameToNum.inc

@@ -90,6 +90,8 @@ SELECT CASE (ObsGroupName)
     OpsFn_ObsGroupNameToNum = ObsGroupMWRI
   CASE ("MWSFY3")
     OpsFn_ObsGroupNameToNum = ObsGroupMWSFY3
+  CASE ("MWSFY3E")
+    OpsFn_ObsGroupNameToNum = ObsGroupMWSFY3E
   CASE ("OceanSound")
     OpsFn_ObsGroupNameToNum = ObsGroupOcean
   CASE ("OceanWinds")

deps/ops/stubs/Ops_Constants/OpsFn_ObsGroupNumToName.inc

@@ -85,6 +85,8 @@ SELECT CASE (ObsGroup)
     OpsFn_ObsGroupNumToName = "MVIRIClr"
   CASE (ObsGroupMWSFY3)
     OpsFn_ObsGroupNumToName = "MWSFY3"
+  CASE (ObsGroupMWSFY3E)
+    OpsFn_ObsGroupNumToName = "MWSFY3E"
   CASE (ObsGroupMWRI)
     OpsFn_ObsGroupNumToName = "MWRI"
   CASE (ObsGroupOcean)

deps/ops/stubs/Ops_Constants/OpsMod_ObsGroupInfo.f90

@@ -80,7 +80,8 @@ MODULE OpsMod_ObsGroupInfo
 INTEGER, PARAMETER :: ObsGroupFCIClr         =  69 ! => FCI clear-sky radiances
 INTEGER, PARAMETER :: ObsGroupSternaMWR      =  70 ! => AWS & Sterna microwave radiances
 INTEGER, PARAMETER :: ObsGroupEPSMWS         =  71 ! => EPS SG MWS 
-INTEGER, PARAMETER :: max_obs_group_num      =  71 ! Number of observation groups
+INTEGER, PARAMETER :: ObsGroupMWSFY3E        =  72 ! => FY3E microwave sounder data
+INTEGER, PARAMETER :: max_obs_group_num      =  72 ! Number of observation groups
 
 CONTAINS