#####	Master file for producing combined images of the front and back side of the sun   #####
#####	Original Author:	Kenny Oslund
#####	Current Author:		Phil Scherrer
#####	Date:	30 Jan 2005
#
#	the program to use this file and produce the combined images may be executed with the 
#		following command line:
#	make_combined_images -front fd_M_96m_01d.2996.0008.fits -back FSI_lowres.2001.03.16_12h.fits -mklist sample_make -v
# pound signs may be used to denote comments


### the first set of instructions makes the whole sun synoptic chart 
### with an aspect ratio of 2.5:1

# first process the front side image from a full resolution magnetogram
!-plate_reproject
-front_fits front.cyleqa.fits -project cyleqa -x_range 4 -naxis1 318 -naxis2 200
-y_range 2 -map_b0 0 -front -gauss_size .025 -los_fix 0.5
# the parameter map_l0 is not specified, so map_l0 will be the same as obs_l0 (by default),
# and obs_l0 will be in the center of the image.
# x_range is in radians, so 4 radians / 318 pixels = PI radians / 250 pixels, which is a 
# 2.5:2 aspect ratio for one hemisphere.

# the open bracket '{' saves the result from that step (by putting it onto
# a stack) so that it will be used in the next step
{
	# next reproject the backside image
	!-reproject
	-back_fits back.cylequ.fits -project cyleqa -x_range 4 -naxis1 318 
	-y_range 2 -map_b0 0 -back
	# map_l0 is not specified, so it is not changed
	{
		# apply a color table to the back and then the front. 
		# different color tables with different bounds will be used.
		!-apply_ct
		-ct both.ct -max -0.40 -min -0.15 -back
		{
			!-apply_ct
			-ct both.ct -max 180 -min 10 -front 
			{ # save clean back and front pair on stack
			    # Add Earthside label to front side image
			    !-label_rgb
			    -weight 0.01 -color 0,0,0 -limb
			    -string "Earthside" -label_y0 188 -label_x0 135 -front
			    # -front_ppm front.cyleqa_lab.ppm
			    { # save Earthside labelled pair
			        # Add Farside label to back side image
			        !-label_rgb
			        -weight 0.01 -color 0,0,0 -limb
			        -string "Farside" -label_y0 188 -label_x0 135 -back
			        # -back_ppm back.cyleqa_lab.ppm
				{ # save Farside and Earthside labelled pair
				    # combine the images into a synoptic chart. 0 deg carrington longitude
				    # is at the left hand side, so 180 deg is at the center of the image,
				    # and is passed (in radians) through the argument map_l0
				    !-combine
				    -naxis1 500 -naxis2 200 -wrap
				    -x_range 6.283185307 -y_range 2 -map_l0 3.14159
				        { # save combined with Earthside labels
				        # now apply some more labels, first grey lat-lon lines
				        !-label_rgb
				        -weight 0.01 -color 200,200,200
				        -lat_space 0.5236 -lon_space 1.0472
					        { # save combined with lat-lon
					        # now add date of observation and save as synoptic3.ppm
					        !-label_rgb
					        -front_ppm synoptic3.ppm -weight 0.01 -color 0,0,0 -limb -date
					        -label_y0 4 -date_format "Date Obs: %.13s:00" -front
					        } # pop date obs, stacktop now with lat-lon
				        } # pop lat-lon, stacktop now combined image
			        } # pop combined, stacktop now pair with Farside
			    } # Farside labelled pair, stacktop now pair with Earthside, will be overwritten
			    # recombine without labels save as synoptic1.ppm
			    !-combine
			    -naxis1 500 -naxis2 200 -file_ppm synoptic1.ppm -wrap
			    -x_range 6.283185307 -y_range 2 -map_l0 3.14159
			} # pop combined, stacktop now clean pair
		# 
		}
# Make a second combined image with labels to test color tables
		# apply a color table to the back and then the front. 
		# different color tables with different bounds will be used.
		!-apply_ct
		-ct linear.ct -min -0.40 -max 0.20 -back
		{
			!-apply_ct
			-ct linear.ct -max 0 -min 200 -front 
			{ # save clean back and front pair on stack
			    # Add Earthside label to front side image
			    !-label_rgb
			    -weight 0.01 -color 0,0,0 -limb
			    -string "Earthside-x" -label_y0 188 -label_x0 135 -front
			    # -front_ppm front.cyleqa_lab.ppm
			    { # save Earthside labelled pair
			        # Add Farside label to back side image
			        !-label_rgb
			        -weight 0.01 -color 0,0,0 -limb
			        -string "Farside-x" -label_y0 188 -label_x0 135 -back
			        # -back_ppm back.cyleqa_lab.ppm
				{ # save Farside and Earthside labelled pair
				    # combine the images into a synoptic chart. 0 deg carrington longitude
				    # is at the left hand side, so 180 deg is at the center of the image,
				    # and is passed (in radians) through the argument map_l0
				    !-combine
				    -naxis1 500 -naxis2 200 -wrap
				    -x_range 6.283185307 -y_range 2 -map_l0 3.14159
				        { # save combined with Earthside labels
				        # now apply some more labels, first grey lat-lon lines
				        !-label_rgb
				        -weight 0.01 -color 200,200,200
				        -lat_space 0.5236 -lon_space 1.0472
					        { # save combined with lat-lon
					        # now add date of observation and save as synoptic3.ppm
					        !-label_rgb
					        -front_ppm synoptic5.ppm -weight 0.01 -color 0,0,0 -limb -date
					        -label_y0 4 -date_format "Date Obs: %.13s:00" -front
					        } # pop date obs, stacktop now with lat-lon
				        } # pop lat-lon, stacktop now combined image
			        } # pop combined, stacktop now pair with Farside
			    } # Farside labelled pair, stacktop now pair with Earthside, will be overwritten
			} # pop combined, stacktop now clean pair
		# now remove rest of the processed images from the stack to make polar views from the input images.
		}
	}
}

### the second set of instructions makes the north polar view

!-plate_reproject
-front_fits front.postel_N.fits -project postels -x_range 3.14159 -y_range 3.14159 
-naxis1 250 -naxis2 250 -map_b0 1.5707 -map_l0 0 -front -gauss_size .03  -los_fix 0.5
{

	!-reproject
	-back_fits back.postel_N.fits -project postels -x_range 3.14159 -y_range 3.14159 
	-naxis1 250 -naxis2 250 -map_b0 1.5707 -map_l0 0 -back
	{
		!-apply_ct
		-ct back_sm.ct -max -0.40 -min -0.15 -back
		{
			!-apply_ct
			-ct both.ct -max 180 -min 10 -front 
			{
				!-combine
				-naxis1 250 -naxis2 250 -file_ppm polar1_N.ppm -wrap
				-x_range 3.14159 -y_range 3.14159 -map_b0 1.5707
				{
					!-label_rgb
					-front_ppm polar2_N.ppm -weight 0.01 -string "North Pole" 
					-lat_space 0.78535 -lon_space 1.047
					
					# because there was no {, polar2.ppm was not saved and polar3.ppm
					# will not have the same labels that polar2.ppm does
					!-label_rgb
					-front_ppm polar3_N.ppm -weight 0.01 -limb
				}
			}
		}
	}
}


### finally, the set of instructions makes the south polar view

!-plate_reproject
-front_fits front.postel_S.fits -project postels -x_range 3.14159 -y_range 3.14159 
-naxis1 250 -naxis2 250 -map_b0 -1.5707 -map_l0 0 -front -gauss_size .03  -los_fix 0.5
{

	!-reproject
	-back_fits back.postel_S.fits -project postels -x_range 3.14159 -y_range 3.14159 
	-naxis1 250 -naxis2 250 -map_b0 -1.5707 -map_l0 0 -back
	{
		!-apply_ct
		-ct back_sm.ct -max 0.40 -min 0.15 -back
		{
			!-apply_ct
			-ct both.ct -max 180 -min 10 -front 
			{
				!-combine
				-naxis1 250 -naxis2 250 -file_ppm polar1_S.ppm -wrap
				-x_range 3.14159 -y_range 3.14159 -map_b0 -1.5707
				{
					!-label_rgb
					-front_ppm polar2_S.ppm -weight 0.01 -string "South Pole" 
					-lat_space 0.78535 -lon_space 1.047
					
					# because there was no {, polar2.ppm was not saved and polar3.ppm
					# will not have the same labels that polar2.ppm does
					!-label_rgb
					-front_ppm polar3_S.ppm -weight 0.01 -limb
				}
			}
		}
	}
}