sesame_tools.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	tb2si1 (dens, pres, ener, nr, nt)
subroutine	tbusr1 (r, p, e, nr, nt, unitab)

Function/Subroutine Documentation

tb2si1()

subroutine tb2si1	(		dens,
			pres,
			ener,
		integer	nr,
		integer	nt )

Definition at line 28 of file sesame_tools.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NR,NT
      my_real
     .   dens(*),pres(nr,*),ener(nr,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J
      my_real
     .   dtab2si,ptab2si,etab2si
C
C---------------------------------------------------------------
C Passage des unites table Sesame type 301-306 aux unites SI
C    1 Mg/m^3 --> 1000 Kg/m^3
C    1 GPa    --> 1E9 Kg/(m.s^2)
C    1 MJ/Kg  --> 1E6 m^2/s^2
C
C DTAB2SI : coef. de passage des masses volumiques tables aux unites SI
C PTAB2SI : coef. de passage des pressions tables aux unites SI
C ETAB2SI : coef. de passage des energies libres tables aux unites SI
C
      dtab2si = 1000.
      ptab2si = 1.e09
      etab2si = 1.e06
C
      DO i=1,nr
          dens(i)=dens(i)*dtab2si
      ENDDO
C
      DO  i=1,nr
        DO  j=1,nt
          pres(i,j)=pres(i,j)*ptab2si
          ener(i,j)=ener(i,j)*etab2si
        ENDDO
      ENDDO
C
      RETURN

tbusr1()

subroutine tbusr1	(		r,
			p,
			e,
		integer	nr,
		integer	nt,
		type(unit_type_), intent(in)	unitab )

Definition at line 82 of file sesame_tools.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
      USE unitab_mod
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NR,NT
      TYPE(UNIT_TYPE_), INTENT(IN) :: UNITAB
      my_real
     .   r(*),p(nr,*),e(nr,*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J
      my_real alscal,amscal,atscal,dscal,pscal,escal
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      amscal=unitab%FAC_M_WORK
      alscal=unitab%FAC_L_WORK
      atscal=unitab%FAC_T_WORK
C----------------------------------------------------------------
C Passage des unites table Sesame type 301-306 aux unites
C utilisateur passer obligatoirement avant l'appel de cette
C routine par la routine TB2SI1
C coefficients de passage des unites SI aux unites utilisateur
C M --> masse
C L --> longueur
C T --> temps
C masse volumique : M.L^-3
C pression        : M.l^-1.T^-2
C energie         : L^2.T^-2
C L'energy array est en fait equivalent a la chaleur specifique 
C de transition de phase. Cf Aide Memoire de Physique, B. YAVORSKI
C ET A.DETLAF, p 927.
C
      dscal = (alscal**3)/amscal
      pscal = alscal*(atscal**2)/amscal
      escal = (atscal**2)/(alscal**2)
C
      DO  i=1,nr
          r(i)=r(i)*dscal
      ENDDO
C
      DO  i=1,nr
        DO  j=1,nt
          p(i,j)=p(i,j)*pscal
          e(i,j)=e(i,j)*escal
        ENDDO
      ENDDO
C
      RETURN