da_satem.f90

1 module da_satem
2 
3    use module_domain, only : domain
4    
5    use da_control, only : obs_qc_pointer,max_ob_levels,missing_r, &
6       v_interp_p, v_interp_h, check_max_iv_print, trace_use_dull, &
7       missing, max_error_uv, max_error_t, rootproc, kts,kte, &
8       max_error_p,max_error_q, check_max_iv_unit,check_max_iv,  &
9       max_stheight_diff,missing_data,max_error_bq,max_error_slp, &
10       max_error_bt, max_error_buv, satem,max_error_thickness
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       maxmin_type
14    use da_par_util, only : da_proc_stats_combine
15    use da_par_util1, only : da_proc_sum_int
16    use da_physics, only : da_tv_profile, da_thickness, da_find_layer, da_thickness_adj, &
17       da_find_layer_adj, da_tv_profile_adj, da_find_layer_tl, &
18       da_thickness_tl, da_tv_profile_tl
19    use da_statistics, only : da_stats_calculate
20    use da_tools, only : da_max_error_qc, da_residual
21    use da_tracing, only : da_trace_entry, da_trace_exit
22 
23    ! The "stats_satem_type" is ONLY used locally in da_satem:
24 
25    type residual_satem1_type
26       real          :: thickness                ! satem thickness
27    end type residual_satem1_type
28 
29    type maxmin_satem_stats_type
30       type (maxmin_type)         :: thickness
31    end type maxmin_satem_stats_type
32 
33    type stats_satem_type
34       type (maxmin_satem_stats_type)  :: maximum, minimum
35       type (residual_satem1_type)     :: average, rms_err
36    end type stats_satem_type
37 
38 
39 contains
40 
41 #include "da_ao_stats_satem.inc"
42 #include "da_jo_and_grady_satem.inc"
43 #include "da_residual_satem.inc"
44 #include "da_oi_stats_satem.inc"
45 #include "da_print_stats_satem.inc"
46 #include "da_transform_xtoy_satem.inc"
47 #include "da_transform_xtoy_satem_adj.inc"
48 #include "da_check_max_iv_satem.inc"
49 #include "da_get_innov_vector_satem.inc"
50 #include "da_calculate_grady_satem.inc"
51 
52 end module da_satem
53