Aim : The nearby TeV blazar 1ES 1959+650 (z=0.047) was reported to be in flaring state during June - July 2016 by Fermi-LAT, FACT, MAGIC and VERITAS collaborations. We studied the spectral energy distributions (SEDs) in different states of the flare during MJD 57530 - 57589 using simultaneous multiwaveband data to understand the possible broadband emission scenario during the flare. Methods : The UV/optical and X-ray data from UVOT and XRT respectively on board Swift and high energy $gamma$-ray data from Fermi-LAT are used to generate multiwaveband lightcurves as well as to obtain high flux states and quiescent state SEDs. The correlation and lag between different energy bands is quantified using discrete correlation function. The synchrotron self Compton (SSC) model was used to reproduce the observed SEDs during flaring and quiescent states of the source. Results : A decent correlation is seen between X-ray and high energy $gamma$-ray fluxes. The spectral hardening with increase in the flux is seen in X-ray band. The powerlaw index vs flux plot in $gamma$-ray band indicates the different emission regions for 0.1 - 3 GeV and 3-300 GeV energy photons. Two zone SSC model satisfactorily fits the observed broadband SEDs. The inner zone is mainly responsible for producing synchrotron peak and high energy $gamma$-ray part of the SED in all states. The second zone is mainly required to produce less variable optical/UV and low energy $gamma$-ray emission. Conclusions : Conventional single zone SSC model does not satisfactorily explain broadband emission during observation period considered. There is an indication of two emission zones in the jet which are responsible for producing broadband emission from optical to high energy $gamma$-rays.