da_test.f90
References to this file elsewhere.
1 module da_test
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
10 #ifdef DM_PARALLEL
11 use module_dm, only : local_communicator, &
12 ntasks_x, ntasks_y, data_order_xyz, mytask, &
13 ntasks, data_order_xy
14 use module_comm_dm, only : halo_psichi_uv_adj_sub, halo_xa_sub, &
15 halo_sfc_xa_sub, halo_ssmi_xa_sub
16 use da_control, only : ips,ipe,jds,jde,jps,jpe,kds,kde,kps,kpe
17 use mpi, only : mpi_sum
18 #endif
19
20 use da_control, only : num_procs
21 use module_domain, only : vp_type, xb_type, x_type, ep_type, &
22 domain
23 use module_state_description, only : dyn_em,dyn_em_tl,dyn_em_ad,p_a_qv
24
25 use da_control, only : trace_use,ims,ime,jms,jme,kms,kme, &
26 ierr, trace_use_dull, comm,global,stdout,rootproc, &
27 its,ite,jts,jte,kts,kte,sfc_assi_options,typical_qrn_rms, &
28 typical_u_rms,typical_v_rms,typical_w_rms,typical_t_rms, typical_p_rms, &
29 typical_q_rms,typical_qcw_rms,print_detail_testing,typical_rh_rms, &
30 fg_format, fg_format_wrf,typical_rf_rms,typical_rv_rms, &
31 typical_thickness_rms, typical_tb19v_rms,typical_tb37h_rms, &
32 typical_tb85h_rms,typical_tb37v_rms,typical_tb85v_rms,typical_tb22v_rms, &
33 typical_tb19h_rms,typical_speed_rms,typical_tpw_rms,typical_ref_rms, &
34 cv_options_hum,inv_typ_vp5_sumsq,inv_typ_vp1_sumsq, &
35 inv_typ_vp3_sumsq,inv_typ_vp2_sumsq,inv_typ_vpalpha_sumsq, &
36 inv_typ_vp4_sumsq,typical_rho_rms,balance_geo,balance_cyc,balance_type, &
37 balance_geocyc, var4d, num_fgat_time,cv_options_hum_specific_humidity, &
38 cv_options_hum_relative_humidity, ids, ide, jds, jde, kds, kde, &
39 sound, synop, profiler, gpsref, gpspw, polaramv, geoamv, ships, metar, &
40 satem, radar, ssmi_rv, ssmi_tb, ssmt1, ssmt2, airsr, pilot, airep, &
41 bogus, buoy, qscat,pseudo, use_radarobs, use_ssmiretrievalobs, &
42 use_gpsrefobs, use_ssmt1obs, use_ssmitbobs, use_ssmt2obs, use_gpspwobs
43 use da_define_structures, only : da_zero_x,da_zero_vp_type,da_allocate_y, &
44 da_deallocate_y,be_type, xbx_type, iv_type, y_type
45 use da_dynamics, only : da_uv_to_divergence,da_uv_to_vorticity, &
46 da_psichi_to_uv, da_psichi_to_uv_adj
47 use da_ffts, only : da_solve_poissoneqn_fct
48 use da_minimisation, only : da_transform_vtoy_adj,da_transform_vtoy
49 use da_obs, only : da_transform_xtoy,da_transform_xtoy_adj
50 use da_par_util, only : da_patch_to_global, da_system
51 #ifdef DM_PARALLEL
52 use da_par_util1, only : true_mpi_real
53 #endif
54 use da_physics, only : da_transform_xtopsfc,da_transform_xtopsfc_adj, &
55 da_pt_to_rho_lin,da_transform_xtotpw,da_transform_xtogpsref_lin, &
56 da_transform_xtowtq, da_transform_xtowtq_adj,da_pt_to_rho_adj, &
57 da_transform_xtotpw_adj
58 use da_reporting, only : da_error, message, da_message
59 use da_spectral, only : da_test_spectral
60 use da_ssmi, only : da_transform_xtoseasfcwind_lin, &
61 da_transform_xtoseasfcwind_adj
62 use da_statistics, only : da_correlation_coeff1d,da_correlation_coeff2d
63 use da_tools_serial, only : da_get_unit,da_free_unit
64 use da_tracing, only : da_trace_entry,da_trace_exit
65 use da_transfer_model, only : da_transfer_wrftltoxa,da_transfer_xatowrftl, &
66 da_transfer_xatowrftl_adj,da_transfer_wrftltoxa_adj
67 ! Don't use, as we pass a 3D array into a 1D one
68 ! use da_wrf_interfaces, only : wrf_dm_bcast_real
69 use da_wrf_interfaces, only : wrf_debug, wrf_shutdown
70 use da_wrfvar_io, only : da_med_initialdata_output,da_med_initialdata_input
71 use da_vtox_transforms, only : da_transform_xtotb_lin, &
72 da_transform_xtotb_adj, da_vertical_transform, da_transform_vptox, &
73 da_transform_xtogpsref_adj,da_transform_vptox_adj,da_transform_vtox, &
74 da_transform_vtox_adj,da_transform_vtovv,da_transform_vtovv_global, &
75 da_transform_vtovv_global_adj, da_transform_vtovv_adj
76
77 implicit none
78
79 contains
80
81 #include "da_check_balance.inc"
82 #include "da_check_cvtovv_adjoint.inc"
83 #include "da_check_vtox_adjoint.inc"
84 #include "da_check_vptox_adjoint.inc"
85 #include "da_check_vp_errors.inc"
86 #include "da_check_vvtovp_adjoint.inc"
87 #include "da_check_xtovptox_errors.inc"
88 #include "da_check_xtoy_adjoint.inc"
89 #include "da_check_xtoy_adjoint_airep.inc"
90 #include "da_check_xtoy_adjoint_gpspw.inc"
91 #include "da_check_xtoy_adjoint_gpsref.inc"
92 #include "da_check_xtoy_adjoint_metar.inc"
93 #include "da_check_xtoy_adjoint_pilot.inc"
94 #include "da_check_xtoy_adjoint_ssmi_rv.inc"
95 #include "da_check_xtoy_adjoint_ssmi_tb.inc"
96 #include "da_check_xtoy_adjoint_satem.inc"
97 #include "da_check_xtoy_adjoint_geoamv.inc"
98 #include "da_check_xtoy_adjoint_polaramv.inc"
99 #include "da_check_xtoy_adjoint_ships.inc"
100 #include "da_check_xtoy_adjoint_radar.inc"
101 #include "da_check_xtoy_adjoint_bogus.inc"
102 #include "da_check_xtoy_adjoint_sound.inc"
103 #include "da_check_xtoy_adjoint_sonde_sfc.inc"
104 #include "da_check_xtoy_adjoint_synop.inc"
105 #include "da_check_xtoy_adjoint_rad.inc"
106 #include "da_transform_xtovp.inc"
107 #include "da_check.inc"
108 #include "da_check_xtoy_adjoint_pseudo.inc"
109 #include "da_check_xtoy_adjoint_qscat.inc"
110 #include "da_check_xtoy_adjoint_ssmt1.inc"
111 #include "da_check_xtoy_adjoint_ssmt2.inc"
112 #include "da_check_xtoy_adjoint_profiler.inc"
113 #include "da_check_xtoy_adjoint_buoy.inc"
114 #include "da_setup_testfield.inc"
115 #include "da_check_sfc_assi.inc"
116 #include "da_check_psfc.inc"
117 #include "da_set_tst_trnsf_fld.inc"
118 #include "da_check_vtoy_adjoint.inc"
119 #include "da_get_y_lhs_value.inc"
120
121 end module da_test