I am new to MPAS-JEDI and am trying to assimilate Nexrad Level II data to a single time in an existing MPAS regional simulation. I have successfully generated the appropriately formatted ioda file (I think), as well as the appropriate yaml file for JEDI. However, when I go to run the mpasjedi_variational.x executable against this yaml I get an error saying
"At line 867 of file /mercer/williams/mpas-bundle/MPAS/src/core_atmosphere/physics/physics_wrf/module_mp_thompson.F (unit = 101, file = ‘MP_THOMPSON_QRacrQG_DATA.DBL’)
Fortran runtime error: End of file
"
I have verified that this table (as well the others generated by mpas_atmosphere_buildtables) are present in my working directory and my previous MPAS run that did not do any data assimilation worked fine with these tables. It is unclear why this is failing with JEDI when it worked with a vanilla MPAS run without any DA. I’ve included my .yaml and namelist/stream files below. Everything has been run in a single directory. I am using a build of mpas-bundle, if that helps. What could be causing this?
namelist.atmosphere:
&nhyd_model
config_time_integration_order = 2
config_dt = 3
config_start_time = ‘2025-01-05_18:00:00’
config_run_duration = ‘0_00:15:00’
config_len_disp=9000.0
config_split_dynamics_transport = true
config_number_of_sub_steps = 2
config_dynamics_split_steps = 3
config_horiz_mixing = ‘2d_smagorinsky’
config_visc4_2dsmag = 0.05
config_scalar_advection = true
config_monotonic = true
config_coef_3rd_order = 0.25
config_epssm = 0.1
config_smdiv = 0.1
/
&limited_area
config_apply_lbcs=true
/
&damping
config_zd = 22000.0
config_xnutr = 0.2
/
&io
config_pio_num_iotasks = 0
config_pio_stride = 1
config_restart_timestamp_name = ‘restart’
/
&decomposition
config_block_decomp_file_prefix = ‘MS_Regional.graph.info.part.’
/
&restart
config_do_restart = false
/
&printout
config_print_global_minmax_vel = true
config_print_detailed_minmax_vel = false
/
&IAU
config_IAU_option = ‘off’
config_IAU_window_length_s = 21600.
/
&physics
config_sst_update = false
config_sstdiurn_update = false
config_deepsoiltemp_update = false
config_radtlw_interval = ‘00:10:00’
config_radtsw_interval = ‘00:10:00’
config_bucket_update = ‘none’
config_gwdo_scheme = ‘off’
config_physics_suite=‘convection_permitting’
/
&soundings
config_sounding_interval = ‘none’
/
streams.atmosphere:
<immutable_stream name=“input”
type=“input”
filename_template=“mpas_ics.nc”
input_interval=“initial_only” />
<immutable_stream name=“restart”
type=“input;output”
filename_template=“restart.$Y-$M-$D_$h.$m.$s.nc”
input_interval = “initial_only”
output_interval=“0:15:00” />
<immutable_stream name=“iau”
type=“input”
filename_template=“x1.40962.AmB.$Y-$M-$D_$h.$m.$s.nc”
filename_interval=“none”
packages=“iau”
input_interval=“initial_only” />
<immutable_stream name=“lbc_in”
type=“input”
filename_template=“lbc.$Y-$M-$D_$h.$m.$s.nc”
filename_interval=“input_interval”
packages=“limited_area”
input_interval=“3:00:00” />
radar.yaml
cost function:
cost type: 3D-Var
analysis variables: [u, v, w, theta, qv, qc, qr]
time window:
begin: 2025-01-05T18:05:00Z
length: PT20M
geometry:
model: MPAS
nml_file: /mercer/williams/MPAS-Model/namelist.atmosphere
streams_file: /mercer/williams/MPAS-Model/streams.atmosphere
background:
filename: restart.2025-01-05_18:15:00.nc
date: 2025-01-05T18:15:00Z
background error:
covariance model: SABER
saber central block:
saber blocks:
- covariance model: BUMP
bump:
io:
data directory: ./bump
files prefix: mpas_b
drivers:
multivariate strategy: univariate
compute covariance: true
modeling:
horizontal length-scale: 500000.0
vertical length-scale: 2000.0
observations: - obs space:
name: KGWX radial velocity
obsdatain:
engine:
type: H5File
obsfile: /mercer/williams/jedi/kgwx_velocity_ioda.nc
obsdataout:
engine:
type: H5File
obsfile: /mercer/williams/jedi/kgwx_velocity_ioda_out.nc
simulated variables: [radialVelocity]
obs operator:
name: RadarRadialVelocity
variables: [wind]
options:
vertical coordinate: height
use_reflectivity: false
obs error:
covariance model: diagonal
diagonal:
value: 25.0
obs filters:
- filter: Domain Check
where: - variable: height
minvalue: 0.0
maxvalue: 12000.0 - filter: Bounds Check
where: - variable: radialVelocity
minvalue: -80.0
maxvalue: 80.0