da_crtm_direct.inc

1 #ifdef CRTM
2 subroutine da_crtm_direct(  nsensor, nchan, nprof, Atmosphere,   &
3                             Surface,      &
4                             GeometryInfo, &
5                             ChannelInfo,  &
6                             RTSolution)
7 
8    integer, intent(in)            :: nsensor, nchan, nprof
9    type (CRTM_RTSolution_type ),  intent(inout)  :: RTSolution(nchan,nprof)
10    type (CRTM_ChannelInfo_type),  intent(in)  :: ChannelInfo(nsensor)
11    type( CRTM_Atmosphere_type ),  intent(in)  :: Atmosphere(nprof)
12    type( CRTM_Surface_type ),     intent(in)  :: Surface(nprof)
13    type( CRTM_GeometryInfo_type ),intent(inout)  :: GeometryInfo(nprof)
14 
15    integer :: Error_Status
16 
17    if (trace_use) call da_trace_entry("da_crtm_direct")
18 
19          Error_Status = CRTM_Forward (Atmosphere,   &
20                             Surface,      &
21                             GeometryInfo, &
22                             ChannelInfo,  &
23                             RTSolution)
24 
25          if ( Error_Status /= 0 ) then
26               call da_error(__FILE__,__LINE__, &
27                  (/"Error in calling CRTM_Forward"/))
28          end if
29 
30   if (trace_use) call da_trace_exit("da_crtm_direct")
31 
32 end subroutine da_crtm_direct
33 #endif