Best fit model determined from the numerical simulations overlaid on the observations at different wavelengths. The corresponding frequency is shown on the right corner in each panel in unit of Hz. The x-axis is the time since trigger in units of seconds. The observed flux density of each instrument is indicated on the y-axis in units of mJy. All data were corrected for MW and host galaxy absorption and extinction effects before modelling. Red solid lines represent the modelled light curves.