Skip to content

Variational Bias Correction of TAMDAR T #904

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
liujake merged 3 commits into from
Nov 22, 2019
Merged

Variational Bias Correction of TAMDAR T #904

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
695cbf3
VarBC of TAMADR temperature observations
FengGaoGrow May 14, 2019
8cabbc0
Merge branch 'develop' into aircraft-varbc
FengGaoGrow May 14, 2019
a483b94
Update da_control.f90
liujake Nov 22, 2019
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
5 changes: 5 additions & 0 deletions Registry/registry.var
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -192,6 +192,11 @@ rconfig logical use_atmsobs namelist,wrfvar4 1 .false. - "use
rconfig logical use_mwtsobs namelist,wrfvar4 1 .false. - "use_mwtsobs" "" ""
rconfig logical use_mwhsobs namelist,wrfvar4 1 .false. - "use_mwhsobs" "" ""
rconfig logical use_mwhs2obs namelist,wrfvar4 1 .false. - "use_mwhs2obs" "" ""
rconfig logical use_varbc_tamdar namelist,wrfvar4 1 .false. - "use_varbc_tamdar" "" ""
rconfig integer varbc_tamdar_bm namelist,wrfvar4 1 1 - "For VarBC: 1-ADR; 2-WMP" "" ""
rconfig integer varbc_tamdar_nbgerr namelist,wrfvar4 1 1000 - "Scaling B-beta" "" ""
rconfig integer varbc_tamdar_nobsmin namelist,wrfvar4 1 6 - "OB count req to run VarBC""" ""
rconfig real varbc_tamdar_pred0 namelist,wrfvar4 1 1.0 - "Predictor(1)" "" ""
rconfig logical check_max_iv namelist,wrfvar5 1 .true. - "check_max_iv" "" ""
rconfig real max_error_t namelist,wrfvar5 1 5.0 - "max_error_t" "" ""
rconfig real max_error_uv namelist,wrfvar5 1 5.0 - "max_error_uv" "" ""
Expand Down
1 change: 1 addition & 0 deletions var/build/da.make
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -22,6 +22,7 @@ WRFVAR_OBJS = \
da_sound.o \
da_mtgirs.o \
da_tamdar.o \
da_varbc_tamdar.o \
da_bogus.o \
da_airep.o \
da_pilot.o \
Expand Down
7 changes: 4 additions & 3 deletions var/build/depend.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -125,7 +125,7 @@ da_interpolation.o : da_interpolation.f90 da_interp_msk_avg_2d_partial.inc da_in
da_lapack.o : da_lapack.f90 dlamch.inc dlarf.inc dlarfg.inc dlarft.inc dlarfb.inc dorg2l.inc dorg2r.inc dsytd2.inc dlatrd.inc dorgqr.inc dorgql.inc dlassq.inc dlapy2.inc dlartg.inc dlasrt.inc dlansy.inc dsytrd.inc dsterf.inc dorgtr.inc dlae2.inc dlasr.inc dlaev2.inc dlascl.inc dlanst.inc dlaset.inc iparmq.inc ieeeck.inc ilaenv.inc dsteqr.inc dsyev.inc da_blas.o
da_mat_cv3.o : da_mat_cv3.f90
da_metar.o : da_metar.f90 da_calculate_grady_metar.inc da_get_innov_vector_metar.inc da_check_max_iv_metar.inc da_transform_xtoy_metar_adj.inc da_transform_xtoy_metar.inc da_print_stats_metar.inc da_oi_stats_metar.inc da_residual_metar.inc da_jo_and_grady_metar.inc da_ao_stats_metar.inc da_tracing.o da_tools.o da_statistics.o da_physics.o da_grid_definitions.o da_par_util1.o da_par_util.o da_interpolation.o da_define_structures.o da_control.o module_domain.o
da_minimisation.o : da_minimisation.f90 da_read_basicstates.inc da_swap_xtraj.inc da_lanczos_io.inc da_kmat_mul.inc da_amat_mul.inc da_sensitivity.inc da_adjoint_sensitivity.inc da_transform_vtoy_adj.inc da_transform_vtoy.inc da_calculate_grady.inc da_minimise_lz.inc da_minimise_cg.inc da_write_diagnostics.inc da_dot_cv.inc da_dot.inc da_get_innov_vector.inc da_get_var_diagnostics.inc da_calculate_residual.inc da_jo_and_grady.inc da_calculate_gradj.inc da_calculate_j.inc da_transform_vtod_wpec.inc da_transform_vtod_wpec_adj.inc module_io_wrf.o da_4dvar.o da_lapack.o module_symbols_util.o da_wrf_interfaces.o da_vtox_transforms.o da_varbc.o da_transfer_model.o da_tracing.o da_tools_serial.o da_statistics.o da_synop.o da_ssmi.o da_sound.o da_ships.o da_satem.o da_reporting.o da_rain.o da_radar.o da_radiance1.o da_radiance.o da_tamdar.o da_mtgirs.o da_qscat.o da_pseudo.o da_profiler.o da_polaramv.o da_par_util1.o da_par_util.o da_pilot.o da_metar.o da_obs_io.o da_gpsref.o da_gpspw.o da_geoamv.o da_obs.o da_define_structures.o da_control.o da_buoy.o da_bogus.o da_airsr.o da_airep.o module_state_description.o module_domain.o module_dm.o module_configure.o da_join_iv_for_multi_inc.o da_wrfvar_io.o da_gpseph.o
da_minimisation.o : da_minimisation.f90 da_read_basicstates.inc da_swap_xtraj.inc da_lanczos_io.inc da_kmat_mul.inc da_amat_mul.inc da_sensitivity.inc da_adjoint_sensitivity.inc da_transform_vtoy_adj.inc da_transform_vtoy.inc da_calculate_grady.inc da_minimise_lz.inc da_minimise_cg.inc da_write_diagnostics.inc da_dot_cv.inc da_dot.inc da_get_innov_vector.inc da_get_var_diagnostics.inc da_calculate_residual.inc da_jo_and_grady.inc da_calculate_gradj.inc da_calculate_j.inc da_transform_vtod_wpec.inc da_transform_vtod_wpec_adj.inc module_io_wrf.o da_4dvar.o da_lapack.o module_symbols_util.o da_wrf_interfaces.o da_vtox_transforms.o da_varbc.o da_transfer_model.o da_tracing.o da_tools_serial.o da_statistics.o da_synop.o da_ssmi.o da_sound.o da_ships.o da_satem.o da_reporting.o da_rain.o da_radar.o da_radiance1.o da_radiance.o da_tamdar.o da_mtgirs.o da_qscat.o da_pseudo.o da_profiler.o da_polaramv.o da_par_util1.o da_par_util.o da_pilot.o da_metar.o da_obs_io.o da_gpsref.o da_gpspw.o da_geoamv.o da_obs.o da_define_structures.o da_control.o da_buoy.o da_bogus.o da_airsr.o da_airep.o module_state_description.o module_domain.o module_dm.o module_configure.o da_join_iv_for_multi_inc.o da_wrfvar_io.o da_gpseph.o da_varbc_tamdar.o
da_module_convert_tool.o : da_module_convert_tool.f90 da_convertor_v_interp.inc
da_module_couple_uv.o : da_module_couple_uv.f90 da_couple.inc da_calc_mu_uv.inc da_couple_uv.inc
da_module_couple_uv_ad.o : da_module_couple_uv_ad.f90 da_couple_ad.inc da_calc_mu_uv_ad.inc da_couple_uv_ad.inc da_module_couple_uv.o
Expand Down Expand Up @@ -160,7 +160,8 @@ da_spectral.o : da_spectral.f90 da_apply_power.inc da_legtra_inv_adj.inc da_vtov
da_ssmi.o : da_ssmi.f90 da_sigma_v_tl.inc da_epsalt_tl.inc da_effang_tl.inc da_spemiss_tl.inc da_roughem_tl.inc da_effht_tl.inc da_tbatmos_tl.inc da_tb_tl.inc da_tbatmos_adj.inc da_spemiss_adj.inc da_roughem_adj.inc da_epsalt_adj.inc da_effht_adj.inc da_effang_adj.inc da_sigma_v_adj.inc da_tb_adj.inc da_calculate_grady_ssmt2.inc da_calculate_grady_ssmt1.inc da_calculate_grady_ssmi_rv.inc da_calculate_grady_ssmi_tb.inc da_transform_xtoy_ssmt2_adj.inc da_transform_xtoy_ssmt2.inc da_transform_xtoy_ssmt1_adj.inc da_transform_xtoy_ssmt1.inc da_print_stats_ssmt2.inc da_print_stats_ssmt1.inc da_oi_stats_ssmt2.inc da_oi_stats_ssmt1.inc da_ao_stats_ssmt2.inc da_ao_stats_ssmt1.inc da_get_innov_vector_ssmt2.inc da_get_innov_vector_ssmt1.inc da_get_innov_vector_ssmi_tb.inc da_get_innov_vector_ssmi_rv.inc da_check_max_iv_ssmt2.inc da_check_max_iv_ssmt1.inc da_check_max_iv_ssmi_tb.inc da_check_max_iv_ssmi_rv.inc da_residual_ssmt2.inc da_residual_ssmt1.inc da_jo_and_grady_ssmt2.inc da_jo_and_grady_ssmt1.inc da_transform_xtozrhoq_adj.inc da_transform_xtozrhoq_lin.inc da_transform_xtozrhoq.inc da_transform_xtoy_ssmi_tb_adj.inc da_transform_xtoy_ssmi_tb.inc da_transform_xtoy_ssmi_rv_adj.inc da_transform_xtoy_ssmi_rv.inc da_transform_xtotb_adj.inc da_transform_xtotb_lin.inc da_transform_xtotb.inc da_transform_xtoseasfcwind_adj.inc da_transform_xtoseasfcwind_lin.inc da_transform_xtoseasfcwind.inc da_transform_xtospeed_adj.inc da_transform_xtospeed_lin.inc da_transform_xtospeed.inc da_oi_stats_ssmi_tb.inc da_oi_stats_ssmi_rv.inc da_residual_ssmi_tb.inc da_residual_ssmi_rv.inc da_jo_and_grady_ssmi_tb.inc da_jo_and_grady_ssmi_rv.inc da_scan_obs_ssmi.inc da_read_obs_ssmi.inc da_ao_stats_ssmi_tb.inc da_ao_stats_ssmi_rv.inc da_tracing.o da_tools_serial.o da_tools.o da_statistics.o da_reporting.o da_par_util1.o da_par_util.o da_interpolation.o da_define_structures.o da_control.o module_ssmi.o module_domain.o
da_statistics.o : da_statistics.f90 da_print_qcstat.inc da_stats_calculate.inc da_data_distribution.inc da_correlation_coeff2d.inc da_correlation_coeff1d.inc da_maxmin_in_field.inc da_analysis_stats.inc da_reporting.o da_tools_serial.o da_tracing.o da_par_util.o da_par_util1.o da_define_structures.o da_control.o module_domain.o
da_synop.o : da_synop.f90 da_check_buddy_synop.inc da_calculate_grady_synop.inc da_check_max_iv_synop.inc da_get_innov_vector_synop.inc da_transform_xtoy_synop_adj.inc da_transform_xtoy_synop.inc da_print_stats_synop.inc da_oi_stats_synop.inc da_residual_synop.inc da_jo_synop_uvtq.inc da_jo_and_grady_synop.inc da_ao_stats_synop.inc da_tracing.o da_tools.o da_statistics.o da_physics.o da_grid_definitions.o da_par_util.o da_par_util1.o da_interpolation.o da_define_structures.o da_control.o module_domain.o
da_tamdar.o : da_tamdar.f90 da_calculate_grady_tamdar_sfc.inc da_check_max_iv_tamdar_sfc.inc da_get_innov_vector_tamdar_sfc.inc da_transform_xtoy_tamdar_sfc_adj.inc da_transform_xtoy_tamdar_sfc.inc da_print_stats_tamdar_sfc.inc da_oi_stats_tamdar_sfc.inc da_residual_tamdar_sfc.inc da_jo_tamdar_sfc_uvtq.inc da_jo_and_grady_tamdar_sfc.inc da_ao_stats_tamdar_sfc.inc da_calculate_grady_tamdar.inc da_get_innov_vector_tamdar.inc da_check_max_iv_tamdar.inc da_transform_xtoy_tamdar_adj.inc da_transform_xtoy_tamdar.inc da_print_stats_tamdar.inc da_oi_stats_tamdar.inc da_residual_tamdar.inc da_jo_tamdar_uvtq.inc da_jo_and_grady_tamdar.inc da_ao_stats_tamdar.inc da_tracing.o da_physics.o da_grid_definitions.o da_par_util1.o da_par_util.o da_tools.o da_statistics.o da_interpolation.o module_domain.o da_define_structures.o da_control.o
da_tamdar.o : da_tamdar.f90 da_calculate_grady_tamdar_sfc.inc da_check_max_iv_tamdar_sfc.inc da_get_innov_vector_tamdar_sfc.inc da_transform_xtoy_tamdar_sfc_adj.inc da_transform_xtoy_tamdar_sfc.inc da_print_stats_tamdar_sfc.inc da_oi_stats_tamdar_sfc.inc da_residual_tamdar_sfc.inc da_jo_tamdar_sfc_uvtq.inc da_jo_and_grady_tamdar_sfc.inc da_ao_stats_tamdar_sfc.inc da_calculate_grady_tamdar.inc da_get_innov_vector_tamdar.inc da_check_max_iv_tamdar.inc da_transform_xtoy_tamdar_adj.inc da_transform_xtoy_tamdar.inc da_print_stats_tamdar.inc da_oi_stats_tamdar.inc da_residual_tamdar.inc da_jo_tamdar_uvtq.inc da_jo_and_grady_tamdar.inc da_ao_stats_tamdar.inc da_tracing.o da_physics.o da_grid_definitions.o da_par_util1.o da_par_util.o da_tools.o da_statistics.o da_interpolation.o module_domain.o da_define_structures.o da_control.o da_varbc_tamdar.o
da_varbc_tamdar.o : da_varbc_tamdar.f90 da_varbc_tamdar_init.inc da_varbc_tamdar_pred.inc da_varbc_tamdar_precond.inc da_varbc_tamdar_direct.inc da_varbc_tamdar_adj.inc da_varbc_tamdar_tl.inc da_varbc_tamdar_update.inc da_tracing.o da_tools_serial.o da_tools.o da_reporting.o da_define_structures.o da_control.o module_dm.o
da_test.o : da_test.f90 da_test_vxtransform.inc da_check_gradient.inc da_get_y_lhs_value.inc da_check_vtoy_adjoint.inc da_set_tst_trnsf_fld.inc da_check_psfc.inc da_check_sfc_assi.inc da_setup_testfield.inc da_check_xtoy_adjoint_buoy.inc da_check_xtoy_adjoint_profiler.inc da_check_xtoy_adjoint_ssmt2.inc da_check_xtoy_adjoint_ssmt1.inc da_check_xtoy_adjoint_qscat.inc da_check_xtoy_adjoint_pseudo.inc da_dot_cv.inc da_dot.inc da_check.inc da_check_gradient.inc da_transform_xtovp.inc da_check_xtoy_adjoint_rad.inc da_check_xtoy_adjoint_synop.inc da_check_xtoy_adjoint_tamdar_sfc.inc da_check_xtoy_adjoint_tamdar.inc da_check_xtoy_adjoint_mtgirs.inc da_check_xtoy_adjoint_sonde_sfc.inc da_check_xtoy_adjoint_sound.inc da_check_xtoy_adjoint_bogus.inc da_check_xtoy_adjoint_rain.inc da_check_xtoy_adjoint_radar.inc da_check_xtoy_adjoint_ships.inc da_check_xtoy_adjoint_polaramv.inc da_check_xtoy_adjoint_geoamv.inc da_check_xtoy_adjoint_satem.inc da_check_xtoy_adjoint_ssmi_tb.inc da_check_xtoy_adjoint_ssmi_rv.inc da_check_xtoy_adjoint_pilot.inc da_check_xtoy_adjoint_metar.inc da_check_xtoy_adjoint_gpsref.inc da_check_xtoy_adjoint_gpspw.inc da_check_xtoy_adjoint_airep.inc da_check_xtoy_adjoint.inc da_check_xtovptox_errors.inc da_check_vvtovp_adjoint.inc da_check_vp_errors.inc da_check_vptox_adjoint.inc da_check_vtox_adjoint.inc da_check_cvtovv_adjoint.inc da_check_balance.inc da_4dvar.o da_vtox_transforms.o da_wrfvar_io.o da_wrf_interfaces.o da_transfer_model.o da_tracing.o da_tools_serial.o da_statistics.o da_ssmi.o da_spectral.o da_reporting.o da_physics.o da_par_util1.o da_par_util.o da_obs.o da_minimisation.o da_ffts.o da_dynamics.o da_define_structures.o module_state_description.o module_domain.o da_control.o module_comm_dm.o module_dm.o module_configure.o da_rain.o da_check_dynamics_adjoint.inc da_check_xtoy_adjoint_gpseph.inc
da_tools.o : da_tools.f90 da_geo2msl1.inc da_msl2geo1.inc da_get_time_slots.inc da_get_julian_time.inc da_get_print_lvl.inc da_get_3d_sum.inc da_get_2d_sum.inc da_set_boundary_3d.inc da_set_boundary_xb.inc da_set_boundary_xa.inc da_ludcmp.inc da_lubksb.inc da_eof_decomposition.inc da_eof_decomposition_test.inc da_buddy_qc.inc da_unifva.inc da_togrid.inc da_togrid_new.inc da_smooth_anl.inc da_openfile.inc da_gaus_noise.inc da_set_randomcv.inc da_random_omb.inc da_max_error_qc.inc da_add_noise_new.inc da_add_noise.inc da_residual_new.inc da_residual.inc da_diff_seconds.inc da_mo_correction.inc da_intpsfc_tem.inc da_intpsfc_prs.inc da_sfcprs.inc da_obs_sfc_correction.inc da_1d_eigendecomposition.inc da_convert_zk.inc da_lc_cone.inc da_set_merc.inc da_map_set.inc da_map_init.inc da_set_ps.inc da_set_lc.inc da_xyll_ps.inc da_xyll_merc.inc da_xyll_lc.inc da_xyll_latlon.inc da_xyll_default.inc da_xyll.inc da_llxy_wrf_new.inc da_llxy_wrf.inc da_llxy_ps_new.inc da_llxy_ps.inc da_llxy_merc_new.inc da_llxy_merc.inc da_llxy_lc_new.inc da_llxy_lc.inc da_llxy_latlon_new.inc da_llxy_latlon.inc da_llxy_rotated_latlon.inc da_llxy_global_new.inc da_llxy_global.inc da_llxy_kma_global_new.inc da_llxy_kma_global.inc da_llxy_default_new.inc da_llxy_default.inc da_llxy_new.inc da_llxy.inc da_map_utils_defines.inc da_lapack.o da_reporting.o da_tracing.o da_tools_serial.o da_define_structures.o da_control.o module_domain.o module_dm.o module_bc.o da_sfc_hori_interp_weights.inc
da_tools_serial.o : da_tools_serial.f90 da_find_fft_trig_funcs.inc da_find_fft_factors.inc da_advance_time.inc da_advance_cymdh.inc da_array_print.inc da_change_date.inc da_free_unit.inc da_get_unit.inc da_reporting.o da_control.o
Expand Down Expand Up @@ -191,7 +192,7 @@ da_wrfvar_esmf.o : da_wrfvar_esmf.f90
da_wrfvar_esmf_super.o : da_wrfvar_esmf_super.f90 da_wrfvar_interface.inc da_esmf_finalize.inc da_esmf_run.inc da_esmf_init.inc
da_wrfvar_io.o : copyfile.c da_wrfvar_io.f90 da_med_initialdata_output_lbc.inc da_med_initialdata_output.inc da_med_initialdata_input.inc da_update_firstguess.inc da_4dvar.o da_tracing.o da_reporting.o da_control.o module_io_domain.o module_domain.o module_date_time.o module_configure.o
da_wrfvar_main.o : da_wrfvar_main.f90 da_4dvar.o da_wrfvar_top.o da_wrf_interfaces.o da_tracing.o da_control.o module_symbols_util.o
da_wrfvar_top.o : da_wrfvar_top.f90 da_solve_init.inc da_solve_dual_res_init.inc da_solve.inc da_wrfvar_finalize.inc da_wrfvar_interface.inc da_wrfvar_run.inc da_wrfvar_init2.inc da_wrfvar_init1.inc da_wrf_interfaces.o da_rain.o da_synop.o da_ssmi.o da_sound.o da_ships.o da_satem.o da_radar.o da_mtgirs.o da_qscat.o da_profiler.o da_polaramv.o da_pilot.o da_metar.o da_gpsref.o da_gpspw.o da_geoamv.o da_buoy.o da_bogus.o da_airsr.o da_airep.o da_crtm.o da_tools.o da_vtox_transforms.o da_transfer_model.o da_tracing.o da_tools_serial.o da_test.o da_setup_structures.o da_reporting.o da_varbc.o da_radiance1.o da_physics.o da_par_util.o da_obs_io.o da_obs.o da_minimisation.o da_define_structures.o da_control.o module_comm_dm.o module_dm.o module_tiles.o module_state_description.o module_radiance.o da_wrfvar_io.o da_4dvar.o module_symbols_util.o module_driver_constants.o module_domain.o module_configure.o module_io_domain.o da_netcdf_interface.o da_gpseph.o
da_wrfvar_top.o : da_wrfvar_top.f90 da_solve_init.inc da_solve_dual_res_init.inc da_solve.inc da_wrfvar_finalize.inc da_wrfvar_interface.inc da_wrfvar_run.inc da_wrfvar_init2.inc da_wrfvar_init1.inc da_wrf_interfaces.o da_rain.o da_synop.o da_ssmi.o da_sound.o da_ships.o da_satem.o da_radar.o da_mtgirs.o da_qscat.o da_profiler.o da_polaramv.o da_pilot.o da_metar.o da_gpsref.o da_gpspw.o da_geoamv.o da_buoy.o da_bogus.o da_airsr.o da_airep.o da_crtm.o da_tools.o da_vtox_transforms.o da_transfer_model.o da_tracing.o da_tools_serial.o da_test.o da_setup_structures.o da_reporting.o da_varbc.o da_radiance1.o da_physics.o da_par_util.o da_obs_io.o da_obs.o da_minimisation.o da_define_structures.o da_control.o module_comm_dm.o module_dm.o module_tiles.o module_state_description.o module_radiance.o da_wrfvar_io.o da_4dvar.o module_symbols_util.o module_driver_constants.o module_domain.o module_configure.o module_io_domain.o da_netcdf_interface.o da_gpseph.o da_varbc_tamdar.o
decode_airs.o : decode_airs.f90 module_read_airs.o
f_qv_from_rh.o : f_qv_from_rh.f90
gamma1.o : gamma1.f90 da_control.o
Expand Down
4 changes: 4 additions & 0 deletions var/da/da_control/da_control.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -278,6 +278,7 @@ module da_control
integer :: check_max_iv_unit, check_buddy_unit, rand_unit, omb_unit, &
filtered_obs_unit
integer :: biasprep_unit, qcstat_conv_unit
integer :: varbc_tamdar_unit

integer,parameter :: filename_len = 200

Expand Down Expand Up @@ -369,6 +370,7 @@ module da_control
integer :: cv_size_domain_jp ! Total jp cv size.
integer :: cv_size_domain_js ! Total js cv size.
integer :: cv_size_domain_jl ! Total jl cv size.
integer :: cv_size_domain_jt ! Total jt cv size.
integer :: cv_size_domain ! Total cv size.

! Hybrid:
Expand Down Expand Up @@ -419,6 +421,8 @@ module da_control
! missing_r
outside_of_domain = -77, & ! Data outside horizontal domain
! or time window, data set to missing_r
fail_varbc_aircraft = -55, & ! Data fail VarBC of aircraft
! => no action
wrong_direction = -15, & ! Wind vector direction <0 or> 360
! => direction set to missing_r
negative_spd = -14, & ! Wind vector norm is negative
Expand Down
23 changes: 23 additions & 0 deletions var/da/da_define_structures/da_define_structures.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -405,6 +405,25 @@ module da_define_structures
type (field_type) , pointer :: q (:) ! q.
end type tamdar_type

type varbc_tamdar_type
character(len=40) :: fmt_param ! Format of parameter table
integer :: nmaxpred ! Max. No. of predictors
integer :: nphase ! No. of flight phases
integer :: nair ! No. of aircrafts in table
integer :: npred ! No. of predictors
integer :: nmaxobs ! Max Obs No.
integer , pointer :: nobs (:,:)! Obs No. in proc
integer , pointer :: nobs_sum(:,:)! Total Obs No.
integer , pointer :: tail_id (:)! Tail ID of aircrafts
integer , pointer :: obs_sn(:,:,:)! Serial No. of Obs in proc
integer , pointer :: ifuse (:,:)! run varbc or not
integer , pointer :: index (:,:,:)! Index in CV
real , pointer :: pred (:,:,:)! Predictors
real , pointer :: param (:,:,:)! Parameters
real , pointer :: bgerr (:,:,:)! Bkg err in Hessian
real , pointer :: vtox(:,:,:,:)! Transformation of CV
end type varbc_tamdar_type

type airsr_type
real , pointer :: h (:) ! Height in m
real , pointer :: p (:) ! pressure.
Expand Down Expand Up @@ -665,6 +684,8 @@ module da_define_structures
type (synop_type) , pointer :: tamdar_sfc(:)
type (rain_type) , pointer :: rain(:)

type (varbc_tamdar_type) :: varbc_tamdar

real :: missing
real :: ptop
end type iv_type
Expand Down Expand Up @@ -953,6 +974,7 @@ module da_define_structures
real :: jl
real :: jd
real :: jm
real :: jt
type (jo_type) :: jo
end type j_type

Expand All @@ -963,6 +985,7 @@ module da_define_structures
integer :: size_jp ! Size of CV array for Jp term.
integer :: size_js ! Size of CV array for Js term.
integer :: size_jl ! Size of CV array for Jl term.
integer :: size_jt ! Size of CV array for Jt term.
integer :: size1c ! Complex size of CV array of 1st variable error.
integer :: size2c ! Complex size of CV array of 2nd variable error.
integer :: size3c ! Complex size of CV array of 3rd variable error.
Expand Down
Loading