docs: def other_profiles(p,ifig=1,xlm=xlm,show=False,xaxis=xaxis_type, figsize2=(10,8)):
      '''Four panels of other profile plots
     
     Parameters
      ----------
  
      p : instance
        mesa_profile instance
  
      xll : tuple
        xlimits: mass_min, mass_max
    
  
      show : Boolean
        False  for batch use
        True   for interactive use
      '''
  
      matplotlib.rc('figure',facecolor='white',figsize=figsize2)
  
      mass = p.get('mass')                      # in units of Msun
      radius = p.get('radius')*ast.rsun_cm/1.e8  # in units of Mm
      if xaxis is "Eulerian":
          xxx = radius
          if xlm[0]==0 and xlm[1] == 0:
              indtop = 0
              indbot = len(mass)-1
          else: 
              indbot = np.where(radius>=xlm[0])[0][-1]
              indtop = np.where(radius<xlm[1])[0][0]
          xll = (radius[indbot],radius[indtop])
          xxlabel = "radius (Mm)"
      elif xaxis is "Lagrangian": 
          xxx = mass
          xll = xlm
          xxlabel = "$M / \mathrm{M_{sun}}$"
      else:
          print("Error: don't understand xaxis choice, must be Lagrangian or Eulerian")
  
  
      # create subplot structure
      t, ([ax1,ax2],[ax3, ax4],[ax5,ax6]) = matplotlib.pyplot.subplots(3, 2, sharex=True, sharey=False)
      
  # panel 1: burns: pp, cno, burn_c
  # panel 2: convection and mixing: entropy, Tgrad
  
  
      # which burning to show
      Enuc = ['pp','cno','tri_alfa','burn_c','burn_o','burn_n','burn_si','burn_mg','burn_na','burn_ne','eps_nuc']
      ax = ax1
      for thing in Enuc:
          ind = Enuc.index(thing)
          ax.plot(xxx, np.log10(p.get(thing)), ls=u.linestylecb(ind,a,b)[0],\
               marker=u.linestylecb(ind,a,b)[1], color=u.linestylecb(ind,a,b)[2],\
               markevery=50,label=thing)
          
      # set x and y lims and labels
      #ax.set_title('Nuclear Energy Production')
      ax.set_ylim(0,15)
      ax.set_xlim(xll)
      ax.legend(loc=1, ncol=2, fontsize='small')
      #ax.set_xlabel(xxlabel)
      ax.set_ylabel('$ \log \epsilon $')
  #--------------------------------------------------------------------------------------------#
  
      # gradients
  
      mix = [['gradr']]
      mix1 = [['grada']]
      
      for i in range(1):
          for thing in mix[i]:
              ind = mix[i].index(thing)
      for i in range(1):
          for thing1 in mix1[i]:
              ind1 = mix1[i].index(thing1)
              ax2.plot(xxx, (np.tanh(np.log10(p.get(thing))-np.log10(p.get(thing1))))\
              ,ls=u.linestylecb(ind,a,b)[0],\
               marker=u.linestylecb(ind,a,b)[1], color=u.linestylecb(ind,a,b)[2],\
               markevery=50,label=thing)
      # set x and y lims and labels
      ax2.axhline(ls='dashed',color='black',label="")
      #ax2.set_title('Mixing Regions')
      ax2.yaxis.tick_right()
      ax2.yaxis.set_label_position("right")
      ax2.set_ylim(-.1,.1)
      ax2.set_xlim(xll)
      ax2.legend(labels='Mixing',loc=1)
      #ax2.set_xlabel(xxlabel)
      ax2.set_ylabel('$\\tanh(\\log(\\frac{\\nabla_{rad}}{\\nabla_{ad}}))$')
  
  #--------------------------------------------------------------------------------------------#
  
      # entropy
      S = ['entropy']
  
      ax = ax5
      
      for thing in S:
          ind = 2
          ax.plot(xxx, p.get(thing), ls=u.linestylecb(ind,a,b)[0],\
           marker=u.linestylecb(ind,a,b)[1], color=u.linestylecb(ind,a,b)[2],\
           markevery=50,label=thing)
      # set x and y lims and labels
      #ax.set_title('Specific Entropy (/A*kerg)')
      ax.set_ylim(0,50)
      ax.set_xlim(xll)
      ax.legend(loc=1)
      ax.set_xlabel(xxlabel)
      ax.set_ylabel(' Specific Entropy')
  
  #--------------------------------------------------------------------------------------------#
      # rho, mu, T
      S = ['logRho','mu','temperature']
      T8 = [False,False,True]
      ax = ax6
  
      for thing in S:
          ind = S.index(thing)
          thisy = p.get(thing)/1.e8 if T8[ind] else p.get(thing)
          ax.plot(xxx, thisy, ls=u.linestylecb(ind,a,b)[0],\
           marker=u.linestylecb(ind,a,b)[1], color=u.linestylecb(ind,a,b)[2],\
           markevery=50,label=thing)
      # set x and y lims and labels                                                              
      #ax.set_title('Rho, mu, T')
      ax.set_ylim(0.,9.)
      ax.set_xlim(xll)
      ax.legend(loc=0)
      ax.set_xlabel(xxlabel)
      ax.set_ylabel('log Rho, mu, T8')
  
  
  #--------------------------------------------------------------------------------------------#
  
      
      # gas pressure fraction and opacity
      S = ['pgas_div_ptotal']
      o = ['log_opacity']
  
      ax = ax4
      axo = ax.twinx()
      
      for thing in S:
          ind = 5
          ax.plot(xxx, p.get(thing), ls=u.linestylecb(ind,a,b)[0],\
           marker=u.linestylecb(ind,a,b)[1], color=u.linestylecb(ind,a,b)[2],\
           markevery=50,label=thing)
  
      for thing in o:
          ind = 3
          axo.plot(xxx, p.get(thing), ls=u.linestylecb(ind,a,b)[0],\
           marker=u.linestylecb(ind,a,b)[1], color=u.linestylecb(ind,a,b)[2],\
           markevery=50,label=thing)
  
     
      # set x and y lims and labels
     # ax.set_title('Pgas fraction + opacity')
     # ax.set_ylim(0,60)
      ax.set_xlim(xll)
      axo.set_xlim(xll)
      ax.legend(loc=0)
      axo.legend(loc=(.15,.85))
      #ax.set_xlabel(xxlabel)
      ax.set_ylabel('$\mathrm{ P_{gas} / P_{tot}}$')
      axo.set_ylabel('$ log(Opacity)$')
  
  #--------------------------------------------------------------------------------------------#
  
      # Diffusion coefficient
      gT = ['log_D_mix','conv_vel_div_csound']
      logy = [False,True]
      
      ax = ax3
      ind = 0
      for thing in gT:
          ind = gT.index(thing)
          thisx = np.log(p.get(thing))+16 if logy[ind] else p.get(thing)
          ax.plot(xxx, thisx, ls=u.linestylecb(ind,a,b)[0],\
           marker=u.linestylecb(ind,a,b)[1], color=u.linestylecb(ind,a,b)[2],\
           markevery=50,label=thing)
  # set x and y lims and labels
      ax.axhline(16,ls='dashed',color='black',label="$\mathrm{Ma}=0$")
     # ax.set_title('Mixing')
      ax.set_ylim(10,17)
      ax.set_xlim(xll)
      ax.legend(loc=0)
     # ax.set_xlabel(xxlabel)
      ax.set_ylabel('$\\log D / [cgs] \\log v_{\mathrm{conv}}/c_s + 16 $ ')
  
      title_str = "Other profiles: "+'t ='+str(title_format%p.header_attr['star_age'])\
                +', dt ='+str(title_format%p.header_attr['time_step'])\
                +', model number ='+str(int(p.header_attr['model_number']))
      t.suptitle(title_str, fontsize=12)
     # t.tight_layout()
      t.subplots_adjust(left=0.1, bottom=0.1, right=0.9, top=0.9, wspace=0.15, hspace=0.1)
      t.savefig('other'+str(int(p.header_attr['model_number'])).zfill(6)+'.png')