Merge remote-tracking branch 'origin/release-v4.5.2'
[WRF.git] / var / da / da_ssmi / da_ssmi.f90
bloba41bb54a92a32b97b3851165f9b3b1f68831e61b
1 module da_ssmi
4 use module_domain, only : xb_type, domain
5 use module_ssmi, only : cal_sigma_v,tb,effht,epsalt,spemiss,tbatmos,&
6 roughem,effang,tbatmos,filter
8 use da_control, only : obs_qc_pointer,max_ob_levels,missing_r, &
9 v_interp_p, v_interp_h, check_max_iv_print, t_roughem, t_kelvin, &
10 missing, max_error_uv, max_error_t, rootproc, pi, trace_use, &
11 max_error_p,max_error_q, check_max_iv_unit,check_max_iv, &
12 max_stheight_diff,missing_data,max_error_bq,max_error_slp, &
13 max_error_bt, max_error_buv, max_error_thickness, mkz, &
14 max_error_rh,max_error_tb, max_error_pw, trace_use_dull, &
15 test_transforms,stdout, use_ssmiretrievalobs, use_ssmitbobs, &
16 global, print_detail_obs, max_ssmi_rv_input, max_ssmi_tb_input, &
17 its,ite,jts,jte,kts,kte,kms,kme,ids,ide,jds,jde,fails_error_max, &
18 ssmi_tb, ssmi_rv, num_ob_indexes, ssmt1, ssmt2, ob_vars,qcstat_conv_unit, &
19 write_rej_obs_conv
20 use da_define_structures, only : maxmin_type, iv_type, y_type, jo_type, &
21 bad_data_type, x_type, number_type, bad_data_type, &
22 maxmin_type,residual_ssmi_rv_type, &
23 residual_ssmi_tb_type, model_loc_type, info_type, field_type, &
24 count_obs_number_type
25 use da_interpolation, only : da_to_zk,da_interp_lin_2d, da_interp_lin_2d_adj, &
26 da_interp_lin_3d,da_interp_lin_3d_adj
27 use da_par_util, only : da_proc_stats_combine
28 use da_par_util1, only : da_proc_sum_int
29 use da_reporting, only : da_warning, message, da_error
30 use da_statistics, only : da_stats_calculate
31 use da_tools, only : da_max_error_qc, da_residual, da_llxy, da_convert_zk,da_get_print_lvl
32 use da_tools_serial, only : da_get_unit, da_free_unit
33 use da_tracing, only : da_trace_entry, da_trace_exit
35 ! The "stats_ssmi_rv_type" is ONLY used locally in da_ssmi_rv:
37 type maxmin_ssmi_rv_stats_type
38 type (maxmin_type) :: tpw ! Toatl precipitable water cm
39 type (maxmin_type) :: Speed ! Wind speed (m/s)
40 end type maxmin_ssmi_rv_stats_type
42 type stats_ssmi_rv_type
43 type (maxmin_ssmi_rv_stats_type) :: maximum, minimum
44 type (residual_ssmi_rv_type) :: average, rms_err
45 end type stats_ssmi_rv_type
47 ! The "stats_ssmi_tb_type" is ONLY used locally in da_ssmi_tb:
49 type maxmin_ssmi_tb_stats_type
50 type (maxmin_type) :: tb19v ! brightness temperature (K)
51 type (maxmin_type) :: tb19h ! brightness temperature (K)
52 type (maxmin_type) :: tb22v ! brightness temperature (K)
53 type (maxmin_type) :: tb37v ! brightness temperature (K)
54 type (maxmin_type) :: tb37h ! brightness temperature (K)
55 type (maxmin_type) :: tb85v ! brightness temperature (K)
56 type (maxmin_type) :: tb85h ! brightness temperature (K)
57 end type maxmin_ssmi_tb_stats_type
59 type stats_ssmi_tb_type
60 type (maxmin_ssmi_tb_stats_type) :: maximum, minimum
61 type (residual_ssmi_tb_type) :: average, rms_err
62 end type stats_ssmi_tb_type
65 contains
67 #include "da_ao_stats_ssmi_rv.inc"
68 #include "da_ao_stats_ssmi_tb.inc"
69 #include "da_read_obs_ssmi.inc"
70 #include "da_scan_obs_ssmi.inc"
71 #include "da_jo_and_grady_ssmi_rv.inc"
72 #include "da_jo_and_grady_ssmi_tb.inc"
73 #include "da_residual_ssmi_rv.inc"
74 #include "da_residual_ssmi_tb.inc"
75 #include "da_oi_stats_ssmi_rv.inc"
76 #include "da_oi_stats_ssmi_tb.inc"
77 #include "da_transform_xtospeed.inc"
78 #include "da_transform_xtospeed_lin.inc"
79 #include "da_transform_xtospeed_adj.inc"
80 #include "da_transform_xtoseasfcwind.inc"
81 #include "da_transform_xtoseasfcwind_lin.inc"
82 #include "da_transform_xtoseasfcwind_adj.inc"
83 #include "da_transform_xtotb.inc"
84 #include "da_transform_xtotb_lin.inc"
85 #include "da_transform_xtotb_adj.inc"
86 #include "da_transform_xtoy_ssmi_rv.inc"
87 #include "da_transform_xtoy_ssmi_rv_adj.inc"
88 #include "da_transform_xtoy_ssmi_tb.inc"
89 #include "da_transform_xtoy_ssmi_tb_adj.inc"
90 #include "da_transform_xtozrhoq.inc"
91 #include "da_transform_xtozrhoq_lin.inc"
92 #include "da_transform_xtozrhoq_adj.inc"
93 #include "da_jo_and_grady_ssmt1.inc"
94 #include "da_jo_and_grady_ssmt2.inc"
95 #include "da_residual_ssmt1.inc"
96 #include "da_residual_ssmt2.inc"
97 #include "da_check_max_iv_ssmi_rv.inc"
98 #include "da_check_max_iv_ssmi_tb.inc"
99 #include "da_check_max_iv_ssmt1.inc"
100 #include "da_check_max_iv_ssmt2.inc"
101 #include "da_get_innov_vector_ssmi_rv.inc"
102 #include "da_get_innov_vector_ssmi_tb.inc"
103 #include "da_get_innov_vector_ssmt1.inc"
104 #include "da_get_innov_vector_ssmt2.inc"
105 #include "da_ao_stats_ssmt1.inc"
106 #include "da_ao_stats_ssmt2.inc"
107 #include "da_oi_stats_ssmt1.inc"
108 #include "da_oi_stats_ssmt2.inc"
109 #include "da_print_stats_ssmt1.inc"
110 #include "da_print_stats_ssmt2.inc"
111 #include "da_transform_xtoy_ssmt1.inc"
112 #include "da_transform_xtoy_ssmt1_adj.inc"
113 #include "da_transform_xtoy_ssmt2.inc"
114 #include "da_transform_xtoy_ssmt2_adj.inc"
115 #include "da_calculate_grady_ssmi_tb.inc"
116 #include "da_calculate_grady_ssmi_rv.inc"
117 #include "da_calculate_grady_ssmt1.inc"
118 #include "da_calculate_grady_ssmt2.inc"
120 #include "da_tb_adj.inc"
121 #include "da_sigma_v_adj.inc"
122 #include "da_effang_adj.inc"
123 #include "da_effht_adj.inc"
124 #include "da_epsalt_adj.inc"
125 #include "da_roughem_adj.inc"
126 #include "da_spemiss_adj.inc"
127 #include "da_tbatmos_adj.inc"
129 #include "da_tb_tl.inc"
130 #include "da_tbatmos_tl.inc"
131 #include "da_effht_tl.inc"
132 #include "da_roughem_tl.inc"
133 #include "da_spemiss_tl.inc"
134 #include "da_effang_tl.inc"
135 #include "da_epsalt_tl.inc"
136 #include "da_sigma_v_tl.inc"
138 end module da_ssmi