From 5db524cc941157d218b4725efbbceadd0ced7b74 Mon Sep 17 00:00:00 2001 From: "jim.jung" Date: Wed, 6 Nov 2024 20:01:58 +0000 Subject: [PATCH] Added/changed comments --- src/gsi/crtm_interface.f90 | 8 ++++++-- src/gsi/read_iasing.f90 | 7 +++++-- 2 files changed, 11 insertions(+), 4 deletions(-) diff --git a/src/gsi/crtm_interface.f90 b/src/gsi/crtm_interface.f90 index 2f2764b82..db8ae6298 100644 --- 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 index 7de2f67bb..4183c5984 100644 --- 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'