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, ep_type, vp_type
10 use module_state_description, only : dyn_em,dyn_em_tl,p_g_qv
11
12 !#ifdef DM_PARALLEL
13 ! use mpi, only : mpi_barrier
14 !#endif
15
16 use da_airep, only : da_calculate_grady_airep, da_ao_stats_airep, &
17 da_oi_stats_airep, da_get_innov_vector_airep, da_residual_airep, &
18 da_jo_and_grady_airep
19 use da_airsr , only : da_calculate_grady_airsr, da_ao_stats_airsr, &
20 da_oi_stats_airsr, da_get_innov_vector_airsr, da_residual_airsr, &
21 da_jo_and_grady_airsr
22 use da_bogus, only : da_calculate_grady_bogus, da_ao_stats_bogus, &
23 da_oi_stats_bogus, da_get_innov_vector_bogus, da_residual_bogus, &
24 da_jo_and_grady_bogus
25 use da_buoy , only : da_calculate_grady_buoy, da_ao_stats_buoy, &
26 da_oi_stats_buoy,da_get_innov_vector_buoy, da_residual_buoy, &
27 da_jo_and_grady_buoy
28 use da_control, only : trace_use,var4d_coupling_disk_simul, &
29 var4d, rootproc,jcdfi_use,var4d_coupling,ierr,comm,num_fgat_time, &
30 stdout, eps, stats_unit, test_dm_exact, global, var4d_multi_inc, &
31 calculate_cg_cost_fn,anal_type_randomcv,cv_size_domain,je_factor, &
32 jb_factor,ntmax,omb_add_noise,write_iv_rad_ascii,use_obs_errfac, &
33 rtm_option,rtm_option_rttov, rtm_option_crtm, anal_type_verify, &
34 write_filtered_rad,omb_set_rand,use_rad,var_scaling2,var_scaling1, &
35 var_scaling4,var_scaling5,var_scaling3, jo_unit, &
36 print_detail_grad,omb_set_rand,grad_unit,cost_unit, &
37 cv_size_domain_je,cv_size_domain_jb, num_pseudo, &
38 sound, sonde_sfc, synop, profiler, gpsref, gpspw, polaramv, geoamv, ships, metar, &
39 satem, radar, ssmi_rv, ssmi_tb, ssmt1, ssmt2, airsr, pilot, airep, &
40 bogus, buoy, qscat,pseudo, radiance, monitor_on
41
42 #ifdef CRTM
43 use da_crtm, only : da_get_innov_vector_crtm
44 !da_get_innov_vector_crtmk
45 #endif
46 use da_define_structures, only : iv_type, y_type, j_type, be_type, &
47 xbx_type, jo_type, da_allocate_y,da_zero_x,da_deallocate_y
48 use da_obs, only : da_transform_xtoy_adj,da_transform_xtoy, &
49 da_add_noise_to_ob,da_random_omb_all
50 use da_geoamv, only : da_calculate_grady_geoamv, da_ao_stats_geoamv, &
51 da_oi_stats_geoamv, da_get_innov_vector_geoamv,da_residual_geoamv, &
52 da_jo_and_grady_geoamv
53 use da_gpspw, only : da_calculate_grady_gpspw, da_ao_stats_gpspw, &
54 da_oi_stats_gpspw, da_get_innov_vector_gpspw, da_residual_gpspw, &
55 da_jo_and_grady_gpspw
56 use da_gpsref, only : da_calculate_grady_gpsref, da_ao_stats_gpsref, &
57 da_oi_stats_gpsref, da_get_innov_vector_gpsref, da_residual_gpsref, &
58 da_jo_and_grady_gpsref
59 use da_obs_io, only : da_final_write_y, da_write_y, da_final_write_obs, &
60 da_write_obs,da_write_obs_etkf,da_write_noise_to_ob, da_use_obs_errfac, &
61 da_write_iv_for_multi_inc, da_read_iv_for_multi_inc
62 use da_metar, only : da_calculate_grady_metar, da_ao_stats_metar, &
63 da_oi_stats_metar, da_get_innov_vector_metar, da_residual_metar, &
64 da_jo_and_grady_metar
65 use da_pilot, only : da_calculate_grady_pilot, da_ao_stats_pilot, &
66 da_oi_stats_pilot, da_get_innov_vector_pilot, da_residual_pilot, &
67 da_jo_and_grady_pilot
68 use da_par_util, only : da_system,da_cv_to_global
69 use da_par_util1, only : da_proc_sum_real,da_proc_sum_ints
70 use da_polaramv, only : da_calculate_grady_polaramv, da_ao_stats_polaramv, &
71 da_oi_stats_polaramv, da_get_innov_vector_polaramv, da_residual_polaramv, &
72 da_jo_and_grady_polaramv
73 use da_profiler, only : da_calculate_grady_profiler, da_ao_stats_profiler, &
74 da_oi_stats_profiler,da_get_innov_vector_profiler, da_residual_profiler, &
75 da_jo_and_grady_profiler
76 use da_pseudo, only : da_calculate_grady_pseudo, da_ao_stats_pseudo, &
77 da_oi_stats_pseudo, da_get_innov_vector_pseudo, da_residual_pseudo, &
78 da_jo_and_grady_pseudo
79 use da_qscat, only : da_calculate_grady_qscat, da_ao_stats_qscat, &
80 da_oi_stats_qscat, da_get_innov_vector_qscat, da_residual_qscat, &
81 da_jo_and_grady_qscat
82 use da_radiance, only : da_calculate_grady_rad, da_write_filtered_rad
83 use da_radiance1, only : da_ao_stats_rad,da_oi_stats_rad, &
84 da_write_iv_rad_ascii,da_residual_rad,da_jo_and_grady_rad
85 use da_radar, only : da_calculate_grady_radar, da_ao_stats_radar, &
86 da_oi_stats_radar, da_get_innov_vector_radar, da_residual_radar, &
87 da_jo_and_grady_radar
88 use da_reporting, only : da_message, da_warning
89 use da_rttov, only : da_get_innov_vector_rttov
90 use da_satem, only : da_calculate_grady_satem, da_ao_stats_satem, &
91 da_oi_stats_satem, da_get_innov_vector_satem, da_residual_satem, &
92 da_jo_and_grady_satem
93 use da_ships, only : da_calculate_grady_ships, da_ao_stats_ships, &
94 da_oi_stats_ships, da_get_innov_vector_ships, da_residual_ships, &
95 da_jo_and_grady_ships
96 use da_sound, only : da_calculate_grady_sound,da_calculate_grady_sonde_sfc, &
97 da_ao_stats_sound, da_oi_stats_sound,da_oi_stats_sound, &
98 da_oi_stats_sonde_sfc,da_ao_stats_sonde_sfc,da_get_innov_vector_sound, &
99 da_get_innov_vector_sonde_sfc,da_jo_and_grady_sound, da_residual_sound, &
100 da_jo_and_grady_sound,da_jo_and_grady_sonde_sfc,da_residual_sonde_sfc
101 use da_ssmi, only : da_calculate_grady_ssmi_tb,da_calculate_grady_ssmi_rv,da_calculate_grady_ssmt1, &
102 da_calculate_grady_ssmt2, da_ao_stats_ssmi_tb ,da_ao_stats_ssmt2, &
103 da_ao_stats_ssmt2, da_oi_stats_ssmt1, da_oi_stats_ssmt2, &
104 da_oi_stats_ssmi_tb,da_oi_stats_ssmi_rv,da_ao_stats_ssmt1,da_get_innov_vector_ssmi_tb, &
105 da_get_innov_vector_ssmi_rv, da_residual_ssmi_rv, da_residual_ssmi_tb, &
106 da_get_innov_vector_ssmt1,da_get_innov_vector_ssmt2, &
107 da_jo_and_grady_ssmt1, da_jo_and_grady_ssmt2,da_jo_and_grady_ssmi_tb, &
108 da_jo_and_grady_ssmi_rv, &
109 da_residual_ssmt1,da_residual_ssmt2, da_ao_stats_ssmi_rv
110 use da_synop, only : da_calculate_grady_synop, da_ao_stats_synop, &
111 da_oi_stats_synop, da_get_innov_vector_synop, da_residual_synop, &
112 da_jo_and_grady_synop
113 use da_statistics, only : da_analysis_stats
114 use da_tools_serial, only : da_get_unit,da_free_unit
115 use da_tracing, only : da_trace_entry, da_trace_exit,da_trace
116 use da_transfer_model, only : da_transfer_wrftltoxa,da_transfer_xatowrftl, &
117 da_transfer_xatowrftl_adj,da_setup_firstguess,da_transfer_wrftltoxa_adj
118 use da_vtox_transforms, only : da_transform_vtox,da_transform_vtox_adj
119 use da_wrf_interfaces, only : wrf_dm_bcast_real, wrf_get_dm_communicator
120 use da_wrfvar_io, only : da_med_initialdata_input
121 use module_symbols_util, only : wrfu_finalize
122
123 implicit none
124
125 private :: da_dot, da_dot_cv
126
127 contains
128
129 #include "da_calculate_j.inc"
130 #include "da_jo_and_grady.inc"
131 #include "da_calculate_residual.inc"
132 #include "da_get_var_diagnostics.inc"
133 #include "da_get_innov_vector.inc"
134 #include "da_dot.inc"
135 #include "da_dot_cv.inc"
136 #include "da_write_diagnostics.inc"
137 #include "da_minimise_cg.inc"
138 #include "da_calculate_grady.inc"
139 #include "da_transform_vtoy.inc"
140 #include "da_transform_vtoy_adj.inc"
141
142 end module da_minimisation