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

nucleardatapy/fig/matter_ENM_fig.py

@@ -1,119 +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.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    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)$')
-    axs[0].set_xlim([0, 0.33])
-    axs[0].set_ylim([0, 30])
-    #
-    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    #axs[1].set_ylabel(r'$e_{sym}(n)$')
-    axs[1].set_xlim([0, 0.33])
-    axs[1].set_ylim([0, 30])
-    axs[1].tick_params('y', labelleft=False)
-    #
-    mb_check = []
-    k = 0
-    #
-    for mb in micro_mbs:
-        #
-        models, models_lower = nuda.matter.micro_models_mb( mb )
-        #
-        for model in models:
-            #
-            enm = nuda.matter.setupMicro( model = model )
-            #
-            if enm.nm_e2a is not None:
-                print('mb:',mb,'model:',model)
-                if mb in mb_check:
-                    if enm.marker:
-                        if enm.err:
-                            axs[0].errorbar( enm.nm_den, enm.nm_e2a, yerr=enm.nm_e2a_err, marker=enm.marker, linestyle=None, errorevery=enm.every, color=nuda.param.col[k] )
-                        else:
-                            axs[0].plot( enm.nm_den, enm.nm_e2a, marker=enm.marker, linestyle=None, markevery=enm.every, color=nuda.param.col[k] )
-                    else:
-                        if enm.err:
-                            axs[0].errorbar( enm.nm_den, enm.nm_e2a, yerr=enm.nm_e2a_err, marker=enm.marker, linestyle=enm.linestyle, errorevery=enm.every, color=nuda.param.col[k] )
-                        else:
-                            axs[0].plot( enm.nm_den, enm.nm_e2a, marker=enm.marker, linestyle=enm.linestyle, markevery=enm.every, color=nuda.param.col[k] )
-                else:
-                    mb_check.append(mb)
-                    k += 1
-                    if enm.marker:
-                        if enm.err:
-                            axs[0].errorbar( enm.nm_den, enm.nm_e2a, yerr=enm.nm_e2a_err, marker=enm.marker, linestyle=None, label=mb, errorevery=enm.every, color=nuda.param.col[k] )
-                        else:
-                            axs[0].plot( enm.nm_den, enm.nm_e2a, marker=enm.marker, linestyle=None, label=mb, markevery=enm.every, color=nuda.param.col[k] )
-                    else:
-                        if enm.err:
-                            axs[0].errorbar( enm.nm_den, enm.nm_e2a, yerr=enm.nm_e2a_err, marker=enm.marker, linestyle=enm.linestyle, label=mb, errorevery=enm.every, color=nuda.param.col[k] )
-                        else:
-                            axs[0].plot( enm.nm_den, enm.nm_e2a, marker=enm.marker, linestyle=enm.linestyle, label=mb, markevery=enm.every, color=nuda.param.col[k] )
-                    #axs[0].plot( esym.den, esym.esym, color=nuda.param.col[k], label=mb )
-            if nuda.env.verb: esm.print_outputs( )
-    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.03,2,'microscopic models',fontsize='10')
-    #axs[0].legend(loc='upper left',fontsize='8', ncol=3)
-    #
-    model_check = []
-    k = 0
-    #
-    for model in pheno_models:
-        #
-        params, params_lower = nuda.matter.pheno_params( model = model )
-        #
-        for param in params:
-            #
-            enm = nuda.matter.setupPheno( model = model, param = param )
-            #
-            if enm.nm_e2a is not None: 
-                print('model:',model,' param:',param)
-                if model in model_check:
-                    axs[1].plot( enm.nm_den, enm.nm_e2a, color=nuda.param.col[k] )
-                else:
-                    model_check.append(model)
-                    k += 1
-                    axs[1].plot( enm.nm_den, enm.nm_e2a, color=nuda.param.col[k], label=model )
-                #pheno.label=None
-                #axs[1].plot( esym.den, esym.esym, label=esym.label )
-            if nuda.env.verb: esym.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.03,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)
-    #
-    #plt.tight_layout()
-    if pname is not None:
-    	plt.savefig(pname, dpi=200)
-    	plt.close()
-    #

nucleardatapy/fig/matter_ESM_fig.py

@@ -1,119 +0,0 @@
-import numpy as np
-import matplotlib.pyplot as plt
-
-import nucleardatapy as nuda
-
-def matter_ESM_fig( pname, micro_mbs, pheno_models, band ):
-    """
-    Plot nucleonic energy per particle E/A in matter.\
-    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.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    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{SM}(n)$')
-    axs[0].set_xlim([0, 0.33])
-    axs[0].set_ylim([-25, 5])
-    #
-    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    #axs[1].set_ylabel(r'$e_{sym}(n)$')
-    axs[1].set_xlim([0, 0.33])
-    axs[1].set_ylim([-25, 5])
-    axs[1].tick_params('y', labelleft=False)
-    #
-    mb_check = []
-    k = 0
-    #
-    for mb in micro_mbs:
-        #
-        models, models_lower = nuda.matter.micro_models_mb( mb )
-        #
-        for model in models:
-            #
-            esm = nuda.matter.setupMicro( model = model )
-            #
-            if esm.sm_e2a is not None:
-                print('mb:',mb,'model:',model)
-                if mb in mb_check:
-                    if esm.marker:
-                        if esm.err:
-                            axs[0].errorbar( esm.sm_den, esm.sm_e2a, yerr=esm.sm_e2a_err, marker=esm.marker, linestyle=None, errorevery=esm.every, color=nuda.param.col[k] )
-                        else:
-                            axs[0].plot( esm.sm_den, esm.sm_e2a, marker=esm.marker, linestyle=None, markevery=esm.every, color=nuda.param.col[k] )
-                    else:
-                        if esm.err:
-                            axs[0].errorbar( esm.sm_den, esm.sm_e2a, yerr=esm.sm_e2a_err, marker=esm.marker, linestyle=esm.linestyle, errorevery=esm.every, color=nuda.param.col[k] )
-                        else:
-                            axs[0].plot( esm.sm_den, esm.sm_e2a, marker=esm.marker, linestyle=esm.linestyle, markevery=esm.every, color=nuda.param.col[k] )
-                else:
-                    mb_check.append(mb)
-                    k += 1
-                    if esm.marker:
-                        if esm.err:
-                            axs[0].errorbar( esm.sm_den, esm.sm_e2a, yerr=esm.sm_e2a_err, marker=esm.marker, linestyle=None, label=mb, errorevery=esm.every, color=nuda.param.col[k] )
-                        else:
-                            axs[0].plot( esm.sm_den, esm.sm_e2a, marker=esm.marker, linestyle=None, label=mb, markevery=esm.every, color=nuda.param.col[k] )
-                    else:
-                        if esm.err:
-                            axs[0].errorbar( esm.sm_den, esm.sm_e2a, yerr=esm.sm_e2a_err, marker=esm.marker, linestyle=esm.linestyle, label=mb, errorevery=esm.every, color=nuda.param.col[k] )
-                        else:
-                            axs[0].plot( esm.sm_den, esm.sm_e2a, marker=esm.marker, linestyle=esm.linestyle, label=mb, markevery=esm.every, color=nuda.param.col[k] )
-                    #axs[0].plot( esym.den, esym.esym, color=nuda.param.col[k], label=mb )
-            if nuda.env.verb: esm.print_outputs( )
-    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.03,2,'microscopic models',fontsize='10')
-    #axs[0].legend(loc='upper left',fontsize='8', ncol=3)
-    #
-    model_check = []
-    k = 0
-    #
-    for model in pheno_models:
-        #
-        params, params_lower = nuda.matter.pheno_params( model = model )
-        #
-        for param in params:
-            #
-            esm = nuda.matter.setupPheno( model = model, param = param )
-            #
-            if esm.sm_e2a is not None: 
-                print('model:',model,' param:',param)
-                if model in model_check:
-                    axs[1].plot( esm.sm_den, esm.sm_e2a, color=nuda.param.col[k] )
-                else:
-                    model_check.append(model)
-                    k += 1
-                    axs[1].plot( esm.sm_den, esm.sm_e2a, color=nuda.param.col[k], label=model )
-                #pheno.label=None
-                #axs[1].plot( esym.den, esym.esym, label=esym.label )
-            if nuda.env.verb: esym.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.03,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=4,frameon=False)
-    #
-    #plt.tight_layout()
-    if pname is not None:
-    	plt.savefig(pname, dpi=200)
-    	plt.close()
-    #

nucleardatapy/fig/matter_Esym_fig.py

@@ -1,122 +0,0 @@
-import numpy as np
-import matplotlib.pyplot as plt
-
-import nucleardatapy as nuda
-
-def matter_Esym_fig( pname, micro_mbs, pheno_models, band ):
-    """
-    Plot nucleonic energy per particle E/A in matter.\
-    The plot is 2x2 with:\
-    [0,0]: E/A versus den.       [0,1]: E/A versus kfn.\
-    [1,0]: E/E_NRFFG versus den. [1,1]: E/E_NRFFG versus kfn.\
-
-    :param pname: name of the figure (*.png)
-    :type pname: str.
-    :param mb: many-body (mb) approach considered.
-    :type mb: str.
-    :param models: models to run on.
-    :type models: array of str.
-    :param band: object instantiated on the reference band.
-    :type band: object.
-    :param matter: can be 'SM' or 'NM'.
-    :type matter: str.
-
-    """
-    #
-    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.9, wspace=0.05, hspace=0.3 )
-    #
-    axs[0].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    axs[0].set_ylabel(r'$e_\text{sym}(n)$')
-    axs[0].set_xlim([0, 0.33])
-    axs[0].set_ylim([0, 60])
-    #
-    axs[1].set_xlabel(r'$n_\text{nuc}$ (fm$^{-3}$)')
-    #axs[1].set_ylabel(r'$e_{sym}(n)$')
-    axs[1].set_xlim([0, 0.33])
-    axs[1].set_ylim([0, 60])
-    axs[1].tick_params('y', labelleft=False)
-    #
-    mb_check = []
-    k = 0
-    #
-    for mb in micro_mbs:
-        #
-        models, models_lower = nuda.matter.micro_esym_models_mb( mb )
-        #
-        for model in models:
-            #
-            esym = nuda.matter.setupMicroEsym( model = model )
-            #
-            if esym.esym is not None:
-                print('mb:',mb,'model:',model)
-                if mb in mb_check:
-                    if esym.marker:
-                        if esym.err:
-                            axs[0].errorbar( esym.den, esym.esym, yerr=esym.esym_err, marker=esym.marker, linestyle=None, errorevery=esym.every, color=nuda.param.col[k] )
-                        else:
-                            axs[0].plot( esym.den, esym.esym, marker=esym.marker, linestyle=None, markevery=esym.every, color=nuda.param.col[k] )
-                    else:
-                        if esym.err:
-                            axs[0].errorbar( esym.den, esym.esym, yerr=esym.esym_err, marker=esym.marker, linestyle=esym.linestyle, errorevery=esym.every, color=nuda.param.col[k] )
-                        else:
-                            axs[0].plot( esym.den, esym.esym, marker=esym.marker, linestyle=esym.linestyle, markevery=esym.every, color=nuda.param.col[k] )
-                else:
-                    mb_check.append(mb)
-                    k += 1
-                    if esym.marker:
-                        if esym.err:
-                            axs[0].errorbar( esym.den, esym.esym, yerr=esym.esym_err, marker=esym.marker, linestyle=None, label=mb, errorevery=esym.every, color=nuda.param.col[k] )
-                        else:
-                            axs[0].plot( esym.den, esym.esym, marker=esym.marker, linestyle=None, label=mb, markevery=esym.every, color=nuda.param.col[k] )
-                    else:
-                        if esym.err:
-                            axs[0].errorbar( esym.den, esym.esym, yerr=esym.esym_err, marker=esym.marker, linestyle=esym.linestyle, label=mb, errorevery=esym.every, color=nuda.param.col[k] )
-                        else:
-                            axs[0].plot( esym.den, esym.esym, marker=esym.marker, linestyle=esym.linestyle, label=mb, markevery=esym.every, color=nuda.param.col[k] )
-                    #axs[0].plot( esym.den, esym.esym, color=nuda.param.col[k], label=mb )
-            if nuda.env.verb: esym.print_outputs( )
-    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.05,5,'microscopic models',fontsize='10')
-    #axs[0].legend(loc='upper left',fontsize='8', ncol=3)
-    #
-    model_check = []
-    k = 0
-    #
-    for model in pheno_models:
-        #
-        params, params_lower = nuda.matter.pheno_params( model = model )
-        #
-        for param in params:
-            #
-            esym = nuda.matter.setupPhenoEsym( model = model, param = param )
-            #
-            if esym.esym is not None: 
-                print('model:',model,' param:',param)
-                if model in model_check:
-                    axs[1].plot( esym.den, esym.esym, color=nuda.param.col[k] )
-                else:
-                    model_check.append(model)
-                    k += 1
-                    axs[1].plot( esym.den, esym.esym, color=nuda.param.col[k], label=model )
-                #pheno.label=None
-                #axs[1].plot( esym.den, esym.esym, label=esym.label )
-            if nuda.env.verb: esym.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.05,5,'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.2,1.0),columnspacing=2,fontsize='8',ncol=4,frameon=False)
-    #
-    #plt.tight_layout()
-    if pname is not None:
-    	plt.savefig(pname, dpi=200)
-    	plt.close()
-    #

nucleardatapy/fig/matter_setupNEPModelDist_fig.py

@@ -1,68 +0,0 @@
-import numpy as np
-import matplotlib.pyplot as plt
-
-import nucleardatapy as nuda
-
-def matter_setupNEPModelDist_fig( pname, models ):
-    """
-    Plot nucleonic energy per particle E/A in matter.\
-    The plot is 5x2 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 models: models to run on.
-    :type models: array of str.
-
-    """
-    #
-    print(f'Plot name: {pname}')
-    #
-    fig, axs = plt.subplots(5,2)
-    fig.tight_layout() # Or equivalently,  "plt.tight_layout()"
-    fig.subplots_adjust(left=0.12, bottom=0.06, right=None, top=0.9, wspace=0.3, hspace=0.4 )
-    #
-    axs[0,0].set_ylabel(r'$E_\text{sat}$')
-    axs[0,0].set_xlim([-16.5, -15])
-    axs[1,0].set_ylabel(r'$n_\text{sat}$')
-    axs[1,0].set_xlim([0.14, 0.18])
-    axs[2,0].set_ylabel(r'$K_\text{sat}$')
-    axs[2,0].set_xlim([180, 360])
-    axs[3,0].set_ylabel(r'$Q_\text{sat}$')
-    axs[3,0].set_xlim([-1000, 1000])
-    axs[4,0].set_ylabel(r'$m_\text{sat}^{*}/m$')
-    axs[4,0].set_xlim([0.35, 1.2])
-    axs[0,1].set_ylabel(r'$E_\text{sym}$')
-    axs[0,1].set_xlim([26, 40])
-    axs[1,1].set_ylabel(r'$L_\text{sym}$')
-    axs[1,1].set_xlim([0, 120])
-    axs[2,1].set_ylabel(r'$K_\text{sym}$')
-    axs[2,1].set_xlim([-400, 220])
-    axs[3,1].set_ylabel(r'$Q_\text{sym}$')
-    axs[3,1].set_xlim([-50, 900])
-    axs[4,1].set_ylabel(r'$\Delta m_\text{sat}^{*}/m$')
-    axs[4,1].set_xlim([-0.5, 1.1])
-    #
-    # Built distribution of NEP
-    #
-    for model in models:
-        #
-        dist = nuda.matter.setupNEPModelDist( model )
-        #
-        axs[0,0].hist( dist.Esat, bins=10, label=model )
-        axs[1,0].hist( dist.nsat, bins=10 )
-        axs[2,0].hist( dist.Ksat, bins=10 )
-        axs[3,0].hist( dist.Qsat, bins=10 )
-        axs[4,0].hist( dist.msat, bins=10 )
-        axs[0,1].hist( dist.Esym, bins=10 )
-        axs[1,1].hist( dist.Lsym, bins=10 )
-        axs[2,1].hist( dist.Ksym, bins=10 )
-        axs[3,1].hist( dist.Qsym, bins=10 )
-        axs[4,1].hist( dist.Dmsat, bins=10 )
-    #
-    #axs[0,0].legend(loc='lower right',fontsize='10',ncol=2)
-    fig.legend(loc='upper left',bbox_to_anchor=(0.2,0.99),columnspacing=2,fontsize='8',ncol=4,frameon=False)
-    #
-    if pname is not None:
-        plt.savefig(pname, dpi=300)
-        plt.close()

nucleardatapy-0.2.0.dist-info/METADATA

@@ -1,115 +0,0 @@
-Metadata-Version: 2.2
-Name: nucleardatapy
-Version: 0.2.0
-Summary: A toolkit for nuclear data processing and analysis.
-Author: Jerome Margueron, Sudhanva Lalit
-Author-email: 
-License: MIT
-Classifier: Intended Audience :: Science/Research
-Classifier: License :: OSI Approved :: MIT License
-Classifier: Programming Language :: Python
-Classifier: Programming Language :: Python :: 3
-Classifier: Operating System :: OS Independent
-Requires-Python: >=3.9
-Description-Content-Type: text/markdown
-License-File: LICENSE
-Requires-Dist: numpy
-Requires-Dist: scipy
-Requires-Dist: pandas
-
-# The toolkit `nucleardatapy`
-
-## Purpose:
-
-The purpose of this toolkit is to simply the access to data, that can be theoretical data or experimental ones. All data are provided with their reference, so when using these data in a scientific paper, reference to data should be provided explicitely. The reference to this toolkit could be given, but it should not mask the reference to data.
-
-This python toolkit is designed to provide: 
-1) microscopic calculations in nuclear matter, 
-2) phenomenological predictions in nuclear matter,
-3) experimental data for finite nuclei.
-
-## Installation of the toolkit:
-
-To install the toolkit, launch:
-```
-$ bash install.sh
-```
-
-In `install.sh`, the default directory where the toolkit is installed is `mylib` in the home directory. There, you have `mylib/nucleardatapy` folder pointing to the version of the toolkit defined in `install.sh` (currently `VER=0.1`). These default options could be changed directly in the header of `install.sh`.
-
-Create an environement variable that will be used by python:
-```
-export NUCLEARDATAPY_TK=/path/to/nucleardatapy
-```
-
-Add this environement variable to the one of python:
-```
-export PYTHONPATH=$NUCLEARDATAPY_TK
-```
-
-In this way, your library will be visible everywhere in your computer.
-
-Put these commands in your `.profile` or `.zprofile` or `.bashrc` for instance (depending on your OS).
-
-The first time, you should run again the `.zprofile` for instance:
-
-```
-$ source .zprofile
-```
-
-Now everything is done about the installation. You can go to the folder `mylib` in your home directory.
-
-## Use nucleardatapy python toolkit
-
-Go to the folder `mylib/nucleardatapy/samples/nucleardatapy_samples/` and try that:
-
-```
-$ python3 sample_SetupMicro.py
-```
-
-## Test the python toolkit
-
-A set of tests can be easily performed. They are stored in tests/ folder.
-
-Launch:
-
-```
-$ bash run_tests.sh
-```
-
-## Get started
-How to obtain microscopic results for APR equation of state:
-
-```Python
-import os
-import sys
-nuda_tk = os.getenv('NUCLEARDATAPY_TK')
-sys.path.insert(0, nuda_tk)
-
-import nucleardatapy as nuda
-
-# Instantiate a microscopic object
-mic = nuda.SetMicroMatter( model = '1998-VAR-AM-APR')
-
-# print outputs
-mic.print_outputs( )
-```
-
-## Contributing
-
-The file `how_to_contribute.md` details how contributors could join our team or share their results.
-
-## License
-
-TBC.
-
-## Report issues
-
-For the current version, we report issues chatting among us. 
-Once this toolkit is released, we should setup a way that users could contact us and report issues or difficulties in installing or using the toolkit.
-
-## Thanks
-
-A special thanks to all contributors who accepted to share their results in this toolkit.
-
-