da_pseudo.f90

1 module da_pseudo
2 
3    use da_control, only : obs_qc_pointer,max_ob_levels,missing_r, &
4       v_interp_p, v_interp_h, &
5       missing, max_error_uv, max_error_t, rootproc, &
6       max_error_p,max_error_q,  &
7       max_stheight_diff,missing_data,max_error_bq,max_error_slp, &
8       max_error_bt, max_error_buv, trace_use,pseudo_var,stdout
9    use da_define_structures, only : maxmin_type, ob_type, y_type, jo_type, &
10       bad_data_type, x_type, number_type, bad_data_type
11    use module_domain, only : xpose_type, xb_type
12    use da_interpolation, only :  &
13       da_interp_obs_lin_2d,da_interp_obs_lin_2d_adj
14    use da_statistics, only : da_stats_calculate
15    use da_tools, only : da_residual
16    use da_par_util, only : da_proc_stats_combine
17    use da_par_util1, only : da_proc_sum_int
18    use da_tracing, only : da_trace_entry, da_trace_exit
19 
20    ! The "stats_pseudo_type" is ONLY used locally in da_pseudo:
21 
22    type residual_pseudo1_type
23       real          :: u                        ! u-wind.
24       real          :: v                        ! v-wind.
25       real          :: t                        ! temperature
26       real          :: p                        ! pressure
27       real          :: q                        ! specific humidity
28    end type residual_pseudo1_type
29 
30    type maxmin_pseudo_stats_type
31       type (maxmin_type)         :: u, v, t, p, q
32    end type maxmin_pseudo_stats_type
33 
34    type stats_pseudo_type
35       type (maxmin_pseudo_stats_type)  :: maximum, minimum
36       type (residual_pseudo1_type)     :: average, rms_err
37    end type stats_pseudo_type
38 
39 contains
40 
41 #include "da_jo_and_grady_pseudo.inc"
42 #include "da_residual_pseudo.inc"
43 #include "da_get_innov_vector_pseudo.inc"
44 #include "da_ao_stats_pseudo.inc"
45 #include "da_oi_stats_pseudo.inc"
46 #include "da_print_stats_pseudo.inc"
47 #include "da_transform_xtoy_pseudo.inc"
48 #include "da_transform_xtoy_pseudo_adj.inc"
49 #include "da_calculate_grady_pseudo.inc"
50 
51 end module da_pseudo
52