analyseMCMC_2dpdf_skymap.f90 File Reference

(AnalyseMCMC version 0.7.1, revision 448, hash 2f54722 (2013-11-19))

Routines to plot a 2D PDF on a sky map. More...

Functions/Subroutines
subroutine	plotthesky (bx10, bx20, by1, by2, raShift)
	Plot a sky map. More...
character function, dimension(10)	getsname (name)
	Get star name from bsc info. More...
real function	getmag (m, mlim)
	Determine size of stellar 'disc'. More...
subroutine	pgimag_project (z, nbx, nby, xb1, xb2, yb1, yb2, z1, z2, clr1, clr2, tr, projection)
	Plot a projected image. More...
subroutine	project_skymap (x, y, raCentre, projection)
	Project a sky map, using projection 'projection'. More...

Detailed Description

Routines to plot a 2D PDF on a sky map.

Function/Subroutine Documentation

subroutine plotthesky	(	real, intent(in)	bx10,
		real, intent(in)	bx20,
		real, intent(in)	by1,
		real, intent(in)	by2,
		real, intent(in)	raShift
	)

Plot a sky map.

Parameters

bx10	RA boundary 1
bx20	RA boundary 2
by1	Dec boundary 1
by2	Dec boundary 2
raShift	Shift in RA needed to centre the PDF

References drev2pi(), getmag(), getsname(), and project_skymap().

Referenced by pdfs2d().

character function, dimension(10) getsname

(

character, dimension(*), intent(in)

name

)

Get star name from bsc info.

Parameters

name	BSC star name

Referenced by plotthesky().

real function getmag	(	real, intent(in)	m,
		real, intent(in)	mlim
	)

Determine size of stellar 'disc'.

Parameters

m	Magnitude of the star
mlim	Limiting magnitude of the map

Referenced by plotthesky().

subroutine pgimag_project	(	real, dimension(nbx,nby), intent(in)	z,
		integer, intent(in)	nbx,
		integer, intent(in)	nby,
		integer, intent(in)	xb1,
		integer, intent(in)	xb2,
		integer, intent(in)	yb1,
		integer, intent(in)	yb2,
		real, intent(in)	z1,
		real, intent(in)	z2,
		integer, intent(in)	clr1,
		integer, intent(in)	clr2,
		real, dimension(6), intent(in)	tr,
		integer, intent(in)	projection
	)

Plot a projected image.

Parameters

z	Binned data array
nbx	Number of bins in the x-direction
nby	Number of bins in the y-direction
xb1	First bin to plot in the x-direction
xb2	Last bin to plot in the x-direction
yb1	First bin to plot in the y-direction
yb2	Last bin to plot in the y-direction
z1	Lowest array value, will appear in first colour (clr1)
z2	Highest array value, will appear in last colour (clr2)
clr1	First colour, use with lowest array valye (z1)
clr2	Last colour, use with lowest array valye (z1)
tr
projection

Note

Clone of pgimag, use projection if projection > 0

References project_skymap().

Referenced by plot_2d_pdf().

subroutine project_skymap	(	real, intent(inout)	x,
		real, intent(inout)	y,
		real, intent(in)	raCentre,
		integer, intent(in)	projection
	)

Project a sky map, using projection 'projection'.

Parameters

x	X coordinate of object (I/O)
y	Y coordinate of object (I/O)
raCentre	Central RA of map
projection	Choice of projection: 1 - Mollweide

Projections:

• 1: Mollweide projection:
  • http://en.wikipedia.org/wiki/Mollweide_projection
  • Newton-Rapson scheme to solve equation: 2*theta + sin(2*theta) = pi*sin(y*rd2r)
  • Convergence is relatively fast, somewhat slower near poles

Referenced by pgimag_project(), plot_values_in_2d_pdf(), and plotthesky().