nucleardatapy 0.2.1__py3-none-any.whl → 1.0.0__py3-none-any.whl

nucleardatapy/fig/eos_setupAM_e_asy_lep_fig.py

@@ -1,125 +0,0 @@
-import numpy as np
-import matplotlib.pyplot as plt
-
-import nucleardatapy as nuda
-
-def eos_setupAM_e_asy_lep_fig( pname, micro_mbs, pheno_models, asy, band ):
-    """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x2 with:\
-    [0]: nuclear chart.
-
-    :param pname: name of the figure (*.png)
-    :type pname: str.
-    :param table: table.
-    :type table: str.
-    :param version: version of table to run on.
-    :type version: str.
-    :param theo_tables: object instantiated on the reference band.
-    :type theo_tables: object.
-
-    """
-    #
-    print(f'Plot name: {pname}')
-    #
-    fig, axs = plt.subplots(1,2)
-    #fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.10, bottom=0.12, right=None, top=0.90, wspace=0.05, hspace=0.3 )
-    #
-    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    axs[0].set_ylabel(r'$E_\text{lep}/A$')
-    axs[0].set_xlim([0, 0.28])
-    axs[0].set_ylim([-10, 35])
-    #
-    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    #axs[1].set_ylabel(r'$E/A$')
-    axs[1].set_xlim([0, 0.28])
-    axs[1].set_ylim([-10, 35])
-    axs[1].tick_params('y', labelleft=False)
-    #
-    mb_check = []
-    #
-    for kmb,mb in enumerate(micro_mbs):
-        #
-        print('mb:',mb,kmb)
-        #
-        models, models_lower = nuda.matter.micro_esym_models_mb( mb )
-        #
-        print('models:',models)
-        #
-        if mb == 'VAR':
-            models.remove('1998-VAR-AM-APR-fit')
-            models_lower.remove('1998-var-am-apr-fit')
-        #
-        for model in models:
-            #
-            micro = nuda.eos.setupAM( model = model, kind = 'micro', asy = asy )
-            if nuda.env.verb_output: micro.print_outputs( )
-            #
-            check = nuda.matter.setupCheck( eos = micro, band = band )
-            #
-            if check.isInside:
-                lstyle = 'solid'
-            else:
-                lstyle = 'dashed'
-                #continue
-            #
-            if micro.esym is not None: 
-                print('model:',model)
-                print('den:',micro.den)
-                print('e2a_lep:',micro.e2a_lep)
-                if mb in mb_check:
-                    axs[0].plot( micro.den, micro.e2a_lep, marker='o', linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
-                else:
-                    mb_check.append(mb)
-                    axs[0].plot( micro.den, micro.e2a_lep, marker='o', linestyle=lstyle, label=mb, markevery=micro.every, color=nuda.param.col[kmb] )
-            # end of model
-        # end of mb
-    axs[0].text(0.02,-8,'microscopic models',fontsize='10')
-    axs[0].text(0.02,-9.5,r'for $\delta=$'+str(asy),fontsize='10')
-    #
-    model_check = []
-    #
-    for kmodel,model in enumerate(pheno_models):
-        #
-        params, params_lower = nuda.matter.pheno_esym_params( model = model )
-        #
-        for param in params:
-            #
-            pheno = nuda.eos.setupAM( model = model, param = param, kind = 'pheno', asy = asy )
-            if nuda.env.verb_output: pheno.print_outputs( )
-            #
-            check = nuda.matter.setupCheck( eos = pheno, band = band )
-            #
-            if check.isInside:
-                lstyle = 'solid'
-            else:
-                lstyle = 'dashed'
-                #continue
-            #
-            if pheno.e2a_lep is not None: 
-                print('model:',model,' param:',param)
-                if model in model_check:
-                    axs[1].plot( pheno.den, pheno.e2a_lep, linestyle=lstyle, markevery=pheno.every, color=nuda.param.col[kmodel] )
-                else:
-                    model_check.append(model)
-                    axs[1].plot( pheno.den, pheno.e2a_lep, linestyle=lstyle, label=model, markevery=pheno.every, color=nuda.param.col[kmodel] )
-            # end of param
-        # end of model
-    #
-    #axs[1].fill_between( band.den, y1=(band.e2a-band.e2a_std), y2=(band.e2a+band.e2a_std), color=band.color, alpha=band.alpha, visible=True )
-    #axs[1].plot( band.den, (band.e2a-band.e2a_std), color='k', linestyle='dashed' )
-    #axs[1].plot( band.den, (band.e2a+band.e2a_std), color='k', linestyle='dashed' )
-    axs[1].text(0.02,-8,'phenomenological models',fontsize='10')
-    axs[1].text(0.02,-9.5,r'for $\delta=$'+str(asy),fontsize='10')
-    #
-    #axs[1].legend(loc='upper left',fontsize='8', ncol=2)
-    #axs[0,1].legend(loc='upper left',fontsize='xx-small', ncol=2)
-    #axs[1].legend(loc='lower center',bbox_to_anchor=(0.5,1.02),mode='expand',columnspacing=0,fontsize='8', ncol=2,frameon=False)
-    #fig.legend(loc='lower center',bbox_to_anchor=(0.5,1.02),mode='expand',columnspacing=0,fontsize='8', ncol=2,frameon=False)
-    fig.legend(loc='upper left',bbox_to_anchor=(0.15,1.0),columnspacing=2,fontsize='8',ncol=5,frameon=False)
-    #
-    if pname is not None: 
-        plt.savefig(pname, dpi=200)
-        plt.close()
-    #

nucleardatapy/fig/eos_setupAM_e_asy_nuc_fig.py

@@ -1,115 +0,0 @@
-import numpy as np
-import matplotlib.pyplot as plt
-
-import nucleardatapy as nuda
-
-def eos_setupAM_e_asy_nuc_fig( pname, micro_mbs, pheno_models, asy, band ):
-    """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x2 with:\
-    [0]: nuclear chart.
-
-    :param pname: name of the figure (*.png)
-    :type pname: str.
-    :param table: table.
-    :type table: str.
-    :param version: version of table to run on.
-    :type version: str.
-    :param theo_tables: object instantiated on the reference band.
-    :type theo_tables: object.
-
-    """
-    #
-    print(f'Plot name: {pname}')
-    #
-    fig, axs = plt.subplots(1,2)
-    #fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.10, bottom=0.12, right=None, top=0.90, wspace=0.05, hspace=0.3 )
-    #
-    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    axs[0].set_ylabel(r'$E_\text{nuc}/A$')
-    axs[0].set_xlim([0, 0.28])
-    axs[0].set_ylim([-10, 35])
-    #
-    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    #axs[1].set_ylabel(r'$E/A$')
-    axs[1].set_xlim([0, 0.28])
-    axs[1].set_ylim([-10, 35])
-    axs[1].tick_params('y', labelleft=False)
-    #
-    mb_check = []
-    #
-    for kmb,mb in enumerate(micro_mbs):
-        #
-        print('mb:',mb,kmb)
-        #
-        models, models_lower = nuda.matter.micro_esym_models_mb( mb )
-        #
-        print('models:',models)
-        #
-        if mb == 'VAR':
-            models.remove('1998-VAR-AM-APR-fit')
-            models_lower.remove('1998-var-am-apr-fit')
-        #
-        for model in models:
-            #
-            micro = nuda.eos.setupAM( model = model, kind = 'micro', asy = asy )
-            if nuda.env.verb_output: micro.print_outputs( )
-            #
-            check = nuda.matter.setupCheck( eos = micro, band = band )
-            #
-            if check.isInside:
-                lstyle = 'solid'
-            else:
-                lstyle = 'dashed'
-                #continue
-            #
-            if micro.e2a_nuc is not None: 
-                if mb in mb_check:
-                    axs[0].plot( micro.den, micro.e2a_nuc, marker='o', linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
-                else:
-                    mb_check.append(mb)
-                    axs[0].plot( micro.den, micro.e2a_nuc, marker='o', linestyle=lstyle, label=mb, markevery=micro.every, color=nuda.param.col[kmb] )
-            # end of model
-        # end of mb
-    axs[0].text(0.02,-8,'microscopic models',fontsize='10')
-    axs[0].text(0.02,-9.5,r'for $\delta=$'+str(asy),fontsize='10')
-    #
-    model_check = []
-    #
-    for kmodel,model in enumerate(pheno_models):
-        #
-        params, params_lower = nuda.matter.pheno_esym_params( model = model )
-        #
-        for param in params:
-            #
-            pheno = nuda.eos.setupAM( model = model, param = param, kind = 'pheno', asy = asy )
-            if nuda.env.verb_output: pheno.print_outputs( )
-            #
-            check = nuda.matter.setupCheck( eos = pheno, band = band )
-            #
-            if check.isInside:
-                lstyle = 'solid'
-            else:
-                lstyle = 'dashed'
-                #continue
-            #
-            if pheno.e2a_nuc is not None: 
-                print('model:',model,' param:',param)
-                if model in model_check:
-                    axs[1].plot( pheno.den, pheno.e2a_nuc, linestyle=lstyle, markevery=pheno.every, color=nuda.param.col[kmodel] )
-                else:
-                    model_check.append(model)
-                    axs[1].plot( pheno.den, pheno.e2a_nuc, linestyle=lstyle, label=model, markevery=pheno.every, color=nuda.param.col[kmodel] )
-            # end of param
-        # end of model
-    #
-    axs[1].text(0.02,-8,'phenomenological models',fontsize='10')
-    axs[1].text(0.02,-9.5,r'for $\delta=$'+str(asy),fontsize='10')
-    #
-    fig.legend(loc='upper left',bbox_to_anchor=(0.15,1.0),columnspacing=2,fontsize='8',ncol=5,frameon=False)
-    #
-    if pname is not None: 
-        plt.savefig(pname, dpi=200)
-        plt.close()
-    #

nucleardatapy/fig/eos_setupAM_e_asy_tot_fig.py

@@ -1,117 +0,0 @@
-import numpy as np
-import matplotlib.pyplot as plt
-
-import nucleardatapy as nuda
-
-def eos_setupAM_e_asy_tot_fig( pname, micro_mbs, pheno_models, asy, band ):
-    """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x2 with:\
-    [0]: nuclear chart.
-
-    :param pname: name of the figure (*.png)
-    :type pname: str.
-    :param table: table.
-    :type table: str.
-    :param version: version of table to run on.
-    :type version: str.
-    :param theo_tables: object instantiated on the reference band.
-    :type theo_tables: object.
-
-    """
-    #
-    print(f'Plot name: {pname}')
-    #
-    fig, axs = plt.subplots(1,2)
-    fig.subplots_adjust(left=0.10, bottom=0.12, right=None, top=0.90, wspace=0.05, hspace=0.3 )
-    #
-    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    axs[0].set_ylabel(r'$E_\text{tot}/A$')
-    axs[0].set_xlim([0, 0.28])
-    axs[0].set_ylim([-10, 35])
-    #
-    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    #axs[1].set_ylabel(r'$E/A$')
-    axs[1].set_xlim([0, 0.28])
-    axs[1].set_ylim([-10, 35])
-    axs[1].tick_params('y', labelleft=False)
-    #
-    mb_check = []
-    #
-    for kmb,mb in enumerate(micro_mbs):
-        #
-        print('mb:',mb,kmb)
-        #
-        models, models_lower = nuda.matter.micro_esym_models_mb( mb )
-        #
-        print('models:',models)
-        #
-        if mb == 'VAR':
-            models.remove('1998-VAR-AM-APR-fit')
-            models_lower.remove('1998-var-am-apr-fit')
-        #
-        for model in models:
-            #
-            micro = nuda.eos.setupAM( model = model, kind = 'micro', asy = asy )
-            if nuda.env.verb_output: micro.print_outputs( )
-            #
-            check = nuda.matter.setupCheck( eos = micro, band = band )
-            #
-            if check.isInside:
-                lstyle = 'solid'
-            else:
-                lstyle = 'dashed'
-                #continue
-            #
-            if micro.esym is not None: 
-                print('model:',model)
-                if mb in mb_check:
-                    axs[0].plot( micro.den, micro.e2a_tot, marker='o', linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
-                else:
-                    mb_check.append(mb)
-                    axs[0].plot( micro.den, micro.e2a_tot, marker='o', linestyle=lstyle, label=mb, markevery=micro.every, color=nuda.param.col[kmb] )
-            # end of model
-        # end of mb
-    #
-    axs[0].text(0.02,-8,'microscopic models',fontsize='10')
-    axs[0].text(0.02,-9.5,r'for $\delta=$'+str(asy),fontsize='10')
-    #
-    model_check = []
-    #
-    for kmodel,model in enumerate(pheno_models):
-        #
-        params, params_lower = nuda.matter.pheno_esym_params( model = model )
-        #
-        for param in params:
-            #
-            pheno = nuda.eos.setupAM( model = model, param = param, kind = 'pheno', asy = asy )
-            if nuda.env.verb_output: pheno.print_outputs( )
-            #
-            check = nuda.matter.setupCheck( eos = pheno, band = band )
-            #
-            if check.isInside:
-                lstyle = 'solid'
-            else:
-                lstyle = 'dashed'
-                #continue
-            #
-            if pheno.e2a_tot is not None: 
-                print('model:',model,' param:',param)
-                #beta.label=None
-                if model in model_check:
-                    axs[1].plot( pheno.den, pheno.e2a_tot, linestyle=lstyle, markevery=pheno.every, color=nuda.param.col[kmodel] )
-                else:
-                    model_check.append(model)
-                    axs[1].plot( pheno.den, pheno.e2a_tot, linestyle=lstyle, label=model, markevery=pheno.every, color=nuda.param.col[kmodel] )
-            # end of param
-        # end of model
-    #
-    axs[1].text(0.02,-8,'phenomenological models',fontsize='10')
-    axs[1].text(0.02,-9.5,r'for $\delta=$'+str(asy),fontsize='10')
-    #
-    fig.legend(loc='upper left',bbox_to_anchor=(0.15,1.0),columnspacing=2,fontsize='8',ncol=5,frameon=False)
-    #
-    if pname is not None: 
-        plt.savefig(pname, dpi=200)
-        plt.close()
-    #

nucleardatapy/fig/eos_setupAM_e_fig.py

@@ -1,173 +0,0 @@
-import numpy as np
-import matplotlib.pyplot as plt
-
-import nucleardatapy as nuda
-
-def eos_setupAM_e_fig( pname, micro_mbs, pheno_models, band ):
-    """
-    Plot nuclear chart (N versus Z).\
-    The plot is 1x2 with:\
-    [0]: nuclear chart.
-
-    :param pname: name of the figure (*.png)
-    :type pname: str.
-    :param table: table.
-    :type table: str.
-    :param version: version of table to run on.
-    :type version: str.
-    :param theo_tables: object instantiated on the reference band.
-    :type theo_tables: object.
-
-    """
-    #
-    print(f'Plot name: {pname}')
-    #
-    fig, axs = plt.subplots(3,2)
-    fig.subplots_adjust(left=0.10, bottom=0.12, right=None, top=0.9, wspace=0.05, hspace=0.05 )
-    #
-    #axs[0,0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    axs[0,0].set_ylabel(r'$E_\text{lep}/A$')
-    axs[0,0].set_xlim([0, 0.28])
-    axs[0,0].set_ylim([-2, 38])
-    axs[0,0].tick_params('x', labelbottom=False)
-    #
-    #axs[0,1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    axs[0,1].set_xlim([0, 0.28])
-    axs[0,1].set_ylim([-2, 38])
-    axs[0,1].tick_params('y', labelleft=False)
-    axs[0,1].tick_params('x', labelbottom=False)
-    #
-    #axs[1,0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    axs[1,0].set_ylabel(r'$E_\text{nuc}/A$')
-    axs[1,0].set_xlim([0, 0.28])
-    axs[1,0].set_ylim([-10, 30])
-    axs[1,0].tick_params('x', labelbottom=False)
-    #
-    #axs[1,1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    axs[1,1].set_xlim([0, 0.28])
-    axs[1,1].set_ylim([-10, 30])
-    axs[1,1].tick_params('y', labelleft=False)
-    axs[1,1].tick_params('x', labelbottom=False)
-    #
-    axs[2,0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    axs[2,0].set_ylabel(r'$E_\text{tot}/A$')
-    axs[2,0].set_xlim([0, 0.28])
-    axs[2,0].set_ylim([-2, 38])
-    #
-    axs[2,1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    axs[2,1].set_xlim([0, 0.28])
-    axs[2,1].set_ylim([-2, 38])
-    axs[2,1].tick_params('y', labelleft=False)
-    #
-    # fix the asymmetry parameters
-    #
-    asys = [ 0.6, 0.8 ]
-    #
-    mb_check = []
-    model_check = []
-    #
-    for asy in asys:
-        #
-        print('asy:',asy)
-        #
-        for kmb,mb in enumerate(micro_mbs):
-            #
-            print('mb:',mb,kmb)
-            #
-            models, models_lower = nuda.matter.micro_esym_models_mb( mb )
-            #models, models_lower = nuda.matter.micro_models_mb( mb )
-            #
-            print('models:',models)
-            #
-            if mb == 'VAR':
-                models.remove('1998-VAR-AM-APR-fit')
-                models_lower.remove('1998-var-am-apr-fit')
-            #
-            for model in models:
-                #
-                micro = nuda.eos.setupAM( model = model, kind = 'micro', asy = asy )
-                if nuda.env.verb_output: micro.print_outputs( )
-                #
-                check = nuda.matter.setupCheck( eos = micro, band = band )
-                #
-                if check.isInside:
-                    lstyle = 'solid'
-                else:
-                    lstyle = 'dashed'
-                    continue
-                #
-                if micro.e2a_lep is not None: 
-                    if mb in mb_check:
-                        print('model:',model)
-                        print('den:',micro.den)
-                        print('e2a_lep:',micro.e2a_lep)
-                        axs[0,0].plot( micro.den, micro.e2a_lep, marker='o', linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
-                        axs[1,0].plot( micro.den, micro.e2a_nuc, marker='o', linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
-                        axs[2,0].plot( micro.den, micro.e2a_tot, marker='o', linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
-                    else:
-                        mb_check.append(mb)
-                        print('mb:',mb)
-                        print('model:',model)
-                        print('den:',micro.den)
-                        print('e2a_lep:',micro.e2a_lep)
-                        axs[0,0].plot( micro.den, micro.e2a_lep, marker='o', linestyle=lstyle, label=mb, markevery=micro.every, color=nuda.param.col[kmb] )
-                        axs[1,0].plot( micro.den, micro.e2a_nuc, marker='o', linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
-                        axs[2,0].plot( micro.den, micro.e2a_tot, marker='o', linestyle=lstyle, markevery=micro.every, color=nuda.param.col[kmb] )
-                # end of model
-            # end of mb
-        #
-        for kmodel,model in enumerate(pheno_models):
-            #
-            params, params_lower = nuda.matter.pheno_esym_params( model = model )
-            #
-            for param in params:
-                #
-                pheno = nuda.eos.setupAM( model = model, param = param, kind = 'pheno', asy = asy )
-                if nuda.env.verb_output: pheno.print_outputs( )
-                #
-                check = nuda.matter.setupCheck( eos = pheno, band = band )
-                #
-                if check.isInside:
-                    lstyle = 'solid'
-                else:
-                    lstyle = 'dashed'
-                    continue
-                #
-                if pheno.e2a_lep is not None: 
-                    print('model:',model,' param:',param)
-                    if model in model_check:
-                        axs[0,1].plot( pheno.den, pheno.e2a_lep, linestyle=lstyle, markevery=pheno.every, color=nuda.param.col[kmodel] )
-                        axs[1,1].plot( pheno.den, pheno.e2a_nuc, linestyle=lstyle, markevery=pheno.every, color=nuda.param.col[kmodel] )
-                        axs[2,1].plot( pheno.den, pheno.e2a_tot, linestyle=lstyle, markevery=pheno.every, color=nuda.param.col[kmodel] )
-                    else:
-                        model_check.append(model)
-                        axs[0,1].plot( pheno.den, pheno.e2a_lep, linestyle=lstyle, label=model, markevery=pheno.every, color=nuda.param.col[kmodel] )
-                        axs[1,1].plot( pheno.den, pheno.e2a_nuc, linestyle=lstyle, markevery=pheno.every, color=nuda.param.col[kmodel] )
-                        axs[2,1].plot( pheno.den, pheno.e2a_tot, linestyle=lstyle, markevery=pheno.every, color=nuda.param.col[kmodel] )
-                # end of param
-            # end of model
-    #
-    axs[0,0].text(0.02,0,'microscopic models',fontsize='10')
-    axs[0,1].text(0.02,0,'phenomenological models',fontsize='10')
-    #
-    axs[0,0].text(0.1,30,r'$\delta=0.6$',fontsize='10')
-    axs[0,1].text(0.1,30,r'$\delta=0.6$',fontsize='10')
-    axs[0,0].text(0.1,13,r'$\delta=0.8$',fontsize='10')
-    axs[0,1].text(0.1,13,r'$\delta=0.8$',fontsize='10')
-    #
-    axs[1,0].text(0.1,-2,r'$\delta=0.6$',fontsize='10')
-    axs[1,1].text(0.1,-2,r'$\delta=0.6$',fontsize='10')
-    axs[1,0].text(0.1,7,r'$\delta=0.8$',fontsize='10')
-    axs[1,1].text(0.1,7,r'$\delta=0.8$',fontsize='10')
-    #
-    axs[2,0].text(0.1,27,r'$\delta=0.6$',fontsize='10')
-    axs[2,1].text(0.1,27,r'$\delta=0.6$',fontsize='10')
-    axs[2,0].text(0.1,15,r'$\delta=0.8$',fontsize='10')
-    axs[2,1].text(0.1,15,r'$\delta=0.8$',fontsize='10')
-    #
-    fig.legend(loc='upper left',bbox_to_anchor=(0.15,1.0),columnspacing=2,fontsize='8',ncol=5,frameon=False)
-    #
-    if pname is not None: 
-        plt.savefig(pname, dpi=200)
-        plt.close()
-    #

nucleardatapy/fig/matter_ENM_fig.py

@@ -1,128 +0,0 @@
-import numpy as np
-import matplotlib.pyplot as plt
-
-import nucleardatapy as nuda
-
-def matter_ENM_fig( pname, micro_mbs, pheno_models, band ):
-    """
-    Plot nucleonic energy per particle E/A in NM.\
-    The plot is 1x2 with:\
-    [0,0]: E/A versus den (micro). [0,1]: E/A versus den (pheno).\
-
-    :param pname: name of the figure (*.png)
-    :type pname: str.
-    :param micro_mbs: many-body (mb) approach considered.
-    :type micro_mbs: str.
-    :param pheno_models: models to run on.
-    :type pheno_models: array of str.
-    :param band: object instantiated on the reference band.
-    :type band: object.
-
-    """
-    #
-    print(f'Plot name: {pname}')
-    #
-    fig, axs = plt.subplots(1,2)
-    fig.subplots_adjust(left=0.10, bottom=0.12, right=None, top=0.9, wspace=0.05, hspace=0.3 )
-    #
-    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    axs[0].set_ylabel(r'$e_\text{NM}(n_\text{nuc})$')
-    axs[0].set_xlim([0, 0.34])
-    axs[0].set_ylim([0, 35])
-    #
-    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    #axs[1].set_ylabel(r'$e_{sym}(n)$')
-    axs[1].set_xlim([0, 0.34])
-    axs[1].set_ylim([0, 35])
-    axs[1].tick_params('y', labelleft=False)
-    #
-    mb_check = []
-    #
-    for kmb,mb in enumerate(micro_mbs):
-        #
-        models, models_lower = nuda.matter.micro_models_mb( mb )
-        #
-        for model in models:
-            #
-            micro = nuda.matter.setupMicro( model = model )
-            if nuda.env.verb: micro.print_outputs( )
-            #
-            check = nuda.matter.setupCheck( eos = micro, band = band )
-            #
-            if check.isInside:
-                lstyle = 'solid'
-            else:
-                lstyle = 'dashed'
-            #
-            if micro.nm_e2a is not None:
-                print('mb:',mb,'model:',model)
-                if mb in mb_check:
-                    if micro.marker:
-                        if micro.e_err:
-                            axs[0].errorbar( micro.nm_den, micro.nm_e2a, yerr=micro.nm_e2a_err, marker=micro.marker, markevery=micro.every, linestyle=lstyle, errorevery=micro.every, color=nuda.param.col[kmb] )
-                        else:
-                            axs[0].plot( micro.nm_den, micro.nm_e2a, marker=micro.marker, markevery=micro.every, linestyle=lstyle, color=nuda.param.col[kmb] )
-                    else:
-                        if micro.e_err:
-                            axs[0].errorbar( micro.nm_den, micro.nm_e2a, yerr=micro.nm_e2a_err, marker=micro.marker, markevery=micro.every, linestyle=lstyle, errorevery=micro.every, color=nuda.param.col[kmb] )
-                        else:
-                            axs[0].plot( micro.nm_den, micro.nm_e2a, marker=micro.marker, markevery=micro.every, linestyle=lstyle, color=nuda.param.col[kmb] )
-                else:
-                    mb_check.append(mb)
-                    if micro.marker:
-                        if micro.e_err:
-                            axs[0].errorbar( micro.nm_den, micro.nm_e2a, yerr=micro.nm_e2a_err, marker=micro.marker, markevery=micro.every, linestyle=lstyle, label=mb, errorevery=micro.every, color=nuda.param.col[kmb] )
-                        else:
-                            axs[0].plot( micro.nm_den, micro.nm_e2a, marker=micro.marker, markevery=micro.every, linestyle=lstyle, label=mb, color=nuda.param.col[kmb] )
-                    else:
-                        if micro.e_err:
-                            axs[0].errorbar( micro.nm_den, micro.nm_e2a, yerr=micro.nm_e2a_err, marker=micro.marker, markevery=micro.every, linestyle=lstyle, label=mb, errorevery=micro.every, color=nuda.param.col[kmb] )
-                        else:
-                            axs[0].plot( micro.nm_den, micro.nm_e2a, marker=micro.marker, markevery=micro.every, linestyle=lstyle, label=mb, color=nuda.param.col[kmb] )
-            # end model
-        # end mb
-    axs[0].fill_between( band.den, y1=(band.e2a-band.e2a_std), y2=(band.e2a+band.e2a_std), color=band.color, alpha=band.alpha, visible=True )
-    axs[0].plot( band.den, (band.e2a-band.e2a_std), color='k', linestyle='dashed' )
-    axs[0].plot( band.den, (band.e2a+band.e2a_std), color='k', linestyle='dashed' )
-    axs[0].text(0.06,2,'microscopic models',fontsize='10')
-    #
-    model_check = []
-    #
-    for kmodel,model in enumerate(pheno_models):
-        #
-        params, params_lower = nuda.matter.pheno_params( model = model )
-        #
-        for param in params:
-            #
-            pheno = nuda.matter.setupPheno( model = model, param = param )
-            if nuda.env.verb: pheno.print_outputs( )
-            #
-            check = nuda.matter.setupCheck( eos = pheno, band = band )
-            #
-            if check.isInside:
-                lstyle = 'solid'
-            else:
-                lstyle = 'dashed'
-            #
-            if pheno.nm_e2a is not None: 
-                print('model:',model,' param:',param)
-                if model in model_check:
-                    axs[1].plot( pheno.nm_den, pheno.nm_e2a, linestyle=lstyle, color=nuda.param.col[kmodel] )
-                else:
-                    model_check.append(model)
-                    axs[1].plot( pheno.nm_den, pheno.nm_e2a, linestyle=lstyle, color=nuda.param.col[kmodel], label=model )
-            # end param
-        # end model
-    axs[1].fill_between( band.den, y1=(band.e2a-band.e2a_std), y2=(band.e2a+band.e2a_std), color=band.color, alpha=band.alpha, visible=True )
-    axs[1].plot( band.den, (band.e2a-band.e2a_std), color='k', linestyle='dashed' )
-    axs[1].plot( band.den, (band.e2a+band.e2a_std), color='k', linestyle='dashed' )
-    axs[1].text(0.06,2,'phenomenological models',fontsize='10')
-    #
-    #axs[1].legend(loc='upper left',fontsize='8', ncol=2)
-    #axs[0,1].legend(loc='upper left',fontsize='xx-small', ncol=2)
-    fig.legend(loc='upper left',bbox_to_anchor=(0.15,1.0),columnspacing=2,fontsize='8',ncol=5,frameon=False)
-    #
-    if pname is not None:
-    	plt.savefig(pname, dpi=200)
-    	plt.close()
-    #