da_minimisation.f90
References to this file elsewhere.
1 module da_minimisation
2
3 !---------------------------------------------------------------------------
4 ! Purpose: Collection of routines associated with minimisation.
5 !---------------------------------------------------------------------------
6
7 use module_configure, only : grid_config_rec_type,nl_set_dyn_opt
8 use module_dm, only : wrf_dm_sum_real
9 use module_domain, only : domain, xpose_type, x_type, ep_type, vp_type, &
10 xb_type
11 use module_state_description, only : dyn_em,dyn_em_tl
12
13 !#ifdef DM_PARALLEL
14 ! use mpi, only : mpi_barrier
15 !#endif
16
17 use da_airep, only : da_calculate_grady_airep, da_ao_stats_airep, &
18 da_oi_stats_airep, da_get_innov_vector_airep, da_residual_airep, &
19 da_jo_and_grady_airep
20 use da_airsr , only : da_calculate_grady_airsr, da_ao_stats_airsr, &
21 da_oi_stats_airsr, da_get_innov_vector_airsr, da_residual_airsr, &
22 da_jo_and_grady_airsr
23 use da_bogus, only : da_calculate_grady_bogus, da_ao_stats_bogus, &
24 da_oi_stats_bogus, da_get_innov_vector_bogus, da_residual_bogus, &
25 da_jo_and_grady_bogus
26 use da_buoy , only : da_calculate_grady_buoy, da_ao_stats_buoy, &
27 da_oi_stats_buoy,da_get_innov_vector_buoy, da_residual_buoy, &
28 da_jo_and_grady_buoy
29 use da_control, only : trace_use,var4d_coupling_disk_simul, &
30 var4d, rootproc,jcdfi_use,var4d_coupling,ierr,comm,num_fgat_time, &
31 stdout, eps, stats_unit, testing_dm_exact, global, &
32 calculate_cg_cost_fn,anal_type_randomcv,cv_size_domain,je_factor, &
33 jb_factor,ntmax,omb_add_noise,write_iv_rad_ascii,use_obs_errfac, &
34 rtm_option,rtm_option_rttov, rtm_option_crtm, use_crtm_kmatrix_fast, &
35 write_filtered_rad,omb_set_rand, num_airsr_tot,jo_unit,num_airep_tot, &
36 num_buoy_tot,num_geoamv_tot,num_gpsref_tot,num_pilot_tot, &
37 num_polaramv_tot,num_radiance_tot,num_qscat_tot,num_profiler_tot, &
38 num_metar_tot,num_gpspw_tot,num_bogus_tot,num_sound_tot,num_ssmi_tot, &
39 num_ships_tot,num_ssmt2_tot,use_radiance,var_scaling2,var_scaling1, &
40 num_synop_tot,var_scaling4,var_scaling5,var_scaling3,num_ssmt1_tot, &
41 num_satem_tot,print_detail_grad,omb_set_rand,grad_unit,cost_unit, &
42 cv_size_domain_je,cv_size_domain_jb
43 use da_crtm, only : da_get_innov_vector_crtmk, da_get_innov_vector_crtm
44 use da_define_structures, only : ob_type, y_type, j_type, be_type, &
45 xbx_type, jo_type, da_allocate_y,da_zero_x,da_deallocate_y
46 use da_obs, only : da_transform_xtoy_adj,da_transform_xtoy, &
47 da_add_noise_to_ob,da_random_omb_all
48 use da_geoamv, only : da_calculate_grady_geoamv, da_ao_stats_geoamv, &
49 da_oi_stats_geoamv, da_get_innov_vector_geoamv,da_residual_geoamv, &
50 da_jo_and_grady_geoamv
51 use da_gpspw, only : da_calculate_grady_gpspw, da_ao_stats_gpspw, &
52 da_oi_stats_gpspw, da_get_innov_vector_gpspw, da_residual_gpspw, &
53 da_jo_and_grady_gpspw
54 use da_gpsref, only : da_calculate_grady_gpsref, da_ao_stats_gpsref, &
55 da_oi_stats_gpsref, da_get_innov_vector_gpsref, da_residual_gpsref, &
56 da_jo_and_grady_gpsref
57 use da_obs_io, only : da_final_write_y, da_write_y, da_final_write_obs, &
58 da_write_obs,da_write_noise_to_ob, da_use_obs_errfac
59 use da_metar, only : da_calculate_grady_metar, da_ao_stats_metar, &
60 da_oi_stats_metar, da_get_innov_vector_metar, da_residual_metar, &
61 da_jo_and_grady_metar
62 use da_pilot, only : da_calculate_grady_pilot, da_ao_stats_pilot, &
63 da_oi_stats_pilot, da_get_innov_vector_pilot, da_residual_pilot, &
64 da_jo_and_grady_pilot
65 use da_par_util, only : da_system,da_cv_to_global
66 use da_par_util1, only : da_proc_sum_real,da_proc_sum_ints
67 use da_polaramv, only : da_calculate_grady_polaramv, da_ao_stats_polaramv, &
68 da_oi_stats_polaramv, da_get_innov_vector_polaramv, da_residual_polaramv, &
69 da_jo_and_grady_polaramv
70 use da_profiler, only : da_calculate_grady_profiler, da_ao_stats_profiler, &
71 da_oi_stats_profiler,da_get_innov_vector_profiler, da_residual_profiler, &
72 da_jo_and_grady_profiler
73 use da_pseudo, only : da_calculate_grady_pseudo, da_ao_stats_pseudo, &
74 da_oi_stats_pseudo, da_get_innov_vector_pseudo, da_residual_pseudo, &
75 da_jo_and_grady_pseudo
76 use da_qscat, only : da_calculate_grady_qscat, da_ao_stats_qscat, &
77 da_oi_stats_qscat, da_get_innov_vector_qscat, da_residual_qscat, &
78 da_jo_and_grady_qscat
79 use da_radiance, only : da_calculate_grady_rad, da_write_filtered_rad
80 use da_radiance1, only : da_oma_stats_rad,da_omb_stats_rad, &
81 da_write_iv_rad_ascii,da_residual_rad,da_jo_and_grady_rad
82 use da_radar, only : da_calculate_grady_radar, da_ao_stats_radar, &
83 da_oi_stats_radar, da_get_innov_vector_radar, da_residual_radar, &
84 da_jo_and_grady_radar
85 use da_reporting, only : da_message, da_warning
86 use da_rttov, only : da_get_innov_vector_rttov
87 use da_satem, only : da_calculate_grady_satem, da_ao_stats_satem, &
88 da_oi_stats_satem, da_get_innov_vector_satem, da_residual_satem, &
89 da_jo_and_grady_satem
90 use da_ships, only : da_calculate_grady_ships, da_ao_stats_ships, &
91 da_oi_stats_ships, da_get_innov_vector_ships, da_residual_ships, &
92 da_jo_and_grady_ships
93 use da_sound, only : da_calculate_grady_sound,da_calculate_grady_sonde_sfc, &
94 da_ao_stats_sound, da_oi_stats_sound,da_oi_stats_sound, &
95 da_oi_stats_sonde_sfc,da_ao_stats_sonde_sfc,da_get_innov_vector_sound, &
96 da_get_innov_vector_sonde_sfc,da_jo_and_grady_sound, da_residual_sound, &
97 da_jo_and_grady_sound,da_jo_and_grady_sonde_sfc,da_residual_sonde_sfc
98 use da_ssmi, only : da_calculate_grady_ssmi,da_calculate_grady_ssmt1, &
99 da_calculate_grady_ssmt2, da_ao_stats_ssmi,da_ao_stats_ssmt2, &
100 da_ao_stats_ssmt2, da_oi_stats_ssmt1, da_oi_stats_ssmt2, &
101 da_oi_stats_ssmi,da_ao_stats_ssmt1,da_get_innov_vector_ssmi, &
102 da_get_innov_vector_ssmt1,da_get_innov_vector_ssmt2, da_residual_ssmi, &
103 da_jo_and_grady_ssmt1, da_jo_and_grady_ssmt2,da_jo_and_grady_ssmi, &
104 da_residual_ssmt1,da_residual_ssmt2
105 use da_synop, only : da_calculate_grady_synop, da_ao_stats_synop, &
106 da_oi_stats_synop, da_get_innov_vector_synop, da_residual_synop, &
107 da_jo_and_grady_synop
108 use da_statistics, only : da_analysis_stats
109 use da_tools1, only : da_get_unit,da_free_unit
110 use da_tracing, only : da_trace_entry, da_trace_exit,da_trace
111 use da_transfer_model, only : da_transfer_wrftltoxa,da_transfer_xatowrftl, &
112 da_transfer_xatowrftl_adj,da_setup_firstguess,da_transfer_wrftltoxa_adj
113 use da_vtox_transforms, only : da_transform_vtox,da_transform_vtox_adj
114 use da_wrf_interfaces, only : wrf_dm_bcast_real
115 use da_wrfvar_io, only : da_med_initialdata_input,da_med_initialdata_output
116
117 implicit none
118
119 private :: da_dot, da_dot_cv
120
121 contains
122
123 #include "da_calculate_j.inc"
124 #include "da_jo_and_grady.inc"
125 #include "da_calculate_residual.inc"
126 #include "da_get_var_diagnostics.inc"
127 #include "da_get_innov_vector.inc"
128 #include "da_dot.inc"
129 #include "da_dot_cv.inc"
130 #include "da_write_diagnostics.inc"
131 #include "da_minimise_cg.inc"
132 #include "da_calculate_grady.inc"
133 #include "da_transform_vtoy.inc"
134 #include "da_transform_vtoy_adj.inc"
135
136 end module da_minimisation