da_crtm.f90

1 module da_crtm
2 
3    !---------------------------------------------------------------------------
4    ! Purpose: module for CRTM radiance data assimilation. 
5    !---------------------------------------------------------------------------
6 
7    use module_domain, only : xpose_type, x_type, xb_type
8 #ifdef CRTM
9    use module_radiance, only : CRTM_RTSolution_type,CRTM_ChannelInfo_type, &
10       CRTM_Atmosphere_type, CRTM_Surface_type,CRTM_GeometryInfo_type, &
11       CRTM_Adjoint,CRTM_Forward,CRTM_K_Matrix,CRTM_Tangent_Linear, &
12       CRTM_Allocate_Atmosphere,H2O_ID,GRAUPEL_CLOUD,ICE_CLOUD,HAIL_CLOUD, &
13       INVALID_WMO_SENSOR_ID,NEW_SNOW,rain_cloud,snow_cloud,O3_ID, GRASS_SOIL, &
14       WMO_AMSRE, WATER_CLOUD, WMO_AMSUB, WMO_AMSUA,WMO_SSMI, Sensor_Descriptor, &
15       crtm_destroy_atmosphere,crtm_sensor_name, &
16       crtm_allocate_surface,crtm_destroy_surface,crtm_assign_atmosphere, &
17       crtm_assign_surface,crtm_zero_surface,CRTM_Zero_Atmosphere, satinfo, &
18       time_slots,crtm_platform_name, crtm_init, &
19       rttov_inst_name,rttov_platform_name
20 #endif
21 
22    use da_control, only : trace_use, crtm_cloud, gravity,stdout, biascorr, &
23       biasprep, qc_rad,missing_r,rtminit_sensor,rtminit_nsensor, filename_len, &
24       use_error_factor_rad,read_biascoef, analysis_date,time_window_max, &
25       time_window_min,num_fgat_time,rtminit_platform, &
26       rtminit_satid, global
27    use da_define_structures, only : y_type, ob_type
28    use da_interpolation, only : da_interp_lin_2d_new, &
29       da_interp_lin_2d,da_interp_lin_2d_adj
30    use da_radiance1, only : da_biasprep,da_detsurtyp,da_biascorr, &
31       da_get_time_slots,da_biasprep,da_read_biascoef
32 #ifdef CRTM
33    use da_radiance1, only : da_qc_crtm
34 #endif
35 
36    use da_reporting, only : da_error,message
37    use da_tools1, only : da_free_unit, da_get_unit
38    use da_tracing, only : da_trace_entry, da_trace_exit
39 
40 #ifdef CRTM
41     TYPE (CRTM_ChannelInfo_type), allocatable, save :: ChannelInfo(:)
42 #endif
43 
44 contains
45 
46 !#include "da_transform_xtoy_crtmk.inc"
47 !#include "da_transform_xtoy_crtmk_f.inc"
48 !#include "da_transform_xtoy_crtmk_f_adj.inc"
49 !#include "da_transform_xtoy_crtmk_adj.inc"
50 #include "da_transform_xtoy_crtm.inc"
51 #include "da_transform_xtoy_crtm_adj.inc"
52 !#include "da_get_innov_vector_crtmk.inc"
53 #include "da_get_innov_vector_crtm.inc"
54 #include "da_crtm_tl.inc"
55 !#include "da_crtm_k.inc"
56 #include "da_crtm_direct.inc"
57 #include "da_crtm_ad.inc"
58 #include "da_crtm_init.inc"
59 #include "da_crtm_sensor_descriptor.inc"
60 
61 end module da_crtm
62