da_crtm_k.inc

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