nucleardatapy 0.2.0__py3-none-any.whl

nucleardatapy/fig/eos_setupAMBeq_fig.py

@@ -0,0 +1,204 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def eos_setupAMBeq_xp_fig( pname, models_micro, models_pheno ):
+    """
+    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.12, bottom=0.12, right=None, top=0.98, wspace=0.05, hspace=0.3 )
+    #
+    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    axs[0].set_ylabel(r'proton fraction $x_p$')
+    axs[0].set_xlim([0, 0.28])
+    axs[0].set_ylim([0, 0.2])
+    #axs[0].set_tick_params('y', right=True)
+    #axs[0].set_tick_params('x', top=True)
+    #
+    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    #axs[1].set_ylabel(r'proton fraction $x_p$')
+    axs[1].set_xlim([0, 0.28])
+    axs[1].set_ylim([0, 0.2])
+    #setp(axs[1].get_yticklabels(), visible=False)
+    axs[1].tick_params('y', labelleft=False)
+    #
+    for model in models_micro:
+        #
+        beta = nuda.eos.setupAMBeq( model = model, kind = 'micro' )
+        if nuda.env.verb_output: beta.print_outputs( )
+        #
+        if beta.esym is not None: 
+            print('model:',model)
+            axs[0].plot( beta.den, beta.x_p, marker='o', linestyle=beta.linestyle, label=beta.label, markevery=beta.every )
+    axs[0].text(0.02,0.18,'microscopic models',fontsize='10')
+    #axs[0].legend(loc='upper left',fontsize='8', ncol=3)
+    #
+    for model in models_pheno:
+        #
+        params, params_lower = nuda.matter.pheno_esym_params( model = model )
+        #
+        for param in params:
+            #
+            beta = nuda.eos.setupAMBeq( model = model, param = param, kind = 'pheno' )
+            if beta.esym is not None: 
+                print('model:',model,' param:',param)
+                #beta.label=None
+                axs[1].plot( beta.den, beta.x_p, linestyle=beta.linestyle, label=beta.label, markevery=beta.every )
+            if nuda.env.verb_output: pheno.print_outputs( )
+    #
+    #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,0.18,'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)
+    #
+    if pname is not None: 
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+
+def eos_setupAMBeq_xe_fig( pname, models_micro, models_pheno ):
+    """
+    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.12, bottom=0.12, right=None, top=0.98, wspace=0.05, hspace=0.3 )
+    #
+    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    axs[0].set_ylabel(r'electron fraction $x_e$')
+    axs[0].set_xlim([0, 0.28])
+    axs[0].set_ylim([0, 0.2])
+    #
+    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    #axs[1].set_ylabel(r'electron fraction $x_e$')
+    axs[1].set_xlim([0, 0.28])
+    axs[1].set_ylim([0, 0.2])
+    axs[1].tick_params('y', labelleft=False)
+    #
+    for model in models_micro:
+        #
+        beta = nuda.eos.setupAMBeq( model = model, kind = 'micro' )
+        if nuda.env.verb_output: beta.print_outputs( )
+        #
+        if beta.esym is not None: 
+            print('model:',model)
+            axs[0].plot( beta.den, beta.x_el, marker='o', linestyle=beta.linestyle, label=beta.label, markevery=beta.every )
+    axs[0].text(0.02,0.18,'microscopic models',fontsize='10')
+    #axs[0].legend(loc='upper left',fontsize='8', ncol=3)
+    #
+    for model in models_pheno:
+        #
+        params, params_lower = nuda.matter.pheno_esym_params( model = model )
+        #
+        for param in params:
+            #
+            beta = nuda.eos.setupAMBeq( model = model, param = param, kind = 'pheno' )
+            if beta.esym is not None: 
+                print('model:',model,' param:',param)
+                #beta.label=None
+                axs[1].plot( beta.den, beta.x_el, linestyle=beta.linestyle, label=beta.label, markevery=beta.every )
+            if nuda.env.verb_output: pheno.print_outputs( )
+    #
+    axs[1].text(0.02,0.18,'phenomenological models',fontsize='10')
+    #
+    if pname is not None: 
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+
+def eos_setupAMBeq_xmu_fig( pname, models_micro, models_pheno ):
+    """
+    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.12, bottom=0.12, right=None, top=0.98, wspace=0.05, hspace=0.3 )
+    #
+    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    axs[0].set_ylabel(r'muon fraction $x_\mu$')
+    axs[0].set_xlim([0, 0.28])
+    axs[0].set_ylim([0, 0.2])
+    #
+    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    #axs[1].set_ylabel(r'muon fraction $x_\mu$')
+    axs[1].set_xlim([0, 0.28])
+    axs[1].set_ylim([0, 0.2])
+    axs[1].tick_params('y', labelleft=False)
+    #
+    for model in models_micro:
+        #
+        beta = nuda.eos.setupAMBeq( model = model, kind = 'micro' )
+        if nuda.env.verb_output: beta.print_outputs( )
+        #
+        if beta.esym is not None: 
+            print('model:',model)
+            axs[0].plot( beta.den, beta.x_mu, marker='o', linestyle=beta.linestyle, label=beta.label, markevery=beta.every )
+    axs[0].text(0.02,0.18,'microscopic models',fontsize='10')
+    #axs[0].legend(loc='upper left',fontsize='8', ncol=3)
+    #
+    for model in models_pheno:
+        #
+        params, params_lower = nuda.matter.pheno_esym_params( model = model )
+        #
+        for param in params:
+            #
+            beta = nuda.eos.setupAMBeq( model = model, param = param, kind = 'pheno' )
+            if beta.esym is not None: 
+                print('model:',model,' param:',param)
+                #beta.label=None
+                axs[1].plot( beta.den, beta.x_mu, linestyle=beta.linestyle, label=beta.label, markevery=beta.every )
+            if nuda.env.verb_output: pheno.print_outputs( )
+    #
+    axs[1].text(0.02,0.18,'phenomenological models',fontsize='10')
+    #
+    if pname is not None: 
+    	plt.savefig(pname, dpi=200)
+    	plt.close()

nucleardatapy/fig/eos_setupAMLeq_fig.py

@@ -0,0 +1,209 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def eos_setupAMLeq_xe_fig( pname, models_micro, models_pheno ):
+    """
+    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}')
+    #
+    # xe at beta-equilibrium
+    #
+    asy = 0.5
+    #
+    fig, axs = plt.subplots(1,2)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.98, wspace=0.05, hspace=0.3 )
+    #
+    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    axs[0].set_ylabel(r'electron fraction $x_e$')
+    axs[0].set_xlim([0, 0.28])
+    axs[0].set_ylim([0.1, 0.3])
+    #
+    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    #axs[1].set_ylabel(r'electron fraction $x_e$')
+    axs[1].set_xlim([0, 0.28])
+    axs[1].set_ylim([0.1, 0.3])
+    axs[1].tick_params('y', labelleft=False)
+    #
+    for model in models_micro:
+        #
+        beta = nuda.eos.setupAMLeq( model = model, kind = 'micro', asy = asy )
+        if nuda.env.verb_output: beta.print_outputs( )
+        #
+        if beta.esym is not None: 
+            print('model:',model)
+            axs[0].plot( beta.den, beta.x_el, linestyle=beta.linestyle, label=beta.label, markevery=beta.every )
+            #axs[0].plot( beta.den, beta.x_el, marker='o', linestyle=beta.linestyle, label=beta.label, markevery=beta.every )
+    axs[0].text(0.08,0.22,'microscopic models',fontsize='10')
+    #axs[0].legend(loc='upper left',fontsize='8', ncol=3)
+    #
+    for model in models_pheno:
+        #
+        params, params_lower = nuda.matter.pheno_esym_params( model = model )
+        #
+        for param in params:
+            #
+            beta = nuda.eos.setupAMLeq( model = model, param = param, kind = 'pheno', asy = asy )
+            if beta.esym is not None: 
+                print('model:',model,' param:',param)
+                #beta.label=None
+                axs[1].plot( beta.den, beta.x_el, linestyle=beta.linestyle, label=beta.label, markevery=beta.every )
+            if nuda.env.verb_output: pheno.print_outputs( )
+    #
+    axs[1].text(0.08,0.22,'phenomenological models',fontsize='10')
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+
+def eos_setupAMLeq_xmu_fig( pname, models_micro, models_pheno ):
+    """
+    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}')
+    #
+    # xmu at beta-equilibrium
+    #
+    asy = 0.5
+    #
+    fig, axs = plt.subplots(1,2)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.98, wspace=0.05, hspace=0.3 )
+    #
+    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    axs[0].set_ylabel(r'muon fraction $x_\mu$')
+    axs[0].set_xlim([0, 0.28])
+    axs[0].set_ylim([0, 0.15])
+    #
+    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    #axs[1].set_ylabel(r'muon fraction $x_\mu$')
+    axs[1].set_xlim([0, 0.28])
+    axs[1].set_ylim([0, 0.15])
+    axs[1].tick_params('y', labelleft=False)
+    #
+    for model in models_micro:
+        #
+        beta = nuda.eos.setupAMLeq( model = model, kind = 'micro', asy = asy )
+        if nuda.env.verb_output: beta.print_outputs( )
+        #
+        if beta.esym is not None: 
+            print('model:',model)
+            axs[0].plot( beta.den, beta.x_mu, marker='o', linestyle=beta.linestyle, label=beta.label, markevery=beta.every )
+    axs[0].text(0.08,0.12,'microscopic models',fontsize='10')
+    #axs[0].legend(loc='upper left',fontsize='8', ncol=3)
+    #
+    for model in models_pheno:
+        #
+        params, params_lower = nuda.matter.pheno_esym_params( model = model )
+        #
+        for param in params:
+            #
+            beta = nuda.eos.setupAMLeq( model = model, param = param, kind = 'pheno', asy = asy )
+            if beta.esym is not None: 
+                print('model:',model,' param:',param)
+                #beta.label=None
+                axs[1].plot( beta.den, beta.x_mu, linestyle=beta.linestyle, label=beta.label, markevery=beta.every )
+            if nuda.env.verb_output: pheno.print_outputs( )
+    #
+    axs[1].text(0.08,0.12,'phenomenological models',fontsize='10')
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+
+def eos_setupAMLeq_xexmu_fig( pname, models_micro, models_pheno ):
+    """
+    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,1)
+    #fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.12, right=None, top=0.95, wspace=0.3, hspace=0.3 )
+    #
+    axs.set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)',fontsize='10')
+    axs.set_ylabel(r'$x_e$, $x_\mu$',fontsize='10')
+    axs.set_xlim([0, 0.28])
+    axs.set_ylim([0, 0.5])
+    #
+    asys = [ 0.1, 0.3, 0.5, 0.7, 0.9 ]
+    #asys = [ 0.5 ]
+    #
+    for inda,asy in enumerate(asys):
+        #
+        for model in models_micro:
+            #
+            continue
+            am = nuda.eos.setupAMLeq( model = model, kind = 'micro', asy = asy )
+            if nuda.env.verb_output: am.print_outputs( )
+            #
+            if am.esym is not None: 
+                print('model:',model)
+                axs.plot( am.den, am.x_mu, marker='o', linestyle=am.linestyle, label=am.label, markevery=am.every )
+        #axs[0].text(0.02,12,'microscopic models',fontsize='10')
+        #axs[0].text(0.02,10,'for $\delta=$'+str(asy),fontsize='10')
+        #
+        for model in models_pheno:
+            #
+            params, params_lower = nuda.matter.pheno_esym_params( model = model )
+            #
+            for param in params:
+                #
+                am = nuda.eos.setupAMLeq( model = model, param = param, kind = 'pheno', asy = asy )
+                if am.esym is not None: 
+                    print('model:',model,' param:',param)
+                    #beta.label=None
+                    axs.plot( am.den, am.x_el, linestyle='solid', label=r'$\delta=$'+str(asy), color=nuda.param.col[inda] )
+                    axs.plot( am.den, am.x_mu, linestyle='dashed', color=nuda.param.col[inda] )
+                if nuda.env.verb_output: pheno.print_outputs( )
+                break
+        #
+        axs.legend(loc='upper right',fontsize='10',ncol=3)
+        axs.text(0.05,0.35,r'$x_e$',fontsize='14')
+        axs.text(0.02,0.10,r'$x_\mu$',fontsize='14')
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()

nucleardatapy/fig/eos_setupAM_fig.py

@@ -0,0 +1,81 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def eos_setupAM_e_fig( pname, models_micro, models_pheno, asy ):
+    """
+    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.95, wspace=0.05, hspace=0.3 )
+    #
+    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
+    axs[0].set_ylabel(r'$E/A$')
+    axs[0].set_xlim([0, 0.28])
+    axs[0].set_ylim([-13, 15])
+    #
+    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([-13, 15])
+    axs[1].tick_params('y', labelleft=False)
+    #
+    for model in models_micro:
+        #
+        am = nuda.eos.setupAM( model = model, kind = 'micro', asy = asy )
+        if nuda.env.verb_output: am.print_outputs( )
+        #
+        if am.esym is not None: 
+            print('model:',model)
+            #print('den:',am.den)
+            #print('e2a:',am.e2a_nuc)
+            axs[0].plot( am.den, am.e2a_nuc, marker='o', linestyle=am.linestyle, label=am.label, markevery=am.every )
+    axs[0].text(0.02,12,'microscopic models',fontsize='10')
+    axs[0].text(0.02,10,r'for $\delta=$'+str(asy),fontsize='10')
+    #axs[0].legend(loc='upper left',fontsize='8', ncol=3)
+    #axs[0].legend(loc='lower center',bbox_to_anchor=(0.5,1.02),mode='expand',columnspacing=0,fontsize='8', ncol=2,frameon=False)
+    #
+    for model in models_pheno:
+        #
+        params, params_lower = nuda.matter.pheno_esym_params( model = model )
+        #
+        for param in params:
+            #
+            am = nuda.eos.setupAM( model = model, param = param, kind = 'pheno', asy = asy )
+            if am.esym is not None: 
+                print('model:',model,' param:',param)
+                #beta.label=None
+                axs[1].plot( am.den, am.e2a_nuc, linestyle=am.linestyle, label=am.label, markevery=am.every )
+            if nuda.env.verb_output: pheno.print_outputs( )
+    #
+    #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,12,'phenomenological models',fontsize='10')
+    axs[1].text(0.02,10,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)
+    #plt.tight_layout(rect=[0,0,1,0.95])
+    #
+    if pname is not None: 
+        plt.savefig(pname, dpi=200)
+        plt.close()

nucleardatapy/fig/hnuc_setupChart_fig.py

@@ -0,0 +1,54 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def hnuc_setupChart_fig( pname, table1L, table2L, table1Xi ):
+    """
+    Plot hyper-nuclear chart (N versus Z).\
+    The plot is 1x1 with:\
+    [0]: nuclear chart.
+
+    :param pname: name of the figure (*.png)
+    :type pname: str.
+    :param table1L: table.
+    :type table1L: str.
+    :param table2L: table.
+    :type table2L: str.
+    :param table1Xi: table.
+    :type table1Xi: str.
+
+    """
+    #
+    print(f'Plot name: {pname}')
+    #
+    hnuc1L = nuda.hnuc.setupRE1LExp( table = table1L )
+    hnuc2L = nuda.hnuc.setupRE2LExp( table = table2L )
+    hnuc1Xi = nuda.hnuc.setupRE1XiExp( table = table1Xi )
+    #
+    # plot
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.15, right=None, top=0.95, wspace=0.3, hspace=0.3)
+    #
+    axs.set_ylabel(r'Proton number $Z$')
+    axs.set_xlabel(r'Neutron number $N$')
+    axs.set_yscale('log')
+    axs.set_xscale('log')
+    axs.set_xlim([0.8, 140])
+    axs.set_ylim([0.8, 90.0])
+    #
+    axs.scatter( hnuc1L.N,      hnuc1L.Z,      marker='s', s=12, color=nuda.param.col[0], label=r'1$\Lambda$ from table '+table1L )
+    axs.scatter( hnuc2L.N+0.05, hnuc2L.Z+0.05, marker='s', s=12, color=nuda.param.col[1], label=r'2$\Lambda$ from table '+table2L )
+    axs.scatter( hnuc1Xi.N+0.1, hnuc1Xi.Z+0.1, marker='s', s=12, color=nuda.param.col[2], label=r'1$\Xi^{-}$ from table '+table1Xi )
+    #
+    #axs.text(0.15,12,r'$K_{sym}$='+str(int(Ksym))+' MeV',fontsize='12')
+    axs.legend(loc='lower right',fontsize='11',ncol=1)
+    #
+    axs.text(1.1,60,'Chart of hypernuclides in nuda toolkit',fontsize='14')
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+    #

nucleardatapy/fig/hnuc_setupRE1LExp_fig.py

@@ -0,0 +1,44 @@
+import numpy as np
+import matplotlib.pyplot as plt
+
+import nucleardatapy as nuda
+
+def hnuc_setupRE1LExp_plot( pname, tables ):
+    #
+    # plot
+    #
+    fig, axs = plt.subplots(1,1)
+    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
+    fig.subplots_adjust(left=0.12, bottom=0.15, right=None, top=0.95, wspace=0.3, hspace=0.3)
+    #
+    axs.set_ylabel(r'Binding Energy (MeV)')
+    axs.set_xlabel(r'$A^{-2/3}$')
+    axs.set_xlim([0.0, 0.28])
+    axs.set_ylim([-5.0, 32.0])
+    #
+    axs.plot( [0.0,0.28], [0.0,0.0], color='k', linestyle='dashed' )
+    for table in tables:
+        #
+        hnuc = nuda.hnuc.setupRE1LExp( table = table )
+        #
+        lab = []
+        for i in range(hnuc.nbdata):
+            if hnuc.label[i] not in lab:
+                lab.append(hnuc.label[i])
+                axs.errorbar( hnuc.A[i]**(-0.6666), hnuc.lre[i], yerr=hnuc.lre_err[i], marker=hnuc.mark[i], color=hnuc.color[i], label=hnuc.label[i], ms=4, ls='none' )
+            else:
+                axs.errorbar( hnuc.A[i]**(-0.6666), hnuc.lre[i], yerr=hnuc.lre_err[i], marker=hnuc.mark[i], color=hnuc.color[i], ms=4, ls='none' )                
+        #
+        axs.text(0.01,26.5,'1s')
+        axs.text(0.01,22,'1p')
+        axs.text(0.01,17,'1d')
+        axs.text(0.01,12,'1f')
+        axs.text(0.01,7,'1g')
+    #
+    #axs.text(0.15,12,r'$K_{sym}$='+str(int(Ksym))+' MeV',fontsize='12')
+    axs.legend(loc='upper right',fontsize='8',ncol=2)
+    #
+    if pname is not None:
+    	plt.savefig(pname, dpi=200)
+    	plt.close()
+    #