var/da/da_radiance/da_radiance.f90

   1 module da_radiance
   2
   3    !---------------------------------------------------------------------------
   4    ! Purpose: module for radiance data assimilation.
   5    !---------------------------------------------------------------------------
   6
   7 #if defined(HDF5)
   8    use hdf5
   9 #endif
  10
  11 #if defined(RTTOV) || defined(CRTM)
  12
  13    use module_domain, only : xb_type, domain
  14    use module_radiance, only : satinfo, &
  15       i_kind,r_kind, r_double, &
  16        one, zero, three,deg2rad,rad2deg, &
  17       q2ppmv, &
  18       init_constants_derived, &
  19       rttov_platform_name, rttov_inst_name, crtm_sensor_name  ! names used by both RTTOV and CRTM
  20 #ifdef RTTOV
  21    use module_radiance, only : coefs, rttov_coefs, rttov_profile, rttov_radiance, &
  22       rttov_transmission,errorstatus_success,gas_id_watervapour,rttov_emissivity
  23 #endif
  24 #ifdef CRTM
  25    use module_radiance, only : crtm_channelinfo_type, crtm_platform_name, crtm_init, &
  26    CRTM_Planck_Radiance, CRTM_Planck_Temperature
  27 #endif
  28
  29 #ifdef DM_PARALLEL
  30 !  use mpi, only : mpi_integer, mpi_status_size, mpi_min, mpi_max, mpi_minloc, &
  31 !      mpi_2double_precision
  32 #endif
  33
  34    use da_control, only : max_ob_levels,missing_r, &
  35       v_interp_p, v_interp_h, trace_use_dull, &
  36       missing, max_error_uv, max_error_t, rootproc, &
  37       max_error_p,max_error_q, radiance, &
  38       max_stheight_diff,missing_data,max_error_bq,max_error_slp, &
  39       max_error_bt, max_error_buv, rtminit_platform,rtminit_satid, &
  40       rtminit_nsensor,rtminit_sensor,filename_len,read_biascoef,analysis_date, &
  41       time_window_max,time_window_min,print_detail_obs,use_hsbobs,use_msuobs, &
  42       use_amsubobs,use_eos_amsuaobs,use_amsuaobs,use_hirs2obs,rtm_option, &
  43       rtm_option_rttov,rtm_option_crtm,use_airsobs,use_kma1dvar,use_hirs3obs, &
  44       use_ssmisobs,use_iasiobs,use_seviriobs,use_filtered_rad,print_detail_rad,stderr, mw_emis_sea, &
  45       rtminit_print, rttov_scatt,comm,root,ierr,biasprep, qc_rad, num_procs, &
  46       tovs_min_transfer,use_error_factor_rad,num_fgat_time,stdout,trace_use, &
  47       qc_good, qc_bad,myproc,biascorr,thinning,thinning_mesh, &
  48       rad_monitoring, monitor_on, kts, kte, kms, kme, calc_weightfunc, &
  49       use_mwtsobs, use_mwhsobs, use_mwhs2obs, use_atmsobs, use_amsr2obs, use_ahiobs,use_gmiobs, &
  50       use_hirs4obs, use_mhsobs,bufr_year, bufr_month,bufr_day,bufr_hour, &
  51       bufr_minute, bufr_second,bufr_solzen, bufr_station_height, &
  52       bufr_landsea_mask,bufr_solazi,tovs_end, max_tovs_input, bufr_satzen, nchan_mhs, &
  53       nchan_msu, nchan_amsua,nchan_hirs2, nchan_hirs3, nchan_hirs4, nchan_airs, &
  54       bufr_lon, bufr_satellite_id, bufr_ifov,nchan_amsub, tovs_start, bufr_lat, &
  55       use_pseudo_rad, pseudo_rad_platid,pseudo_rad_satid, pseudo_rad_senid, &
  56       pseudo_rad_ichan, pseudo_rad_inv, pseudo_rad_lat,pseudo_rad_lon, &
  57       pseudo_rad_err, use_simulated_rad,use_rttov_kmatrix, use_crtm_kmatrix , &
  58       use_rad,crtm_cloud, DT_cloud_model, global, use_varbc, freeze_varbc, &
  59       airs_warmest_fov, time_slots, interp_option, ids, ide, jds, jde, &
  60       ips, ipe, jps, jpe, simulated_rad_ngrid, obs_qc_pointer, use_blacklist_rad, use_satcv, &
  61       use_goesimgobs, pi, earth_radius, satellite_height,use_clddet_zz, ahi_superob_halfwidth, ahi_apply_clrsky_bias
  62
  63 #ifdef CRTM
  64    use da_crtm, only : da_crtm_init, da_get_innov_vector_crtm
  65 #endif
  66    use da_define_structures, only : maxmin_type, iv_type, y_type, jo_type, j_type, &
  67       bad_data_type, x_type, number_type, bad_data_type, &
  68       airsr_type,info_type, model_loc_type, varbc_info_type, varbc_type,clddet_geoir_type, superob_type
  69    use da_interpolation, only : da_to_zk, da_to_zk_new
  70    use da_tools_serial, only : da_get_unit, da_free_unit
  71    use da_par_util1, only : da_proc_sum_int,da_proc_sum_ints
  72 #ifdef DM_PARALLEL
  73    use da_par_util, only :  da_proc_stats_combine, true_mpi_real
  74 #else
  75    use da_par_util, only :  da_proc_stats_combine
  76 #endif
  77    use da_physics, only : da_sfc_pre, da_transform_xtopsfc, &
  78       da_transform_xtopsfc_adj,da_tpq_to_slp_lin,da_tpq_to_slp_adj
  79    use da_radiance1, only : num_tovs_before,num_tovs_after,tovs_copy_count, &
  80       tovs_send_pe, tovs_recv_pe, tovs_send_start, tovs_send_count, &
  81       tovs_recv_start,con_vars_type,aux_vars_type, datalink_type,da_qc_amsub, &
  82       da_qc_amsua,da_biascorr, da_detsurtyp,da_biasprep, &
  83       da_qc_rad, da_cld_eff_radius, da_read_biascoef
  84    use da_reporting, only : da_message, da_warning, message, da_error
  85 #ifdef RTTOV
  86    use da_rttov, only : da_rttov_init, da_get_innov_vector_rttov
  87 #endif
  88    use da_statistics, only : da_stats_calculate
  89    use da_tools, only : da_residual, da_obs_sfc_correction, &
  90       da_llxy, da_llxy_new, da_togrid_new, da_get_julian_time, da_get_time_slots, &
  91       da_xyll, map_info
  92    use da_tracing, only : da_trace_entry, da_trace_exit, da_trace, &
  93       da_trace_int_sort
  94    use da_varbc, only : da_varbc_direct,da_varbc_coldstart,da_varbc_precond, &
  95       da_varbc_pred
  96    use da_wrf_interfaces, only : wrf_dm_bcast_real
  97    use gsi_thinning, only : r999,r360,rlat_min,rlat_max,rlon_min,rlon_max, &
  98                             dlat_grid,dlon_grid,thinning_grid, &
  99                             makegrids,map2grids, &
 100                             destroygrids
 101    use mod_clddet_geoir, only: qc_SDob,qc_RTCT,qc_rtmt,qc_cwvt,qc_tit, &
 102                     find_std,calc_rtmt,calc_correlation
 103
 104    implicit none
 105
 106    include 'netcdf.inc'
 107
 108 #ifdef DM_PARALLEL
 109    include 'mpif.h'
 110 #endif
 111
 112 contains
 113
 114 #include "da_calculate_grady_rad.inc"
 115 #include "da_read_filtered_rad.inc"
 116 #include "da_read_simulated_rad.inc"
 117 #include "da_write_filtered_rad.inc"
 118 #include "da_read_obs_bufrtovs.inc"
 119 #include "da_read_obs_fy3.inc"
 120 #include "da_read_obs_hdf5mwhs2.inc"
 121 #include "da_read_obs_bufratms.inc"
 122 #include "ATMS_Spatial_Average.inc"
 123 #include "da_read_obs_bufrairs.inc"
 124 #include "da_read_obs_bufrssmis.inc"
 125 #include "da_read_obs_bufriasi.inc"
 126 #include "da_read_obs_bufrseviri.inc"
 127 #include "da_read_obs_hdf5amsr2.inc"
 128 #include "da_read_obs_hdf5ahi.inc"
 129 #include "da_read_obs_netcdf4ahi_geocat.inc"
 130 #include "da_read_obs_netcdf4ahi_jaxa.inc"
 131 #include "da_read_obs_ncgoesimg.inc"
 132 #include "da_read_obs_hdf5gmi.inc"
 133 #include "da_get_satzen.inc"
 134 #include "da_allocate_rad_iv.inc"
 135 #include "da_initialize_rad_iv.inc"
 136 #include "da_read_kma1dvar.inc"
 137 #include "da_sort_rad.inc"
 138 #include "da_setup_radiance_structures.inc"
 139 #include "da_radiance_init.inc"
 140 #include "da_get_innov_vector_radiance.inc"
 141 #include "da_read_pseudo_rad.inc"
 142 #include "da_blacklist_rad.inc"
 143 #include "da_deallocate_radiance.inc"
 144
 145 #endif
 146
 147 end module da_radiance
 148