pro mk_square2, dummy, qdebug=qdebug
;
;uses normalized coordinates
;
;
data_dx = abs(!cxmax - !cxmin)			;in data coordinates
data_dy = abs(!cymax - !cymin)
sc_dx = !d.x_size
sc_dy = !d.y_size
;
x_ppdu = sc_dx/data_dx		;pixels per data unit
y_ppdu = sc_dy/data_dy
if (x_ppdu gt y_ppdu) then begin
    ;need to shrink x axis
    !p.position(0) = 0.15
    !p.position(1) = 0.15
    !p.position(2) = !p.position(0)+0.75*y_ppdu/x_ppdu	;data_dx/data_dy * (sc_dy/sc_dx)
    !p.position(3) = 0.90
end else begin
    !p.position(0) = 0.15
    !p.position(1) = 0.15
    !p.position(2) = 0.90
    !p.position(3) = !p.position(1)+0.75*x_ppdu/y_ppdu	;data_dy/data_dx * (sc_dx/sc_dy)
end
;
if (keyword_set(qdebug)) then begin
    print, 'data_dx,dy     ', data_dx, data_dy
    print, 'sc_dx,dy       ', sc_dx, sc_dy
    print, 'x_ppdu, y_ppdu ', x_ppdu, y_ppdu
    print, '!p.region      ', !p.region
    print, '!p.position    ', !p.position
    print, '!d.x,y_size    ', !d.x_size, !d.y_size
    print, '!x,y.window    ', !x.window, !y.window
end
;stop
end