Esope source files

   1 : C MODMAX    SOURCE    CHAT      05/01/13    01:50:03     5004           
   2 :       SUBROUTINE MODMAX(TP,EU,M,CODU,CODL,COD,BETJEF,BETFLU)
   3 : C
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8(A-H,O-Z) 
   6 :       DIMENSION CODU(9,9),CODL(8,8),COD(8)
   7 : C
   8 : *
   9 :       SEGMENT BETJEF
  10 :        REAL*8 AA,BETA,FC,PALF,YOUN,XNU,GFC,GFT,CAR,ETA,TDEF,
  11 :      &          TCON,DPSTF1,DPSTF2,TETA,PDT,TP00
  12 :        INTEGER ICT,ICC,IMOD,IVIS,ITR,
  13 :      &                ISIM,IBB,IGAU,IZON
  14 :       ENDSEGMENT
  15 :       SEGMENT BETFLU
  16 :        REAL*8 DATCOU,DATCUR,DATSEC,E28,PGTZO,PGDUR,TAU1,TAU2,
  17 :      &         TP0,TZER
  18 :        INTEGER ITYPE,IMD,NBRC,NCOE,NTZERO,NTPS,IFOR
  19 :       ENDSEGMENT
  20 : C  
  21 :       NC = NCOE+1       
  22 : C      
  23 : C----------------------------------------------------------------
  24 :       IF (IFOR.EQ.1) THEN
  25 : C     Modele visco elasto plastique solide de MAXWELL      
  26 : C----------------------------------------------------------------        
  27 :       EU=0.D0          
  28 : C
  29 :       DO 20 J=1,M
  30 : C      
  31 :       IF (J.EQ.M) THEN          
  32 : C
  33 : C      
  34 :       DO 30 N=1,NC     
  35 :         IF(N.EQ.1) THEN
  36 : C
  37 :           EU=EU+CODU(J,N)
  38 : C
  39 :         ELSE
  40 : C
  41 :           EUD=0.D0
  42 :           EUD=CODU(J,N)*EXP((DATCUR-(TP/86400-DATCOU))
  43 :      *        *COD(N-1))
  44 :           EU=EU+EUD
  45 : C
  46 :         ENDIF
  47 :    30 CONTINUE
  48 : C     
  49 : C      ELSE      
  50 : C      WRITE(*,*)'PAS DE VALEURS TROUVE POUR LE MODULE'
  51 : C      
  52 :       ENDIF
  53 : C      
  54 :    20 CONTINUE
  55 : C
  56 : C
  57 : C MODIFICATION DU 19/01/00
  58 : C********************************************************************
  59 : C
  60 : C----------------------------------------------------------
  61 :       ELSE IF (IFOR.EQ.2) THEN
  62 : C     Modele visco elasto plastique liquide de MAXWELL
  63 : C----------------------------------------------------------
  64 : C  
  65 :       EU=0.D0          
  66 : C
  67 :       DO 40 J=1,M
  68 : C      
  69 :       IF (J.EQ.M) THEN          
  70 : C      
  71 :         DO 50 N=1,NCOE     
  72 :           EUD=0.D0
  73 :           EUD=CODL(J,N)*EXP((DATCUR-(TP/86400-DATCOU))
  74 :      *        *COD(N))
  75 :           EU=EU+EUD
  76 :    50 CONTINUE
  77 : C     
  78 : C      ELSE     
  79 : C      WRITE(*,*)'PAS DE VALEURS TROUVE POUR LE MODULE'
  80 : C      
  81 :       ENDIF
  82 : C      
  83 :    40 CONTINUE
  84 : C 
  85 : C------------------------------------------------------------
  86 :       ENDIF
  87 : C------------------------------------------------------------
  88 : C
  89 :       RETURN
  90 :       END
  91 :  
  92 :