head	8.7;
access;
symbols;
locks; strict;
comment	@c @;


8.7
date	2026.05.11.12.01.11;	author mcalabre;	state Exp;
branches;
next	8.6;

8.6
date	2026.03.29.13.53.57;	author mcalabre;	state Exp;
branches
	8.6.1.1;
next	8.5;

8.5
date	2025.12.06.13.47.42;	author mcalabre;	state Exp;
branches
	8.5.1.1;
next	8.4;

8.4
date	2024.10.28.13.56.18;	author mcalabre;	state Exp;
branches
	8.4.1.1;
next	8.3;

8.3
date	2024.05.13.16.33.01;	author mcalabre;	state Exp;
branches
	8.3.1.1;
next	8.2;

8.2
date	2023.11.16.10.00.31;	author mcalabre;	state Exp;
branches
	8.2.1.1;
next	8.1;

8.1
date	2023.07.05.17.12.07;	author mcalabre;	state Exp;
branches
	8.1.1.1;
next	8.0;

8.0
date	2023.07.01.06.15.05;	author mcalabre;	state Exp;
branches
	8.0.1.1;
next	7.13;

7.13
date	2022.10.07.03.43.35;	author mcalabre;	state Exp;
branches
	7.13.1.1;
next	7.12;

7.12
date	2022.09.09.04.57.58;	author mcalabre;	state Exp;
branches
	7.12.1.1;
next	7.11;

7.11
date	2022.04.26.06.13.52;	author mcalabre;	state Exp;
branches
	7.11.1.1;
next	7.10;

7.10
date	2022.04.23.14.25.15;	author mcalabre;	state Exp;
branches
	7.10.1.1;
next	7.9;

7.9
date	2022.03.25.15.14.49;	author mcalabre;	state Exp;
branches
	7.9.1.1;
next	7.8;

7.8
date	2022.03.25.06.11.50;	author mcalabre;	state Exp;
branches
	7.8.1.1;
next	7.7;

7.7
date	2021.07.12.06.36.49;	author mcalabre;	state Exp;
branches
	7.7.1.1;
next	7.6;

7.6
date	2021.04.13.12.57.02;	author mcalabre;	state Exp;
branches
	7.6.1.1;
next	7.5;

7.5
date	2021.03.20.05.54.59;	author mcalabre;	state Exp;
branches
	7.5.1.1;
next	7.4;

7.4
date	2021.01.31.02.24.52;	author mcalabre;	state Exp;
branches
	7.4.1.1;
next	7.3;

7.3
date	2020.06.03.03.37.03;	author mcalabre;	state Exp;
branches
	7.3.1.1;
next	7.2;

7.2
date	2020.03.09.07.31.24;	author mcalabre;	state Exp;
branches
	7.2.1.1;
next	7.1;

7.1
date	2019.12.31.13.25.20;	author mcalabre;	state Exp;
branches
	7.1.1.1;
next	6.4;

6.4
date	2019.08.15.09.30.18;	author mcalabre;	state Exp;
branches
	6.4.1.1;
next	6.3;

6.3
date	2019.07.12.07.33.40;	author mcalabre;	state Exp;
branches
	6.3.1.1;
next	6.2;

6.2
date	2018.10.20.10.03.14;	author mcalabre;	state Exp;
branches
	6.2.1.1;
next	6.1;

6.1
date	2018.10.19.08.35.17;	author mcalabre;	state Exp;
branches
	6.1.1.1;
next	5.20;

5.20
date	2018.10.05.05.32.42;	author mcalabre;	state Exp;
branches
	5.20.1.1;
next	5.19;

5.19
date	2018.07.26.15.37.05;	author mcalabre;	state Exp;
branches
	5.19.1.1;
next	5.18;

5.18
date	2018.01.10.08.32.14;	author mcalabre;	state Exp;
branches
	5.18.1.1;
next	5.17;

5.17
date	2017.09.18.08.44.25;	author mcalabre;	state Exp;
branches
	5.17.1.1;
next	5.16;

5.16
date	2017.01.15.04.25.03;	author mcalabre;	state Exp;
branches
	5.16.1.1;
next	5.15;

5.15
date	2016.04.05.12.55.13;	author mcalabre;	state Exp;
branches
	5.15.1.1;
next	5.14;

5.14
date	2016.02.07.10.49.33;	author mcalabre;	state Exp;
branches
	5.14.1.1;
next	5.13;

5.13
date	2016.01.26.08.51.16;	author mcalabre;	state Exp;
branches
	5.13.1.1;
next	5.12;

5.12
date	2015.11.15.02.13.03;	author mcalabre;	state Exp;
branches
	5.12.1.1;
next	5.11;

5.11
date	2015.10.18.09.13.08;	author mcalabre;	state Exp;
branches
	5.11.1.1;
next	5.10;

5.10
date	2015.10.09.08.19.18;	author mcalabre;	state Exp;
branches
	5.10.1.1;
next	5.9;

5.9
date	2015.07.21.09.20.04;	author mcalabre;	state Exp;
branches
	5.9.1.1;
next	5.8;

5.8
date	2015.07.08.11.04.03;	author mcalabre;	state Exp;
branches
	5.8.1.1;
next	5.7;

5.7
date	2015.06.29.02.44.19;	author mcalabre;	state Exp;
branches
	5.7.1.1;
next	5.6;

5.6
date	2015.06.14.07.11.26;	author mcalabre;	state Exp;
branches
	5.6.1.1;
next	5.5;

5.5
date	2015.05.05.13.16.34;	author mcalabre;	state Exp;
branches
	5.5.1.1;
next	5.4;

5.4
date	2015.04.21.14.41.08;	author mcalabre;	state Exp;
branches
	5.4.1.1;
next	5.3;

5.3
date	2015.04.21.02.50.54;	author mcalabre;	state Exp;
branches
	5.3.1.1;
next	5.2;

5.2
date	2015.04.15.12.35.10;	author mcalabre;	state Exp;
branches
	5.2.1.1;
next	5.1;

5.1
date	2015.04.07.10.56.39;	author mcalabre;	state Exp;
branches
	5.1.1.1;
next	5.0;

5.0
date	2015.04.05.12.25.03;	author mcalabre;	state Exp;
branches
	5.0.1.1;
next	4.25;

4.25
date	2014.12.14.14.29.38;	author mcalabre;	state Exp;
branches
	4.25.1.1;
next	4.24;

4.24
date	2014.09.18.15.25.02;	author mcalabre;	state Exp;
branches
	4.24.1.1;
next	4.23;

4.23
date	2014.05.13.05.50.53;	author mcalabre;	state Exp;
branches
	4.23.1.1;
next	4.22;

4.22
date	2014.04.12.15.03.54;	author mcalabre;	state Exp;
branches
	4.22.1.1;
next	4.21;

4.21
date	2014.03.24.05.12.29;	author mcalabre;	state Exp;
branches
	4.21.1.1;
next	4.20;

4.20
date	2013.12.18.05.42.51;	author mcalabre;	state Exp;
branches
	4.20.1.1;
next	4.19;

4.19
date	2013.09.29.14.17.53;	author mcalabre;	state Exp;
branches
	4.19.1.1;
next	4.18;

4.18
date	2013.07.13.10.00.06;	author mcalabre;	state Exp;
branches
	4.18.1.1;
next	4.17;

4.17
date	2013.01.29.05.29.22;	author cal103;	state Exp;
branches
	4.17.1.1;
next	4.16;

4.16
date	2012.11.07.04.42.45;	author cal103;	state Exp;
branches
	4.16.1.1;
next	4.15;

4.15
date	2012.09.26.14.26.06;	author cal103;	state Exp;
branches
	4.15.1.1;
next	4.14;

4.14
date	2012.07.13.10.02.47;	author cal103;	state Exp;
branches
	4.14.1.1;
next	4.13;

4.13
date	2012.03.14.07.39.31;	author cal103;	state Exp;
branches
	4.13.1.1;
next	4.12;

4.12
date	2012.02.29.07.57.57;	author cal103;	state Exp;
branches
	4.12.1.1;
next	4.11;

4.11
date	2012.02.20.23.30.03;	author cal103;	state Exp;
branches
	4.11.1.1;
next	4.10;

4.10
date	2012.02.05.23.41.45;	author cal103;	state Exp;
branches
	4.10.1.1;
next	4.9;

4.9
date	2012.01.24.05.40.30;	author cal103;	state Exp;
branches
	4.9.1.1;
next	4.8;

4.8
date	2011.08.15.08.05.55;	author cal103;	state Exp;
branches
	4.8.1.1;
next	4.7;

4.7
date	2011.02.07.07.03.43;	author cal103;	state Exp;
branches
	4.7.1.1;
next	4.6;

4.6
date	2010.11.16.06.15.32;	author cal103;	state Exp;
branches
	4.6.1.1;
next	4.5;

4.5
date	2010.07.16.07.01.27;	author cal103;	state Exp;
branches
	4.5.1.1;
next	4.4;

4.4
date	2009.08.06.07.54.40;	author cal103;	state Exp;
branches
	4.4.1.1;
next	4.3;

4.3
date	2007.12.27.05.49.14;	author cal103;	state Exp;
branches
	4.3.1.1;
next	4.2;

4.2
date	2005.09.21.13.27.15;	author cal103;	state Exp;
branches
	4.2.1.1;
next	4.1;

4.1
date	2005.08.26.12.54.06;	author cal103;	state Exp;
branches;
next	4.0;

4.0
date	2005.02.04.12.22.47;	author cal103;	state Exp;
branches;
next	3.6;

3.6
date	2004.08.24.08.42.14;	author cal103;	state Exp;
branches;
next	3.5;

3.5
date	2004.06.28.05.37.31;	author mcalabre;	state Exp;
branches;
next	3.4;

3.4
date	2004.02.10.23.56.37;	author mcalabre;	state Exp;
branches;
next	3.3;

3.3
date	2003.10.22.08.55.52;	author mcalabre;	state Exp;
branches;
next	3.2;

3.2
date	2003.09.09.08.04.14;	author mcalabre;	state Exp;
branches;
next	3.1;

3.1
date	2003.04.28.08.53.24;	author mcalabre;	state Exp;
branches;
next	1.22;

1.22
date	2002.11.06.00.40.53;	author mcalabre;	state Exp;
branches;
next	1.21;

1.21
date	2002.04.03.06.21.31;	author mcalabre;	state Exp;
branches;
next	1.20;

1.20
date	2001.11.20.11.39.28;	author mcalabre;	state Exp;
branches;
next	1.19;

1.19
date	2001.02.19.02.18.18;	author mcalabre;	state Exp;
branches;
next	1.18;

1.18
date	2000.12.07.00.04.49;	author mcalabre;	state Exp;
branches;
next	1.17;

1.17
date	2000.05.24.05.24.06;	author mcalabre;	state Exp;
branches;
next	1.16;

1.16
date	2000.05.10.04.54.48;	author mcalabre;	state Exp;
branches;
next	1.15;

1.15
date	99.12.01.04.10.42;	author mcalabre;	state Exp;
branches;
next	1.14;

1.14
date	99.07.22.02.04.16;	author mcalabre;	state Exp;
branches;
next	1.13;

1.13
date	97.04.11.02.15.17;	author mcalabre;	state Exp;
branches;
next	1.12;

1.12
date	97.04.09.01.57.19;	author mcalabre;	state Exp;
branches;
next	1.11;

1.11
date	97.04.03.04.14.08;	author mcalabre;	state Exp;
branches;
next	1.10;

1.10
date	97.03.19.00.29.05;	author mcalabre;	state Exp;
branches;
next	1.9;

1.9
date	97.03.18.07.30.16;	author mcalabre;	state Exp;
branches;
next	1.8;

1.8
date	97.03.14.06.58.29;	author mcalabre;	state Exp;
branches;
next	1.7;

1.7
date	97.03.13.06.57.17;	author mcalabre;	state Exp;
branches;
next	1.6;

1.6
date	97.03.11.06.53.24;	author mcalabre;	state Exp;
branches;
next	1.5;

1.5
date	97.03.10.22.31.50;	author mcalabre;	state Exp;
branches;
next	1.4;

1.4
date	97.02.19.06.58.56;	author mcalabre;	state Exp;
branches;
next	1.3;

1.3
date	97.02.18.03.39.53;	author mcalabre;	state Exp;
branches;
next	1.2;

1.2
date	97.02.07.05.48.22;	author mcalabre;	state Exp;
branches;
next	1.1;

1.1
date	97.02.07.03.57.14;	author mcalabre;	state Exp;
branches;
next	;

4.2.1.1
date	2005.12.05.03.23.34;	author cal103;	state Exp;
branches;
next	4.2.1.2;

4.2.1.2
date	2005.12.05.03.32.04;	author cal103;	state Exp;
branches;
next	4.2.1.3;

4.2.1.3
date	2005.12.16.03.42.59;	author cal103;	state Exp;
branches;
next	4.2.1.4;

4.2.1.4
date	2007.12.03.05.41.29;	author cal103;	state Exp;
branches;
next	;

4.3.1.1
date	2007.12.27.13.26.31;	author cal103;	state Exp;
branches;
next	4.3.1.2;

4.3.1.2
date	2009.07.23.04.17.17;	author cal103;	state Exp;
branches;
next	;

4.4.1.1
date	2009.08.10.08.55.45;	author cal103;	state Exp;
branches;
next	4.4.1.2;

4.4.1.2
date	2010.01.19.03.45.24;	author cal103;	state Exp;
branches;
next	;

4.5.1.1
date	2010.07.16.07.05.33;	author cal103;	state Exp;
branches;
next	;

4.6.1.1
date	2010.11.16.06.16.20;	author cal103;	state Exp;
branches;
next	4.6.1.2;

4.6.1.2
date	2011.01.14.06.35.20;	author cal103;	state Exp;
branches;
next	;

4.7.1.1
date	2011.02.07.07.04.23;	author cal103;	state Exp;
branches;
next	;

4.8.1.1
date	2011.08.15.08.07.07;	author cal103;	state Exp;
branches;
next	4.8.1.2;

4.8.1.2
date	2012.01.24.05.24.50;	author cal103;	state Exp;
branches;
next	;

4.9.1.1
date	2012.01.24.05.41.46;	author cal103;	state Exp;
branches;
next	;

4.10.1.1
date	2012.02.05.23.42.47;	author cal103;	state Exp;
branches;
next	;

4.11.1.1
date	2012.02.20.23.31.14;	author cal103;	state Exp;
branches;
next	;

4.12.1.1
date	2012.02.29.07.59.11;	author cal103;	state Exp;
branches;
next	;

4.13.1.1
date	2012.03.14.07.40.39;	author cal103;	state Exp;
branches;
next	4.13.1.2;

4.13.1.2
date	2012.07.13.03.43.52;	author cal103;	state Exp;
branches;
next	;

4.14.1.1
date	2012.07.13.10.05.05;	author cal103;	state Exp;
branches;
next	;

4.15.1.1
date	2012.09.26.14.27.53;	author cal103;	state Exp;
branches;
next	;

4.16.1.1
date	2012.11.07.04.45.58;	author cal103;	state Exp;
branches;
next	4.16.1.2;

4.16.1.2
date	2013.01.29.05.11.32;	author cal103;	state Exp;
branches;
next	;

4.17.1.1
date	2013.01.29.05.32.33;	author cal103;	state Exp;
branches;
next	;

4.18.1.1
date	2013.07.13.10.02.40;	author mcalabre;	state Exp;
branches;
next	;

4.19.1.1
date	2013.09.29.14.20.41;	author mcalabre;	state Exp;
branches;
next	;

4.20.1.1
date	2013.12.18.05.46.57;	author mcalabre;	state Exp;
branches;
next	4.20.1.2;

4.20.1.2
date	2014.03.21.06.45.26;	author mcalabre;	state Exp;
branches;
next	;

4.21.1.1
date	2014.03.24.05.15.31;	author mcalabre;	state Exp;
branches;
next	;

4.22.1.1
date	2014.04.12.15.06.34;	author mcalabre;	state Exp;
branches;
next	;

4.23.1.1
date	2014.05.13.05.53.52;	author mcalabre;	state Exp;
branches;
next	;

4.24.1.1
date	2014.09.18.15.28.12;	author mcalabre;	state Exp;
branches;
next	;

4.25.1.1
date	2014.12.14.14.32.48;	author mcalabre;	state Exp;
branches;
next	4.25.1.2;

4.25.1.2
date	2015.01.06.01.03.14;	author mcalabre;	state Exp;
branches;
next	;

5.0.1.1
date	2015.04.05.12.28.04;	author mcalabre;	state Exp;
branches;
next	;

5.1.1.1
date	2015.04.07.11.01.13;	author mcalabre;	state Exp;
branches;
next	;

5.2.1.1
date	2015.04.15.12.39.18;	author mcalabre;	state Exp;
branches;
next	;

5.3.1.1
date	2015.04.21.02.54.06;	author mcalabre;	state Exp;
branches;
next	;

5.4.1.1
date	2015.04.21.14.44.32;	author mcalabre;	state Exp;
branches;
next	;

5.5.1.1
date	2015.05.05.13.21.04;	author mcalabre;	state Exp;
branches;
next	;

5.6.1.1
date	2015.06.14.07.14.48;	author mcalabre;	state Exp;
branches;
next	;

5.7.1.1
date	2015.06.29.02.47.27;	author mcalabre;	state Exp;
branches;
next	;

5.8.1.1
date	2015.07.08.11.07.22;	author mcalabre;	state Exp;
branches;
next	;

5.9.1.1
date	2015.07.21.09.28.12;	author mcalabre;	state Exp;
branches;
next	;

5.10.1.1
date	2015.10.09.08.22.40;	author mcalabre;	state Exp;
branches;
next	;

5.11.1.1
date	2015.10.18.09.16.31;	author mcalabre;	state Exp;
branches;
next	;

5.12.1.1
date	2015.11.15.02.16.18;	author mcalabre;	state Exp;
branches;
next	5.12.1.2;

5.12.1.2
date	2016.01.25.05.22.41;	author mcalabre;	state Exp;
branches;
next	;

5.13.1.1
date	2016.01.26.08.54.36;	author mcalabre;	state Exp;
branches;
next	;

5.14.1.1
date	2016.02.07.10.52.48;	author mcalabre;	state Exp;
branches;
next	;

5.15.1.1
date	2016.04.05.12.58.39;	author mcalabre;	state Exp;
branches;
next	5.15.1.2;

5.15.1.2
date	2017.01.15.04.08.24;	author mcalabre;	state Exp;
branches;
next	;

5.16.1.1
date	2017.01.15.04.29.04;	author mcalabre;	state Exp;
branches;
next	;

5.17.1.1
date	2017.09.18.08.46.01;	author mcalabre;	state Exp;
branches;
next	5.17.1.2;

5.17.1.2
date	2018.01.10.08.27.33;	author mcalabre;	state Exp;
branches;
next	;

5.18.1.1
date	2018.01.10.08.33.04;	author mcalabre;	state Exp;
branches;
next	;

5.19.1.1
date	2018.07.26.15.41.42;	author mcalabre;	state Exp;
branches;
next	;

5.20.1.1
date	2018.10.05.05.33.33;	author mcalabre;	state Exp;
branches;
next	;

6.1.1.1
date	2018.10.19.08.36.15;	author mcalabre;	state Exp;
branches;
next	;

6.2.1.1
date	2018.10.20.10.04.15;	author mcalabre;	state Exp;
branches;
next	6.2.1.2;

6.2.1.2
date	2019.05.10.02.13.12;	author mcalabre;	state Exp;
branches;
next	;

6.3.1.1
date	2019.07.12.07.34.37;	author mcalabre;	state Exp;
branches;
next	;

6.4.1.1
date	2019.08.15.09.31.29;	author mcalabre;	state Exp;
branches;
next	6.4.1.2;

6.4.1.2
date	2019.12.31.03.29.47;	author mcalabre;	state Exp;
branches;
next	;

7.1.1.1
date	2019.12.31.13.28.19;	author mcalabre;	state Exp;
branches;
next	;

7.2.1.1
date	2020.03.09.07.34.50;	author mcalabre;	state Exp;
branches;
next	7.2.1.2;

7.2.1.2
date	2020.05.07.04.56.37;	author mcalabre;	state Exp;
branches;
next	;

7.3.1.1
date	2020.06.03.03.38.11;	author mcalabre;	state Exp;
branches;
next	7.3.1.2;

7.3.1.2
date	2021.01.19.14.12.19;	author mcalabre;	state Exp;
branches;
next	7.3.1.3;

7.3.1.3
date	2021.01.28.12.32.51;	author mcalabre;	state Exp;
branches;
next	;

7.4.1.1
date	2021.01.31.02.26.24;	author mcalabre;	state Exp;
branches;
next	;

7.5.1.1
date	2021.03.20.05.56.42;	author mcalabre;	state Exp;
branches;
next	;

7.6.1.1
date	2021.04.13.13.00.37;	author mcalabre;	state Exp;
branches;
next	7.6.1.2;

7.6.1.2
date	2021.07.02.02.31.54;	author mcalabre;	state Exp;
branches;
next	7.6.1.3;

7.6.1.3
date	2021.07.02.02.50.23;	author mcalabre;	state Exp;
branches;
next	7.6.1.4;

7.6.1.4
date	2021.07.02.02.56.19;	author mcalabre;	state Exp;
branches;
next	;

7.7.1.1
date	2021.07.12.06.38.28;	author mcalabre;	state Exp;
branches;
next	7.7.1.2;

7.7.1.2
date	2022.03.25.06.06.01;	author mcalabre;	state Exp;
branches;
next	;

7.8.1.1
date	2022.03.25.06.14.15;	author mcalabre;	state Exp;
branches;
next	;

7.9.1.1
date	2022.03.25.15.16.40;	author mcalabre;	state Exp;
branches;
next	;

7.10.1.1
date	2022.04.23.14.26.55;	author mcalabre;	state Exp;
branches;
next	;

7.11.1.1
date	2022.04.26.06.15.33;	author mcalabre;	state Exp;
branches;
next	;

7.12.1.1
date	2022.09.09.05.00.36;	author mcalabre;	state Exp;
branches;
next	;

7.13.1.1
date	2022.10.07.03.45.47;	author mcalabre;	state Exp;
branches;
next	7.13.1.2;

7.13.1.2
date	2023.06.15.04.17.22;	author mcalabre;	state Exp;
branches;
next	;

8.0.1.1
date	2023.07.01.06.22.45;	author mcalabre;	state Exp;
branches;
next	8.0.1.2;

8.0.1.2
date	2023.07.04.12.59.10;	author mcalabre;	state Exp;
branches;
next	8.0.1.3;

8.0.1.3
date	2023.07.04.15.15.23;	author mcalabre;	state Exp;
branches;
next	;

8.1.1.1
date	2023.07.05.17.15.01;	author mcalabre;	state Exp;
branches;
next	;

8.2.1.1
date	2023.11.16.10.05.59;	author mcalabre;	state Exp;
branches;
next	8.2.1.2;

8.2.1.2
date	2024.04.02.13.43.15;	author mcalabre;	state Exp;
branches;
next	;

8.3.1.1
date	2024.05.13.16.35.40;	author mcalabre;	state Exp;
branches;
next	;

8.4.1.1
date	2024.10.28.13.58.42;	author mcalabre;	state Exp;
branches;
next	8.4.1.2;

8.4.1.2
date	2025.12.06.13.37.20;	author mcalabre;	state Exp;
branches;
next	;

8.5.1.1
date	2025.12.06.14.05.17;	author mcalabre;	state Exp;
branches;
next	8.5.1.2;

8.5.1.2
date	2026.03.29.13.36.20;	author mcalabre;	state Exp;
branches;
next	;

8.6.1.1
date	2026.03.29.14.19.28;	author mcalabre;	state Exp;
branches;
next	;


desc
@FORTRAN test program for PGSBOX.
@


8.7
log
@WCSLIB 8.7
@
text
@*=======================================================================
*
* PGSBOX 8.7 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2026, Mark Calabretta
*
* This file is part of PGSBOX.
*
* PGSBOX is free software: you can redistribute it and/or modify it
* under the terms of the GNU Lesser General Public License as published
* by the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* PGSBOX is distributed in the hope that it will be useful, but WITHOUT
* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
* FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
* License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* License along with PGSBOX.  If not, see http://www.gnu.org/licenses.
*
* Author: Mark Calabretta, Australia Telescope National Facility, CSIRO.
* http://www.atnf.csiro.au/computing/software/wcs
* $Id: pgtest.f,v 8.6.1.1 2026/03/29 14:19:28 mcalabre Exp mcalabre $
*=======================================================================
      PROGRAM PGTEST
*=======================================================================
      LOGICAL   LARGE
      INTEGER   C0(7), CI(7), GCODE(2), IC, IERR, J, NAXIS(2),
     :          NLIPRM(2), STATUS
      REAL      BLC(2), SCL, TRC(2)
      DOUBLE PRECISION CACHE(4,0:256), DEC0, GRID1(0:8), GRID2(0:8),
     :          NLDPRM(8), ROTN, TIKLEN
      CHARACTER CTYPE(2)*72, DEVTYP*16, ESCAPE*1, FCODE(2)*4,
     :          IDENTS(3)*80, NLCPRM(1)*1, OPT(2)*1

      DOUBLE PRECISION D2R, PI
      PARAMETER (PI = 3.141592653589793238462643D0)
      PARAMETER (D2R = PI/180D0)

      INCLUDE 'wcs.inc'
      INTEGER   WCS(WCSLEN)

      EXTERNAL LNGVEL
      EXTERNAL FSCAN
      EXTERNAL PGCRFN
      EXTERNAL PGWCSL
*-----------------------------------------------------------------------
*     Setup.
      NAXIS(1) = 512
      NAXIS(2) = 512

      BLC(1) = 0.5
      BLC(2) = 0.5
      TRC(1) = NAXIS(1) + 0.5
      TRC(2) = NAXIS(2) + 0.5

      ESCAPE = CHAR(92)

*     CALL PGBEG (0, '?', 1, 1)
      call pgbeg (0, '/xw', 1, 1)

      CALL PGQINF ('TYPE', DEVTYP, J)
      IF (DEVTYP.EQ.'PS'  .OR.
     :    DEVTYP.EQ.'VPS' .OR.
     :    DEVTYP.EQ.'CPS' .OR.
     :    DEVTYP.EQ.'VCPS') THEN
*       Switch black and white.
        CALL PGSCR (0, 1.0, 1.0, 1.0)
        CALL PGSCR (1, 0.0, 0.0, 0.0)
      END IF

      LARGE = DEVTYP.EQ.'XWINDOW'
      IF (LARGE) THEN
        SCL = 1.0
        CALL PGVSTD ()
      ELSE
        SCL = 0.7
        CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
      END IF

*     Yellow.
      CALL PGSCR (2, 1.0, 1.0, 0.0)
*     White.
      CALL PGSCR (3, 1.0, 1.0, 1.0)
*     Pale blue.
      CALL PGSCR (4, 0.5, 0.5, 0.8)
*     Pale red.
      CALL PGSCR (5, 0.8, 0.5, 0.5)
*     Grey.
      CALL PGSCR (6, 0.7, 0.7, 0.7)
*     Dark green.
      CALL PGSCR (7, 0.3, 0.5, 0.3)

      C0(1) = -1
      C0(2) = -1
      C0(3) = -1
      C0(4) = -1
      C0(5) = -1
      C0(6) = -1
      C0(7) = -1

      CALL PGWNAD (0.0, 1.0, 0.0, 1.0)
      CALL PGASK (.TRUE.)
      CALL PGPAGE ()

      STATUS = WCSPTI (WCS, WCS_FLAG, -1, 0, 0)

*-----------------------------------------------------------------------

*     Longitude-velocity map; the y-axis is regularly spaced in
*     frequency but is to be labelled as a true relativistic velocity.
*       - PGSBOX uses subroutine LNGVEL.
*       - Separable (i.e. orthogonal), non-linear coordinate system.
*       - Automatic choice of coordinate increments.
*       - Extraction of a common scaling factor.
*       - Automatic choice of what edges to label.
*       - Request for tickmarks (internal) for one coordinate and grid
*         lines for the other.
*       - Simple two-colour grid using two calls with deferred
*         labelling on the first call.
*       - Degree labelling.
*       - Suppression of zero arcmin and arcsec fields in sexagesimal
*         degree format.

      WRITE (*, '(/,A)') 'Longitude-velocity map'

*     Reference pixel coordinates.
      NLDPRM(1) =   1D0
      NLDPRM(2) = 256D0

*     Reference pixel values.
      NLDPRM(3) = 0D0
      NLDPRM(4) = 1.420D9

*     Coordinate increments.
      NLDPRM(5) = 360D0/(NAXIS(1)-1)
      NLDPRM(6) = 4D6

*     Rest frequency.
      NLDPRM(7) = 1.420D9

*     Annotation.
      IDENTS(1) = 'galactic longitude'
      IDENTS(2) = 'velocity (m/s)'
      IDENTS(3) = 'HI line'

      OPT(1) = 'F'
      OPT(2) = ' '

*     Normal size lettering.
      CALL PGSCH (1.0*SCL)

*     Yellow tick marks for longitude and grid lines for velocity.
      CALL PGSCI (2)
      GCODE(1) = 1
      GCODE(2) = 2

*     Defer labelling.
      IC = -1
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 2D0,
     :  0, GRID1, 0, GRID2, .FALSE., LNGVEL, 1, 1, 7, ' ', NLIPRM,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Draw fiducial grid lines in white and do labelling.
      CALL PGSCI (1)
      GCODE(1) = 2
      GCODE(2) = 2
      GRID1(1) = 180D0
      GRID2(1) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 0, C0, GCODE, 0D0,
     :  1, GRID1, 1, GRID2, .FALSE., LNGVEL, 1, 1, 7, ' ', NLIPRM,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Draw the frame.
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL WCSLEEP (250)
      CALL PGPAGE ()

*-----------------------------------------------------------------------

*     Azimuth-frequency scan; this sort of output might be obtained from
*     an antenna that scans in azimuth with a receiver that scans
*     simultaneously in frequency.
*       - PGSBOX uses subroutine FSCAN.
*       - Non-separable (i.e. non-orthogonal) non-linear coordinate
*         system.
*       - Automatic choice of what edges to label; results in labelling
*         the bottom, left and right sides of the plot.
*       - Cyclic labelling.  FSCAN returns the azimuth in the range
*         0 - 720 degrees but PGSBOX is set to normalize this to two
*         cycles of 0 - 360 degrees.
*       - Logarithmic labelling.
*       - Automatic choice of coordinate increments but with request
*         for all grid lines for the logarithmic coordinate.
*       - Degree labelling.
*       - Suppression of common zero arcmin and arcsec fields in
*         sexagesimal degree format.

      WRITE (*, '(/,A)') 'Azimuth-frequency scan'

*     Reference pixel coordinates.
      NLDPRM(1) = 0.5D0
      NLDPRM(2) = 0.5D0

*     Reference pixel values.
      NLDPRM(3) = 0D0
      NLDPRM(4) = 8.5D0

*     Coordinate increments.
      NLDPRM(5) = 720D0/(NAXIS(1)+1)
      NLDPRM(6) = 0.002D0

*     Rate of change of NLDPRM(4) with x-pixel.
      NLDPRM(7) = -0.002D0

*     Annotation.
      IDENTS(1) = 'azimuth'
      IDENTS(2) = ESCAPE // 'gn/Hz'
      IDENTS(3) = 'Frequency/azimuth scan'

      OPT(1) = 'D'
      OPT(2) = 'L'

*     Normal size lettering.
      CALL PGSCH (1.0*SCL)

*     Draw full grid lines.
      CALL PGSCI (1)
      GCODE(1) = 2
      GCODE(2) = 2

*     Setting LABDEN = 9900 forces all logarithmic grid lines to be
*     drawn.
      IC = -1
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 9900, C0, GCODE, 2D0,
     :  0, GRID1, 0, GRID2, .FALSE., FSCAN, 1, 1, 7, ' ', NLIPRM,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Draw the frame.
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL WCSLEEP (250)
      CALL PGPAGE ()

*-----------------------------------------------------------------------

*     Z versus time plot.
*       - PGSBOX uses subroutine PGCRFN.
*       - Separable (i.e. orthogonal), non-linear coordinate system.
*       - Use of function PGCRFN for separable axis types.
*       - Automatic choice of what edges to label; results in labelling
*         the bottom and left sides of the plot.
*       - Automatic choice of coordinate increments.
*       - Logarithmic labelling over many orders of magnitude.
*       - Single-character annotation on a vertical axis is upright.

      WRITE (*, '(/,A)') 'Z versus time plot'

*     Function types.
      FCODE(1) = 'Lin '
      FCODE(2) = 'Log '

*     Reference pixel coordinates.
      NLDPRM(1) = 0.5D0
      NLDPRM(2) = -50D0

*     Coordinate increments.
      NLDPRM(3) = 0.04D0
      NLDPRM(4) = 0.02D0

*     Reference pixel values.
      NLDPRM(5) = -3.0D0
      NLDPRM(6) =  1.0D0

*     Annotation.
      IDENTS(1) = 'Age of universe (sec)'
      IDENTS(2) = 'Y'
      IDENTS(3) = ' '

      OPT(1) = 'L'
      OPT(2) = ' '

*     Normal size lettering.
      CALL PGSCH (1.0*SCL)

*     Draw ticks for the first coordinate, grid lines for the second.
      CALL PGSCI (1)
      GCODE(1) = 1
      GCODE(2) = 2

      IC = -1
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 0, C0, GCODE, 2D0,
     :  0, GRID1, 0, GRID2, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Draw the frame.
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL WCSLEEP (250)
      CALL PGPAGE ()

*-----------------------------------------------------------------------

*     Simple SIN projection near the south celestial pole.
*       - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*       - Non-separable (i.e. non-orthogonal) curvilinear coordinate
*         system.
*       - Demonstrate parameter definition for PGWCSL.
*       - Discovery of grid lines that do not cross any axis.
*       - Automatic choice of what edges to label; results in labelling
*         all sides of the plot.
*       - Automatic choice of coordinate increments but with request
*         for increased grid density for each coordinate.
*       - Double precision accuracy.
*       - Cyclic coordinates.  PGWCSL returns the right ascension in
*         the range -180 to +180 degrees, i.e. with a discontinuity
*         at +/- 180 degrees.
*       - Labelling of degrees as time in the range 0 - 24h.
*       - Suppression of labels that would overlap one another.
*       - Sexagesimal degree labelling with automatically determined
*         precision.
*       - Suppression of common zero minute and second fields in
*         sexagesimal time format.

      WRITE (*, '(/,A)') 'Simple SIN projection'

      STATUS = WCSINI (2, WCS)

*     Set projection type to SIN.
      CTYPE(1) = 'RA---SIN'
      CTYPE(2) = 'DEC--SIN'
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)

*     Reference pixel coordinates.
      STATUS = WCSPTD (WCS, WCS_CRPIX, 384D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX, 256D0, 2, 0)

*     Coordinate increments.
      STATUS = WCSPTD (WCS, WCS_CDELT, -1D0/3600000D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT,  1D0/3600000D0, 2, 0)

*     Spherical coordinate references.
      DEC0 = -89.99995D0
      STATUS = WCSPTD (WCS, WCS_CRVAL, 25D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, DEC0, 2, 0)

*     Set parameters for an NCP projection.
      DEC0 = DEC0*D2R
      STATUS = WCSPTD (WCS, WCS_PV, 0D0, 2, 1)
      STATUS = WCSPTD (WCS, WCS_PV, COS(DEC0)/SIN(DEC0), 2, 2)

*     Annotation.
      IDENTS(1) = 'Right ascension'
      IDENTS(2) = 'Declination'
      IDENTS(3) = 'WCS SIN projection'

      OPT(1) = 'G'
      OPT(2) = 'E'

*     Reduced size lettering.
      CALL PGSCH (0.8*SCL)

*     Draw full grid lines.
      CALL PGSCI (1)
      GCODE(1) = 2
      GCODE(2) = 2

*     Draw the celestial grid.  The grid density is set for each world
*     coordinate by specifying LABDEN = 1224.
      IC = -1
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 1224, C0, GCODE, 0D0,
     :  0, GRID1, 0, GRID2, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Draw the frame.
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL WCSLEEP (250)
      CALL PGPAGE ()

*-----------------------------------------------------------------------

*     Conic equal area projection.
*       - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*       - Non-separable (i.e. non-orthogonal) curvilinear coordinate
*         system.
*       - Coordinate system undefined in areas of the plot.
*       - Demonstrate parameter definition for PGWCSL.
*       - Discontinuous grid lines handled by PGWCSL.
*       - Discovery of grid lines that do not cross any axis.
*       - Colour control for grid and labelling.
*       - Reduced size lettering.
*       - Automatic choice of what edges to label; results in labelling
*         all sides of the plot.
*       - Automatic choice of coordinate increments.
*       - Cyclic coordinates.  PGWCSL returns the longitude in the
*         range -180 to +180 degrees, i.e. with a discontinuity at
*         +/- 180 degrees.
*       - Suppression of labels that would overlap one another.
*       - Suppression of common zero arcmin and arcsec fields in
*         sexagesimal degree format.

      WRITE (*, '(/,A)') 'Conic equal area projection'

      STATUS = WCSINI (2, WCS)

*     Set projection type to conic equal-area.
      CTYPE(1) = 'RA---COE'
      CTYPE(2) = 'DEC--COE'
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)

*     Reference pixel coordinates.
      STATUS = WCSPTD (WCS, WCS_CRPIX, 256D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX, 256D0, 2, 0)

*     Coordinate increments.
      STATUS = WCSPTD (WCS, WCS_CDELT, -1D0/3D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT,  1D0/3D0, 2, 0)

*     Spherical coordinate references.
      STATUS = WCSPTD (WCS, WCS_CRVAL, 90D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, 30D0, 2, 0)
      STATUS = WCSPTD (WCS, WCS_LONPOLE, 150D0, 0, 0)

*     Middle latitude and offset from standard parallels.
      STATUS = WCSPTD (WCS, WCS_PV, 60D0, 2, 1)
      STATUS = WCSPTD (WCS, WCS_PV, 15D0, 2, 2)

*     Annotation.
      IDENTS(1) = 'longitude'
      IDENTS(2) = 'latitude'
      IDENTS(3) = 'WCS conic equal area projection'

      OPT(1) = 'E'
      OPT(2) = 'E'

*     Reduced size lettering.
      CALL PGSCH (0.8*SCL)

*     Draw full grid lines.
      GCODE(1) = 2
      GCODE(2) = 2

*     Use colour to associate grid lines and labels.
*     Meridians in red.
      CALL PGSCR (10, 0.5, 0.0, 0.0)
*     Parallels in blue.
      CALL PGSCR (11, 0.0, 0.2, 0.5)
*     Longitudes in red.
      CALL PGSCR (12, 0.8, 0.3, 0.0)
*     Latitudes in blue.
      CALL PGSCR (13, 0.0, 0.4, 0.7)
*     Longitude labels in red.
      CALL PGSCR (14, 0.8, 0.3, 0.0)
*     Latitude labels in blue.
      CALL PGSCR (15, 0.0, 0.4, 0.7)
*     Title in cyan.
      CALL PGSCR (16, 0.3, 1.0, 1.0)

      CI(1) = 10
      CI(2) = 11
      CI(3) = 12
      CI(4) = 13
      CI(5) = 14
      CI(6) = 15
      CI(7) = 16

*     Draw the celestial grid letting PGSBOX choose the increments.
      IC = -1
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 0, CI, GCODE, 0D0,
     :  0, GRID1, 0, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Set parameters to draw the native grid.
      STATUS = WCSPTD (WCS, WCS_CRVAL,  0D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, 60D0, 2, 0)
      STATUS = WCSPTD (WCS, WCS_LONPOLE, 999D0, 0, 0)
      STATUS = WCSPTD (WCS, WCS_LATPOLE, 999D0, 0, 0)

*     We just want to delineate the boundary, in green.
      CALL PGSCI (7)
      GRID1(1) = -180D0
      GRID1(2) =  180D0
      GRID2(1) =  -90D0
      GRID2(2) =   90D0

      IC = -1
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0,
     :  2, GRID1, 2, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Draw the frame.
      CALL PGSCI (1)
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL WCSLEEP (250)
      CALL PGPAGE ()

*-----------------------------------------------------------------------

*     Polyconic projection with colour-coded grid.
*       - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*       - Non-separable (i.e. non-orthogonal) curvilinear coordinate
*         system.
*       - Coordinate system undefined in areas of the plot.
*       - Demonstrate parameter definition for PGWCSL.
*       - Discontinuous grid lines handled by PGWCSL.
*       - Colour coded labelling.
*       - Colour coded grid implemented by the caller.
*       - Basic management of the axis-crossing table (see code).
*       - Reduced size lettering.
*       - Tick marks external to the frame.
*       - User selection of what edges to label with request for both
*         coordinates to be labelled on bottom, left and top edges.
*       - User selection of grid lines to plot.
*       - Concatenation of annotation at bottom and left; automatically
*         suppressed at the top since only one coordinate is labelled
*         there.
*       - Suppression of labels that would overlap one another.
*       - Degree labelling.
*       - Labelling of degrees as time in the range -12 - +12h.
*       - Suppression of common zero minute and second fields in
*         sexagesimal time format.

      WRITE (*, '(/,A)') 'Polyconic projection with colour-coded grid'

      STATUS = WCSINI (2, WCS)

*     Set projection type to polyconic.
      CTYPE(1) = 'RA---PCO'
      CTYPE(2) = 'DEC--PCO'
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)

*     Reference pixel coordinates.
      STATUS = WCSPTD (WCS, WCS_CRPIX, 192D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX, 640D0, 2, 0)

*     Rotate 30 degrees.
      ROTN = 30D0*D2R
      STATUS = WCSPTD (WCS, WCS_PC,  COS(ROTN), 1, 1)
      STATUS = WCSPTD (WCS, WCS_PC,  SIN(ROTN), 1, 2)
      STATUS = WCSPTD (WCS, WCS_PC, -SIN(ROTN), 2, 1)
      STATUS = WCSPTD (WCS, WCS_PC,  COS(ROTN), 2, 2)

*     Coordinate increments.
      STATUS = WCSPTD (WCS, WCS_CDELT, -1D0/5D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT,  1D0/5D0, 2, 0)

*     Spherical coordinate references.
      STATUS = WCSPTD (WCS, WCS_CRVAL, 332D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL,  40D0, 2, 0)
      STATUS = WCSPTD (WCS, WCS_LONPOLE, -30D0, 0, 0)

*     Annotation.
      IDENTS(1) = 'Hour angle'
      IDENTS(2) = 'Declination'
      IDENTS(3) = 'WCS polyconic projection'

      OPT(1) = 'H'
      OPT(2) = 'B'

*     Reduced size lettering.
      CALL PGSCH (0.9*SCL)

*     Draw external (TIKLEN < 0) tick marks every 5 degrees.
      CALL PGSCI (6)
      GCODE(1) = 1
      GCODE(2) = 1
      TIKLEN = -2D0

      IC = -1
      GRID1(0) = 5D0
      GRID2(0) = 5D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, TIKLEN,
     :  0, GRID1, 0, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Resetting the table index to zero causes information about the
*     tick marks to be discarded.
      IC = 0

*     Draw full grid lines in yellow rather than tick marks.
      CALL PGSCI (2)
      GCODE(1) = 2
      GCODE(2) = 2

*     Draw the primary meridian and equator.
      GRID1(1) = 0D0
      GRID2(1) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0,
     :  1, GRID1, 1, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     At this point the axis-crossing table will have entries for the
*     primary meridian and equator.  Labelling was deferred in the
*     previous call, and the table is passed intact on the second call
*     to accumulate further axis-crossings.

*     Draw 90 degree meridians and poles in white.
      CALL PGSCI (3)
      GRID1(1) =  90D0
      GRID1(2) = 180D0
      GRID1(3) = 270D0
      GRID2(1) = -90D0
      GRID2(2) =  90D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0,
     :  3, GRID1, 2, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Draw the first set of 15 degree meridians and parallels in blue.
      CALL PGSCI (4)
      GRID1(1) =  15D0
      GRID1(2) =  60D0
      GRID1(3) = 105D0
      GRID1(4) = 150D0
      GRID1(5) = 195D0
      GRID1(6) = 240D0
      GRID1(7) = 285D0
      GRID1(8) = 330D0
      GRID2(1) = -75D0
      GRID2(2) = -30D0
      GRID2(3) =  15D0
      GRID2(4) =  60D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0,
     :  8, GRID1, 4, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Draw the second set of 15 degree meridians and parallels in red.
      CALL PGSCI (5)
      GRID1(1) =  30D0
      GRID1(2) =  75D0
      GRID1(3) = 120D0
      GRID1(4) = 165D0
      GRID1(5) = 210D0
      GRID1(6) = 255D0
      GRID1(7) = 300D0
      GRID1(8) = 345D0
      GRID2(1) = -60D0
      GRID2(2) = -15D0
      GRID2(3) =  30D0
      GRID2(4) =  75D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0,
     :  8, GRID1, 4, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     The axis-crossing table has now accumulated information for all of
*     the preceding meridians and parallels but no labels have been
*     produced.  It will acquire information for the the next set of
*     meridians and parallels before being processed by this call to
*     PGSBOX which finally produces the labels.

*     Draw the 45 degree meridians and parallels in grey and use colour
*     to differentiate grid labels.
*     Meridians and parallels in grey.
      CALL PGSCR (10, 0.7, 0.7, 0.7)
      CALL PGSCR (11, 0.7, 0.7, 0.7)
*     Longitudes tinged red.
      CALL PGSCR (12, 1.0, 0.9, 0.6)
*     Latitudes tinged green.
      CALL PGSCR (13, 0.8, 1.0, 0.9)
*     Longitude labels tinged red.
      CALL PGSCR (14, 1.0, 0.9, 0.6)
*     Latitude labels tinged green.
      CALL PGSCR (15, 0.8, 1.0, 0.9)
*     Title in white.
      CALL PGSCR (16, 1.0, 1.0, 1.0)

      CI(1) = 10
      CI(2) = 11
      CI(3) = 12
      CI(4) = 13
      CI(5) = 14
      CI(6) = 15
      CI(7) = 16

      CALL PGSCI (6)

*     Tell PGSBOX what edges to label.
      GRID1(1) =  45D0
      GRID1(2) = 135D0
      GRID1(3) = 225D0
      GRID1(4) = 315D0
      GRID2(1) = -45D0
      GRID2(2) =  45D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 2333, 0, CI, GCODE, 0D0,
     :  4, GRID1, 2, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Native grid in green (delineates boundary).
      CALL PGSCI (7)
      GRID1(1) = -180D0
      GRID1(2) =  180D0
      GRID2(1) = -999D0

      STATUS = WCSPTD (WCS, WCS_CRVAL, 0D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, 0D0, 2, 0)
      STATUS = WCSPTD (WCS, WCS_LONPOLE, 999D0, 0, 0)

      IC = -1
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0,
     :  2, GRID1, 1, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Draw the frame.
      CALL PGSCI (1)
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL WCSLEEP (250)
      CALL PGPAGE ()

*-----------------------------------------------------------------------

*     Plate Carree projection.
*       - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*       - Rectangular image.
*       - Dual coordinate grids.
*       - Non-separable (i.e. non-orthogonal) curvilinear coordinate
*         system.
*       - Demonstrate parameter definition for PGWCSL.
*       - Discontinuous grid lines handled by PGWCSL.
*       - Colour coding of grid and labelling.
*       - Reduced size lettering.
*       - Manual labelling control.
*       - Manual and automatic choice of coordinate increments.
*       - Cyclic coordinates.  PGWCSL returns the longitude in the
*         range -180 to +180 degrees, i.e. with a discontinuity at
*         +/- 180 degrees.
*       - Suppression of labels that would overlap one another.

      WRITE (*, '(/,A)') 'Plate Carree projection'

      NAXIS(1) = 181
      NAXIS(2) =  91

      BLC(1) = 0.5
      BLC(2) = 0.5
      TRC(1) = NAXIS(1) + 0.5
      TRC(2) = NAXIS(2) + 0.5

*     Reset viewport for rectangular image.

      IF (LARGE) THEN
        CALL PGVSTD ()
      ELSE
        CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
      END IF
      CALL PGWNAD (0.0, 1.0, 0.0, REAL(NAXIS(2))/REAL(NAXIS(1)))

      STATUS = WCSINI (2, WCS)

*     Set projection type to plate carree.
      CTYPE(1) = 'GLON-CAR'
      CTYPE(2) = 'GLAT-CAR'
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)

*     Reference pixel coordinates.
      STATUS = WCSPTD (WCS, WCS_CRPIX, 226D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX,  46D0, 2, 0)

*     Linear transformation matrix.
      ROTN = 15D0*D2R
      STATUS = WCSPTD (WCS, WCS_PC,  COS(ROTN), 1, 1)
      STATUS = WCSPTD (WCS, WCS_PC,  SIN(ROTN), 1, 2)
      STATUS = WCSPTD (WCS, WCS_PC, -SIN(ROTN), 2, 1)
      STATUS = WCSPTD (WCS, WCS_PC,  COS(ROTN), 2, 2)

*     Coordinate increments.
      STATUS = WCSPTD (WCS, WCS_CDELT, -1D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT,  1D0, 2, 0)

*     Set parameters to draw the native grid.
      STATUS = WCSPTD (WCS, WCS_CRVAL, 0D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, 0D0, 2, 0)

*     The reference pixel was defined so that the native longitude runs
*     from 225 deg to 45 deg and this will cause the grid to be
*     truncated at the 180 deg boundary.  However, being a cylindrical
*     projection it is possible to recentre it in longitude.  CYLFIX
*     will modify CRPIX, CRVAL, and LONPOLE to suit.

      CALL CYLFIX (NAXIS, WCS)

*     Annotation.
      IDENTS(1) = ' '
      IDENTS(2) = ' '
      IDENTS(3) = 'WCS plate caree projection'

      OPT(1) = 'C'
      OPT(2) = 'C'

*     Reduced size lettering.
      CALL PGSCH (0.8*SCL)

*     Draw full grid lines.
      GCODE(1) = 2
      GCODE(2) = 2

*     Draw native grid in green.
      CALL PGSCR (16, 0.0, 0.2, 0.0)
*     Title in cyan.
      CALL PGSCR (17, 0.3, 1.0, 1.0)

      CI(1) = 16
      CI(2) = 16
      CI(3) =  7
      CI(4) =  7
      CI(5) = -1
      CI(6) = -1
      CI(7) = 17

      IC = -1
      GRID1(0) = 15D0
      GRID2(0) = 15D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 2100, 0, CI, GCODE, 0D0,
     :  0, GRID1, 0, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Reset CRPIX previously modified by CYLFIX.
      STATUS = WCSPTD (WCS, WCS_CRPIX, 226D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX,  46D0, 2, 0)

*     Galactic reference coordinates.
      STATUS = WCSPTD (WCS, WCS_CRVAL, 30D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, 35D0, 2, 0)
      STATUS = WCSPTD (WCS, WCS_LONPOLE, 999D0, 0, 0)

      CALL CYLFIX (NAXIS, WCS)

*     Annotation.
      IDENTS(1) = 'longitude'
      IDENTS(2) = 'latitude'
      IDENTS(3) = ' '

      OPT(1) = 'E'
      OPT(2) = 'E'

*     Use colour to associate grid lines and labels.
*     Meridians in red.
      CALL PGSCR (10, 0.5, 0.0, 0.0)
*     Parallels in blue.
      CALL PGSCR (11, 0.0, 0.2, 0.5)
*     Longitudes in red.
      CALL PGSCR (12, 0.8, 0.3, 0.0)
*     Latitudes in blue.
      CALL PGSCR (13, 0.0, 0.4, 0.7)
*     Longitude labels in red.
      CALL PGSCR (14, 0.8, 0.3, 0.0)
*     Latitude labels in blue.
      CALL PGSCR (15, 0.0, 0.4, 0.7)

      CI(1) = 10
      CI(2) = 11
      CI(3) = 12
      CI(4) = 13
      CI(5) = 14
      CI(6) = 15
      CI(7) = -1

*     Draw the celestial grid letting PGSBOX choose the increments.
      IC = -1
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 21, 0, CI, GCODE, 0D0,
     :  0, GRID1, 0, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Draw the frame.
      CALL PGSCI (1)
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL WCSLEEP (250)
      CALL PGPAGE ()

*-----------------------------------------------------------------------

*     Plate Carree projection.
*       - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*       - BLC, TRC unrelated to pixel coordinates.
*       - Demonstrate parameter definition for PGWCSL.
*       - Poles and 180 meridian projected along edges of the frame.
*       - Reduced size lettering.
*       - Manual and automatic choice of coordinate increments.
*       - Suppression of common zero minute and second fields in
*         sexagesimal time format.

      WRITE (*, '(/,A)') 'Plate Carree projection'

      BLC(1) = -180.0
      BLC(2) =  -90.0
      TRC(1) =  180.0
      TRC(2) =  +90.0

*     Reset viewport for rectangular image.

      IF (LARGE) THEN
        CALL PGVSTD ()
      ELSE
        CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
      END IF
      CALL PGWNAD (BLC(1), TRC(1), BLC(2), TRC(2))

      STATUS = WCSINI (2, WCS)

*     Set projection type to plate carree.
      CTYPE(1) = 'RA---CAR'
      CTYPE(2) = 'DEC--CAR'
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)

*     Reference pixel coordinates.
      STATUS = WCSPTD (WCS, WCS_CRPIX, 0D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX, 0D0, 2, 0)

*     Coordinate increments.
      STATUS = WCSPTD (WCS, WCS_CDELT, -1D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT,  1D0, 2, 0)

*     Set parameters to draw the native grid.
      STATUS = WCSPTD (WCS, WCS_CRVAL, 0D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, 0D0, 2, 0)

*     Annotation.
      IDENTS(1) = 'Right ascension'
      IDENTS(2) = 'Declination'
      IDENTS(3) = 'WCS plate caree projection'

      OPT(1) = 'G'
      OPT(2) = 'E'

*     Reduced size lettering.
      CALL PGSCH (0.7*SCL)

*     Draw full grid lines.
      GCODE(1) = 2
      GCODE(2) = 2

      CALL PGSCI (1)

      IC = -1
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 2121, 1212, C0, GCODE, 0D0,
     :  0, GRID1, 0, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Draw the frame.
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL WCSLEEP (250)
      CALL PGPAGE ()

*-----------------------------------------------------------------------

*     Cylindrical perspective projection.
*       - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*       - BLC, TRC unrelated to pixel coordinates.
*       - Demonstrate parameter definition for PGWCSL.
*       - Reduced size lettering.
*       - Manual and automatic choice of coordinate increments.
*       - Suppression of common zero minute and second fields in
*         sexagesimal time format.

      WRITE (*, '(/,A)') 'Cylindrical perspective projection'

      BLC(1) = -180.0
      BLC(2) =  -90.0
      TRC(1) =  180.0
      TRC(2) =  +90.0

*     Reset viewport for rectangular image.

      IF (LARGE) THEN
        CALL PGVSTD ()
      ELSE
        CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
      END IF
      CALL PGWNAD (BLC(1), TRC(1), BLC(2), TRC(2))

      STATUS = WCSINI (2, WCS)

*     Set projection type to cylindrical perspective.
      CTYPE(1) = 'RA---CYP'
      CTYPE(2) = 'DEC--CYP'
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)

*     Reference pixel coordinates.
      STATUS = WCSPTD (WCS, WCS_CRPIX, 0D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX, 0D0, 2, 0)

*     Coordinate increments.
      STATUS = WCSPTD (WCS, WCS_CDELT, -1D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT,  1D0, 2, 0)

*     Set parameters to draw the native grid.
      STATUS = WCSPTD (WCS, WCS_CRVAL,  45D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, -90D0, 2, 0)
      STATUS = WCSPTD (WCS, WCS_LONPOLE, 999D0, 0, 0)

*     mu and lambda projection parameters.
      STATUS = WCSPTD (WCS, WCS_PV, 0D0, 2, 1)
      STATUS = WCSPTD (WCS, WCS_PV, 1D0, 2, 2)

*     Annotation.
      IDENTS(1) = 'Right ascension'
      IDENTS(2) = 'Declination'
      IDENTS(3) = 'WCS cylindrical perspective projection'

      OPT(1) = 'G'
      OPT(2) = 'E'

*     Reduced size lettering.
      CALL PGSCH (0.7*SCL)

*     Draw full grid lines.
      GCODE(1) = 2
      GCODE(2) = 2

      CALL PGSCI (1)

      IC = -1
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 2121, 1212, C0, GCODE, 0D0,
     :  0, GRID1, 0, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Draw the frame.
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL WCSLEEP (250)
      CALL PGPAGE ()

*-----------------------------------------------------------------------
*     Gnomonic projection.
*       - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*       - Demonstrate parameter definition for PGWCSL.
*       - Reduced size lettering.
*       - Manual and automatic choice of coordinate increments.
*       - Suppression of common zero minute and second fields in
*         sexagesimal time format.

      WRITE (*, '(/,A)') 'TAN projection'

      NAXIS(1) = 100
      NAXIS(2) = 100

      BLC(1) = 0.5
      BLC(2) = 0.5
      TRC(1) = NAXIS(1) + 0.5
      TRC(2) = NAXIS(2) + 0.5

*     Reset viewport for rectangular image.
      IF (LARGE) THEN
        CALL PGVSTD ()
      ELSE
        CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
      END IF
      CALL PGWNAD (0.0, 1.0, 0.0, REAL(NAXIS(2))/REAL(NAXIS(1)))

      STATUS = WCSINI (2, WCS)

*     Set projection type to gnomonic.
      CTYPE(1) = 'RA---TAN'
      CTYPE(2) = 'DEC--TAN'
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)

*     Reference pixel coordinates.
      STATUS = WCSPTD (WCS, WCS_CRPIX, 50.5D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX,  1.0D0, 2, 0)

*     Coordinate increments.
      STATUS = WCSPTD (WCS, WCS_CDELT, 1D-3, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT, 1D-3, 2, 0)

*     Set parameters to draw the native grid.
      STATUS = WCSPTD (WCS, WCS_CRVAL, -45.0D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, -89.7D0, 2, 0)
      STATUS = WCSPTD (WCS, WCS_LONPOLE, 999D0, 0, 0)

*     Annotation.
      IDENTS(1) = 'Right ascension'
      IDENTS(2) = 'Declination'
      IDENTS(3) = 'WCS TAN projection'

      OPT(1) = 'E'
      OPT(2) = 'E'

*     Reduced size lettering.
      CALL PGSCH (0.7*SCL)

*     Draw full grid lines.
      GCODE(1) = 2
      GCODE(2) = 2

      CALL PGSCI (1)

      IC = -1
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 1212, C0, GCODE, 0D0,
     :  0, GRID1, 0, GRID2, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Draw the frame.
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL WCSLEEP (250)
      CALL PGPAGE ()

*-----------------------------------------------------------------------

*     Linear-linear plot with two types of alternative labelling.
*       - PGSBOX uses subroutine PGCRFN.
*       - Separable (i.e. orthogonal), linear coordinate system.
*       - Use of function PGCRFN for separable axis types.
*       - Alternative labelling and axis annotation.
*       - Direct manipulation of the axis-crossing table.
*       - Tick mark and grid line control.
*       - User selection of what edges to label.
*       - Automatic choice of coordinate increments.

      WRITE (*, '(/,A)') 'Linear plot with alternative labelling'

      IF (LARGE) THEN
        CALL PGVSTD ()
      ELSE
        CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
      END IF
      CALL PGWNAD (0.0, 1.0, 0.0, 1.0)

      NAXIS(1) = 512
      NAXIS(2) = 512

      BLC(1) = 0.5
      BLC(2) = 0.5
      TRC(1) = NAXIS(1) + 0.5
      TRC(2) = NAXIS(2) + 0.5

*     Function types.
      FCODE(1) = 'Lin '
      FCODE(2) = 'Lin '

*     Reference pixel coordinates.
      NLDPRM(1) = 0.5D0
      NLDPRM(2) = 0.5D0

*     Coordinate increments.
      NLDPRM(3) = 0.03D0
      NLDPRM(4) = 0.03D0

*     Reference pixel values.
      NLDPRM(5) = 20D0
      NLDPRM(6) =  0D0

*     Annotation.
      IDENTS(1) = 'temperature of frog (' // ESCAPE // 'uo' //
     :            ESCAPE // 'dC)'
      IDENTS(2) = 'distance hopped (m)'
      IDENTS(3) = ' '

      OPT(1) = ' '
      OPT(2) = ' '

*     Reduced size lettering.
      CALL PGSCH (0.8*SCL)

*     Draw tick marks at the bottom for the first coordinate, grid lines
*     for the second.  Setting GCODE(1) = -1 inhibits information being
*     stored for labels on the top edge while GCODE(2) = 2 causes
*     information to be stored for labels on the right edge even if
*     those labels are not actually produced.
      CALL PGSCI (1)
      GCODE(1) = -1
      GCODE(2) = 2

*     Set LABCTL = 21 to label the bottom and left edges only.
      IC = -1
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 21, 0, C0, GCODE, 2D0,
     :  0, GRID1, 0, GRID2, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Information for labels on the right edge was stored in the
*     crossing table on the first call to PGSBOX.  We now want to
*     manipulate it to convert metres to feet.  Note that while it's
*     a simple matter to draw alternative sets of tick marks on opposite
*     edges of the frame, as with the two temperature scales, we have
*     the slightly more difficult requirement of labelling grid lines
*     with different values at each end.
      DO 10 J = 1, IC
*       Look for entries associated with the right edge of the frame.
        IF (CACHE(1,J).EQ.4D0) THEN
*         Convert to feet, rounding to the nearest 0.1.
          CACHE(4,J) = CACHE(4,J) * 1D3/(25.4*12D0)
          CACHE(4,J) = AINT(CACHE(4,J)*10D0 + 0.5D0)/10D0
        END IF
 10   CONTINUE

*     Annotation for the right edge.
      IDENTS(1) = ' '
      IDENTS(2) = '(feet)'

*     Set LABCTL = 12000 to label the right edge with the second
*     coordinate without redrawing the grid lines.
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 12000, 0, C0, GCODE, 2D0,
     :  0, GRID1, 0, GRID2, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     The alternative temperature scale in Fahrenheit is to be
*     constructed with a new set of tick marks.
      NLDPRM(3) = NLDPRM(3)*1.8D0
      NLDPRM(5) = NLDPRM(5)*1.8D0 + 32D0

*     Draw tick marks at the top for the first coordinate, don't redo
*     grid lines for the second.
      GCODE(1) = -100
      GCODE(2) = 0

*     Annotation for the top edge.
      IDENTS(1) = '(' // ESCAPE // 'uo' // ESCAPE // 'dF)'
      IDENTS(2) = ' '

*     Set LABCTL = 100 to label the top edge; Set IC = -1 to redetermine
*     the coordinate extrema.
      IC = -1
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 100, 0, C0, GCODE, 2D0,
     :  0, GRID1, 0, GRID2, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Draw the frame.
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL WCSLEEP (250)
      CALL PGPAGE ()

*-----------------------------------------------------------------------

*     Calendar axes using subroutine PGLBOX.
*       - Separable (i.e. orthogonal), linear coordinate system.
*       - Use of PGLBOX for simple linear axis types.
*       - Automatic choice of what edges to label; results in labelling
*         the bottom and left sides of the plot.
*       - Automatic choice of coordinate increments.
*       - Calendar date axis labelling.
*       - Single-character annotation on a vertical axis is upright.

      WRITE (*, '(/,A)') 'Calendar axes using subroutine PGLBOX'

      CALL PGSWIN (51900.0, 52412.0, 51900.0, 57020.0)

*     Annotation.
      IDENTS(1) = 'Date started'
      IDENTS(2) = 'Date finished'
      IDENTS(3) = 'Calendar axes using subroutine PGLBOX'

      OPT(1) = 'Y'
      OPT(2) = 'Y'

*     Reduced size lettering.
      CALL PGSCH (0.7*SCL)

*     Draw tick marks on each axis.
      CALL PGSCI (1)
      GCODE(1) = 1
      GCODE(2) = 1

      IC = -1
      CALL PGLBOX (IDENTS, OPT, 0, 0, C0, GCODE, 2D0, 0, 0D0, 0, 0D0,
     :  .FALSE., 256, IC, CACHE, IERR)

*     Draw the frame.
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL WCSLEEP (250)
      CALL PGPAGE ()

*-----------------------------------------------------------------------

*     Simple linear axes handled by PGWCSL.
*       - Separable (i.e. orthogonal), linear coordinate system.
*       - Automatic choice of what edges to label; results in labelling
*         the bottom and left sides of the plot.
*       - Automatic choice of coordinate increments.
*       - Tick marks and labels at the edges of the frame.
*       - Single-character annotation on a vertical axis is upright.

      WRITE (*, '(/,A)') 'Simple linear axes handled by pgwcsl()'

      NAXIS(1) = 3
      NAXIS(2) = 3

      BLC(1) = 0.5
      BLC(2) = 0.5
      TRC(1) = NAXIS(1) + 0.5
      TRC(2) = NAXIS(2) + 0.5

      STATUS = WCSINI (2, WCS)

      CTYPE(1) = 'X'
      CTYPE(2) = 'Y'
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)

*     Reference pixel coordinates.
      STATUS = WCSPTD (WCS, WCS_CRPIX, 2D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX, 2D0, 2, 0)

*     Coordinate increments.
      STATUS = WCSPTD (WCS, WCS_CDELT, 1D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT, 1D0, 2, 0)

*     Spherical coordinate references.
      STATUS = WCSPTD (WCS, WCS_CRVAL, 2D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, 2D0, 2, 0)

*     Annotation.
      IDENTS(1) = 'X'
      IDENTS(2) = 'Y'
      IDENTS(3) = 'Simple linear axes handled by pgwcsl()'

      OPT(1) = ' '
      OPT(2) = ' '

*     Reduced size lettering.
      CALL PGSCH (0.8*SCL)

*     Draw tick marks only.
      CALL PGSCI (1)
      GCODE(1) = 1
      GCODE(2) = 1

      IC = -1
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 0, C0, GCODE, 2D0,
     :  0, GRID1, 0, GRID2, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)

*     Draw the frame.
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL WCSLEEP (250)
      CALL PGPAGE ()

*-----------------------------------------------------------------------

      STATUS = WCSFREE (WCS)

      CALL PGASK (.FALSE.)
      CALL PGEND ()

      END
@


8.6
log
@WCSLIB 8.6
@
text
@d3 1
a3 1
* PGSBOX 8.6 - draw curvilinear coordinate axes for PGPLOT.
d23 1
a23 1
* $Id: pgtest.f,v 8.5.1.2 2026/03/29 13:36:20 mcalabre Exp mcalabre $
@


8.6.1.1
log
@WCSLIB 8.6 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 8.6 2026/03/29 13:53:57 mcalabre Exp $
@


8.5
log
@WCSLIB 8.5
@
text
@d3 2
a4 2
* PGSBOX 8.5 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2025, Mark Calabretta
d22 2
a23 2
* http://www.atnf.csiro.au/people/Mark.Calabretta
* $Id: pgtest.f,v 8.4.1.2 2025/12/06 13:37:20 mcalabre Exp mcalabre $
@


8.5.1.1
log
@WCSLIB 8.5 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 8.5 2025/12/06 13:47:42 mcalabre Exp $
@


8.5.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2026, Mark Calabretta
d22 2
a23 2
* http://www.atnf.csiro.au/computing/software/wcs
* $Id: pgtest.f,v 8.5.1.1 2025/12/06 14:05:17 mcalabre Exp mcalabre $
@


8.4
log
@WCSLIB 8.4
@
text
@d3 2
a4 2
* PGSBOX 8.4 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2024, Mark Calabretta
d23 1
a23 1
* $Id: pgtest.f,v 8.3.1.1 2024/05/13 16:35:40 mcalabre Exp mcalabre $
@


8.4.1.1
log
@WCSLIB 8.4 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 8.4 2024/10/28 13:56:18 mcalabre Exp $
@


8.4.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2025, Mark Calabretta
d23 1
a23 1
* $Id: pgtest.f,v 8.4.1.1 2024/10/28 13:58:42 mcalabre Exp mcalabre $
@


8.3
log
@WCSLIB 8.3
@
text
@d3 1
a3 1
* PGSBOX 8.3 - draw curvilinear coordinate axes for PGPLOT.
d23 1
a23 1
* $Id: pgtest.f,v 8.2.1.2 2024/04/02 13:43:15 mcalabre Exp mcalabre $
@


8.3.1.1
log
@WCSLIB 8.3 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 8.3 2024/05/13 16:33:01 mcalabre Exp $
@


8.2
log
@WCSLIB 8.2
@
text
@d3 2
a4 2
* PGSBOX 8.2 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2023, Mark Calabretta
d23 1
a23 1
* $Id: pgtest.f,v 8.1.1.1 2023/07/05 17:15:01 mcalabre Exp mcalabre $
@


8.2.1.1
log
@WCSLIB 8.2 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 8.2 2023/11/16 10:00:31 mcalabre Exp $
@


8.2.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2024, Mark Calabretta
d23 1
a23 1
* $Id: pgtest.f,v 8.2.1.1 2023/11/16 10:05:59 mcalabre Exp mcalabre $
@


8.1
log
@WCSLIB 8.1
@
text
@d3 1
a3 1
* PGSBOX 8.1 - draw curvilinear coordinate axes for PGPLOT.
d23 1
a23 1
* $Id: pgtest.f,v 8.0.1.3 2023/07/04 15:15:23 mcalabre Exp mcalabre $
@


8.1.1.1
log
@WCSLIB 8.1 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 8.1 2023/07/05 17:12:07 mcalabre Exp $
@


8.0
log
@WCSLIB 8.0
@
text
@d3 1
a3 1
* PGSBOX 8.0 - draw curvilinear coordinate axes for PGPLOT.
d23 1
a23 1
* $Id: pgtest.f,v 7.13.1.2 2023/06/15 04:17:22 mcalabre Exp mcalabre $
d163 1
a163 1
     :  0, GRID1, 0, GRID2, .FALSE., LNGVEL, 1, 1, 7, ' ', 0,
d173 1
a173 1
     :  1, GRID1, 1, GRID2, .FALSE., LNGVEL, 1, 1, 7, ' ', 0,
d179 1
a179 1

d241 1
a241 1
     :  0, GRID1, 0, GRID2, .FALSE., FSCAN, 1, 1, 7, ' ', 0,
d247 1
d306 1
d388 1
d509 1
d722 1
d886 1
d964 1
d1046 1
d1124 1
d1255 1
d1296 1
d1364 1
@


8.0.1.1
log
@WCSLIB 8.0 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 8.0 2023/07/01 06:15:05 mcalabre Exp $
@


8.0.1.2
log
@Quell compiler warnings from gfortran 12.1.0.  Added SLEEP(1) calls to
slow down the plots, otherwise they were just a blur.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 8.0.1.1 2023/07/01 06:22:45 mcalabre Exp mcalabre $
d163 1
a163 1
     :  0, GRID1, 0, GRID2, .FALSE., LNGVEL, 1, 1, 7, ' ', NLIPRM,
d173 1
a173 1
     :  1, GRID1, 1, GRID2, .FALSE., LNGVEL, 1, 1, 7, ' ', NLIPRM,
d179 1
a179 1
      CALL SLEEP (1)
d241 1
a241 1
     :  0, GRID1, 0, GRID2, .FALSE., FSCAN, 1, 1, 7, ' ', NLIPRM,
a246 1
      CALL SLEEP (1)
a304 1
      CALL SLEEP (1)
a385 1
      CALL SLEEP (1)
a505 1
      CALL SLEEP (1)
a717 1
      CALL SLEEP (1)
a880 1
      CALL SLEEP (1)
a957 1
      CALL SLEEP (1)
a1038 1
      CALL SLEEP (1)
a1115 1
      CALL SLEEP (1)
a1245 1
      CALL SLEEP (1)
a1285 1
      CALL SLEEP (1)
a1352 1
      CALL SLEEP (1)
@


8.0.1.3
log
@Replaced the SLEEP(1) calls, which were too slow, with calls to the new
Fortran test utility function, WCSLEEP(250).
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 8.0.1.2 2023/07/04 12:59:10 mcalabre Exp mcalabre $
d179 1
a179 1
      CALL WCSLEEP (250)
d247 1
a247 1
      CALL WCSLEEP (250)
d306 1
a306 1
      CALL WCSLEEP (250)
d388 1
a388 1
      CALL WCSLEEP (250)
d509 1
a509 1
      CALL WCSLEEP (250)
d722 1
a722 1
      CALL WCSLEEP (250)
d886 1
a886 1
      CALL WCSLEEP (250)
d964 1
a964 1
      CALL WCSLEEP (250)
d1046 1
a1046 1
      CALL WCSLEEP (250)
d1124 1
a1124 1
      CALL WCSLEEP (250)
d1255 1
a1255 1
      CALL WCSLEEP (250)
d1296 1
a1296 1
      CALL WCSLEEP (250)
d1364 1
a1364 1
      CALL WCSLEEP (250)
@


7.13
log
@WCSLIB 7.13
@
text
@d3 2
a4 2
* PGSBOX 7.13 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2022, Mark Calabretta
d23 1
a23 1
* $Id: pgtest.f,v 7.12.1.1 2022/09/09 05:00:36 mcalabre Exp mcalabre $
@


7.13.1.1
log
@WCSLIB 7.13 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 7.13 2022/10/07 03:43:35 mcalabre Exp $
@


7.13.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2023, Mark Calabretta
d23 1
a23 1
* $Id: pgtest.f,v 7.13.1.1 2022/10/07 03:45:47 mcalabre Exp mcalabre $
@


7.12
log
@WCSLIB 7.12
@
text
@d3 1
a3 1
* PGSBOX 7.12 - draw curvilinear coordinate axes for PGPLOT.
d23 1
a23 1
* $Id: pgtest.f,v 7.11.1.1 2022/04/26 06:15:33 mcalabre Exp mcalabre $
@


7.12.1.1
log
@WCSLIB 7.12 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 7.12 2022/09/09 04:57:58 mcalabre Exp $
@


7.11
log
@WCSLIB 7.11
@
text
@d3 1
a3 1
* PGSBOX 7.11 - draw curvilinear coordinate axes for PGPLOT.
d23 1
a23 1
* $Id: pgtest.f,v 7.10.1.1 2022/04/23 14:26:55 mcalabre Exp mcalabre $
@


7.11.1.1
log
@WCSLIB 7.11 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 7.11 2022/04/26 06:13:52 mcalabre Exp $
@


7.10
log
@WCSLIB 7.10
@
text
@d3 1
a3 1
* PGSBOX 7.10 - draw curvilinear coordinate axes for PGPLOT.
d23 1
a23 1
* $Id: pgtest.f,v 7.9.1.1 2022/03/25 15:16:40 mcalabre Exp mcalabre $
@


7.10.1.1
log
@WCSLIB 7.10 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 7.10 2022/04/23 14:25:15 mcalabre Exp $
@


7.9
log
@WCSLIB 7.9
@
text
@d3 1
a3 1
* PGSBOX 7.9 - draw curvilinear coordinate axes for PGPLOT.
d23 1
a23 1
* $Id: pgtest.f,v 7.8.1.1 2022/03/25 06:14:15 mcalabre Exp mcalabre $
@


7.9.1.1
log
@WCSLIB 7.9 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 7.9 2022/03/25 15:14:49 mcalabre Exp $
@


7.8
log
@WCSLIB 7.8
@
text
@d3 1
a3 1
* PGSBOX 7.8 - draw curvilinear coordinate axes for PGPLOT.
d23 1
a23 1
* $Id: pgtest.f,v 7.7.1.2 2022/03/25 06:06:01 mcalabre Exp mcalabre $
@


7.8.1.1
log
@WCSLIB 7.8 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 7.8 2022/03/25 06:11:50 mcalabre Exp $
@


7.7
log
@WCSLIB 7.7
@
text
@d3 2
a4 2
* PGSBOX 7.7 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2021, Mark Calabretta
d23 1
a23 1
* $Id: pgtest.f,v 7.6.1.4 2021/07/02 02:56:19 mcalabre Exp mcalabre $
@


7.7.1.1
log
@WCSLIB 7.7 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 7.7 2021/07/12 06:36:49 mcalabre Exp $
@


7.7.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2022, Mark Calabretta
d23 1
a23 1
* $Id: pgtest.f,v 7.7.1.1 2021/07/12 06:38:28 mcalabre Exp mcalabre $
@


7.6
log
@WCSLIB 7.6
@
text
@d3 1
a3 1
* PGSBOX 7.6 - draw curvilinear coordinate axes for PGPLOT.
d23 1
a23 1
* $Id: pgtest.f,v 7.5.1.1 2021/03/20 05:56:42 mcalabre Exp mcalabre $
a46 3

*     A portability fix - does '\' itself need to be escaped?
      DATA ESCAPE /'\\'/
d57 3
a59 1
*      CALL PGBEG (0, '?', 1, 1)
d160 5
a164 3
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 2D0, 0, 0D0,
     :  0, 0D0, .FALSE., LNGVEL, 1, 1, 7, ' ', 0, NLDPRM, 256, IC,
     :  CACHE, IERR)
d172 3
a174 3
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 0, C0, GCODE, 0D0, 1,
     :  GRID1, 1, GRID2, .FALSE., LNGVEL, 1, 1, 7, ' ', 0, NLDPRM, 256,
     :  IC, CACHE, IERR)
d238 5
a242 3
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 9900, C0, GCODE, 2D0, 0,
     :  0D0, 0, 0D0, .FALSE., FSCAN, 1, 1, 7, ' ', 0, NLDPRM, 256, IC,
     :  CACHE, IERR)
d296 5
a300 3
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 0, C0, GCODE, 2D0, 0, 0D0,
     :  0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM, 256,
     :  IC, CACHE, IERR)
d377 4
a380 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 1224, C0, GCODE, 0D0, 0,
     :  0D0, 0, 0D0, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d478 5
a482 3
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 0, CI, GCODE, 0D0, 0, 0D0,
     :  0, 0D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM, 256,
     :  IC, CACHE, IERR)
d498 2
a499 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0, 2,
     :  GRID1, 2, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d582 5
a586 3
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, TIKLEN, 0,
     :  5D0, 0, 5D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM,
     :  256, IC, CACHE, IERR)
d600 2
a601 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0, 1,
     :  GRID1, 1, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d616 2
a617 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0, 3,
     :  GRID1, 2, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d634 2
a635 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0, 8,
     :  GRID1, 4, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d652 2
a653 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0, 8,
     :  GRID1, 4, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d695 2
a696 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 2333, 0, CI, GCODE, 0D0, 4,
     :  GRID1, 2, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d710 2
a711 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0, 2,
     :  GRID1, 1, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d822 4
a825 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 2100, 0, CI, GCODE, 0D0, 0,
     :  15D0, 0, 15D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d871 5
a875 3
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 21, 0, CI, GCODE, 0D0, 0, 0D0,
     :  0, 0D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM, 256,
     :  IC, CACHE, IERR)
d949 5
a953 3
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 2121, 1212, C0, GCODE, 0D0, 0,
     :  0D0, 0, 0D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM,
     :  256, IC, CACHE, IERR)
d1030 5
a1034 3
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 2121, 1212, C0, GCODE, 0D0, 0,
     :  0D0, 0, 0D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM,
     :  256, IC, CACHE, IERR)
d1107 4
a1110 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 1212, C0, GCODE, 0D0, 0,
     :  0D0, 0, 0D0, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d1186 5
a1190 3
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 21, 0, C0, GCODE, 2D0, 0, 0D0,
     :  0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM, 256,
     :  IC, CACHE, IERR)
d1214 5
a1218 3
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 12000, 0, C0, GCODE, 2D0, 0,
     :  0D0, 0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM,
     :  256, IC, CACHE, IERR)
d1237 5
a1241 3
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 100, 0, C0, GCODE, 2D0, 0,
     :  0D0, 0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM,
     :  256, IC, CACHE, IERR)
d1344 4
a1347 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 0, C0, GCODE, 2D0, 0,
     :  0D0, 0, 0D0, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d1359 1
a1359 1
      CALL PGASK (0)
@


7.6.1.1
log
@WCSLIB 7.6 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 7.6 2021/04/13 12:57:02 mcalabre Exp $
@


7.6.1.2
log
@Define ESCAPE using the CHAR() function to avoid the confusion that arises
by gfortran interpreting it as a C-style escape character.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 7.6.1.1 2021/04/13 13:00:37 mcalabre Exp mcalabre $
d47 3
d60 1
a60 3
      ESCAPE = CHAR(92)

*     CALL PGBEG (0, '?', 1, 1)
@


7.6.1.3
log
@Replaced a call to PGASK(0) with the formally required PGASK(.FALSE.).
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 7.6.1.2 2021/07/02 02:31:54 mcalabre Exp mcalabre $
d1329 1
a1329 1
      CALL PGASK (.FALSE.)
@


7.6.1.4
log
@gfortran 11.1.0 now compares the argument list in all calls to the same
subroutine or function that occur anywhere in a source file.  A compilation
error results if any arguments do not match in rank (scalar, vector, array,
etc.).  As there are no compiler options to suppress this annoying behaviour,
it is necessary to make the subroutine or function calls match, or at least
appear so.  This particularly affected calls to PGSBOX in relation to the
GRID1 and GRID2 arguments.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 7.6.1.3 2021/07/02 02:50:23 mcalabre Exp mcalabre $
d160 3
a162 5
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 2D0,
     :  0, GRID1, 0, GRID2, .FALSE., LNGVEL, 1, 1, 7, ' ', 0,
     :  NLDPRM, 256, IC, CACHE, IERR)
d170 3
a172 3
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 0, C0, GCODE, 0D0,
     :  1, GRID1, 1, GRID2, .FALSE., LNGVEL, 1, 1, 7, ' ', 0,
     :  NLDPRM, 256, IC, CACHE, IERR)
d236 3
a238 5
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 9900, C0, GCODE, 2D0,
     :  0, GRID1, 0, GRID2, .FALSE., FSCAN, 1, 1, 7, ' ', 0,
     :  NLDPRM, 256, IC, CACHE, IERR)
d292 3
a294 5
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 0, C0, GCODE, 2D0,
     :  0, GRID1, 0, GRID2, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM,
     :  NLDPRM, 256, IC, CACHE, IERR)
d371 2
a372 4
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 1224, C0, GCODE, 0D0,
     :  0, GRID1, 0, GRID2, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d470 3
a472 5
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 0, CI, GCODE, 0D0,
     :  0, GRID1, 0, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
d488 2
a489 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0,
     :  2, GRID1, 2, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d572 3
a574 5
      GRID1(0) = 5D0
      GRID2(0) = 5D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, TIKLEN,
     :  0, GRID1, 0, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
d588 2
a589 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0,
     :  1, GRID1, 1, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d604 2
a605 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0,
     :  3, GRID1, 2, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d622 2
a623 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0,
     :  8, GRID1, 4, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d640 2
a641 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0,
     :  8, GRID1, 4, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d683 2
a684 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 2333, 0, CI, GCODE, 0D0,
     :  4, GRID1, 2, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d698 2
a699 2
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, -1, 0, C0, GCODE, 0D0,
     :  2, GRID1, 1, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d810 2
a811 4
      GRID1(0) = 15D0
      GRID2(0) = 15D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 2100, 0, CI, GCODE, 0D0,
     :  0, GRID1, 0, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d857 3
a859 5
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 21, 0, CI, GCODE, 0D0,
     :  0, GRID1, 0, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
d933 3
a935 5
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 2121, 1212, C0, GCODE, 0D0,
     :  0, GRID1, 0, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
d1012 3
a1014 5
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 2121, 1212, C0, GCODE, 0D0,
     :  0, GRID1, 0, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
d1087 2
a1088 4
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 1212, C0, GCODE, 0D0,
     :  0, GRID1, 0, GRID2, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
d1164 3
a1166 5
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 21, 0, C0, GCODE, 2D0,
     :  0, GRID1, 0, GRID2, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM,
     :  NLDPRM, 256, IC, CACHE, IERR)
d1190 3
a1192 5
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 12000, 0, C0, GCODE, 2D0,
     :  0, GRID1, 0, GRID2, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM,
     :  NLDPRM, 256, IC, CACHE, IERR)
d1211 3
a1213 5
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 100, 0, C0, GCODE, 2D0,
     :  0, GRID1, 0, GRID2, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM,
     :  NLDPRM, 256, IC, CACHE, IERR)
d1316 2
a1317 4
      GRID1(0) = 0D0
      GRID2(0) = 0D0
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 0, C0, GCODE, 2D0,
     :  0, GRID1, 0, GRID2, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
@


7.5
log
@WCSLIB 7.5
@
text
@d3 1
a3 1
* PGSBOX 7.5 - draw curvilinear coordinate axes for PGPLOT.
d23 1
a23 1
* $Id: pgtest.f,v 7.4.1.1 2021/01/31 02:26:24 mcalabre Exp mcalabre $
@


7.5.1.1
log
@WCSLIB 7.5 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 7.5 2021/03/20 05:54:59 mcalabre Exp $
@


7.4
log
@WCSLIB 7.4
@
text
@d3 1
a3 1
* PGSBOX 7.4 - draw curvilinear coordinate axes for PGPLOT.
d23 1
a23 1
* $Id: pgtest.f,v 7.3.1.3 2021/01/28 12:32:51 mcalabre Exp mcalabre $
@


7.4.1.1
log
@WCSLIB 7.4 development branch.
@
text
@d23 1
a23 1
* $Id: pgtest.f,v 7.4 2021/01/31 02:24:52 mcalabre Exp $
@


7.3
log
@WCSLIB 7.3
@
text
@d3 2
a4 2
* PGSBOX 7.3 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2020, Mark Calabretta
a20 2
* Direct correspondence concerning PGSBOX to mark@@calabretta.id.au
*
d23 1
a23 1
* $Id: pgtest.f,v 7.2.1.2 2020/05/07 04:56:37 mcalabre Exp mcalabre $
@


7.3.1.1
log
@WCSLIB 7.3 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 7.3 2020/06/03 03:37:03 mcalabre Exp $
@


7.3.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2021, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 7.3.1.1 2020/06/03 03:38:11 mcalabre Exp mcalabre $
@


7.3.1.3
log
@Removed my email address from the copyright notice.
@
text
@d21 2
d25 1
a25 1
* $Id: pgtest.f,v 7.3.1.2 2021/01/19 14:12:19 mcalabre Exp mcalabre $
@


7.2
log
@WCSLIB 7.2
@
text
@d3 1
a3 1
* PGSBOX 7.2 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 7.1.1.1 2019/12/31 13:28:19 mcalabre Exp mcalabre $
d109 1
a109 1
      STATUS = WCSPUT (WCS, WCS_FLAG, -1, 0, 0)
d334 2
a335 2
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d338 2
a339 2
      STATUS = WCSPUT (WCS, WCS_CRPIX, 384D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRPIX, 256D0, 2, 0)
d342 2
a343 2
      STATUS = WCSPUT (WCS, WCS_CDELT, -1D0/3600000D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CDELT,  1D0/3600000D0, 2, 0)
d347 2
a348 2
      STATUS = WCSPUT (WCS, WCS_CRVAL, 25D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRVAL, DEC0, 2, 0)
d352 2
a353 2
      STATUS = WCSPUT (WCS, WCS_PV, 0D0, 2, 1)
      STATUS = WCSPUT (WCS, WCS_PV, COS(DEC0)/SIN(DEC0), 2, 2)
d412 2
a413 2
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d416 2
a417 2
      STATUS = WCSPUT (WCS, WCS_CRPIX, 256D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRPIX, 256D0, 2, 0)
d420 2
a421 2
      STATUS = WCSPUT (WCS, WCS_CDELT, -1D0/3D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CDELT,  1D0/3D0, 2, 0)
d424 3
a426 3
      STATUS = WCSPUT (WCS, WCS_CRVAL, 90D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRVAL, 30D0, 2, 0)
      STATUS = WCSPUT (WCS, WCS_LONPOLE, 150D0, 0, 0)
d429 2
a430 2
      STATUS = WCSPUT (WCS, WCS_PV, 60D0, 2, 1)
      STATUS = WCSPUT (WCS, WCS_PV, 15D0, 2, 2)
d478 4
a481 4
      STATUS = WCSPUT (WCS, WCS_CRVAL,  0D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRVAL, 60D0, 2, 0)
      STATUS = WCSPUT (WCS, WCS_LONPOLE, 999D0, 0, 0)
      STATUS = WCSPUT (WCS, WCS_LATPOLE, 999D0, 0, 0)
d534 2
a535 2
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d538 2
a539 2
      STATUS = WCSPUT (WCS, WCS_CRPIX, 192D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRPIX, 640D0, 2, 0)
d543 4
a546 4
      STATUS = WCSPUT (WCS, WCS_PC,  COS(ROTN), 1, 1)
      STATUS = WCSPUT (WCS, WCS_PC,  SIN(ROTN), 1, 2)
      STATUS = WCSPUT (WCS, WCS_PC, -SIN(ROTN), 2, 1)
      STATUS = WCSPUT (WCS, WCS_PC,  COS(ROTN), 2, 2)
d549 2
a550 2
      STATUS = WCSPUT (WCS, WCS_CDELT, -1D0/5D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CDELT,  1D0/5D0, 2, 0)
d553 3
a555 3
      STATUS = WCSPUT (WCS, WCS_CRVAL, 332D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRVAL,  40D0, 2, 0)
      STATUS = WCSPUT (WCS, WCS_LONPOLE, -30D0, 0, 0)
d696 3
a698 3
      STATUS = WCSPUT (WCS, WCS_CRVAL, 0D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRVAL, 0D0, 2, 0)
      STATUS = WCSPUT (WCS, WCS_LONPOLE, 999D0, 0, 0)
d754 2
a755 2
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d758 2
a759 2
      STATUS = WCSPUT (WCS, WCS_CRPIX, 226D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRPIX,  46D0, 2, 0)
d763 4
a766 4
      STATUS = WCSPUT (WCS, WCS_PC,  COS(ROTN), 1, 1)
      STATUS = WCSPUT (WCS, WCS_PC,  SIN(ROTN), 1, 2)
      STATUS = WCSPUT (WCS, WCS_PC, -SIN(ROTN), 2, 1)
      STATUS = WCSPUT (WCS, WCS_PC,  COS(ROTN), 2, 2)
d769 2
a770 2
      STATUS = WCSPUT (WCS, WCS_CDELT, -1D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CDELT,  1D0, 2, 0)
d773 2
a774 2
      STATUS = WCSPUT (WCS, WCS_CRVAL, 0D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRVAL, 0D0, 2, 0)
d818 2
a819 2
      STATUS = WCSPUT (WCS, WCS_CRPIX, 226D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRPIX,  46D0, 2, 0)
d822 3
a824 3
      STATUS = WCSPUT (WCS, WCS_CRVAL, 30D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRVAL, 35D0, 2, 0)
      STATUS = WCSPUT (WCS, WCS_LONPOLE, 999D0, 0, 0)
d903 2
a904 2
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d907 2
a908 2
      STATUS = WCSPUT (WCS, WCS_CRPIX, 0D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRPIX, 0D0, 2, 0)
d911 2
a912 2
      STATUS = WCSPUT (WCS, WCS_CDELT, -1D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CDELT,  1D0, 2, 0)
d915 2
a916 2
      STATUS = WCSPUT (WCS, WCS_CRVAL, 0D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRVAL, 0D0, 2, 0)
d977 2
a978 2
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d981 2
a982 2
      STATUS = WCSPUT (WCS, WCS_CRPIX, 0D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRPIX, 0D0, 2, 0)
d985 2
a986 2
      STATUS = WCSPUT (WCS, WCS_CDELT, -1D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CDELT,  1D0, 2, 0)
d989 3
a991 3
      STATUS = WCSPUT (WCS, WCS_CRVAL,  45D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRVAL, -90D0, 2, 0)
      STATUS = WCSPUT (WCS, WCS_LONPOLE, 999D0, 0, 0)
d994 2
a995 2
      STATUS = WCSPUT (WCS, WCS_PV, 0D0, 2, 1)
      STATUS = WCSPUT (WCS, WCS_PV, 1D0, 2, 2)
d1056 2
a1057 2
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d1060 2
a1061 2
      STATUS = WCSPUT (WCS, WCS_CRPIX, 50.5D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRPIX,  1.0D0, 2, 0)
d1064 2
a1065 2
      STATUS = WCSPUT (WCS, WCS_CDELT, 1D-3, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CDELT, 1D-3, 2, 0)
d1068 3
a1070 3
      STATUS = WCSPUT (WCS, WCS_CRVAL, -45.0D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRVAL, -89.7D0, 2, 0)
      STATUS = WCSPUT (WCS, WCS_LONPOLE, 999D0, 0, 0)
d1287 2
a1288 2
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d1291 2
a1292 2
      STATUS = WCSPUT (WCS, WCS_CRPIX, 2D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRPIX, 2D0, 2, 0)
d1295 2
a1296 2
      STATUS = WCSPUT (WCS, WCS_CDELT, 1D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CDELT, 1D0, 2, 0)
d1299 2
a1300 2
      STATUS = WCSPUT (WCS, WCS_CRVAL, 2D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRVAL, 2D0, 2, 0)
@


7.2.1.1
log
@WCSLIB 7.2 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 7.2 2020/03/09 07:31:24 mcalabre Exp $
@


7.2.1.2
log
@Replaced calls to WCSPUT with their type-specific equivalents.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 7.2.1.1 2020/03/09 07:34:50 mcalabre Exp mcalabre $
d109 1
a109 1
      STATUS = WCSPTI (WCS, WCS_FLAG, -1, 0, 0)
d334 2
a335 2
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d338 2
a339 2
      STATUS = WCSPTD (WCS, WCS_CRPIX, 384D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX, 256D0, 2, 0)
d342 2
a343 2
      STATUS = WCSPTD (WCS, WCS_CDELT, -1D0/3600000D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT,  1D0/3600000D0, 2, 0)
d347 2
a348 2
      STATUS = WCSPTD (WCS, WCS_CRVAL, 25D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, DEC0, 2, 0)
d352 2
a353 2
      STATUS = WCSPTD (WCS, WCS_PV, 0D0, 2, 1)
      STATUS = WCSPTD (WCS, WCS_PV, COS(DEC0)/SIN(DEC0), 2, 2)
d412 2
a413 2
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d416 2
a417 2
      STATUS = WCSPTD (WCS, WCS_CRPIX, 256D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX, 256D0, 2, 0)
d420 2
a421 2
      STATUS = WCSPTD (WCS, WCS_CDELT, -1D0/3D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT,  1D0/3D0, 2, 0)
d424 3
a426 3
      STATUS = WCSPTD (WCS, WCS_CRVAL, 90D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, 30D0, 2, 0)
      STATUS = WCSPTD (WCS, WCS_LONPOLE, 150D0, 0, 0)
d429 2
a430 2
      STATUS = WCSPTD (WCS, WCS_PV, 60D0, 2, 1)
      STATUS = WCSPTD (WCS, WCS_PV, 15D0, 2, 2)
d478 4
a481 4
      STATUS = WCSPTD (WCS, WCS_CRVAL,  0D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, 60D0, 2, 0)
      STATUS = WCSPTD (WCS, WCS_LONPOLE, 999D0, 0, 0)
      STATUS = WCSPTD (WCS, WCS_LATPOLE, 999D0, 0, 0)
d534 2
a535 2
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d538 2
a539 2
      STATUS = WCSPTD (WCS, WCS_CRPIX, 192D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX, 640D0, 2, 0)
d543 4
a546 4
      STATUS = WCSPTD (WCS, WCS_PC,  COS(ROTN), 1, 1)
      STATUS = WCSPTD (WCS, WCS_PC,  SIN(ROTN), 1, 2)
      STATUS = WCSPTD (WCS, WCS_PC, -SIN(ROTN), 2, 1)
      STATUS = WCSPTD (WCS, WCS_PC,  COS(ROTN), 2, 2)
d549 2
a550 2
      STATUS = WCSPTD (WCS, WCS_CDELT, -1D0/5D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT,  1D0/5D0, 2, 0)
d553 3
a555 3
      STATUS = WCSPTD (WCS, WCS_CRVAL, 332D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL,  40D0, 2, 0)
      STATUS = WCSPTD (WCS, WCS_LONPOLE, -30D0, 0, 0)
d696 3
a698 3
      STATUS = WCSPTD (WCS, WCS_CRVAL, 0D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, 0D0, 2, 0)
      STATUS = WCSPTD (WCS, WCS_LONPOLE, 999D0, 0, 0)
d754 2
a755 2
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d758 2
a759 2
      STATUS = WCSPTD (WCS, WCS_CRPIX, 226D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX,  46D0, 2, 0)
d763 4
a766 4
      STATUS = WCSPTD (WCS, WCS_PC,  COS(ROTN), 1, 1)
      STATUS = WCSPTD (WCS, WCS_PC,  SIN(ROTN), 1, 2)
      STATUS = WCSPTD (WCS, WCS_PC, -SIN(ROTN), 2, 1)
      STATUS = WCSPTD (WCS, WCS_PC,  COS(ROTN), 2, 2)
d769 2
a770 2
      STATUS = WCSPTD (WCS, WCS_CDELT, -1D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT,  1D0, 2, 0)
d773 2
a774 2
      STATUS = WCSPTD (WCS, WCS_CRVAL, 0D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, 0D0, 2, 0)
d818 2
a819 2
      STATUS = WCSPTD (WCS, WCS_CRPIX, 226D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX,  46D0, 2, 0)
d822 3
a824 3
      STATUS = WCSPTD (WCS, WCS_CRVAL, 30D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, 35D0, 2, 0)
      STATUS = WCSPTD (WCS, WCS_LONPOLE, 999D0, 0, 0)
d903 2
a904 2
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d907 2
a908 2
      STATUS = WCSPTD (WCS, WCS_CRPIX, 0D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX, 0D0, 2, 0)
d911 2
a912 2
      STATUS = WCSPTD (WCS, WCS_CDELT, -1D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT,  1D0, 2, 0)
d915 2
a916 2
      STATUS = WCSPTD (WCS, WCS_CRVAL, 0D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, 0D0, 2, 0)
d977 2
a978 2
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d981 2
a982 2
      STATUS = WCSPTD (WCS, WCS_CRPIX, 0D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX, 0D0, 2, 0)
d985 2
a986 2
      STATUS = WCSPTD (WCS, WCS_CDELT, -1D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT,  1D0, 2, 0)
d989 3
a991 3
      STATUS = WCSPTD (WCS, WCS_CRVAL,  45D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, -90D0, 2, 0)
      STATUS = WCSPTD (WCS, WCS_LONPOLE, 999D0, 0, 0)
d994 2
a995 2
      STATUS = WCSPTD (WCS, WCS_PV, 0D0, 2, 1)
      STATUS = WCSPTD (WCS, WCS_PV, 1D0, 2, 2)
d1056 2
a1057 2
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d1060 2
a1061 2
      STATUS = WCSPTD (WCS, WCS_CRPIX, 50.5D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX,  1.0D0, 2, 0)
d1064 2
a1065 2
      STATUS = WCSPTD (WCS, WCS_CDELT, 1D-3, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT, 1D-3, 2, 0)
d1068 3
a1070 3
      STATUS = WCSPTD (WCS, WCS_CRVAL, -45.0D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, -89.7D0, 2, 0)
      STATUS = WCSPTD (WCS, WCS_LONPOLE, 999D0, 0, 0)
d1287 2
a1288 2
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPTC (WCS, WCS_CTYPE, CTYPE(2), 2, 0)
d1291 2
a1292 2
      STATUS = WCSPTD (WCS, WCS_CRPIX, 2D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRPIX, 2D0, 2, 0)
d1295 2
a1296 2
      STATUS = WCSPTD (WCS, WCS_CDELT, 1D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CDELT, 1D0, 2, 0)
d1299 2
a1300 2
      STATUS = WCSPTD (WCS, WCS_CRVAL, 2D0, 1, 0)
      STATUS = WCSPTD (WCS, WCS_CRVAL, 2D0, 2, 0)
@


7.1
log
@WCSLIB 7.1
@
text
@d3 1
a3 1
* PGSBOX 7.1 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 6.4.1.2 2019/12/31 03:29:47 mcalabre Exp mcalabre $
@


7.1.1.1
log
@WCSLIB 7.1 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 7.1 2019/12/31 13:25:20 mcalabre Exp $
@


6.4
log
@WCSLIB 6.4
@
text
@d3 2
a4 2
* PGSBOX 6.4 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2019, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 6.3.1.1 2019/07/12 07:34:37 mcalabre Exp mcalabre $
@


6.4.1.1
log
@WCSLIB 6.4 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 6.4 2019/08/15 09:30:18 mcalabre Exp $
@


6.4.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2020, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 6.4.1.1 2019/08/15 09:31:29 mcalabre Exp mcalabre $
@


6.3
log
@WCSLIB 6.3
@
text
@d3 1
a3 1
* PGSBOX 6.3 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 6.2.1.2 2019/05/10 02:13:12 mcalabre Exp mcalabre $
@


6.3.1.1
log
@WCSLIB 6.3 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 6.3 2019/07/12 07:33:40 mcalabre Exp $
@


6.2
log
@WCSLIB 6.2
@
text
@d3 2
a4 2
* PGSBOX 6.2 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2018, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 6.1.1.1 2018/10/19 08:36:15 mcalabre Exp mcalabre $
@


6.2.1.1
log
@WCSLIB 6.2 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 6.2 2018/10/20 10:03:14 mcalabre Exp $
@


6.2.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2019, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 6.2.1.1 2018/10/20 10:04:15 mcalabre Exp mcalabre $
@


6.1
log
@WCSLIB 6.1
@
text
@d3 1
a3 1
* PGSBOX 6.1 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.20.1.1 2018/10/05 05:33:33 mcalabre Exp mcalabre $
@


6.1.1.1
log
@WCSLIB 6.1 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 6.1 2018/10/19 08:35:17 mcalabre Exp $
@


5.20
log
@WCSLIB 5.20
@
text
@d3 1
a3 1
* PGSBOX 5.20 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.19.1.1 2018/07/26 15:41:42 mcalabre Exp mcalabre $
@


5.20.1.1
log
@WCSLIB 5.20 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.20 2018/10/05 05:32:42 mcalabre Exp $
@


5.19
log
@WCSLIB 5.19
@
text
@d3 1
a3 1
* PGSBOX 5.19 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.18.1.1 2018/01/10 08:33:04 mcalabre Exp mcalabre $
@


5.19.1.1
log
@WCSLIB 5.19 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.19 2018/07/26 15:37:05 mcalabre Exp $
@


5.18
log
@WCSLIB 5.18
@
text
@d3 1
a3 1
* PGSBOX 5.18 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.17.1.2 2018/01/10 08:27:33 mcalabre Exp mcalabre $
@


5.18.1.1
log
@WCSLIB 5.18 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.18 2018/01/10 08:32:14 mcalabre Exp $
@


5.17
log
@WCSLIB 5.17
@
text
@d3 2
a4 2
* PGSBOX 5.17 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2017, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 5.16.1.1 2017/01/15 04:29:04 mcalabre Exp mcalabre $
@


5.17.1.1
log
@WCSLIB 5.17 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.17 2017/09/18 08:44:25 mcalabre Exp $
@


5.17.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2018, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 5.17.1.1 2017/09/18 08:46:01 mcalabre Exp mcalabre $
@


5.16
log
@WCSLIB 5.16
@
text
@d3 1
a3 1
* PGSBOX 5.16 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.15.1.2 2017/01/15 04:08:24 mcalabre Exp mcalabre $
@


5.16.1.1
log
@WCSLIB 5.16 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.16 2017/01/15 04:25:03 mcalabre Exp $
@


5.15
log
@WCSLIB 5.15
@
text
@d3 2
a4 2
* PGSBOX 5.15 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2016, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 5.14.1.1 2016/02/07 10:52:48 mcalabre Exp mcalabre $
@


5.15.1.1
log
@WCSLIB 5.15 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.15 2016/04/05 12:55:13 mcalabre Exp $
@


5.15.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2017, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 5.15.1.1 2016/04/05 12:58:39 mcalabre Exp mcalabre $
@


5.14
log
@WCSLIB 5.14
@
text
@d3 1
a3 1
* PGSBOX 5.14 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.13.1.1 2016/01/26 08:54:36 mcalabre Exp mcalabre $
@


5.14.1.1
log
@WCSLIB 5.14 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.14 2016/02/07 10:49:33 mcalabre Exp $
@


5.13
log
@WCSLIB 5.13
@
text
@d3 1
a3 1
* PGSBOX 5.13 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.12.1.2 2016/01/25 05:22:41 mcalabre Exp mcalabre $
@


5.13.1.1
log
@WCSLIB 5.13 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.13 2016/01/26 08:51:16 mcalabre Exp $
@


5.12
log
@WCSLIB 5.12
@
text
@d3 2
a4 2
* PGSBOX 5.12 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2015, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 5.11.1.1 2015/10/18 09:16:31 mcalabre Exp mcalabre $
@


5.12.1.1
log
@WCSLIB 5.12 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.12 2015/11/15 02:13:03 mcalabre Exp $
@


5.12.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2016, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 5.12.1.1 2015/11/15 02:16:18 mcalabre Exp mcalabre $
@


5.11
log
@WCSLIB 5.11
@
text
@d3 1
a3 1
* PGSBOX 5.11 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.10.1.1 2015/10/09 08:22:40 mcalabre Exp mcalabre $
@


5.11.1.1
log
@WCSLIB 5.11 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.11 2015/10/18 09:13:08 mcalabre Exp $
@


5.10
log
@WCSLIB 5.10
@
text
@d3 1
a3 1
* PGSBOX 5.10 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.9.1.1 2015/07/21 09:28:12 mcalabre Exp mcalabre $
@


5.10.1.1
log
@WCSLIB 5.10 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.10 2015/10/09 08:19:18 mcalabre Exp $
@


5.9
log
@WCSLIB 5.9
@
text
@d3 1
a3 1
* PGSBOX 5.9 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.8.1.1 2015/07/08 11:07:22 mcalabre Exp mcalabre $
@


5.9.1.1
log
@WCSLIB 5.9 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.9 2015/07/21 09:20:04 mcalabre Exp $
@


5.8
log
@WCSLIB 5.8
@
text
@d3 1
a3 1
* PGSBOX 5.8 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.7.1.1 2015/06/29 02:47:27 mcalabre Exp mcalabre $
@


5.8.1.1
log
@WCSLIB 5.8 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.8 2015/07/08 11:04:03 mcalabre Exp $
@


5.7
log
@WCSLIB 5.7
@
text
@d3 1
a3 1
* PGSBOX 5.7 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.6.1.1 2015/06/14 07:14:48 mcalabre Exp mcalabre $
@


5.7.1.1
log
@WCSLIB 5.7 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.7 2015/06/29 02:44:19 mcalabre Exp $
@


5.6
log
@WCSLIB 5.6
@
text
@d3 1
a3 1
* PGSBOX 5.6 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.5.1.1 2015/05/05 13:21:04 mcalabre Exp mcalabre $
@


5.6.1.1
log
@WCSLIB 5.6 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.6 2015/06/14 07:11:26 mcalabre Exp $
@


5.5
log
@WCSLIB 5.5
@
text
@d3 1
a3 1
* PGSBOX 5.5 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.4.1.1 2015/04/21 14:44:32 mcalabre Exp mcalabre $
@


5.5.1.1
log
@WCSLIB 5.5 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.5 2015/05/05 13:16:34 mcalabre Exp $
@


5.4
log
@WCSLIB 5.4
@
text
@d3 1
a3 1
* PGSBOX 5.4 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.3.1.1 2015/04/21 02:54:06 mcalabre Exp mcalabre $
@


5.4.1.1
log
@WCSLIB 5.4 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.4 2015/04/21 14:41:08 mcalabre Exp $
@


5.3
log
@WCSLIB 5.3
@
text
@d3 1
a3 1
* PGSBOX 5.3 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.2.1.1 2015/04/15 12:39:18 mcalabre Exp mcalabre $
@


5.3.1.1
log
@WCSLIB 5.3 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.3 2015/04/21 02:50:54 mcalabre Exp $
@


5.2
log
@WCSLIB 5.2
@
text
@d3 1
a3 1
* PGSBOX 5.2 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.1.1.1 2015/04/07 11:01:13 mcalabre Exp mcalabre $
@


5.2.1.1
log
@WCSLIB 5.2 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.2 2015/04/15 12:35:10 mcalabre Exp $
@


5.1
log
@WCSLIB 5.1
@
text
@d3 1
a3 1
* PGSBOX 5.1 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 5.0.1.1 2015/04/05 12:28:04 mcalabre Exp mcalabre $
@


5.1.1.1
log
@WCSLIB 5.1 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.1 2015/04/07 10:56:39 mcalabre Exp $
@


5.0
log
@WCSLIB 5.0
@
text
@d3 1
a3 1
* PGSBOX 5.0 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 4.25.1.2 2015/01/06 01:03:14 mcalabre Exp mcalabre $
@


5.0.1.1
log
@WCSLIB 5.0 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 5.0 2015/04/05 12:25:03 mcalabre Exp $
@


4.25
log
@WCSLIB 4.25
@
text
@d3 2
a4 2
* PGSBOX 4.25 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2014, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 4.24.1.1 2014/09/18 15:28:12 mcalabre Exp mcalabre $
@


4.25.1.1
log
@WCSLIB 4.25 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 4.25 2014/12/14 14:29:38 mcalabre Exp $
@


4.25.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2015, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 4.25.1.1 2014/12/14 14:32:48 mcalabre Exp mcalabre $
@


4.24
log
@WCSLIB 4.24
@
text
@d3 1
a3 1
* PGSBOX 4.24 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 4.23.1.1 2014/05/13 05:53:52 mcalabre Exp mcalabre $
@


4.24.1.1
log
@WCSLIB 4.24 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 4.24 2014/09/18 15:25:02 mcalabre Exp $
@


4.23
log
@WCSLIB 4.23
@
text
@d3 1
a3 1
* PGSBOX 4.23 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 4.22.1.1 2014/04/12 15:06:34 mcalabre Exp mcalabre $
@


4.23.1.1
log
@WCSLIB 4.23 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 4.23 2014/05/13 05:50:53 mcalabre Exp $
@


4.22
log
@WCSLIB 4.22
@
text
@d3 1
a3 1
* PGSBOX 4.22 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 4.21.1.1 2014/03/24 05:15:31 mcalabre Exp mcalabre $
@


4.22.1.1
log
@WCSLIB 4.22 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 4.22 2014/04/12 15:03:54 mcalabre Exp $
@


4.21
log
@WCSLIB 4.21
@
text
@d3 1
a3 1
* PGSBOX 4.21 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 4.20.1.2 2014/03/21 06:45:26 mcalabre Exp mcalabre $
@


4.21.1.1
log
@WCSLIB 4.21 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 4.21 2014/03/24 05:12:29 mcalabre Exp $
@


4.20
log
@WCSLIB 4.20
@
text
@d3 2
a4 2
* PGSBOX 4.20 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2013, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 4.19.1.1 2013/09/29 14:20:41 mcalabre Exp mcalabre $
@


4.20.1.1
log
@WCSLIB 4.20 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 4.20 2013/12/18 05:42:51 mcalabre Exp $
@


4.20.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2014, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 4.20.1.1 2013/12/18 05:46:57 mcalabre Exp mcalabre $
@


4.19
log
@WCSLIB 4.19
@
text
@d3 1
a3 1
* PGSBOX 4.19 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 4.18.1.1 2013/07/13 10:02:40 mcalabre Exp mcalabre $
@


4.19.1.1
log
@WCSLIB 4.19 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 4.19 2013/09/29 14:17:53 mcalabre Exp $
@


4.18
log
@WCSLIB 4.18
@
text
@d3 1
a3 1
* PGSBOX 4.18 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 4.17.1.1 2013/01/29 05:32:33 cal103 Exp cal103 $
@


4.18.1.1
log
@WCSLIB 4.18 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 4.18 2013/07/13 10:00:06 mcalabre Exp $
@


4.17
log
@WCSLIB 4.17
@
text
@d3 1
a3 1
* PGSBOX 4.17 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 4.16.1.2 2013/01/29 05:11:32 cal103 Exp cal103 $
@


4.17.1.1
log
@WCSLIB 4.17 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 4.17 2013/01/29 05:29:22 cal103 Exp $
@


4.16
log
@WCSLIB 4.16
@
text
@d3 2
a4 2
* PGSBOX 4.16 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2012, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 4.15.1.1 2012/09/26 14:27:53 cal103 Exp cal103 $
@


4.16.1.1
log
@WCSLIB 4.16 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 4.16 2012/11/07 04:42:45 cal103 Exp $
@


4.16.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2013, Mark Calabretta
d25 1
a25 1
* $Id: pgtest.f,v 4.16.1.1 2012/11/07 04:45:58 cal103 Exp cal103 $
@


4.15
log
@WCSLIB 4.15
@
text
@d3 1
a3 1
* PGSBOX 4.15 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 4.14.1.1 2012/07/13 10:05:05 cal103 Exp cal103 $
@


4.15.1.1
log
@WCSLIB 4.15 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 4.15 2012/09/26 14:26:06 cal103 Exp $
@


4.14
log
@WCSLIB 4.14
@
text
@d3 1
a3 1
* PGSBOX 4.14 - draw curvilinear coordinate axes for PGPLOT.
d25 1
a25 1
* $Id: pgtest.f,v 4.13.1.2 2012/07/13 03:43:52 cal103 Exp cal103 $
@


4.14.1.1
log
@WCSLIB 4.14 development branch.
@
text
@d25 1
a25 1
* $Id: pgtest.f,v 4.14 2012/07/13 10:02:47 cal103 Exp $
@


4.13
log
@WCSLIB 4.13
@
text
@d3 1
a3 1
* PGSBOX 4.13 - draw curvilinear coordinate axes for PGPLOT.
d21 1
a21 7
* Correspondence concerning PGSBOX may be directed to:
*   Internet email: mcalabre@@atnf.csiro.au
*   Postal address: Dr. Mark Calabretta
*                   Australia Telescope National Facility, CSIRO
*                   PO Box 76
*                   Epping NSW 1710
*                   AUSTRALIA
d23 3
a25 3
* Author: Mark Calabretta, Australia Telescope National Facility
* http://www.atnf.csiro.au/~mcalabre/index.html
* $Id: pgtest.f,v 4.12.1.1 2012/02/29 07:59:11 cal103 Exp cal103 $
@


4.13.1.1
log
@WCSLIB 4.13 development branch.
@
text
@d31 1
a31 1
* $Id: pgtest.f,v 4.13 2012/03/14 07:39:31 cal103 Exp $
@


4.13.1.2
log
@Changed contact details.
@
text
@d21 7
a27 1
* Direct correspondence concerning PGSBOX to mark@@calabretta.id.au
d29 3
a31 3
* Author: Mark Calabretta, Australia Telescope National Facility, CSIRO.
* http://www.atnf.csiro.au/people/Mark.Calabretta
* $Id: pgtest.f,v 4.13.1.1 2012/03/14 07:40:39 cal103 Exp cal103 $
@


4.12
log
@WCSLIB 4.12
@
text
@d3 1
a3 1
* PGSBOX 4.12 - draw curvilinear coordinate axes for PGPLOT.
d31 1
a31 1
* $Id: pgtest.f,v 4.11.1.1 2012/02/20 23:31:14 cal103 Exp cal103 $
@


4.12.1.1
log
@WCSLIB 4.12 development branch.
@
text
@d31 1
a31 1
* $Id: pgtest.f,v 4.12 2012/02/29 07:57:57 cal103 Exp $
@


4.11
log
@WCSLIB 4.11
@
text
@d3 1
a3 1
* PGSBOX 4.11 - draw curvilinear coordinate axes for PGPLOT.
d31 1
a31 1
* $Id: pgtest.f,v 4.10.1.1 2012/02/05 23:42:47 cal103 Exp cal103 $
@


4.11.1.1
log
@WCSLIB 4.11 development branch.
@
text
@d31 1
a31 1
* $Id: pgtest.f,v 4.11 2012/02/20 23:30:03 cal103 Exp $
@


4.10
log
@WCSLIB 4.10
@
text
@d3 1
a3 1
* PGSBOX 4.10 - draw curvilinear coordinate axes for PGPLOT.
d31 1
a31 1
* $Id: pgtest.f,v 4.9.1.1 2012/01/24 05:41:46 cal103 Exp cal103 $
@


4.10.1.1
log
@WCSLIB 4.10 development branch.
@
text
@d31 1
a31 1
* $Id: pgtest.f,v 4.10 2012/02/05 23:41:45 cal103 Exp $
@


4.9
log
@WCSLIB 4.9
@
text
@d3 1
a3 1
* PGSBOX 4.9 - draw curvilinear coordinate axes for PGPLOT.
d31 1
a31 1
* $Id: pgtest.f,v 4.8.1.2 2012/01/24 05:24:50 cal103 Exp cal103 $
@


4.9.1.1
log
@WCSLIB 4.9 development branch.
@
text
@d31 1
a31 1
* $Id: pgtest.f,v 4.9 2012/01/24 05:40:30 cal103 Exp $
@


4.8
log
@WCSLIB 4.8
@
text
@d3 2
a4 2
* PGSBOX 4.8 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2011, Mark Calabretta
d31 1
a31 1
* $Id: pgtest.f,v 4.7.1.1 2011/02/07 07:04:23 cal103 Exp cal103 $
@


4.8.1.1
log
@WCSLIB 4.8 development branch.
@
text
@d31 1
a31 1
* $Id: pgtest.f,v 4.8 2011/08/15 08:05:55 cal103 Exp $
@


4.8.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2012, Mark Calabretta
d31 1
a31 1
* $Id: pgtest.f,v 4.8.1.1 2011/08/15 08:07:07 cal103 Exp cal103 $
@


4.7
log
@WCSLIB 4.7
@
text
@d3 1
a3 1
* PGSBOX 4.7 - draw curvilinear coordinate axes for PGPLOT.
d31 1
a31 1
* $Id: pgtest.f,v 4.6.1.2 2011/01/14 06:35:20 cal103 Exp cal103 $
@


4.7.1.1
log
@WCSLIB 4.7 development branch.
@
text
@d31 1
a31 1
* $Id: pgtest.f,v 4.7 2011/02/07 07:03:43 cal103 Exp $
@


4.6
log
@WCSLIB 4.6
@
text
@d3 2
a4 2
* PGSBOX 4.6 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2010, Mark Calabretta
d31 1
a31 1
* $Id: pgtest.f,v 4.5.1.1 2010/07/16 07:05:33 cal103 Exp cal103 $
@


4.6.1.1
log
@WCSLIB 4.6 development branch.
@
text
@d31 1
a31 1
* $Id: pgtest.f,v 4.6 2010/11/16 06:15:32 cal103 Exp $
@


4.6.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2011, Mark Calabretta
d31 1
a31 1
* $Id: pgtest.f,v 4.6.1.1 2010/11/16 06:16:20 cal103 Exp cal103 $
@


4.5
log
@WCSLIB 4.5
@
text
@d3 1
a3 1
* PGSBOX 4.5 - draw curvilinear coordinate axes for PGPLOT.
d31 1
a31 1
* $Id: pgtest.f,v 4.4.1.2 2010/01/19 03:45:24 cal103 Exp cal103 $
@


4.5.1.1
log
@WCSLIB 4.5 development branch.
@
text
@d31 1
a31 1
* $Id: pgtest.f,v 4.5 2010/07/16 07:01:27 cal103 Exp $
@


4.4
log
@WCSLIB 4.4
@
text
@d3 2
a4 2
* PGSBOX 4.4 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2009, Mark Calabretta
d31 1
a31 1
* $Id: pgtest.f,v 4.3.1.2 2009/07/23 04:17:17 cal103 Exp cal103 $
@


4.4.1.1
log
@WCSLIB 4.4 development branch.
@
text
@d31 1
a31 1
* $Id: pgtest.f,v 4.4 2009/08/06 07:54:40 cal103 Exp $
@


4.4.1.2
log
@Updated date in copyright notice.
@
text
@d4 1
a4 1
* Copyright (C) 1997-2010, Mark Calabretta
d31 1
a31 1
* $Id: pgtest.f,v 4.4.1.1 2009/08/10 08:55:45 cal103 Exp cal103 $
@


4.3
log
@PGSBOX 4.3
@
text
@d3 2
a4 2
*   PGSBOX 4.3 - an implementation of the FITS WCS standard.
*   Copyright (C) 1997-2007, Mark Calabretta
d6 1
a6 1
*   This file is part of PGSBOX.
d8 4
a11 4
*   PGSBOX is free software: you can redistribute it and/or modify
*   it under the terms of the GNU Lesser General Public License as
*   published by the Free Software Foundation, either version 3 of
*   the License, or (at your option) any later version.
d13 4
a16 4
*   PGSBOX is distributed in the hope that it will be useful, but
*   WITHOUT ANY WARRANTY; without even the implied warranty of
*   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
*   GNU Lesser General Public License for more details.
d18 2
a19 2
*   You should have received a copy of the GNU Lesser General Public
*   License along with PGSBOX.  If not, see http://www.gnu.org/licenses.
d21 7
a27 7
*   Correspondence concerning PGSBOX may be directed to:
*      Internet email: mcalabre@@atnf.csiro.au
*      Postal address: Dr. Mark Calabretta
*                      Australia Telescope National Facility, CSIRO
*                      PO Box 76
*                      Epping NSW 1710
*                      AUSTRALIA
d29 3
a31 3
*   Author: Mark Calabretta, Australia Telescope National Facility
*   http://www.atnf.csiro.au/~mcalabre/index.html
*   $Id: pgtest.f,v 4.2.1.4 2007/12/03 05:41:29 cal103 Exp cal103 $
d76 3
a78 3
*        Switch black and white.
         CALL PGSCR (0, 1.0, 1.0, 1.0)
         CALL PGSCR (1, 0.0, 0.0, 0.0)
d83 2
a84 2
         SCL = 1.0
         CALL PGVSTD ()
d86 2
a87 2
         SCL = 0.7
         CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
d121 12
a132 12
*        - PGSBOX uses subroutine LNGVEL.
*        - Separable (i.e. orthogonal), non-linear coordinate system.
*        - Automatic choice of coordinate increments.
*        - Extraction of a common scaling factor.
*        - Automatic choice of what edges to label.
*        - Request for tickmarks (internal) for one coordinate and grid
*          lines for the other.
*        - Simple two-colour grid using two calls with deferred
*          labelling on the first call.
*        - Degree labelling.
*        - Suppression of zero arcmin and arcsec fields in sexagesimal
*          degree format.
d170 2
a171 2
     :   0, 0D0, .FALSE., LNGVEL, 1, 1, 7, ' ', 0, NLDPRM, 256, IC,
     :   CACHE, IERR)
d180 2
a181 2
     :   GRID1, 1, GRID2, .FALSE., LNGVEL, 1, 1, 7, ' ', 0, NLDPRM, 256,
     :   IC, CACHE, IERR)
d194 14
a207 14
*        - PGSBOX uses subroutine FSCAN.
*        - Non-separable (i.e. non-orthogonal) non-linear coordinate
*          system.
*        - Automatic choice of what edges to label; results in labelling
*          the bottom, left and right sides of the plot.
*        - Cyclic labelling.  FSCAN returns the azimuth in the range
*          0 - 720 degrees but PGSBOX is set to normalize this to two
*          cycles of 0 - 360 degrees.
*        - Logarithmic labelling.
*        - Automatic choice of coordinate increments but with request
*          for all grid lines for the logarithmic coordinate.
*        - Degree labelling.
*        - Suppression of common zero arcmin and arcsec fields in
*          sexagesimal degree format.
d246 2
a247 2
     :   0D0, 0, 0D0, .FALSE., FSCAN, 1, 1, 7, ' ', 0, NLDPRM, 256, IC,
     :   CACHE, IERR)
d257 8
a264 8
*        - PGSBOX uses subroutine PGCRFN.
*        - Separable (i.e. orthogonal), non-linear coordinate system.
*        - Use of function PGCRFN for separable axis types.
*        - Automatic choice of what edges to label; results in labelling
*          the bottom and left sides of the plot.
*        - Automatic choice of coordinate increments.
*        - Logarithmic labelling over many orders of magnitude.
*        - Single-character annotation on a vertical axis is upright.
d302 2
a303 2
     :   0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM, 256,
     :   IC, CACHE, IERR)
d313 19
a331 19
*        - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*        - Non-separable (i.e. non-orthogonal) curvilinear coordinate
*          system.
*        - Demonstrate parameter definition for PGWCSL.
*        - Discovery of grid lines that do not cross any axis.
*        - Automatic choice of what edges to label; results in labelling
*          all sides of the plot.
*        - Automatic choice of coordinate increments but with request
*          for increased grid density for each coordinate.
*        - Double precision accuracy.
*        - Cyclic coordinates.  PGWCSL returns the right ascension in
*          the range -180 to +180 degrees, i.e. with a discontinuity
*          at +/- 180 degrees.
*        - Labelling of degrees as time in the range 0 - 24h.
*        - Suppression of labels that would overlap one another.
*        - Sexagesimal degree labelling with automatically determined
*          precision.
*        - Suppression of common zero minute and second fields in
*          sexagesimal time format.
d381 2
a382 2
     :   0D0, 0, 0D0, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :   NLDPRM, 256, IC, CACHE, IERR)
d392 18
a409 18
*        - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*        - Non-separable (i.e. non-orthogonal) curvilinear coordinate
*          system.
*        - Coordinate system undefined in areas of the plot.
*        - Demonstrate parameter definition for PGWCSL.
*        - Discontinuous grid lines handled by PGWCSL.
*        - Discovery of grid lines that do not cross any axis.
*        - Colour control for grid and labelling.
*        - Reduced size lettering.
*        - Automatic choice of what edges to label; results in labelling
*          all sides of the plot.
*        - Automatic choice of coordinate increments.
*        - Cyclic coordinates.  PGWCSL returns the longitude in the
*          range -180 to +180 degrees, i.e. with a discontinuity at
*          +/- 180 degrees.
*        - Suppression of labels that would overlap one another.
*        - Suppression of common zero arcmin and arcsec fields in
*          sexagesimal degree format.
d480 2
a481 2
     :   0, 0D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM, 256,
     :   IC, CACHE, IERR)
d498 2
a499 2
     :   GRID1, 2, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :   NLDPRM, 256, IC, CACHE, IERR)
d510 22
a531 22
*        - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*        - Non-separable (i.e. non-orthogonal) curvilinear coordinate
*          system.
*        - Coordinate system undefined in areas of the plot.
*        - Demonstrate parameter definition for PGWCSL.
*        - Discontinuous grid lines handled by PGWCSL.
*        - Colour coded labelling.
*        - Colour coded grid implemented by the caller.
*        - Basic management of the axis-crossing table (see code).
*        - Reduced size lettering.
*        - Tick marks external to the frame.
*        - User selection of what edges to label with request for both
*          coordinates to be labelled on bottom, left and top edges.
*        - User selection of grid lines to plot.
*        - Concatenation of annotation at bottom and left; automatically
*          suppressed at the top since only one coordinate is labelled
*          there.
*        - Suppression of labels that would overlap one another.
*        - Degree labelling.
*        - Labelling of degrees as time in the range -12 - +12h.
*        - Suppression of common zero minute and second fields in
*          sexagesimal time format.
d582 2
a583 2
     :   5D0, 0, 5D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM,
     :   256, IC, CACHE, IERR)
d598 2
a599 2
     :   GRID1, 1, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :   NLDPRM, 256, IC, CACHE, IERR)
d614 2
a615 2
     :   GRID1, 2, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :   NLDPRM, 256, IC, CACHE, IERR)
d632 2
a633 2
     :   GRID1, 4, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :   NLDPRM, 256, IC, CACHE, IERR)
d650 2
a651 2
     :   GRID1, 4, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :   NLDPRM, 256, IC, CACHE, IERR)
d693 2
a694 2
     :   GRID1, 2, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :   NLDPRM, 256, IC, CACHE, IERR)
d708 2
a709 2
     :   GRID1, 1, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :   NLDPRM, 256, IC, CACHE, IERR)
d720 15
a734 15
*        - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*        - Rectangular image.
*        - Dual coordinate grids.
*        - Non-separable (i.e. non-orthogonal) curvilinear coordinate
*          system.
*        - Demonstrate parameter definition for PGWCSL.
*        - Discontinuous grid lines handled by PGWCSL.
*        - Colour coding of grid and labelling.
*        - Reduced size lettering.
*        - Manual labelling control.
*        - Manual and automatic choice of coordinate increments.
*        - Cyclic coordinates.  PGWCSL returns the longitude in the
*          range -180 to +180 degrees, i.e. with a discontinuity at
*          +/- 180 degrees.
*        - Suppression of labels that would overlap one another.
d749 1
a749 1
         CALL PGVSTD ()
d751 1
a751 1
         CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
d820 2
a821 2
     :   15D0, 0, 15D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :   NLDPRM, 256, IC, CACHE, IERR)
d867 2
a868 2
     :   0, 0D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM, 256,
     :   IC, CACHE, IERR)
d879 8
a886 8
*        - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*        - BLC, TRC unrelated to pixel coordinates.
*        - Demonstrate parameter definition for PGWCSL.
*        - Poles and 180 meridian projected along edges of the frame.
*        - Reduced size lettering.
*        - Manual and automatic choice of coordinate increments.
*        - Suppression of common zero minute and second fields in
*          sexagesimal time format.
d898 1
a898 1
         CALL PGVSTD ()
d900 1
a900 1
         CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
d943 2
a944 2
     :   0D0, 0, 0D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM,
     :   256, IC, CACHE, IERR)
d954 7
a960 7
*        - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*        - BLC, TRC unrelated to pixel coordinates.
*        - Demonstrate parameter definition for PGWCSL.
*        - Reduced size lettering.
*        - Manual and automatic choice of coordinate increments.
*        - Suppression of common zero minute and second fields in
*          sexagesimal time format.
d972 1
a972 1
         CALL PGVSTD ()
d974 1
a974 1
         CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
d1022 2
a1023 2
     :   0D0, 0, 0D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM,
     :   256, IC, CACHE, IERR)
d1032 6
a1037 6
*        - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*        - Demonstrate parameter definition for PGWCSL.
*        - Reduced size lettering.
*        - Manual and automatic choice of coordinate increments.
*        - Suppression of common zero minute and second fields in
*          sexagesimal time format.
d1051 1
a1051 1
         CALL PGVSTD ()
d1053 1
a1053 1
         CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
d1097 2
a1098 2
     :   0D0, 0, 0D0, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :   NLDPRM, 256, IC, CACHE, IERR)
d1108 8
a1115 8
*        - PGSBOX uses subroutine PGCRFN.
*        - Separable (i.e. orthogonal), linear coordinate system.
*        - Use of function PGCRFN for separable axis types.
*        - Alternative labelling and axis annotation.
*        - Direct manipulation of the axis-crossing table.
*        - Tick mark and grid line control.
*        - User selection of what edges to label.
*        - Automatic choice of coordinate increments.
d1120 1
a1120 1
         CALL PGVSTD ()
d1122 1
a1122 1
         CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
d1174 2
a1175 2
     :   0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM, 256,
     :   IC, CACHE, IERR)
d1185 6
a1190 6
*        Look for entries associated with the right edge of the frame.
         IF (CACHE(1,J).EQ.4D0) THEN
*           Convert to feet, rounding to the nearest 0.1.
            CACHE(4,J) = CACHE(4,J) * 1D3/(25.4*12D0)
            CACHE(4,J) = AINT(CACHE(4,J)*10D0 + 0.5D0)/10D0
         END IF
d1200 2
a1201 2
     :   0D0, 0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM,
     :   256, IC, CACHE, IERR)
d1221 2
a1222 2
     :   0D0, 0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM,
     :   256, IC, CACHE, IERR)
d1232 7
a1238 7
*        - Separable (i.e. orthogonal), linear coordinate system.
*        - Use of PGLBOX for simple linear axis types.
*        - Automatic choice of what edges to label; results in labelling
*          the bottom and left sides of the plot.
*        - Automatic choice of coordinate increments.
*        - Calendar date axis labelling.
*        - Single-character annotation on a vertical axis is upright.
d1262 1
a1262 1
     :   .FALSE., 256, IC, CACHE, IERR)
d1272 6
a1277 6
*        - Separable (i.e. orthogonal), linear coordinate system.
*        - Automatic choice of what edges to label; results in labelling
*          the bottom and left sides of the plot.
*        - Automatic choice of coordinate increments.
*        - Tick marks and labels at the edges of the frame.
*        - Single-character annotation on a vertical axis is upright.
d1326 2
a1327 2
     :   0D0, 0, 0D0, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :   NLDPRM, 256, IC, CACHE, IERR)
@


4.3.1.1
log
@PGSBOX 4.3 development branch.
@
text
@d3 2
a4 2
*   PGSBOX 4.4 - an implementation of the FITS WCS standard.
*   Copyright (C) 1997-2008, Mark Calabretta
d31 1
a31 1
*   $Id: pgtest.f,v 4.3 2007/12/27 05:49:14 cal103 Exp cal103 $
@


4.3.1.2
log
@Changed indentation.
@
text
@d3 2
a4 2
* PGSBOX 4.4 - draw curvilinear coordinate axes for PGPLOT.
* Copyright (C) 1997-2009, Mark Calabretta
d6 1
a6 1
* This file is part of PGSBOX.
d8 4
a11 4
* PGSBOX is free software: you can redistribute it and/or modify it
* under the terms of the GNU Lesser General Public License as published
* by the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
d13 4
a16 4
* PGSBOX is distributed in the hope that it will be useful, but WITHOUT
* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
* FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
* License for more details.
d18 2
a19 2
* You should have received a copy of the GNU Lesser General Public
* License along with PGSBOX.  If not, see http://www.gnu.org/licenses.
d21 7
a27 7
* Correspondence concerning PGSBOX may be directed to:
*   Internet email: mcalabre@@atnf.csiro.au
*   Postal address: Dr. Mark Calabretta
*                   Australia Telescope National Facility, CSIRO
*                   PO Box 76
*                   Epping NSW 1710
*                   AUSTRALIA
d29 3
a31 3
* Author: Mark Calabretta, Australia Telescope National Facility
* http://www.atnf.csiro.au/~mcalabre/index.html
* $Id: pgtest.f,v 4.3.1.1 2007/12/27 13:26:31 cal103 Exp cal103 $
d76 3
a78 3
*       Switch black and white.
        CALL PGSCR (0, 1.0, 1.0, 1.0)
        CALL PGSCR (1, 0.0, 0.0, 0.0)
d83 2
a84 2
        SCL = 1.0
        CALL PGVSTD ()
d86 2
a87 2
        SCL = 0.7
        CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
d121 12
a132 12
*       - PGSBOX uses subroutine LNGVEL.
*       - Separable (i.e. orthogonal), non-linear coordinate system.
*       - Automatic choice of coordinate increments.
*       - Extraction of a common scaling factor.
*       - Automatic choice of what edges to label.
*       - Request for tickmarks (internal) for one coordinate and grid
*         lines for the other.
*       - Simple two-colour grid using two calls with deferred
*         labelling on the first call.
*       - Degree labelling.
*       - Suppression of zero arcmin and arcsec fields in sexagesimal
*         degree format.
d170 2
a171 2
     :  0, 0D0, .FALSE., LNGVEL, 1, 1, 7, ' ', 0, NLDPRM, 256, IC,
     :  CACHE, IERR)
d180 2
a181 2
     :  GRID1, 1, GRID2, .FALSE., LNGVEL, 1, 1, 7, ' ', 0, NLDPRM, 256,
     :  IC, CACHE, IERR)
d194 14
a207 14
*       - PGSBOX uses subroutine FSCAN.
*       - Non-separable (i.e. non-orthogonal) non-linear coordinate
*         system.
*       - Automatic choice of what edges to label; results in labelling
*         the bottom, left and right sides of the plot.
*       - Cyclic labelling.  FSCAN returns the azimuth in the range
*         0 - 720 degrees but PGSBOX is set to normalize this to two
*         cycles of 0 - 360 degrees.
*       - Logarithmic labelling.
*       - Automatic choice of coordinate increments but with request
*         for all grid lines for the logarithmic coordinate.
*       - Degree labelling.
*       - Suppression of common zero arcmin and arcsec fields in
*         sexagesimal degree format.
d246 2
a247 2
     :  0D0, 0, 0D0, .FALSE., FSCAN, 1, 1, 7, ' ', 0, NLDPRM, 256, IC,
     :  CACHE, IERR)
d257 8
a264 8
*       - PGSBOX uses subroutine PGCRFN.
*       - Separable (i.e. orthogonal), non-linear coordinate system.
*       - Use of function PGCRFN for separable axis types.
*       - Automatic choice of what edges to label; results in labelling
*         the bottom and left sides of the plot.
*       - Automatic choice of coordinate increments.
*       - Logarithmic labelling over many orders of magnitude.
*       - Single-character annotation on a vertical axis is upright.
d302 2
a303 2
     :  0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM, 256,
     :  IC, CACHE, IERR)
d313 19
a331 19
*       - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*       - Non-separable (i.e. non-orthogonal) curvilinear coordinate
*         system.
*       - Demonstrate parameter definition for PGWCSL.
*       - Discovery of grid lines that do not cross any axis.
*       - Automatic choice of what edges to label; results in labelling
*         all sides of the plot.
*       - Automatic choice of coordinate increments but with request
*         for increased grid density for each coordinate.
*       - Double precision accuracy.
*       - Cyclic coordinates.  PGWCSL returns the right ascension in
*         the range -180 to +180 degrees, i.e. with a discontinuity
*         at +/- 180 degrees.
*       - Labelling of degrees as time in the range 0 - 24h.
*       - Suppression of labels that would overlap one another.
*       - Sexagesimal degree labelling with automatically determined
*         precision.
*       - Suppression of common zero minute and second fields in
*         sexagesimal time format.
d381 2
a382 2
     :  0D0, 0, 0D0, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
d392 18
a409 18
*       - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*       - Non-separable (i.e. non-orthogonal) curvilinear coordinate
*         system.
*       - Coordinate system undefined in areas of the plot.
*       - Demonstrate parameter definition for PGWCSL.
*       - Discontinuous grid lines handled by PGWCSL.
*       - Discovery of grid lines that do not cross any axis.
*       - Colour control for grid and labelling.
*       - Reduced size lettering.
*       - Automatic choice of what edges to label; results in labelling
*         all sides of the plot.
*       - Automatic choice of coordinate increments.
*       - Cyclic coordinates.  PGWCSL returns the longitude in the
*         range -180 to +180 degrees, i.e. with a discontinuity at
*         +/- 180 degrees.
*       - Suppression of labels that would overlap one another.
*       - Suppression of common zero arcmin and arcsec fields in
*         sexagesimal degree format.
d480 2
a481 2
     :  0, 0D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM, 256,
     :  IC, CACHE, IERR)
d498 2
a499 2
     :  GRID1, 2, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
d510 22
a531 22
*       - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*       - Non-separable (i.e. non-orthogonal) curvilinear coordinate
*         system.
*       - Coordinate system undefined in areas of the plot.
*       - Demonstrate parameter definition for PGWCSL.
*       - Discontinuous grid lines handled by PGWCSL.
*       - Colour coded labelling.
*       - Colour coded grid implemented by the caller.
*       - Basic management of the axis-crossing table (see code).
*       - Reduced size lettering.
*       - Tick marks external to the frame.
*       - User selection of what edges to label with request for both
*         coordinates to be labelled on bottom, left and top edges.
*       - User selection of grid lines to plot.
*       - Concatenation of annotation at bottom and left; automatically
*         suppressed at the top since only one coordinate is labelled
*         there.
*       - Suppression of labels that would overlap one another.
*       - Degree labelling.
*       - Labelling of degrees as time in the range -12 - +12h.
*       - Suppression of common zero minute and second fields in
*         sexagesimal time format.
d582 2
a583 2
     :  5D0, 0, 5D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM,
     :  256, IC, CACHE, IERR)
d598 2
a599 2
     :  GRID1, 1, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
d614 2
a615 2
     :  GRID1, 2, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
d632 2
a633 2
     :  GRID1, 4, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
d650 2
a651 2
     :  GRID1, 4, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
d693 2
a694 2
     :  GRID1, 2, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
d708 2
a709 2
     :  GRID1, 1, GRID2, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
d720 15
a734 15
*       - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*       - Rectangular image.
*       - Dual coordinate grids.
*       - Non-separable (i.e. non-orthogonal) curvilinear coordinate
*         system.
*       - Demonstrate parameter definition for PGWCSL.
*       - Discontinuous grid lines handled by PGWCSL.
*       - Colour coding of grid and labelling.
*       - Reduced size lettering.
*       - Manual labelling control.
*       - Manual and automatic choice of coordinate increments.
*       - Cyclic coordinates.  PGWCSL returns the longitude in the
*         range -180 to +180 degrees, i.e. with a discontinuity at
*         +/- 180 degrees.
*       - Suppression of labels that would overlap one another.
d749 1
a749 1
        CALL PGVSTD ()
d751 1
a751 1
        CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
d820 2
a821 2
     :  15D0, 0, 15D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
d867 2
a868 2
     :  0, 0D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM, 256,
     :  IC, CACHE, IERR)
d879 8
a886 8
*       - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*       - BLC, TRC unrelated to pixel coordinates.
*       - Demonstrate parameter definition for PGWCSL.
*       - Poles and 180 meridian projected along edges of the frame.
*       - Reduced size lettering.
*       - Manual and automatic choice of coordinate increments.
*       - Suppression of common zero minute and second fields in
*         sexagesimal time format.
d898 1
a898 1
        CALL PGVSTD ()
d900 1
a900 1
        CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
d943 2
a944 2
     :  0D0, 0, 0D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM,
     :  256, IC, CACHE, IERR)
d954 7
a960 7
*       - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*       - BLC, TRC unrelated to pixel coordinates.
*       - Demonstrate parameter definition for PGWCSL.
*       - Reduced size lettering.
*       - Manual and automatic choice of coordinate increments.
*       - Suppression of common zero minute and second fields in
*         sexagesimal time format.
d972 1
a972 1
        CALL PGVSTD ()
d974 1
a974 1
        CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
d1022 2
a1023 2
     :  0D0, 0, 0D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM,
     :  256, IC, CACHE, IERR)
d1032 6
a1037 6
*       - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*       - Demonstrate parameter definition for PGWCSL.
*       - Reduced size lettering.
*       - Manual and automatic choice of coordinate increments.
*       - Suppression of common zero minute and second fields in
*         sexagesimal time format.
d1051 1
a1051 1
        CALL PGVSTD ()
d1053 1
a1053 1
        CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
d1097 2
a1098 2
     :  0D0, 0, 0D0, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
d1108 8
a1115 8
*       - PGSBOX uses subroutine PGCRFN.
*       - Separable (i.e. orthogonal), linear coordinate system.
*       - Use of function PGCRFN for separable axis types.
*       - Alternative labelling and axis annotation.
*       - Direct manipulation of the axis-crossing table.
*       - Tick mark and grid line control.
*       - User selection of what edges to label.
*       - Automatic choice of coordinate increments.
d1120 1
a1120 1
        CALL PGVSTD ()
d1122 1
a1122 1
        CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
d1174 2
a1175 2
     :  0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM, 256,
     :  IC, CACHE, IERR)
d1185 6
a1190 6
*       Look for entries associated with the right edge of the frame.
        IF (CACHE(1,J).EQ.4D0) THEN
*         Convert to feet, rounding to the nearest 0.1.
          CACHE(4,J) = CACHE(4,J) * 1D3/(25.4*12D0)
          CACHE(4,J) = AINT(CACHE(4,J)*10D0 + 0.5D0)/10D0
        END IF
d1200 2
a1201 2
     :  0D0, 0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM,
     :  256, IC, CACHE, IERR)
d1221 2
a1222 2
     :  0D0, 0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM,
     :  256, IC, CACHE, IERR)
d1232 7
a1238 7
*       - Separable (i.e. orthogonal), linear coordinate system.
*       - Use of PGLBOX for simple linear axis types.
*       - Automatic choice of what edges to label; results in labelling
*         the bottom and left sides of the plot.
*       - Automatic choice of coordinate increments.
*       - Calendar date axis labelling.
*       - Single-character annotation on a vertical axis is upright.
d1262 1
a1262 1
     :  .FALSE., 256, IC, CACHE, IERR)
d1272 6
a1277 6
*       - Separable (i.e. orthogonal), linear coordinate system.
*       - Automatic choice of what edges to label; results in labelling
*         the bottom and left sides of the plot.
*       - Automatic choice of coordinate increments.
*       - Tick marks and labels at the edges of the frame.
*       - Single-character annotation on a vertical axis is upright.
d1326 2
a1327 2
     :  0D0, 0, 0D0, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :  NLDPRM, 256, IC, CACHE, IERR)
@


4.2
log
@PGSBOX 4.2
@
text
@d3 2
a4 2
*   PGSBOX 4.2 - a non-linear coordinate axis plotter for PGPLOT.
*   Copyright (C) 1997-2005, Mark Calabretta
d6 1
a6 4
*   PGSBOX is free software; you can redistribute it and/or modify it under
*   the terms of the GNU General Public License as published by the Free
*   Software Foundation; either version 2 of the License, or (at your option)
*   any later version.
d8 4
a11 4
*   PGSBOX is distributed in the hope that it will be useful, but WITHOUT ANY
*   WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
*   FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
*   details.
d13 7
a19 3
*   You should have received a copy of the GNU General Public License along
*   with PGSBOX; if not, write to the Free Software Foundation, Inc.,
*   59 Temple Place, Suite 330, Boston, MA  02111-1307, USA
d31 1
a31 1
*   $Id: pgtest.f,v 4.1 2005/08/26 12:54:06 cal103 Exp cal103 $
a58 6
*     Uncomment the following two lines to raise SIGFPE on floating
*     point exceptions for the Sun FORTRAN compiler.  This signal can
*     be caught within 'dbx' by issuing the command "catch FPE".
*#include <f77_floatingpoint.h>
*      call ieee_handler ('set', 'common', SIGFPE_ABORT)
*-----------------------------------------------------------------------
d953 154
@


4.2.1.1
log
@PGSBOX 4.2 development branch.
@
text
@d3 1
a3 1
*   PGSBOX 4.3 - a non-linear coordinate axis plotter for PGPLOT.
d30 1
a30 1
*   $Id: pgtest.f,v 4.2 2005/09/21 13:27:15 cal103 Exp cal103 $
@


4.2.1.2
log
@Added two new tests, oblique CYP and TAN graticules arising from checking
example headers in William Thompson's solar coordinates paper.
@
text
@d30 1
a30 1
*   $Id: pgtest.f,v 4.2.1.1 2005/12/05 03:23:34 cal103 Exp cal103 $
a949 155

*     Draw the frame.
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL PGPAGE ()

*-----------------------------------------------------------------------

*     Cylindrical perspective projection.
*        - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*        - BLC, TRC unrelated to pixel coordinates.
*        - Demonstrate parameter definition for PGWCSL.
*        - Reduced size lettering.
*        - Manual and automatic choice of coordinate increments.
*        - Suppression of common zero minute and second fields in
*          sexagesimal time format.

      WRITE (*, '(/,A)') 'Cylindrical perspective projection'

      BLC(1) = -180.0
      BLC(2) =  -90.0
      TRC(1) =  180.0
      TRC(2) =  +90.0

*     Reset viewport for rectangular image.

      IF (LARGE) THEN
         CALL PGVSTD ()
      ELSE
         CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
      END IF
      CALL PGWNAD (BLC(1), TRC(1), BLC(2), TRC(2))

      STATUS = WCSINI (2, WCS)

*     Set projection type to cylindrical perspective.
      CTYPE(1) = 'RA---CYP'
      CTYPE(2) = 'DEC--CYP'
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(2), 2, 0)

*     Reference pixel coordinates.
      STATUS = WCSPUT (WCS, WCS_CRPIX, 0D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRPIX, 0D0, 2, 0)

*     Coordinate increments.
      STATUS = WCSPUT (WCS, WCS_CDELT, -1D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CDELT,  1D0, 2, 0)

*     Set parameters to draw the native grid.
      STATUS = WCSPUT (WCS, WCS_CRVAL,  45D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRVAL, -90D0, 2, 0)
      STATUS = WCSPUT (WCS, WCS_LONPOLE, 999D0, 0, 0)

*     mu and lambda projection parameters.
      STATUS = WCSPUT (WCS, WCS_PV, 0D0, 2, 1)
      STATUS = WCSPUT (WCS, WCS_PV, 1D0, 2, 2)

*     Annotation.
      IDENTS(1) = 'Right ascension'
      IDENTS(2) = 'Declination'
      IDENTS(3) = 'WCS cylindrical perspective projection'

      OPT(1) = 'G'
      OPT(2) = 'E'

*     Reduced size lettering.
      CALL PGSCH (0.7*SCL)

*     Draw full grid lines.
      GCODE(1) = 2
      GCODE(2) = 2

      CALL PGSCI (1)

      IC = -1
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 2121, 1212, C0, GCODE, 0D0, 0,
     :   0D0, 0, 0D0, .TRUE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS, NLDPRM,
     :   256, IC, CACHE, IERR)

*     Draw the frame.
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL PGPAGE ()

*-----------------------------------------------------------------------
*     Gnomonic projection.
*        - PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*        - Demonstrate parameter definition for PGWCSL.
*        - Reduced size lettering.
*        - Manual and automatic choice of coordinate increments.
*        - Suppression of common zero minute and second fields in
*          sexagesimal time format.

      WRITE (*, '(/,A)') 'TAN projection'

      NAXIS(1) = 100
      NAXIS(2) = 100

      BLC(1) = 0.5
      BLC(2) = 0.5
      TRC(1) = NAXIS(1) + 0.5
      TRC(2) = NAXIS(2) + 0.5

*     Reset viewport for rectangular image.

      IF (LARGE) THEN
         CALL PGVSTD ()
      ELSE
         CALL PGVSIZ (1.0, 3.0, 1.0, 3.0)
      END IF
      CALL PGWNAD (0.0, 1.0, 0.0, REAL(NAXIS(2))/REAL(NAXIS(1)))

      STATUS = WCSINI (2, WCS)

*     Set projection type to gnomonic.
      CTYPE(1) = 'RA---TAN'
      CTYPE(2) = 'DEC--TAN'
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(1), 1, 0)
      STATUS = WCSPUT (WCS, WCS_CTYPE, CTYPE(2), 2, 0)

*     Reference pixel coordinates.
      STATUS = WCSPUT (WCS, WCS_CRPIX, 50.5D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRPIX,  1.0D0, 2, 0)

*     Coordinate increments.
      STATUS = WCSPUT (WCS, WCS_CDELT, 1D-3, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CDELT, 1D-3, 2, 0)

*     Set parameters to draw the native grid.
      STATUS = WCSPUT (WCS, WCS_CRVAL, -45.0D0, 1, 0)
      STATUS = WCSPUT (WCS, WCS_CRVAL, -89.7D0, 2, 0)
      STATUS = WCSPUT (WCS, WCS_LONPOLE, 999D0, 0, 0)

*     Annotation.
      IDENTS(1) = 'Right ascension'
      IDENTS(2) = 'Declination'
      IDENTS(3) = 'WCS TAN projection'

      OPT(1) = 'E'
      OPT(2) = 'E'

*     Reduced size lettering.
      CALL PGSCH (0.7*SCL)

*     Draw full grid lines.
      GCODE(1) = 2
      GCODE(2) = 2

      CALL PGSCI (1)

      IC = -1
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 1212, C0, GCODE, 0D0, 0,
     :   0D0, 0, 0D0, .FALSE., PGWCSL, 1, WCSLEN, 1, NLCPRM, WCS,
     :   NLDPRM, 256, IC, CACHE, IERR)
@


4.2.1.3
log
@Removed comment re catching floating point exceptions in SunOS.
@
text
@d30 1
a30 1
*   $Id: pgtest.f,v 4.2.1.2 2005/12/05 03:32:04 cal103 Exp cal103 $
d58 6
d1055 1
@


4.2.1.4
log
@Switched to LGPL 3.0.
@
text
@d3 2
a4 2
*   PGSBOX 4.3 - an implementation of the FITS WCS standard.
*   Copyright (C) 1997-2007, Mark Calabretta
d6 4
a9 1
*   This file is part of PGSBOX.
d11 4
a14 4
*   PGSBOX is free software: you can redistribute it and/or modify
*   it under the terms of the GNU Lesser General Public License as
*   published by the Free Software Foundation, either version 3 of
*   the License, or (at your option) any later version.
d16 3
a18 7
*   PGSBOX is distributed in the hope that it will be useful, but
*   WITHOUT ANY WARRANTY; without even the implied warranty of
*   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
*   GNU Lesser General Public License for more details.
*
*   You should have received a copy of the GNU Lesser General Public
*   License along with PGSBOX.  If not, see http://www.gnu.org/licenses.
d30 1
a30 1
*   $Id: pgtest.f,v 4.2.1.3 2005/12/16 03:42:59 cal103 Exp cal103 $
@


4.1
log
@PGSBOX 4.1: in the COE example, LATPOLE must be reset when drawing the
second (native) grid because it will have been set to a non-default value
by WCSSET when the first grid was drawn; set wcs.flag to -1 before WCSINIT
and call WCSFREE at the end.
@
text
@d3 1
a3 1
*   PGSBOX 4.1 - a non-linear coordinate axis plotter for PGPLOT.
d30 1
a30 1
*   $Id: pgtest.f,v 4.0 2005/02/04 12:22:47 cal103 Exp cal103 $
@


4.0
log
@PGSBOX 4.0: changed from LGPL to GPL.
@
text
@d3 1
a3 1
*   PGSBOX 4.0 - a non-linear coordinate axis plotter for PGPLOT.
d30 1
a30 1
*   $Id: pgtest.f,v 3.6 2004/08/24 08:42:14 cal103 Exp cal103 $
d492 1
d1187 2
a1190 1

@


3.6
log
@PGSBOX 3.6: Changed "WRITE (6, ...)" to "WRITE (*, ...)" as being standard
FORTRAN and notionally more portable.
@
text
@d3 2
a4 2
*   PGSBOX 3.6 - a non-linear coordinate axis plotter for PGPLOT.
*   Copyright (C) 1997-2004, Mark Calabretta
d6 4
a9 4
*   This program is free software; you can redistribute it and/or modify
*   it under the terms of the GNU General Public License as published by
*   the Free Software Foundation; either version 2 of the License, or
*   (at your option) any later version.
d11 4
a14 4
*   This program is distributed in the hope that it will be useful, but
*   WITHOUT ANY WARRANTY; without even the implied warranty of
*   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
*   General Public License for more details.
d16 3
a18 3
*   You should have received a copy of the GNU General Public License
*   along with this library; if not, write to the Free Software
*   Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
d30 1
a30 1
*   $Id: pgtest.f,v 3.5 2004/06/28 05:37:31 mcalabre Exp mcalabre $
d197 1
a197 1
*     an antenna which scans in azimuth with a receiver which scans
d322 1
a322 1
*        - Discovery of grid lines which do not cross any axis.
d332 1
a332 1
*        - Suppression of labels which would overlap one another.
d403 1
a403 1
*        - Discovery of grid lines which do not cross any axis.
d412 1
a412 1
*        - Suppression of labels which would overlap one another.
d527 1
a527 1
*        - User selection of which grid lines to plot.
d531 1
a531 1
*        - Suppression of labels which would overlap one another.
d738 1
a738 1
*        - Suppression of labels which would overlap one another.
@


3.5
log
@PGSBOX 3.5: print the name of each example as it is run; added titles to the
calendar axes and the PGWCSL linear axes plots.
@
text
@d3 1
a3 1
*   PGSBOX 3.5 - a non-linear coordinate axis plotter for PGPLOT.
d30 1
a30 1
*   $Id: pgtest.f,v 3.4 2004/02/10 23:56:37 mcalabre Exp mcalabre $
d139 1
a139 1
      WRITE (6, '(/,A)') 'Longitude-velocity map'
d214 1
a214 1
      WRITE (6, '(/,A)') 'Azimuth-frequency scan'
d271 1
a271 1
      WRITE (6, '(/,A)') 'Z versus time plot'
d338 1
a338 1
      WRITE (6, '(/,A)') 'Simple SIN projection'
d416 1
a416 1
      WRITE (6, '(/,A)') 'Conic equal area projection'
d537 1
a537 1
      WRITE (6, '(/,A)') 'Polyconic projection with colour-coded grid'
d740 1
a740 1
      WRITE (6, '(/,A)') 'Plate Carree projection'
d892 1
a892 1
      WRITE (6, '(/,A)') 'Plate Carree projection'
d967 1
a967 1
      WRITE (6, '(/,A)') 'Linear plot with alternative labelling'
d1090 1
a1090 1
      WRITE (6, '(/,A)') 'Calendar axes using subroutine PGLBOX'
d1129 1
a1129 1
      WRITE (6, '(/,A)') 'Simple linear axes handled by pgwcsl()'
@


3.4
log
@PGSBOX 3.4
@
text
@d3 1
a3 1
*   PGSBOX 3.4 - a non-linear coordinate axis plotter for PGPLOT.
d29 2
a30 1
*   $Id: pgtest.f,v 3.3 2003/10/22 08:55:52 mcalabre Exp mcalabre $
d39 1
a39 1
     :           NLDPRM(8), ROTN, TIKLEN
d139 2
d214 2
d271 2
d338 2
d416 2
d537 2
d740 2
d892 2
d967 2
d1090 2
d1097 1
a1097 1
      IDENTS(3) = ' '
d1121 1
a1121 2
*     Simple linear axes.
*        - Simple linear axes handled by PGWCSL.
d1129 2
d1161 1
a1161 1
      IDENTS(3) = ' '
@


3.3
log
@PGSBOX 3.3: fixed a bug in setting CRVAL1 for the SIN example; added two
further examples - a plate carree projection with the poles and +/- 180 deg
meridians coincident with the frame, and use of pgwcsl_() for a simple
linear plot.
@
text
@d3 2
a4 2
*   PGSBOX 3.3 - a non-linear coordinate axis plotter for PGPLOT.
*   Copyright (C) 1997-2003, Mark Calabretta
d29 1
a29 1
*   $Id: pgtest.f,v 3.2 2003/09/09 08:04:14 mcalabre Exp mcalabre $
@


3.2
log
@PGSBOX 3.2: set WCS_FLAG = -1 to initialize memory management in WCSINI;
changed the ninth example to use PGLBOX.
@
text
@d3 1
a3 1
*   PGSBOX 3.2 - a non-linear coordinate axis plotter for PGPLOT.
d29 1
a29 1
*   $Id: pgtest.f,v 3.1 2003/04/28 08:53:24 mcalabre Exp mcalabre $
d119 2
d123 3
a125 3
*     First example: longitude-velocity map.  The y-axis is regularly
*     spaced in frequency but is to be labelled as a true relativistic
*     velocity.  PGSBOX uses subroutine LNGVEL.  Features:
d193 4
a196 4
*     Second example: azimuth-frequency scan.  This sort of output
*     might be obtained from an antenna which scans in azimuth with a
*     receiver which scans simultaneously in frequency.  PGSBOX uses
*     subroutine FSCAN.  Features:
d256 2
a257 2
*     Third example: Z versus time plot.  PGSBOX uses subroutine PGCRFN.
*     Features:
d310 2
a311 3
*     Fourth example: A simple SIN projection near the south celestial
*     pole.  PGSBOX uses subroutine PGWCSL to interface to WCSLIB.
*     Features:
a330 1
      STATUS = WCSPUT (WCS, WCS_FLAG, -1, 0, 0)
d349 1
a349 1
      STATUS = WCSPUT (WCS, WCS_CRVAL, 25D0 1, 0)
d387 2
a388 2
*     Fifth example: Conic equal area.  PGSBOX uses subroutine PGWCSL
*     to interface to WCSLIB.  Features:
d502 2
a503 2
*     Sixth example: Polyconic with colour-coded grid.  PGSBOX uses
*     subroutine PGWCSL to interface to WCSLIB.  Features:
d710 2
a711 2
*     Seventh example: Plate Carree projection.  PGSBOX uses subroutine
*     PGWCSL to interface to WCSLIB.  Features:
d867 75
a941 2
*     Eighth example: linear-linear plot with two types of alternative
*     labelling.  PGSBOX uses subroutine PGCRFN.  Features:
d990 2
a991 2
*     Normal size lettering.
      CALL PGSCH (1.0*SCL)
d1062 1
a1062 2
*     Ninth example: Calendar axes using subroutine PGLBOX.
*     Features:
d1096 66
d1164 1
@


3.1
log
@PGSBOX 3.1
Use the WCSLIB 3.x FORTRAN wrappers to set up the wcsprm data structure for
PGWCSL; eliminated use of trigd functions; added example 9 for calendar date
axes; amended copyright notice.
@
text
@d3 1
a3 1
*   PGSBOX 3.1 - a non-linear coordinate axis plotter for PGPLOT.
d29 1
a29 1
*   $Id: pgtest.f,v 3.1 2003/04/28 05:24:34 mcalabre Exp $
d330 1
d989 1
a989 1
*     Ninth example: Calendar axes.  PGSBOX uses subroutine PGCRFN.
d992 1
a992 1
*        - Use of function PGCRFN for separable axis types.
d999 1
a999 15
*     Function types.
      FCODE(1) = 'Lin '
      FCODE(2) = 'Lin '

*     Reference pixel coordinates.
      NLDPRM(1) = 0.5D0
      NLDPRM(2) = 0.5D0

*     Coordinate increments.
      NLDPRM(3) = 1.0D0
      NLDPRM(4) = 10.0D0

*     Reference pixel values.
      NLDPRM(5) = 51900D0
      NLDPRM(6) = 51900D0
d1018 2
a1019 3
      CALL PGSBOX (BLC, TRC, IDENTS, OPT, 0, 0, C0, GCODE, 2D0, 0, 0D0,
     :   0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM, 256,
     :   IC, CACHE, IERR)
@


1.22
log
@Use CYLFIX to repair the header of example 7; use EQUIVALENCEs for CEL and
PRJ in constructing NLDPRM; renamed AXEN -> NAXIS, LNGPOL -> LONPOL.
@
text
@d3 2
a4 2
*   PGSBOX - a non-linear coordinate axis plotter for PGPLOT.
*   Copyright (C) 1997-2001, Mark Calabretta
d6 4
a9 4
*   This library is free software; you can redistribute it and/or
*   modify it under the terms of the GNU Library General Public
*   License as published by the Free Software Foundation; either
*   version 2 of the License, or (at your option) any later version.
d11 2
a12 2
*   This library is distributed in the hope that it will be useful,
*   but WITHOUT ANY WARRANTY; without even the implied warranty of
d14 1
a14 1
*   Library General Public License for more details.
d16 3
a18 3
*   You should have received a copy of the GNU Library General Public
*   License along with this library; if not, write to the Free
*   Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
d22 4
a25 4
*      Postal address: Dr. Mark Calabretta,
*                      Australia Telescope National Facility,
*                      P.O. Box 76,
*                      Epping, NSW, 2121,
d29 1
a29 1
*   $Id: pgtest.f,v 1.21 2002/04/03 06:21:31 mcalabre Exp mcalabre $
d35 1
a35 1
     :          NLIPRM(2), WCS(4)
d37 8
a44 21
      DOUBLE PRECISION CDELT(2), CACHE(4,0:256), CEL(10), CRPIX(2),
     :          CRVAL(2), GRID1(0:10), GRID2(0:10), LATPOL, LIN(19),
     :          LONPOL, NLDPRM(50), PC(2,2), PRJ(0:20), PROJP(0:9),
     :          ROTN, TIKLEN
      CHARACTER CTYPE(2)*8, DEVTYP*16, ESCAPE*1, FCODE(2)*4,
     :          IDENTS(3)*80, OPT(2)*1

*     From the prologue of PGWCSL.
      EQUIVALENCE (LIN(4),  CRPIX)
      EQUIVALENCE (LIN(6),  PC)
      EQUIVALENCE (LIN(10), CDELT)

      EQUIVALENCE (CEL(1), CRVAL)
      EQUIVALENCE (CEL(3), LONPOL)
      EQUIVALENCE (CEL(4), LATPOL)

      EQUIVALENCE (PRJ(0), PROJP)

      EQUIVALENCE (NLDPRM(1),  LIN)
      EQUIVALENCE (NLDPRM(20), CEL)
      EQUIVALENCE (NLDPRM(30), PRJ)
d46 3
a55 2

      DOUBLE PRECISION COSD, SIND
d137 1
a137 1
      NLDPRM(1) = 1D0
d162 1
a162 2
*     Do tick marks for longitude and grid lines for velocity in yellow.
*     Defer labelling.
d166 2
d232 3
d236 1
a239 3
*     Normal size lettering.
      CALL PGSCH (1.0*SCL)

a241 1
      CALL PGSCI (1)
d283 1
a283 1
      IDENTS(3) = ''
d288 3
d292 1
a295 4
*     Normal size lettering.
      CALL PGSCH (1.0*SCL)

      CALL PGSCI (1)
d313 1
a313 1
*        - Use of equivalences to define parameters for PGWCSL.
d330 2
d335 6
d342 13
a354 23
*     Reference pixel coordinates (EQUIVALENCE'd into NLDPRM above).
      CRPIX(1) = 384D0
      CRPIX(2) = 256D0

*     Linear transformation matrix (EQUIVALENCE'd into NLDPRM).
      PC(1,1) = 1D0
      PC(1,2) = 0D0
      PC(2,1) = 0D0
      PC(2,2) = 1D0

*     Coordinate increments (EQUIVALENCE'd into NLDPRM).
      CDELT(1) = -1D0/3600000D0
      CDELT(2) =  1D0/3600000D0

*     Spherical coordinate references (EQUIVALENCE'd into NLDPRM).
      CRVAL(1) =  25.000D0
      CRVAL(2) = -89.99995D0
      LONPOL = 999D0
      LATPOL = 999D0

*     Set parameters for an NCP projection (EQUIVALENCE'd into NLDPRM).
      PROJP(1)  = 0D0
      PROJP(2)  = COSD(CRVAL(2))/SIND(CRVAL(2))
d364 3
d368 1
a371 3
*     Compact lettering.
      CALL PGSCH (0.8*SCL)

a373 1
      CALL PGSCI (1)
d376 2
a377 2
     :   0D0, 0, 0D0, .FALSE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM,
     :   256, IC, CACHE, IERR)
d391 1
a391 1
*        - Use of equivalences to define parameters for PGWCSL.
d406 2
d411 6
d418 8
a425 13
*     Reference pixel coordinates (EQUIVALENCE'd into NLDPRM above).
      CRPIX(1) = 256D0
      CRPIX(2) = 256D0

*     Linear transformation matrix (EQUIVALENCE'd into NLDPRM).
      PC(1,1) = 1D0
      PC(1,2) = 0D0
      PC(2,1) = 0D0
      PC(2,2) = 1D0

*     Coordinate increments (EQUIVALENCE'd into NLDPRM).
      CDELT(1) = -1D0/3D0
      CDELT(2) =  1D0/3D0
d428 2
a429 8
      PROJP(1)  = 60D0
      PROJP(2)  = 15D0

*     Spherical coordinate references (EQUIVALENCE'd into NLDPRM).
      CRVAL(1) = 90D0
      CRVAL(2) = 30D0
      LONPOL = 150D0
      LATPOL = 999D0
d439 3
a445 3
*     Reduced size lettering.
      CALL PGSCH (0.8*SCL)

d461 1
a461 1
 
d473 2
a474 2
     :   0, 0D0, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
     :   CACHE, IERR)
d477 3
a479 4
      CRVAL(1) = 0D0
      CRVAL(2) = PROJP(1)
      LONPOL = 999D0
      LATPOL = 999D0
d482 1
d487 1
a487 1
      CALL PGSCI (7)
d490 2
a491 2
     :   GRID1, 2, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM,
     :   256, IC, CACHE, IERR)
d506 1
a506 1
*        - Use of equivalences to define parameters for PGWCSL.
d525 2
d530 2
d533 3
a535 3
*     Reference pixel coordinates (EQUIVALENCE'd into NLDPRM above).
      CRPIX(1) = 192D0
      CRPIX(2) = 640D0
d538 14
a551 14
      PC(1,1) =  COSD(30D0)
      PC(1,2) =  SIND(30D0)
      PC(2,1) = -SIND(30D0)
      PC(2,2) =  COSD(30D0)

*     Coordinate increments (EQUIVALENCE'd into NLDPRM).
      CDELT(1) = -1D0/5D0
      CDELT(2) =  1D0/5D0

*     Spherical coordinate references (EQUIVALENCE'd into NLDPRM).
      CRVAL(1) = 332D0
      CRVAL(2) =  40D0
      LONPOL = -30D0
      LATPOL = 999D0
d565 1
d569 1
a569 1
      CALL PGSCI (6)
d572 1
a572 1
     :   5D0, 0, 5D0, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM,
d579 2
a580 1
*     Now reset to draw full grid lines rather than tick marks.
d584 1
a584 1
*     Draw the primary meridian and equator in yellow.
a586 1
      CALL PGSCI (2)
d588 2
a589 2
     :   GRID1, 1, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM,
     :   256, IC, CACHE, IERR)
d597 1
a602 1
      CALL PGSCI (3)
d604 2
a605 2
     :   GRID1, 2, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM,
     :   256, IC, CACHE, IERR)
d608 1
a620 1
      CALL PGSCI (4)
d622 2
a623 2
     :   GRID1, 4, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM,
     :   256, IC, CACHE, IERR)
d626 1
a638 1
      CALL PGSCI (5)
d640 2
a641 2
     :   GRID1, 4, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM,
     :   256, IC, CACHE, IERR)
d673 2
a681 1
      CALL PGSCI (6)
d683 2
a684 2
     :   GRID1, 2, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM,
     :   256, IC, CACHE, IERR)
d687 1
d691 5
a695 5
      CRVAL(1) = 0D0
      CRVAL(2) = 0D0
      LONPOL = 999D0
      LATPOL = 999D0
      CALL PGSCI (7)
d698 2
a699 2
     :   GRID1, 1, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM,
     :   256, IC, CACHE, IERR)
d715 1
a715 1
*        - Use of equivalences to define parameters for PGWCSL.
d743 2
d748 13
d762 3
a764 14
*     Reference pixel coordinates (EQUIVALENCE'd into NLDPRM above).
      CRPIX(1) = 226.0D0
      CRPIX(2) =  46.0D0

*     Linear transformation matrix (EQUIVALENCE'd into NLDPRM).
      ROTN = 15D0
      PC(1,1) =  COSD(ROTN)
      PC(1,2) =  SIND(ROTN)
      PC(2,1) = -SIND(ROTN)
      PC(2,2) =  COSD(ROTN)

*     Coordinate increments (EQUIVALENCE'd into NLDPRM).
      CDELT(1) = -1D0
      CDELT(2) =  1D0
d767 2
a768 4
      CRVAL(1) = 0D0
      CRVAL(2) = 0D0
      LONPOL = 999D0
      LATPOL = 999D0
d773 2
a774 3
*     projection it is possible to recentre it in longitude.  This will
*     modify CRPIX (via the LIN equivalence), and CRVAL and LONPOLE (via
*     the CEL equivalence).
d776 1
a776 2
      CALL CYLFIX (NAXIS, CRPIX, PC, CDELT, CTYPE, PROJP, CRVAL, LONPOL,
     :   LATPOL, IERR)
d779 2
a780 2
      IDENTS(1) = ''
      IDENTS(2) = ''
d786 3
a792 3
*     Reduced size lettering.
      CALL PGSCH (0.8*SCL)

d808 2
a809 2
     :   15D0, 0, 15D0, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM,
     :   256, IC, CACHE, IERR)
d812 2
a813 2
      CRPIX(1) = 226.0D0
      CRPIX(2) =  46.0D0
d815 4
a818 5
*     Galactic reference coordinates (EQUIVALENCE'd into NLDPRM).
      CRVAL(1) = 30D0
      CRVAL(2) = 35D0
      LONPOL = 999D0
      LATPOL = 999D0
d820 1
a820 2
      CALL CYLFIX (NAXIS, CRPIX, PC, CDELT, CTYPE, PROJP, CRVAL, LONPOL,
     :   LATPOL, IERR)
d825 1
a825 1
      IDENTS(3) = ''
d855 2
a856 2
     :   0, 0D0, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
     :   CACHE, IERR)
d911 1
a911 1
      IDENTS(3) = ''
d916 3
d924 1
a927 3
*     Normal size lettering.
      CALL PGSCH (1.0*SCL)

a928 1
      CALL PGSCI (1)
d951 1
a951 1
      IDENTS(1) = ''
d964 1
a964 1
 
d972 2
a973 2
      IDENTS(2) = ''
 
d980 54
@


1.21
log
@Handle output to PostScript devices.
@
text
@d29 1
a29 1
*   $Id: pgtest.f,v 1.20 2001/11/20 11:39:28 mcalabre Exp $
d34 2
a35 2
      INTEGER   AXEN(2), C0(7), CI(7), GCODE(2), IC, IERR, J, NLIPRM(2),
     :          WCS(4)
d37 4
a40 3
      DOUBLE PRECISION CDELT(2), CACHE(4,0:256), CRPIX(2), CRVAL(2),
     :          GRID1(0:10), GRID2(0:10), LATPOL, LIN(19), LNGPOL,
     :          NLDPRM(50), PC(2,2), PROJP(0:9), TIKLEN
d48 7
d56 2
a57 4
      EQUIVALENCE (NLDPRM(20), CRVAL)
      EQUIVALENCE (NLDPRM(22), LNGPOL)
      EQUIVALENCE (NLDPRM(23), LATPOL)
      EQUIVALENCE (NLDPRM(30), PROJP)
d76 2
a77 2
      AXEN(1) = 512
      AXEN(2) = 512
d81 2
a82 2
      TRC(1) = AXEN(1) + 0.5
      TRC(2) = AXEN(2) + 0.5
d84 2
a85 1
      CALL PGBEG (0, '?', 1, 1)
d97 1
a97 1
      LARGE = DEVTYP.EQ.'XWINDOWS'
d157 1
a157 1
      NLDPRM(5) = 360D0/(AXEN(1)-1)
d229 1
a229 1
      NLDPRM(5) = 720D0/(AXEN(1)+1)
d362 1
a362 1
      LNGPOL = 999D0
d444 1
a444 1
      LNGPOL = 150D0
d495 1
a495 1
      LNGPOL = 999D0
d562 1
a562 1
      LNGPOL = -30D0
d702 1
a702 1
      LNGPOL = 999D0
d735 2
a736 2
      AXEN(1) = 181
      AXEN(2) =  91
d740 2
a741 2
      TRC(1) = AXEN(1) + 0.5
      TRC(2) = AXEN(2) + 0.5
d750 1
a750 1
      CALL PGWNAD (0.0, 1.0, 0.0, REAL(AXEN(2))/REAL(AXEN(1)))
a759 2
      crpix(1) =  crpix(1) - 180d0

d761 5
a765 4
      PC(1,1) = 1D0
      PC(1,2) = 0D0
      PC(2,1) = 0D0
      PC(2,2) = 1D0
d774 1
a774 1
      LNGPOL = 999D0
d777 9
a785 2
      crval(1) =  crval(1) + 180d0
      crval(2) = -crval(2)
d820 3
d824 1
a824 1
*     Spherical coordinate references (EQUIVALENCE'd into NLDPRM).
d827 1
a827 1
      LNGPOL = 999D0
d830 2
a831 2
      crval(1) =  crval(1) + 180d0
      crval(2) = -crval(2)
d894 2
a895 2
      AXEN(1) = 512
      AXEN(2) = 512
d899 2
a900 2
      TRC(1) = AXEN(1) + 0.5
      TRC(2) = AXEN(2) + 0.5
@


1.20
log
@Declare COSD and SIND for linux; fixed an argument in a call to PGSBOX.
@
text
@d29 1
a29 1
*   $Id: pgtest.f,v 1.19 2001/02/19 02:18:18 mcalabre Exp mcalabre $
d33 1
d36 1
a36 1
      REAL      BLC(2), TRC(2)
d40 2
a41 1
      CHARACTER CTYPE(2)*8, ESCAPE*1, FCODE(2)*4, IDENTS(3)*80, OPT(2)*1
d78 21
a98 5
      CALL PGBEG (0, '/xwindow', 1, 1)
*     Black.
      CALL PGSCR (0, 0.0, 0.0, 0.0)
*     White.
      CALL PGSCR (1, 1.0, 1.0, 1.0)
d120 1
a121 2
      CALL PGVSTD ()
      CALL PGWNAD (0.0, 1.0, 0.0, 1.0)
d165 1
a165 1
      CALL PGSCH (1.0)
d241 1
a241 1
      CALL PGSCH (1.0)
d297 1
a297 1
      CALL PGSCH (1.0)
d375 1
a375 1
      CALL PGSCH (0.8)
d453 1
a453 1
      CALL PGSCH (0.8)
d567 1
a567 1
      CALL PGSCH (0.9)
d737 6
a742 1
      CALL PGVSTD
d787 1
a787 1
      CALL PGSCH (0.8)
d871 5
a875 1
      CALL PGVSTD
d920 1
a920 1
      CALL PGSCH (1.0)
@


1.19
log
@PGCRVL -> PGSBOX.
@
text
@d4 1
a4 1
*   Copyright (C) 1997-2000, Mark Calabretta
d29 1
a29 1
*   $Id: pgtest.f,v 1.18 2000/12/07 00:04:49 mcalabre Exp mcalabre $
d58 2
d157 1
a157 1
     :   0, 0D0, .FALSE., LNGVEL, 1, 1, 7, ' ', 0D0, NLDPRM, 256, IC,
d740 2
a741 2
      CDELT(1) = -1.0D0
      CDELT(2) =  1.0D0
@


1.18
log
@Added new plate carree example.
@
text
@d3 1
a3 1
*   PGCRVL - a non-linear coordinate axis plotter for PGPLOT.
d20 1
a20 1
*   Correspondence concerning PGCRVL may be directed to:
d29 1
a29 1
*   $Id: pgtest.f,v 1.17 2000/05/24 05:24:06 mcalabre Exp mcalabre $
d35 1
d69 5
d109 1
a109 1
*     velocity.  PGCRVL uses subroutine LNGVEL.  Features:
d148 1
a148 1
*     Draw tick marks for longitude and grid lines for velocity in yellow.
d154 2
a155 2
      CALL PGCRVL (AXEN, IDENTS, OPT, -1, 0, C0, GCODE, 2D0, 0, 0D0, 0,
     :   0D0, .FALSE., LNGVEL, 1, 1, 7, ' ', 0D0, NLDPRM, 256, IC,
d164 3
a166 3
      CALL PGCRVL (AXEN, IDENTS, OPT, 0, 0, C0, GCODE, 0D0, 1, GRID1,
     :   1, GRID2, .FALSE., LNGVEL, 1, 1, 7, ' ', 0, NLDPRM, 256, IC,
     :   CACHE, IERR)
d178 1
a178 1
*     receiver which scans simultaneously in frequency.  PGCRVL uses
d185 1
a185 1
*          0 - 720 degrees but PGCRVL is set to normalize this to two
d228 2
a229 2
      CALL PGCRVL (AXEN, IDENTS, OPT, 0, 9900, C0, GCODE, 2D0, 0, 0D0,
     :   0, 0D0, .FALSE., FSCAN, 1, 1, 7, ' ', 0, NLDPRM, 256, IC,
d239 1
a239 1
*     Third example: Z versus time plot.  PGCRVL uses subroutine PGCRFN.
d282 1
a282 1
      CALL PGCRVL (AXEN, IDENTS, OPT, 0, 0, C0, GCODE, 2D0, 0, 0D0,
d294 1
a294 1
*     pole.  PGCRVL uses subroutine PGWCSL to interface to WCSLIB.
d362 3
a364 3
      CALL PGCRVL (AXEN, IDENTS, OPT, 0, 1224, C0, GCODE, 0D0, 0, 0D0,
     :   0, 0D0, .FALSE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     :   IC, CACHE, IERR)
d373 1
a373 1
*     Fifth example: Conic equal area.  PGCRVL uses subroutine PGWCSL
d460 1
a460 1
*     Draw the celestial grid letting PGCRVL choose the increments.
d462 2
a463 2
      CALL PGCRVL (AXEN, IDENTS, OPT, 0, 0, CI, GCODE, 0D0, 0, 0D0, 0,
     :   0D0, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
d479 3
a481 3
      CALL PGCRVL (AXEN, IDENTS, OPT, -1, 0, C0, GCODE, 0D0, 2, GRID1,
     :   2, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     :   IC, CACHE, IERR)
d491 1
a491 1
*     Sixth example: Polyconic with colour-coded grid.  PGCRVL uses
d556 3
a558 3
      CALL PGCRVL (AXEN, IDENTS, OPT, -1, 0, C0, GCODE, TIKLEN, 0, 5D0,
     :   0, 5D0, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
     :   CACHE, IERR)
d572 3
a574 3
      CALL PGCRVL (AXEN, IDENTS, OPT, -1, 0, C0, GCODE, 0D0, 1, GRID1,
     :   1, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     :   IC, CACHE, IERR)
d588 3
a590 3
      CALL PGCRVL (AXEN, IDENTS, OPT, -1, 0, C0, GCODE, 0D0, 3, GRID1,
     :   2, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     :   IC, CACHE, IERR)
d606 3
a608 3
      CALL PGCRVL (AXEN, IDENTS, OPT, -1, 0, C0, GCODE, 0D0, 8, GRID1,
     :   4, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     :   IC, CACHE, IERR)
d624 3
a626 3
      CALL PGCRVL (AXEN, IDENTS, OPT, -1, 0, C0, GCODE, 0D0, 8, GRID1,
     :   4, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     :   IC, CACHE, IERR)
d632 1
a632 1
*     PGCRVL which finally produces the labels.
d658 1
a658 1
*     Tell PGCRVL what edges to label.
d666 1
a666 1
      CALL PGCRVL (AXEN, IDENTS, OPT, 2333, 0, CI, GCODE, 0D0, 4,
d680 3
a682 3
      CALL PGCRVL (AXEN, IDENTS, OPT, -1, 0, C0, GCODE, 0D0, 2, GRID1,
     :   1, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     :   IC, CACHE, IERR)
d692 1
a692 1
*     Seventh example: Plate Carree projection.  PGCRVL uses subroutine
d712 5
d779 3
a781 3
      CALL PGCRVL (AXEN, IDENTS, OPT, 2100, 0, CI, GCODE, 0D0, 0, 15D0,
     :   0, 15D0, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     :   IC, CACHE, IERR)
d823 1
a823 1
*     Draw the celestial grid letting PGCRVL choose the increments.
d825 2
a826 2
      CALL PGCRVL (AXEN, IDENTS, OPT, 21, 0, CI, GCODE, 0D0, 0, 0D0, 0,
     :   0D0, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
d838 1
a838 1
*     labelling.  PGCRVL uses subroutine PGCRFN.  Features:
d853 5
d897 1
a897 1
      CALL PGCRVL (AXEN, IDENTS, OPT, 21, 0, C0, GCODE, 2D0, 0, 0D0,
d902 1
a902 1
*     crossing table on the first call to PGCRVL.  We now want to
d923 1
a923 1
      CALL PGCRVL (AXEN, IDENTS, OPT, 12000, 0, C0, GCODE, 2D0, 0,
d944 3
a946 3
      CALL PGCRVL (AXEN, IDENTS, OPT, 100, 0, C0, GCODE, 2D0, 0, 0D0,
     :   0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM, 256,
     :   IC, CACHE, IERR)
@


1.17
log
@Fix for a portability problem, some Fortran compilers require the backslash
character to be escaped (with a backslash - hence '\\') and others don't.
@
text
@d29 1
a29 1
*   $Id: pgtest.f,v 1.16 2000/05/10 04:54:48 mcalabre Exp mcalabre $
a376 1
*        - Discovery of grid lines which do not cross any axis.
d686 141
a826 1
*     Seventh example: linear-linear plot with two types of alternative
d836 6
a858 1
*-----------------------------------------------------------------------
@


1.16
log
@Updated date in copyright notice.
@
text
@d29 1
a29 1
*   $Id: pgwcsl.f,v 1.8 1999/12/01 04:18:09 mcalabre Exp mcalabre $
d38 1
a38 1
      CHARACTER CTYPE(2)*8, FCODE(2)*4, IDENTS(3)*80, OPT(2)*1
d54 3
d205 1
a205 1
      IDENTS(2) = '\\gn/Hz'
d714 3
a716 1
      IDENTS(1) = 'temperature of frog (\\uo\\dC)'
d778 1
a778 1
      IDENTS(1) = '(\\uo\\dF)'
@


1.15
log
@Amended copyright notice; converted continuation character from "*" to ":".
@
text
@d4 1
a4 1
*   Copyright (C) 1997-1999, Mark Calabretta
d28 2
@


1.14
log
@Added copyright notice.
@
text
@d20 1
a20 1
*   Correspondence concerning WCSLIB may be directed to:
d32 1
a32 1
     *          WCS(4)
d34 2
a35 2
     *          GRID1(0:10), GRID2(0:10), LATPOL, LIN(19), LNGPOL,
     *          NLDPRM(50), PC(2,2), PROJP(0:9), TIKLEN
d144 2
a145 2
     *   0D0, .FALSE., LNGVEL, 1, 1, 7, ' ', 0D0, NLDPRM, 256, IC,
     *   CACHE, IERR)
d154 2
a155 2
     *   1, GRID2, .FALSE., LNGVEL, 1, 1, 7, ' ', 0, NLDPRM, 256, IC,
     *   CACHE, IERR)
d218 2
a219 2
     *   0, 0D0, .FALSE., FSCAN, 1, 1, 7, ' ', 0, NLDPRM, 256, IC,
     *   CACHE, IERR)
d272 2
a273 2
     *   0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM, 256,
     *   IC, CACHE, IERR)
d352 2
a353 2
     *   0, 0D0, .FALSE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     *   IC, CACHE, IERR)
d453 2
a454 2
     *   0D0, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
     *   CACHE, IERR)
d470 2
a471 2
     *   2, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     *   IC, CACHE, IERR)
d547 2
a548 2
     *   0, 5D0, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
     *   CACHE, IERR)
d563 2
a564 2
     *   1, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     *   IC, CACHE, IERR)
d579 2
a580 2
     *   2, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     *   IC, CACHE, IERR)
d597 2
a598 2
     *   4, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     *   IC, CACHE, IERR)
d615 2
a616 2
     *   4, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     *   IC, CACHE, IERR)
d657 2
a658 2
     *   GRID1, 2, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM,
     *   256, IC, CACHE, IERR)
d671 2
a672 2
     *   1, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     *   IC, CACHE, IERR)
d731 2
a732 2
     *   0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM, 256,
     *   IC, CACHE, IERR)
d757 2
a758 2
     *   0D0, 0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM,
     *   256, IC, CACHE, IERR)
d778 2
a779 2
     *   0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM, 256,
     *   IC, CACHE, IERR)
@


1.13
log
@*** empty log message ***
@
text
@d1 28
@


1.12
log
@*** empty log message ***
@
text
@d3 2
a4 1
      INTEGER   AXEN(2), GCODE(2), IC, IERR, J, NLIPRM(2), WCS(4)
d8 1
a8 1
      CHARACTER CTYPE(2)*8, FCODE(2)*4, IDENTS(3)*80
d53 8
d103 3
d115 1
a115 1
      CALL PGCRVL (AXEN, IDENTS, 'F ', -1, 0, -1, GCODE, 2D0, 0, 0D0, 0,
d125 1
a125 1
      CALL PGCRVL (AXEN, IDENTS, 'F ', 0, 0, -1, GCODE, 0D0, 1, GRID1,
d175 3
d189 1
a189 1
      CALL PGCRVL (AXEN, IDENTS, 'DL', 0, 9900, -1, GCODE, 0D0, 0, 0D0,
d231 3
d243 1
a243 1
      CALL PGCRVL (AXEN, IDENTS, 'L ', 0, 0, -1, GCODE, 2D0, 0, 0D0,
d309 3
d323 1
a323 1
      CALL PGCRVL (AXEN, IDENTS, 'GE', 0, 1224, -1, GCODE, 0D0, 0, 0D0,
d388 3
d413 8
d424 1
a424 1
      CALL PGCRVL (AXEN, IDENTS, 'EE', 0, 0, 10, GCODE, 0D0, 0, 0D0, 0,
d441 1
a441 1
      CALL PGCRVL (AXEN, IDENTS, 'EE', -1, 0, -1, GCODE, 0D0, 2, GRID1,
d506 3
d518 1
a518 1
      CALL PGCRVL (AXEN, IDENTS, 'HB', -1, 0, -1, GCODE, TIKLEN, 0, 5D0,
d534 1
a534 1
      CALL PGCRVL (AXEN, IDENTS, 'HB', -1, 0, -1, GCODE, 0D0, 1, GRID1,
d550 1
a550 1
      CALL PGCRVL (AXEN, IDENTS, 'HB', -1, 0, -1, GCODE, 0D0, 3, GRID1,
d568 1
a568 1
      CALL PGCRVL (AXEN, IDENTS, 'HB', -1, 0, -1, GCODE, 0D0, 8, GRID1,
d586 1
a586 1
      CALL PGCRVL (AXEN, IDENTS, 'HB', -1, 0, -1, GCODE, 0D0, 8, GRID1,
d612 8
d628 1
a628 1
      CALL PGCRVL (AXEN, IDENTS, 'HB', 2333, 0, 10, GCODE, 0D0, 4,
d642 1
a642 1
      CALL PGCRVL (AXEN, IDENTS, 'HB', -1, 0, -1, GCODE, 0D0, 2, GRID1,
d685 3
d702 1
a702 1
      CALL PGCRVL (AXEN, IDENTS, '  ', 21, 0, -1, GCODE, 2D0, 0, 0D0,
d728 1
a728 1
      CALL PGCRVL (AXEN, IDENTS, '  ', 12000, 0, -1, GCODE, 2D0, 0,
d749 1
a749 1
      CALL PGCRVL (AXEN, IDENTS, '  ', 100, 0, -1, GCODE, 2D0, 0, 0D0,
@


1.11
log
@*** empty log message ***
@
text
@d4 1
a4 1
      DOUBLE PRECISION CDELT(2), CROSS(4,256), CRPIX(2), CRVAL(2),
d102 1
a102 1
      IC = 0
d105 1
a105 1
     *   CROSS, IERR)
d115 1
a115 1
     *   CROSS, IERR)
d173 1
a173 1
      IC = 0
d176 1
a176 1
     *   CROSS, IERR)
d224 1
a224 1
      IC = 0
d227 1
a227 1
     *   IC, CROSS, IERR)
d301 1
a301 1
      IC = 0
d304 1
a304 1
     *   IC, CROSS, IERR)
d391 1
a391 1
      IC = 0
d394 1
a394 1
     *   CROSS, IERR)
d408 1
a408 1
      IC = 0
d411 1
a411 1
     *   IC, CROSS, IERR)
d482 1
a482 1
      IC = 0
d485 1
a485 1
     *   CROSS, IERR)
d501 1
a501 1
     *   IC, CROSS, IERR)
d517 1
a517 1
     *   IC, CROSS, IERR)
d535 1
a535 1
     *   IC, CROSS, IERR)
d553 1
a553 1
     *   IC, CROSS, IERR)
d587 1
a587 1
     *   256, IC, CROSS, IERR)
d598 1
d601 1
a601 1
     *   IC, CROSS, IERR)
d644 1
a644 1
*     stored for labels on top edge while GCODE(2) = 2 causes
d655 1
a655 1
      IC = 0
d658 1
a658 17
     *   IC, CROSS, IERR)

*     The alternative temperature scale in Fahrenheit is to be
*     constructed with a new set of tick marks but the distance is
*     to be rescaled to feet for the existing grid lines.
      NLDPRM(3) = NLDPRM(3)*1.8D0
      NLDPRM(5) = NLDPRM(5)*1.8D0 + 32D0

*     Draw tick marks at the top for the first coordinate, don't redo
*     grid lines for the second.
      GCODE(1) = -100
      GCODE(2) = 0

*     Set LABCTL = -1 to defer labelling.
      CALL PGCRVL (AXEN, IDENTS, '  ', -1, 0, -1, GCODE, 2D0, 0, 0D0,
     *   0, 0D0, .FALSE., PGCRFN, 8, 2, 4, FCODE, NLIPRM, NLDPRM, 256,
     *   IC, CROSS, IERR)
d669 1
a669 1
         IF (CROSS(1,J).EQ.4D0) THEN
d671 2
a672 2
            CROSS(4,J) = CROSS(4,J) * 1D3/(25.4*12D0)
            CROSS(4,J) = AINT(CROSS(4,J)*10D0 + 0.5D0)/10D0
d676 2
a677 2
*     New annotation.
      IDENTS(1) = '(\\uo\\dF)'
d680 3
a682 3
*     Set LABCTL = 12100 to label the top and right edges without
*     redoing the ticks or grid lines.
      CALL PGCRVL (AXEN, IDENTS, '  ', 12100, -1, -1, GCODE, 2D0, 0,
d684 22
a705 1
     *   256, IC, CROSS, IERR)
@


1.10
log
@*** empty log message ***
@
text
@d320 1
@


1.9
log
@*** empty log message ***
@
text
@d14 4
a17 4
      EQUIVALENCE (NLDPRM(20), PROJP)
      EQUIVALENCE (NLDPRM(41), CRVAL)
      EQUIVALENCE (NLDPRM(43), LNGPOL)
      EQUIVALENCE (NLDPRM(44), LATPOL)
@


1.8
log
@*** empty log message ***
@
text
@d3 1
a3 1
      INTEGER   AXEN(2), GCODE(2), IC, IERR, WCS(4)
d7 1
a7 1
      CHARACTER CTYPE(2)*8, IDENTS(3)*80
d21 1
d61 12
a72 3
*     velocity.  PGCRVL uses subroutine LNGVEL and chooses what axes to
*     label and how to label them, extracting a common scaling factor
*     from the velocity.
d103 1
a103 1
      CALL PGCRVL (AXEN, IDENTS, 'F ', -1, -1, GCODE, 2D0, 0, 0D0, 0,
d113 2
a114 2
      CALL PGCRVL (AXEN, IDENTS, 'F ', 0, -1, GCODE, 0D0, 1, GRID1, 1,
     *   GRID2, .FALSE., LNGVEL, 1, 1, 7, ' ', 0, NLDPRM, 256, IC,
d124 1
a124 1
 
d127 15
a141 5
*     receiver which scans simultaneously in frequency.  It illustrates
*     the use of cyclic and logarithmic axes.  Subroutine FSCAN returns
*     the azimuth in the range 0 - 720 degrees but PGCRVL is set to
*     normalize this to two cycles of 0 - 360 degrees.  It chooses what
*     axes to label and how to label them.
d146 1
a146 1
 
d149 2
a150 2
      NLDPRM(4) = 8.5
 
d157 1
a157 1
 
d166 1
a166 1
 
d170 1
a170 1
*     Setting LABCTL = -9900 forces all logarithmic grid lines to be
d174 54
a227 3
      CALL PGCRVL (AXEN, IDENTS, 'DL', -9900, -1, GCODE, 0D0, 0, 0D0, 0,
     *   0D0, .FALSE., FSCAN, 1, 1, 7, ' ', 0, NLDPRM, 256, IC, CROSS,
     *   IERR)
d236 21
a256 5
*     Third example: A simple SIN projection near the south celestial
*     pole using PGWCSL and WCSLIB to demonstrate sexagesimal axis
*     labelling and angles cycling through 360 degrees.  PGCRVL chooses
*     the coordinate increments and determines what axes to label using
*     a higher grid density than normal.
d261 1
a261 1
 
d265 1
a265 1
 
d271 1
a271 1
 
d294 1
a294 1
 
d299 1
a299 1
*     coordinate by specifying LABCTL = -1224.
d302 4
a305 4
      CALL PGCRVL (AXEN, IDENTS, 'GE', -1224, -1, GCODE, 0D0, 0, 0D0, 0,
     *   0D0, .FALSE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
     *   CROSS, IERR)
 
d308 1
a308 1
 
d313 19
a331 4
*     Fourth example: Conic equal area.  Uses PGWCSL and WCSLIB.
*     This demonstrates some of the nasty things that can happen with
*     general curvilinear coordinate systems.  PGCRVL chooses the
*     coordinate increments and determines what axes to label.
d336 1
a336 1
 
d340 1
a340 1
 
d346 1
a346 1
 
d350 1
a350 1
 
d369 1
a369 1
 
d388 1
a388 1
 
d391 1
a391 1
      CALL PGCRVL (AXEN, IDENTS, 'EE', 0, 10, GCODE, 0D0, 0, 0D0, 0,
d400 1
a400 1
 
d408 3
a410 3
      CALL PGCRVL (AXEN, IDENTS, 'EE', -1, -1, GCODE, 0D0, 2, GRID1, 2,
     *   GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
     *   CROSS, IERR)
d415 1
a415 1
 
d420 23
a442 3
*     Fifth example: Polyconic with colour-coded grid.  Uses PGWCSL and
*     WCSLIB for a virtuoso performance.  This example illustrates
*     basic management of the axis-crossing table.
d452 1
a452 1
*     Rotate 30 degrees, just because it's Friday.
d472 1
a472 1
 
d482 2
a483 2
      CALL PGCRVL (AXEN, IDENTS, 'HB', -1, -1, GCODE, TIKLEN, 0, 5D0, 0,
     *   5D0, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
d498 3
a500 3
      CALL PGCRVL (AXEN, IDENTS, 'HB', -1, -1, GCODE, 0D0, 1, GRID1, 1,
     *   GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
     *   CROSS, IERR)
d514 3
a516 3
      CALL PGCRVL (AXEN, IDENTS, 'HB', -1, -1, GCODE, 0D0, 3, GRID1, 2,
     *   GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
     *   CROSS, IERR)
d532 3
a534 3
      CALL PGCRVL (AXEN, IDENTS, 'HB', -1, -1, GCODE, 0D0, 8, GRID1, 4,
     *   GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
     *   CROSS, IERR)
d550 3
a552 3
      CALL PGCRVL (AXEN, IDENTS, 'HB', -1, -1, GCODE, 0D0, 8, GRID1, 4,
     *   GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
     *   CROSS, IERR)
d558 1
a558 2
*     PGCRVL which finally produces the labels.  It will zero the table
*     index on return.
d576 1
a576 1
*     Tell PGCRVL what axes to label.
d584 3
a586 3
      CALL PGCRVL (AXEN, IDENTS, 'HB', 2333, 10, GCODE, 0D0, 4, GRID1,
     *   2, GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256,
     *   IC, CROSS, IERR)
d597 3
a599 3
      CALL PGCRVL (AXEN, IDENTS, 'HB', -1, -1, GCODE, 0D0, 2, GRID1, 1,
     *   GRID2, .TRUE., PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
     *   CROSS, IERR)
d603 98
@


1.7
log
@*** empty log message ***
@
text
@d20 1
a20 1
      EXTERNAL RADAR
d118 4
a121 4
*     the use of cyclic axes.  Subroutine RADAR returns the azimuth in
*     the range 0 - 720 degrees but PGCRVL is set to normalize this to
*     two cycles of 0 - 360 degrees.  It chooses what axes to label and
*     how to label them.
d129 1
a129 1
      NLDPRM(4) = 1420D0
d133 1
a133 1
      NLDPRM(6) = 1D0
d136 1
a136 1
      NLDPRM(7) = -1D0
d140 2
a141 2
      IDENTS(2) = 'frequency (MHz)'
      IDENTS(3) = 'Radar scan'
d150 2
d154 2
a155 2
      CALL PGCRVL (AXEN, IDENTS, 'D ', 0, -1, GCODE, 0D0, 0, 0D0, 0,
     *   0D0, .FALSE., RADAR, 1, 1, 7, ' ', 0, NLDPRM, 256, IC, CROSS,
@


1.6
log
@*** empty log message ***
@
text
@d61 2
a62 1
*     label.
d80 1
a80 1
      IDENTS(1) = 'longitude (deg)'
d93 1
a93 1
      CALL PGCRVL (AXEN, IDENTS, 'C ', -1, -1, GCODE, 2D0, 0, 0D0, 0,
d103 1
a103 1
      CALL PGCRVL (AXEN, IDENTS, 'C ', 0, -1, GCODE, 0D0, 1, GRID1, 1,
d120 2
a121 1
*     two cycles of 0 - 360 degrees.  It chooses what axes to label.
d139 1
a139 1
      IDENTS(1) = 'azimuth (deg)'
d152 1
a152 1
      CALL PGCRVL (AXEN, IDENTS, 'A ', 0, -1, GCODE, 0D0, 0, 0D0, 0,
d163 5
a167 3
*     Third example: A simple SIN projection using PGWCSL and WCSLIB
*     to demonstrate sexagesimal axis labelling.  PGCRVL chooses the
*     coordinate increments and determines what axes to label.
d174 1
a174 1
      CRPIX(1) = 256D0
d184 2
a185 2
      CDELT(1) = -1D0/360D0
      CDELT(2) =  1D0/360D0
d188 2
a189 2
      CRVAL(1) =   0D0
      CRVAL(2) = -85D0
d209 2
a210 1
*     Draw the celestial grid.
d213 1
a213 1
      CALL PGCRVL (AXEN, IDENTS, 'AB', 0, -1, GCODE, 0D0, 0, 0D0, 0,
d258 2
a259 2
      IDENTS(1) = 'longitude (deg)'
      IDENTS(2) = 'latitude (deg)'
d287 1
a287 1
      CALL PGCRVL (AXEN, IDENTS, 'AB', 0, 10, GCODE, 0D0, 0, 0D0, 0,
d304 1
a304 1
      CALL PGCRVL (AXEN, IDENTS, 'AB', -1, -1, GCODE, 0D0, 2, GRID1, 2,
d345 2
a346 2
      IDENTS(1) = 'longitude (deg)'
      IDENTS(2) = 'latitude (deg)'
d358 1
a358 1
      CALL PGCRVL (AXEN, IDENTS, 'AB', -1, -1, GCODE, TIKLEN, 0, 5D0, 0,
d374 1
a374 1
      CALL PGCRVL (AXEN, IDENTS, 'AB', -1, -1, GCODE, 0D0, 1, GRID1, 1,
d390 1
a390 1
      CALL PGCRVL (AXEN, IDENTS, 'AB', -1, -1, GCODE, 0D0, 3, GRID1, 2,
d408 1
a408 1
      CALL PGCRVL (AXEN, IDENTS, 'AB', -1, -1, GCODE, 0D0, 8, GRID1, 4,
d426 1
a426 1
      CALL PGCRVL (AXEN, IDENTS, 'AB', -1, -1, GCODE, 0D0, 8, GRID1, 4,
d461 1
a461 1
      CALL PGCRVL (AXEN, IDENTS, 'AB', 2333, 10, GCODE, 0D0, 4, GRID1,
d474 1
a474 1
      CALL PGCRVL (AXEN, IDENTS, 'AB', -1, -1, GCODE, 0D0, 2, GRID1, 1,
@


1.5
log
@*** empty log message ***
@
text
@d3 1
a3 1
      INTEGER   AXEN(2), GCODE(2), IC, IERR, j, k, WCS(4)
d6 1
a6 1
     *          NLDPRM(50), PC(2,2), PROJP(0:9)
d34 16
a49 8
      CALL PGSCR (0, 0.00, 0.00, 0.00)
      CALL PGSCR (1, 1.00, 1.00, 1.00)
      CALL PGSCR (2, 1.00, 1.00, 0.00)
      CALL PGSCR (3, 1.00, 1.00, 1.00)
      CALL PGSCR (4, 0.50, 0.50, 0.80)
      CALL PGSCR (5, 0.80, 0.50, 0.50)
      CALL PGSCR (6, 0.70, 0.70, 0.70)
      CALL PGSCR (7, 0.30, 0.50, 0.30)
d60 2
a61 1
*     velocity.  Uses subroutine LNGVEL.
d78 1
a78 6
      CALL PGSCH (1.0)
      CALL PGSCI (1)
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

      CALL PGSCI (1)

d82 8
a89 3
      GRID1(1) = 180D0
      GRID2(1) = 0D0
      GCODE(1) = 2
d92 2
a93 2
      CALL PGCRVL (AXEN, IDENTS, 'C ', GCODE, 0D0, 1, GRID1, 1, GRID2,
     *   .FALSE.,  -1, LNGVEL, 1, 1, 7, ' ', 0, NLDPRM, 256, IC,
d96 3
a98 2
      CALL PGSCI (2)
      GCODE(1) = 1
d100 8
a107 3
      CALL PGCRVL (AXEN, IDENTS, 'C ', GCODE, 2D0, 0, 0D0, 0, 0D0,
     *   .FALSE., 0, LNGVEL, 1, 1, 7, ' ', 0D0, NLDPRM, 256, IC, CROSS,
     *   IERR)
a108 4
*     write (6, 10) ((cross(j,k), j=1,4), k=1,ic)
*10   format (f5.1,f12.5,f5.1,1p,g15.5,0p)
*     write (6, 20) cross(1,ic+1), cross(2,ic+1)
*20   format ('Separability indices: ',2f10.4)
d119 1
a119 1
*     two cycles of 0 - 360 degrees.
d136 1
a136 6
      CALL PGSCH (1.0)
      CALL PGSCI (1)
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)
 
      CALL PGSCI (1)
 
d140 2
a141 2
      GRID1(1) = 0D0
      GRID2(1) = 0D0
d144 5
d150 2
a151 2
      CALL PGCRVL (AXEN, IDENTS, 'A ', GCODE, 0D0, 0, GRID1, 0, GRID2,
     *   .FALSE., 0, RADAR, 1, 1, 7, ' ', 0, NLDPRM, 256, IC, CROSS,
a152 1
*     2021
d154 3
d162 4
a165 1
*     to demonstrate sexagesimal axis labelling.
d169 1
d173 1
d179 1
d183 1
d189 1
a189 1
*     Set parameters for an NCP projection.
d193 1
d197 2
d202 1
a202 2
*     Draw celestial grid letting PGCRVL choose the increments.
      CALL PGSCI (7)
d204 3
d208 5
a213 3
      CALL PGCRVL (AXEN, IDENTS, 'AB', GCODE, 0D0, 0, 0D0, 0, 0D0,
     *   .FALSE., 0, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
     *   CROSS, IERR)
d221 4
a224 1
*     general curvilinear coordinate systems.
d228 1
d232 1
d238 1
d246 5
a250 3
      CALL PGSCH (1.0)
      CALL PGSCI (1)
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)
d252 1
d256 2
d260 3
d264 15
a278 5
*     Draw celestial grid letting PGCRVL choose the increments.
      CRVAL(1) = 90D0
      CRVAL(2) = 30D0
      LNGPOL = 150D0
      LATPOL = 999D0
d280 1
a280 2
      CALL PGSCI (2)
      CALL PGSCH (0.8)
d282 2
a283 2
      CALL PGCRVL (AXEN, IDENTS, 'AB', GCODE, 0D0, 0, 0D0, 0, 0D0,
     *   .TRUE., 0, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
d286 3
a288 3
*     Native grid in green (delineates boundary).
      CRVAL(1) =  0D0
      CRVAL(2) = 60D0
d292 1
d299 2
a300 2
      CALL PGCRVL (AXEN, IDENTS, 'AB', GCODE, 0D0, 2, GRID1, 2, GRID2,
     *   .TRUE., -1, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
d302 4
d314 2
d319 1
d329 1
d333 1
a333 5
      CALL PGSCH (1.0)
      CALL PGSCI (1)
      CALL PGBOX ('BC', 0.0, 0, 'BC', 0.0, 0)

*     Draw celestial grid.
d339 1
d343 3
d347 1
a347 2
*     Tick marks every 5 degrees.
      CALL PGSCI (6)
d350 2
d353 2
a354 2
      CALL PGCRVL (AXEN, IDENTS, 'AB', GCODE, -2D0, 0, 5D0, 0, 5D0,
     *   .TRUE., -1, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
d361 5
a365 1
*     Primary meridian and equator in yellow.
d369 2
a370 4
      GCODE(1) = 2
      GCODE(2) = 2
      CALL PGCRVL (AXEN, IDENTS, 'AB', GCODE, 0D0, 1, GRID1, 1, GRID2,
     *   .TRUE., -1, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
d378 1
a378 1
*     90 degree meridians and poles in white.
d385 2
a386 2
      CALL PGCRVL (AXEN, IDENTS, 'AB', GCODE, 0D0, 3, GRID1, 2, GRID2,
     *   .TRUE., -1, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
d389 1
a389 1
*     First set of 15 degree meridians and parallels in blue.
d403 2
a404 2
      CALL PGCRVL (AXEN, IDENTS, 'AB', GCODE, 0D0, 8, GRID1, 4, GRID2,
     *   .TRUE., -1, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
d407 1
a407 1
*     Second set of 15 degree meridians and parallels in red.
a419 1
      CALL PGSCH (0.9)
d421 2
a422 2
      CALL PGCRVL (AXEN, IDENTS, 'AB', GCODE, 0D0, 8, GRID1, 4, GRID2,
     *   .TRUE., -1, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
d432 17
a448 1
*     45 degree meridians and parallels in grey.
d456 3
a458 3
      CALL PGCRVL (AXEN, IDENTS, 'AB', GCODE, 0D0, 4, GRID1, 2, GRID2,
     *   .TRUE., 2333, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
     *   CROSS, IERR)
d469 2
a470 2
      CALL PGCRVL (AXEN, IDENTS, 'AB', GCODE, 0D0, 2, GRID1, 1, GRID2,
     *   .TRUE., -1, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC,
d472 4
@


1.4
log
@*** empty log message ***
@
text
@d7 1
a7 1
      CHARACTER CTYPE(2)*8
d19 2
a21 1
      EXTERNAL LNGVEL
a71 1
      CALL PGLAB ('longitude (deg)', 'velocity (m/s)', 'HI line')
d74 4
d83 3
a85 2
      CALL PGCRVL (AXEN, 'A ', GCODE, 0D0, 1, GRID1, 1, GRID2, .FALSE., 
     *   -1, LNGVEL, 1, 1, 7, ' ', 0D0, NLDPRM, 256, IC, CROSS, IERR)
d90 73
a162 2
      CALL PGCRVL (AXEN, 'A ', GCODE, 2D0, 0, 0D0, 0, 0D0, .FALSE.,
     *   21, LNGVEL, 1, 1, 7, ' ', 0D0, NLDPRM, 256, IC, CROSS, IERR)
d164 4
a167 4
      write (6, 10) ((cross(j,k), j=1,4), k=1,ic)
 10   format (f5.1,f12.5,f5.1,1p,g15.5,0p)
      write (6, 20) cross(1,ic+1), cross(2,ic+1)
 20   format ('Separability indices: ',2f10.4)
d169 19
d192 1
a192 1
*     Second example: Conic equal area.  Uses PGWCSL and WCSLIB.
a215 1
      CALL PGLAB ('', '', 'WCS conic equal area projection')
d217 3
d230 1
a230 1
      CALL PGSCH (0.5)
d232 3
a234 6
      CALL PGCRVL (AXEN, 'AA', GCODE, 0D0, 0, 0D0, 0, 0D0, .TRUE.,
     *   1221, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS,
     *   IERR)

      write (6, 10) ((cross(j,k), j=1,4), k=1,ic)
      write (6, 20) cross(1,ic+1), cross(2,ic+1)
d248 3
a250 3
      CALL PGCRVL (AXEN, 'AA', GCODE, 0D0, 2, GRID1, 2, GRID2, .TRUE.,
     *   -1, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS,
     *   IERR)
d256 3
a258 2
*     Third example: Polyconic with colour-coded grid.  Uses PGWCSL and
*     WCSLIB for a virtuoso performance.
a273 1
*     These axis labels are contentious.
a276 1
      CALL PGLAB ('longitude', 'latitude', 'WCS polyconic projection')
d279 1
a279 1
      CRVAL(1) = -28D0
d284 4
d293 7
a299 3
      CALL PGCRVL (AXEN, 'AA', GCODE, 2D0, 0, 5D0, 0, 5D0, .TRUE.,
     *   -1, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS,
     *   IERR)
d307 8
a314 4
      IC = 0
      CALL PGCRVL (AXEN, 'AA', GCODE, 0D0, 1, GRID1, 1, GRID2, .TRUE.,
     *   -1, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS,
     *   IERR)
d323 3
a325 3
      CALL PGCRVL (AXEN, 'AA', GCODE, 0D0, 3, GRID1, 2, GRID2, .TRUE.,
     *   -1, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS,
     *   IERR)
d341 3
a343 3
      CALL PGCRVL (AXEN, 'AA', GCODE, 0D0, 8, GRID1, 4, GRID2, .TRUE.,
     *   -1, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS,
     *   IERR)
d358 1
a358 1
      CALL PGSCH (0.7)
d360 10
a369 3
      CALL PGCRVL (AXEN, 'AA', GCODE, 0D0, 8, GRID1, 4, GRID2, .TRUE.,
     *   -1, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS,
     *   IERR)
d379 3
a381 6
      CALL PGCRVL (AXEN, 'AA', GCODE, 0D0, 4, GRID1, 2, GRID2, .TRUE.,
     *   2121, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS,
     *   IERR)

      write (6, 10) ((cross(j,k), j=1,4), k=1,ic)
      write (6, 20) cross(1,ic+1), cross(2,ic+1)
d392 3
a394 4
      IC = 0
      CALL PGCRVL (AXEN, 'AA', GCODE, 0D0, 2, GRID1, 1, GRID2, .TRUE.,
     *   -1, PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS,
     *   IERR)
@


1.3
log
@*** empty log message ***
@
text
@d54 1
a54 1
      NLDPRM(1) = 256D0
d58 1
a58 1
      NLDPRM(3) = 180D0
d62 1
a62 1
      NLDPRM(5) = 360D0/AXEN(1)
d68 1
d80 1
a80 1
     *   LNGVEL, 1, 1, 7, ' ', 0D0, NLDPRM, 256, IC, CROSS, IERR)
d86 1
a86 1
     *   LNGVEL, 1, 1, 7, ' ', 0D0, NLDPRM, 256, IC, CROSS, IERR)
d89 3
a91 1
 10   format (1p,g12.5,0p,f5.1,2f12.5)
d118 1
d133 1
d136 2
a137 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS, IERR)
d140 1
d142 1
a142 1
*     Native grid in green (defines boundary).
d153 1
d155 2
a156 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS, IERR)
d180 1
d186 1
a186 1
      CRVAL(1) = -30D0
d197 2
a198 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS, IERR)
d208 2
a209 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS, IERR)
d219 2
a220 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS, IERR)
d237 2
a238 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS, IERR)
d253 1
d256 2
a257 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS, IERR)
d268 2
a269 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS, IERR)
d272 1
d274 1
a274 1
*     Native grid in green (defines boundary).
d283 1
d285 2
a286 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, 256, IC, CROSS, IERR)
@


1.2
log
@*** empty log message ***
@
text
@d3 4
a6 4
      INTEGER   AXEN(2), GCODE(2), IERR, WCS(4)
      DOUBLE PRECISION CDELT(2), CRPIX(2), CRVAL(2), GRID1(0:10),
     *          GRID2(0:10), LATPOL, LIN(19), LNGPOL, NLDPRM(50),
     *          PC(2,2), PROJP(0:9)
d22 6
d77 1
d79 1
a79 1
     *   LNGVEL, 1, 1, 7, ' ', 0D0, NLDPRM, IERR)
d85 4
a88 1
     *   LNGVEL, 1, 1, 7, ' ', 0D0, NLDPRM, IERR)
d129 1
d131 1
a131 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, IERR)
d133 2
d147 1
a147 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, IERR)
d185 1
d187 1
a187 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, IERR)
d195 1
d197 1
a197 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, IERR)
d207 1
a207 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, IERR)
d224 1
a224 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, IERR)
d241 1
a241 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, IERR)
d252 3
a254 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, IERR)
d266 1
a266 1
     *   PGWCSL, 16, 4, 50, CTYPE, WCS, NLDPRM, IERR)
@


1.1
log
@Initial revision
@
text
@d3 1
a3 1
      INTEGER   AXEN(2), IERR, WCS(4)
d69 3
a71 1
      CALL PGCRVL (AXEN, 'A ', 0, 1, GRID1, 1, GRID2, .FALSE., 
d75 3
a77 1
      CALL PGCRVL (AXEN, 'A ', 2, 0, 0D0, 0, 0D0, .FALSE.,
d109 3
d119 1
a119 1
      CALL PGCRVL (AXEN, 'AA', 0, 0, 0D0, 0, 0D0, .TRUE.,
d133 1
a133 1
      CALL PGCRVL (AXEN, 'AA', 0, 2, GRID1, 2, GRID2, .TRUE.,
d170 3
a172 1
      CALL PGCRVL (AXEN, 'AA', 2, 0, 5D0, 0, 5D0, .TRUE.,
d179 3
a181 1
      CALL PGCRVL (AXEN, 'AA', 0, 1, GRID1, 1, GRID2, .TRUE.,
d191 1
a191 1
      CALL PGCRVL (AXEN, 'AA', 0, 3, GRID1, 2, GRID2, .TRUE.,
d208 1
a208 1
      CALL PGCRVL (AXEN, 'AA', 0, 8, GRID1, 4, GRID2, .TRUE.,
d225 1
a225 1
      CALL PGCRVL (AXEN, 'AA', 0, 8, GRID1, 4, GRID2, .TRUE.,
d236 1
a236 1
      CALL PGCRVL (AXEN, 'AA', 0, 4, GRID1, 2, GRID2, .TRUE.,
d248 1
a248 1
      CALL PGCRVL (AXEN, 'AA', 0, 2, GRID1, 1, GRID2, .TRUE.,
@
