The variational method for vortex flow (VF) analyses, called VF-Var (formulated in Part I), is applied to the 20 May 2013 Newcastle–Moore tornadic mesocyclone observed from the operational KTLX radar and an experimental phased array radar. The dual-Doppler analyzed vortex flow (VF) field reveals the following features: The axisymmetric part of the VF is a well-defined slantwise two-cell vortex in which the maximum tangential velocity is nearly 40 ms−1 at the edge of the vortex core (0.6 km from the vortex center), the central downdraft velocity reaches -35 ms−1 at 3 km height, and the surrounding-updraft velocity reaches 26 ms−1 at 5 km height. The total VF field is a loosely-defined slantwise two-cell vortex consisting of a nearly-axisymmetric vortex core (in which the ground-relative surface wind speed reaches 50 ms−1 on the southeast edge), a strong non-axisymmetric slantwise downdraft in the vortex core, and a main updraft in a banana-shape area southeast of the vortex core which extends slantwise upward and spirals cyclonically around the vortex core. The single-Doppler analysis with observations from the KTLX radar only exhibits roughly the same features as the dual-Doppler analysis but contains spurious vertical-motion structures in and around the vortex core. Analysis errors are assessed by leveraging the findings from Parts II and III, which indicate that the dual-Doppler analyzed VF is accurate enough to represent the true VF but the single-Doppler analyzed VF is not (especially for non-axisymmetric vertical motions in and around the vortex core), so the dual-Doppler analyzed VF should be useful for initializing/verifying high-resolution tornado simulations.