27 use crtm_channelinfo_define
, only: crtm_channelinfo_type
30 type(crtm_channelinfo_type
),
intent(inout) :: channelinfo
31 integer,
intent(in) :: nchannels,channels(nchannels)
33 integer :: temp(nchannels)
35 if(nchannels>channelinfo%n_channels)
then
36 write(6,*)
'ERROR*** tried to use more channels than sensor has'
37 write(6,*)
' ',nchannels,
' > ',channelinfo%n_channels
41 check:
do i=1,nchannels
42 if(channels(i)<1 .or. channels(i)>channelinfo%n_channels)
then
43 write(6,*)
'ERROR*** invalid channel id: ',channels(i)
44 write(6,*)
' in SELECT_CHANNELS at index ',i
47 temp(i)=channelinfo%Channel_Index(channels(i))
50 channelinfo%n_channels=nchannels
51 channelinfo%Channel_Index(1:nchannels)=temp
55 subroutine select_channels_l(channelinfo,nchannels,channels,L,igot)
60 use crtm_channelinfo_define
, only: crtm_channelinfo_type
63 type(crtm_channelinfo_type
),
intent(inout) :: channelinfo
64 integer,
intent(in) :: nchannels,channels(nchannels)
66 integer :: temp(nchannels)
67 integer :: l(nchannels)
70 if(nchannels>channelinfo%n_channels)
then
71 write(6,*)
'ERROR*** tried to use more channels than sensor has'
72 write(6,*)
' ',nchannels,
' > ',channelinfo%n_channels
78 if(channels(i)<1 .or. channels(i)>channelinfo%n_channels)
then
79 write(6,*)
'ERROR*** invalid channel id: ',channels(i)
80 write(6,*)
' in SELECT_CHANNELS at index ',i
85 temp(k)=channelinfo%Channel_Index(channels(i))
88 channelinfo%Process_Channel(channels(i))=.false.
101 end subroutine select_channels_l
subroutine select_channels(channelinfo, nchannels, channels)
select_channels() verifies channel information.