da_crtm_k.inc

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