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    if (trace_use) call da_trace_entry("da_crtm_sensor_descriptor")  
11 
12   do i=1,nsensor
13      
14      platform_id  = rtminit_platform(i)
15      satellite_id = rtminit_satid(i)
16      sensor_id    = rtminit_sensor(i)
17 
18      if (trim(crtm_platform_name(platform_id)) == 'eos') then
19         if (satellite_id == 2) crtm_sat='aqua'
20         if (satellite_id == 1) crtm_sat='terra'
21      else
22         write(crtm_sat, '(a,i2.2)')  &
23              trim( crtm_platform_name(platform_id) ),satellite_id
24      end if
25 
26      if ( trim(crtm_sensor_name(sensor_id)) == 'airs' ) then
27         crtm_sensor='airs281SUBSET'
28      elseif ( trim(crtm_sensor_name(sensor_id)) == 'hirs' ) then
29         if (satellite_id <= 14) crtm_sensor='hirs2'
30         if (satellite_id >= 15 .and. satellite_id <= 17) crtm_sensor='hirs3'
31         if (satellite_id == 18) crtm_sensor='hirs4'
32      elseif ( trim(crtm_sensor_name(sensor_id)) == 'avhrr' ) then
33         if (satellite_id <= 14) crtm_sensor='avhrr2'
34         if (satellite_id >= 15 .and. satellite_id <= 17) crtm_sensor='avhrr3'
35         if (satellite_id == 18) crtm_sensor='avhrr4'
36      else
37         crtm_sensor=crtm_sensor_name(sensor_id)
38      end if
39 
40      sensor_descriptor(i)=trim(crtm_sensor)//'_'//trim(crtm_sat)
41 
42   end do
43 
44    if (trace_use) call da_trace_exit("da_crtm_sensor_descriptor")  
45 
46 end subroutine da_crtm_sensor_descriptor
47 #endif