define oplot_common_HHe (xmin, xmax, ymin, ymax, ilabel, redshift)
{
   % oplot the  H-like and He-triplet lines for elements
   
   variable lams, ilam, lamclrs, elstr, thisval;
   %        H       He-R   He-I    He-F  
   lams = [1.530 , 1.588, 1.595,  1.604,     % Ni
	   1.7781, 1.851, 1.857,  1.868,     % Fe
	   3.018,  3.177, 3.191,  3.211,     % Ca
	   3.731,  3.949, 3.9675, 3.994,     % Ar
	   4.7274, 5.0387, 5.0665, 5.1015,     % S
	   6.180,  6.6477, 6.6879, 6.740,      % Si
	   8.4190, 9.1685, 9.2310, 9.314,     % Mg
	   12.132, 13.447, 13.553, 13.700,     % Ne
	   18.967, 21.602, 21.804, 22.101,     % O
	   24.779, 28.787, 29.082,  29.534,    % N
   	   33.737, 40.268]   ;                  % C

   % Adjust for the redshift:
   lams = lams*(1.0+redshift);

   lamclrs = [10, 10, 10, 10,                        % Ni 
	      5, 5, 5, 5,                        % Fe
	      15, 15, 15, 15,                    % Ca
	      12, 12, 12, 12,                          % Ar
	      8, 8, 8, 8,                           % S
	      14, 14, 14, 14,                    % Si
	      10, 10, 10, 10,                        % Mg
	      13, 13, 13, 13,                          % Ne
	      5, 5, 5, 5,                         % O
	      15, 15, 15, 15,                       % N
   	      1, 1];                                  % C

   elstr = ["Ni", "Ni", "Ni", "Ni",
	    "Fe", "Fe", "Fe", "Fe", 
	      "Ca", "Ca", "Ca", "Ca",  
	      "Ar", "Ar", "Ar", "Ar", 
	      "S", "S", "S", "S", 
	      "Si", "Si", "Si", "Si", 
	      "Mg", "Mg", "Mg", "Mg", 
	      "Ne", "Ne", "Ne", "Ne", 
	      "O", "O", "O", "O",
	      "N", "N", "N", "N",
	      "C", "C"] ;

   _for ilam (0,length(lams)-1,1)
     {
	thisval = Const_hc/lams[ilam];
	if (thisval > xmin and thisval < xmax)
	  {
	     oplot(thisval*[1.0,1.0], [ymin,ymax], lamclrs[ilam] );
	  }
     }


   % Add the element string
   if (ilabel == 1)
     {	
	_for ilam (0,length(lams)-1,4)
	  {
	     color(lamclrs[ilam]);	  
	     thisval = Const_hc*2.0/(lams[ilam]+lams[ilam+1]);
	     if (thisval > xmin and thisval < xmax)
	       {
		  xylabel(thisval, ymax + 0.1*(ymax-ymin), elstr[ilam], 0.0, 0.5);
	       }
	  }
     }
   
}