da_crtm_tl.inc

1 #ifdef CRTM
2 subroutine da_crtm_tl(  nchan, Atmosphere,   &
3                             Surface,      &
4                             Atmosphere_TL,&
5                             Surface_TL,   &
6                             GeometryInfo, &
7                             ChannelInfo,  &
8                             RTSolution,   &
9                             RTSolution_TL)
10 
11    integer, intent(in)            :: nchan
12    TYPE (CRTM_RTSolution_type ),  intent(inout)  :: RTSolution(nchan),RTSolution_TL(nchan)
13    TYPE (CRTM_ChannelInfo_type),  intent(in)  :: ChannelInfo
14    TYPE( CRTM_Atmosphere_type ),  intent(inout)  :: Atmosphere, Atmosphere_TL
15    TYPE( CRTM_Surface_type ),     intent(in)  :: Surface,Surface_TL
16    TYPE( CRTM_GeometryInfo_type ),intent(inout)  :: GeometryInfo
17 
18    integer :: Error_Status
19 
20    if (trace_use) call da_trace_entry("da_crtm_tl")
21 
22          Error_Status = CRTM_Tangent_Linear(Atmosphere,   &
23                             Surface,      &
24                             Atmosphere_TL,&
25                             Surface_TL,   &
26                             GeometryInfo, &
27                             ChannelInfo,  &
28                             RTSolution,   &
29                             RTSolution_TL) 
30 
31          if ( Error_Status /= 0 ) THEN
32               call da_error(__FILE__,__LINE__, &
33                  (/"Error in calling CRTM_Tangent_Linear"/))
34          endif
35 
36 
37    if (trace_use) call da_trace_exit("da_crtm_tl")
38 
39 end subroutine da_crtm_tl
40 #endif