docs: def gtd7d(Input, flags, output): # pragma: no cover
      '''A separate subroutine (GTD7D) computes the ‘‘effec-
      tive’’ mass density by summing the thermospheric mass
      density and the mass density of the anomalous oxygen
      component. Below 500 km, the effective mass density is
      equivalent to the thermospheric mass density, and we drop
      the distinction.
      '''
      gtd7(Input, flags, output)
      output.d[5] = 1.66E-24 * (4.0 * output.d[0] + 16.0 * output.d[1] + 28.0 * output.d[2] + 32.0 * output.d[3] + 40.0 * output.d[4] + output.d[6] + 14.0 * output.d[7] + 16.0 * output.d[8]);
      if (flags.sw[0]):
          output.d[5]=output.d[5]/1000;
      return