Transmission spectra data of periodically modulated graphene structure

All files are transmission through periodically modulated graphene structure. The transmission is calculated using python by the ratio of normal component of Poynting vector after and before the periodic structure, making use of scattering matrix theory and quantum mechanical graphene conductivity. All file names including "10" are for turbostratic graphene, where the conductivity model in calculation is multiplied by a factor of 10, otherwise they are for single graphene layer. File names containing "changeQ" have columns with different wavenumber with constant temperature. Files names containing "change" have columns with different temperatures but constant wavenumber. The first line in all data files declares values of parameters. The second line in all data files are the header names. The parameter “mubeta” is the unitless product \mu \beta occurring in the formula of graphene conductivity, whilst "b/d" is the filling ratio of grating structure, and "num_sheets" is the factor multiplied by the conductivity of graphene.  The columns headered "W" are the unitless frequencies normalised by the chemical potential in graphene conductivity. The columns headered by the form "Q=#G_1/2" are transmission values of the grating structure, at in-plane wavenumber equal to # times the wavenumber at edge of the first Brillouin zone. The datas in each file are occurring in blocks, separated by empty (whitespaced) lines. The parameter lists "crossing_ws", and "diff_orders" correspond to the periodicity of grating (so that the value is the reciprocal wavenumber unit), and diffraction order of transmission, respectively. From one block to the next, crossing_ws and diff_orders are cycled through in a nested manner, where diff_orders are cycled through first.