Numérotation des lignes :

C PRITH1    SOURCE    BECC      09/12/07    21:15:57     6579      subroutine prith1(nordpo, acv, Tmaxcv, T, ether, cvtot)CC     This subroutine creates ether and cvtot (T and acv given.CC     INPUTC     *****CC     nordpo           = order of polynomial for cp and cv (see alsoC                        Tmaxcv)CC     acv              = vector such thatC                        cv = \sum_{i=1,nordpo+1} acv(i) T^{i-1}C                        ether = \sum_{i=1,nordpo+1} acv(i) T^{i} / (i)CC     Tmaxcv           = maximum temperature for cv polynomial expansionC                        cv(T) = cv(Tmaxcv) if T > TmaxcvCC     T                = temperatureCC     OUTPUTC     ******CC     cvtot     = cvCC     ether     = sensible energy (0 at 0K)CC************************************************************************CC     HISTORY AND MODIFICATIONSCc     07/12/2009      Created.CC************************************************************************CC      implicit none      integer nordpo      real*8 acv(1:(nordpo + 1)), Tmaxcv, T, ether, cvtot     &     , t2, funt, funcvC     t2 = min (t, tmaxcv)      t2 = 0.5d0 * ((t + tmaxcv) - abs (t - tmaxcv))      cvtot = acv(1)      ether = t2 * cvtot      funt = 1.0d0      do i1 = 2, nordpo + 1, 1         funt = funt * t2         funcv = acv(i1) * funt         cvtot = cvtot + funcv         ether = ether + (funcv * t2/ i1)      enddoc     t2 = max ((t - tmaxcv),0)      t2 = 0.5d0 * ((t - tmaxcv) + abs (t - tmaxcv))      ether = ether + (cvtot * t2)      return      endC

© Cast3M 2003 - Tous droits réservés.
Mentions légales