|
|
DataMuseum.dkPresents historical artifacts from the history of: DKUUG/EUUG Conference tapes |
This is an automatic "excavation" of a thematic subset of
See our Wiki for more about DKUUG/EUUG Conference tapes Excavated with: AutoArchaeologist - Free & Open Source Software. |
top - metrics - downloadIndex: S T
Length: 2227 (0x8b3)
Types: TextFile
Names: »S00005.f«
└─⟦db229ac7e⟧ Bits:30007240 EUUGD20: SSBA 1.2 / AFW Benchmarks
└─⟦this⟧ »EUUGD20/AFUU-ssba1.21/ssba1.21E/doduc/S00005.f«
└─⟦this⟧ »EUUGD20/AFUU-ssba1.21/ssba1.21F/doduc/S00005.f«
SUBROUTINE S00005(IORG)
IMPLICITDOUBLEPRECISION(A-H,O-Z)
COMMON/AAA11/V0011N(20),V0012X(20),V0011,V0012,
1TLIM,ZLIMIN,ZLIMAX,DTVOI,TIMIMP(20),DIFREF(20),DIFMIC(20),
3V00001,V22202,TEMIMP,TEMMIC,DELM,DELV,DELVI,DELVS,
6V22201,V22203,DTPRE,DTPHY,DTTEMP,DTCHAN,DTDIS,DTBEF,DTORG
7,NIMIC,NMULTR,IZONE,N00011(20),ISORT,NECRIT,NSTAR,IPLOT,NZONE
+,IFRE(20),ICO,IIMP,NPAS1(20),NPAS2(20)
COMMON/AAA13/TEM,P,DT,HNIV(22),VITESS
1 ,ZNIV,RI(2),ISS,NITERA,NIV,ISLNI,IBTR,IHTR,ITYP(22)
COMMON/AAA55/KPRESS,KDEB,KCON,KTGAI,KT876,ICAT(3),KIMP,KGLISS
KIMP=0
I=IZONE
IF(I.GE.NZONE)GOTO 1000
IF(TEM.GE.TIMIMP(I+1))I=I+1
1000 DFRE=DIFREF(I)
IFR=IFRE(I)
DMIC=DIFMIC(I)
NP1=NPAS1(I)
NP2=NPAS2(I)
V0011=V0011N(I)
V0012=V0012X(I)
IF(I.EQ.IZONE)GOTO 1200
KIMP=1
ICO=0
IZONE=I
IF(IFR.LT.10)GOTO 3000
GOTO 2000
1200 ECT=TEM-TEMIMP
IF(ECT.LT.DFRE)GOTO 1400
KIMP=1
IF(IFR.LT.10)GOTO 3000
GOTO 2000
1400 IF(IFR.LT.10)GOTO 3000
IF(IIMP.LT.NP1)GOTO 1600
KIMP=1
GOTO 2000
1600 IF(ICO.LT.NP2)GOTO 1800
KIMP=1
ICO=0
GOTO 2000
1800 IF(IFR.LT.20)GOTO 2000
IF(IORG.EQ.1)KIMP=1
2000 DTMI=V0011+1.D-6
IF(DT.LT.DTMI)GOTO 2100
ICO=0
GOTO 3000
2100 IF(ICO.LE.1)GOTO 2200
ICO=ICO+1
GOTO 3000
2200 IF(IORG.EQ.1)GOTO 3000
KIMP=1
ICO=ICO+1
3000 IIMP=IIMP+1
IF(KIMP.EQ.0)GOTO 4000
IOK=0
IF(ISORT.EQ.0)GOTO 3500
NSTAR=NSTAR+1
IF(NSTAR.GE.N00011(I))GOTO 3200
GOTO 3500
3200 NSTAR=0
NECRIT=NECRIT+1
IOK=1
3500 CONTINUE
IIMP=0
X=TEM/DFRE
IX=IDINT(X)
TEMIMP=IX*DFRE
NUM=NITERA-NIMIC
IF(NUM.LT.1)NUM=1
DTMOY=(TEM-TEMMIC)/NUM
TEMMIC=TEMIMP
NIMIC=NITERA
K=IFR-10*(IFR/10)
PRINT12345,TEM,NITERA
12345 FORMAT(9H TEMPS = ,F15.8,12H , NITERA : ,I10)
RETURN
4000 ECT=TEM-TEMMIC
IF(ECT.LT.DMIC)RETURN
X=TEM/DMIC
IX=IDINT(X)
NUM=NITERA-NIMIC
IF(NUM.LT.1)NUM=1
DTMOY=(TEM-TEMMIC)/NUM
TEMMIC=IX*DMIC
NIMIC=NITERA
RETURN
END