Updated doc/make
Added function calcgmst()

* Added function calcgmst():
  * to astro.f90
  * computes Greenwich Mean Siderial time for a given Julian day