da_qscat.f90

1 module da_qscat 
2 
3    use module_domain, only : domain
4    
5    use da_control, only : obs_qc_pointer,max_ob_levels,missing_r, &
6       check_max_iv_print, check_max_iv_unit, v_interp_p, v_interp_h, &
7       check_max_iv, missing, max_error_uv, max_error_t, rootproc, &
8       qscat, max_error_p,max_error_q, trace_use_dull, &
9       max_stheight_diff,missing_data,max_error_bq,max_error_slp, &
10       max_error_bt, max_error_buv, anal_type_verify, kms,kme,kts,kte
11    use da_define_structures, only : maxmin_type, iv_type, y_type, jo_type, &
12       bad_data_type, x_type, number_type, bad_data_type
13    use da_interpolation, only : da_to_zk, &
14       da_interp_lin_3d,da_interp_lin_3d_adj
15    use da_par_util, only : da_proc_stats_combine
16    use da_par_util1, only : da_proc_sum_int
17    use da_statistics, only : da_stats_calculate
18    use da_tools, only : da_max_error_qc, da_residual, da_convert_zk
19    use da_tracing, only : da_trace_entry, da_trace_exit
20 
21    ! The "stats_qscat_type" is ONLY used locally in da_qscat:
22 
23    type residual_qscat1_type
24       real          :: u                        ! u-wind.
25       real          :: v                        ! v-wind.
26    end type residual_qscat1_type
27 
28    type maxmin_qscat_stats_type
29       type (maxmin_type)         :: u, v
30    end type maxmin_qscat_stats_type
31 
32    type stats_qscat_type
33       type (maxmin_qscat_stats_type)  :: maximum, minimum
34       type (residual_qscat1_type)     :: average, rms_err
35    end type stats_qscat_type
36 
37 contains
38 
39 #include "da_jo_and_grady_qscat.inc"
40 #include "da_residual_qscat.inc"
41 #include "da_check_max_iv_qscat.inc"
42 #include "da_get_innov_vector_qscat.inc"
43 #include "da_ao_stats_qscat.inc"
44 #include "da_oi_stats_qscat.inc"
45 #include "da_print_stats_qscat.inc"
46 #include "da_transform_xtoy_qscat.inc"
47 #include "da_transform_xtoy_qscat_adj.inc"
48 #include "da_calculate_grady_qscat.inc"
49 
50 end module da_qscat