da_crtm_sensor_descriptor.inc

1 #ifdef CRTM
2 subroutine da_crtm_sensor_descriptor(nsensor,sensor_descriptor)
3 
4  integer,        intent(in)  :: nsensor
5  character(len=80), intent(inout) :: sensor_descriptor(nsensor)
6 
7  integer :: i,platform_id,satellite_id,sensor_id
8  character (len=80) :: crtm_sat, crtm_sensor
9 
10   do i=1,nsensor
11      
12      platform_id  = rtminit_platform(i)
13      satellite_id = rtminit_satid(i)
14      sensor_id    = rtminit_sensor(i)
15 
16      if (trim(crtm_platform_name(platform_id)) == 'eos') then
17         if (satellite_id == 2) crtm_sat='aqua'
18         if (satellite_id == 1) crtm_sat='terra'
19      else
20         write(crtm_sat, '(a,i2.2)')  &
21              trim( crtm_platform_name(platform_id) ),satellite_id
22      endif
23 
24      if ( trim(crtm_sensor_name(sensor_id)) == 'airs' ) then
25         crtm_sensor='airs281SUBSET'
26      elseif ( trim(crtm_sensor_name(sensor_id)) == 'hirs' ) then
27         if (satellite_id <= 14) crtm_sensor='hirs2'
28         if (satellite_id >= 15 .and. satellite_id <= 17) crtm_sensor='hirs3'
29         if (satellite_id == 18) crtm_sensor='hirs4'
30      elseif ( trim(crtm_sensor_name(sensor_id)) == 'avhrr' ) then
31         if (satellite_id <= 14) crtm_sensor='avhrr2'
32         if (satellite_id >= 15 .and. satellite_id <= 17) crtm_sensor='avhrr3'
33         if (satellite_id == 18) crtm_sensor='avhrr4'
34      else
35         crtm_sensor=crtm_sensor_name(sensor_id)
36      endif
37 
38      sensor_descriptor(i)=trim(crtm_sensor)//'_'//trim(crtm_sat)
39 
40   enddo
41      write(6,*) sensor_descriptor
42 
43 end subroutine da_crtm_sensor_descriptor
44 #endif