Merge branch 'IASI-NG' of https://github.com/wx20jjung/GSI into IASI-NG

NOAA-EMC · Nov 7, 2024 · 7518ed8 · 7518ed8
2 parents 493e8e6 + 5db524c
commit 7518ed8
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 4 deletions.
diff --git a/src/gsi/crtm_interface.f90 b/src/gsi/crtm_interface.f90
@@ -667,8 +667,12 @@ subroutine init_crtm(init_pass,mype_diaghdr,mype,nchanl,nreal,isis,obstype,radmo
    error_status = crtm_channelinfo_subset(channelinfo(1), &
         channel_subset = nuchan(subset_start:subset_end))
 
-!TODO JAJ the correct one    else if (channelinfo(1)%sensor_id(1:7) == 'iasi-ng' .AND. isis(1:7) == 'iasi-ng') then
-else if (channelinfo(1)%sensor_id(1:3) == '999' .AND. isis(1:7) == 'iasi-ng') then  ! TODO To be removed.
+! TODO The CRTM spectral coefficient files have the instrument name in the beginning of the file.  The current iasi-ng coefficient
+! TODO file contains '999' instead of the instrument name.  When the final coefficient file is built, it will have 'iasi-ng'.
+! TODO  else if (channelinfo(1)%sensor_id(1:7) == 'iasi-ng' .AND. isis(1:7) == 'iasi-ng') then
+! TODO when this file exists, use the above line.
+else if (channelinfo(1)%sensor_id(1:3) == '999' .AND. isis(1:7) == 'iasi-ng') then
+! TODO and remove the above line.
    sensorindex = 1
    subset_start = 0
    subset_end = 0

diff --git a/src/gsi/read_iasing.f90 b/src/gsi/read_iasing.f90
@@ -297,9 +297,12 @@ subroutine read_iasing(mype,val_iasing,ithin,isfcalc,rmesh,jsatid,gstime,&
 
 !  find IASI-NG sensorindex
   sensorindex_iasing = 0
+!TODO The final iasi-ng spectral coefficient file is not available yet.  The current file has '999' in place of iasi-ng. 
+!TODO The line below is a temporary fix to work around this problem in the CRTM spectral coefficient file. 
   if ( sc(1)%sensor_id(1:7) == 'iasi-ng' .or. sc(1)%sensor_id == '999') then
-!TODO This is the correct code once the CRTM spectral coefficient is fixed.          
-!TODO JAJ  if ( sc(1)%sensor_id(1:7) == 'iasi-ng' ) then   
+!TODO Remove the above line when the correct CRTM file is available
+!TODO   if ( sc(1)%sensor_id(1:7) == 'iasi-ng' ) then
+!TODO Use the above line once the CRTM spectral coefficient file is fixed 
      sensorindex_iasing = 1
   else
      write(6,*)'READ_IASI-NG: ***ERROR*** sensorindex_iasi-ng not set  NO IASI-NG DATA USED'