A three-dimensional subsonic aerodynamic panel code (VSAERO) was used to predict the effects of upper and lower external nozzle flap geometry on the external afterbody/nozzle pressure coefficient distributions and external nozzle drag of nonaxisymmetric convergent-divergent exhaust nozzles having parallel external sidewalls installed on a generic twin-engine high performance aircraft model. Nozzle static pressure coefficient distributions along the upper and lower surfaces near the model centerline and near the outer edges (corner) of the two surfaces were calculated, and nozzle drag was predicted using these surface pressure distributions. A comparison between the theoretical predictions and experimental wind tunnel data is made to evaluate the utility of the code in calculating the flow about these types of non-axisymmetric afterbody configurations. For free-stream Mach numbers of 0.60 and 0.90, the conditions where the flows were attached on the boattails yielded the best comparison between the theoretical predictions and the experimental data. For the Boattail terminal angles of greater than 15 deg., the experimental data for M = 0.60 and 0.90 indicated areas of separated flow, so the theoretical predictions failed to match the experimental data. Even though calculations of regions of separated flows are within the capabilities of the theoretical method, acceptable solutions were not obtained.