From: Apmsj@SLZ.NL
Subject: [NMusers] exceeding max. number of INPUT-variables
Date: Wed, 24 Aug 2005 11:36:06 +0200
dear colleagues
the attempt to code fortran (GET_SUBS.FOR) for excessive (>20)
INPUT-variables in a PD-model on a doublet anticancer drug treatment
creates the attached error-messages when running on DOS (NONMEM.ErrorMessages)
GET_SUBS.FOR was added to the control-stread as follows,
where the CSV-file was reduced to items not-constant within IDs,
and items constant within IDs were added via TEXT-file DATA.PRN:
$INPUT C=DROP ID TIME AMT RATE DV CMT MDV EVID
BSL=DROP V1=DROP VM1=DROP KM2=DROP VM2=DROP V3=DROP
Q=DROP K13=DROP K31=DROP CL=DROP V4=DROP QCA=DROP
V5=DROP
$DATA P59.CSV IGNORE=C
$SUBROUTINES OTHER=GET_SUBS.FOR ADVAN6 TOL=5
$MODEL COMP=(PAC)
COMP=(PER1)
COMP=(PER2)
COMP=(CARBO)
COMP=(PER3)
COMP=(AUC1.P)
COMP=(AUC2.C)
COMP=(PROL)
COMP=(TRANS1)
COMP=(TRANS2)
COMP=(TRANS3)
COMP=(ANC)
$PK
" FIRST
" COMMON/PRCOMG/IDUM1,IDUM2,IMAX,IDUM4,IDUM5
" INTEGER IDUM1,IDUM2,IMAX,IDUM4,IDUM5
" IMAX=1000000
" REAL CURID,ICOV1,ICOV2,ICOV3,ICOV4,ICOV5,ICOV6,ICOV7,ICOV8,ICOV9,ICOV10,ICOV11,ICOV12,ICOV13,ICOV14,ICOV15,ICOV16,MID
" LOGICAL READ
" IF(.NOT.READ)THEN
" CALL READER()
" CURID=1
" READ=.TRUE.
" ENDIF
" IF(NEWIND.LT.2)CALL GET_SUB(NEWIND,ID,CURID,ICOV1,ICOV2,ICOV3,ICOV4,ICOV5,ICOV6,ICOV7,ICOV8,ICOV9,ICOV10,ICOV11,ICOV12,ICOV13,ICOV14,ICOV15,MID)
CID=MID
IF(CID.NE.ID)THEN
"PRINT*,'ERROR CHECKING FAILED,CID=',CID,'ID=',ID
ENDIF
COV1=ICOV1
COV2=ICOV2
COV3=ICOV3
COV4=ICOV4
COV5=ICOV5
COV6=ICOV6
COV7=ICOV7
COV8=ICOV8
COV9=ICOV9
COV10=ICOV10
COV11=ICOV11
COV12=ICOV12
COV13=ICOV13
COV14=ICOV14
COV15=ICOV15
BSL=COV1
V1=COV2
KM1=COV3
VM1=COV4
KM2=COV5
VM2=COV6
K21=COV7
V3=COV8
Q=COV9
K13=COV10
K31=COV11
CL=COV12
V4=COV13
QCA=COV14
V5=COV15
I suspect errors within GET.SUBS.FOR, and any help would be appreciated.
with best regards
markus
Markus Joerger MD
Division of Experimental Oncology
The Netherlands Cancer Institute
Plesmanlaan 121
1066 CX AMSTERDAM
THE NETHERLANDS
+31/20-512-4657
+31/6-15-496290
apmsj@slz.nl
postal adres:
Department of Pharmacy&Pharmacology
Slotervaart Ziekenhuis
Louwesweg 6
1066 EC AMSTERDAM
THE NETHERLANDS
c:\data\tessa>call set path=c:\nmv\util;"C:\Perl\bin\;C:\Program Files\Microsoft
Visual Studio\Common\Tools;C:\Program Files\Microsoft Visual Studio\Common\Msde
v98\BIN;C:\Program Files\Microsoft Visual Studio\DF98\BIN;C:\Program Files\Micro
soft Visual Studio\VC98\BIN;C:\WINDOWS\system32;C:\WINDOWS;C:\WINDOWS\System32\W
bem"
Setting environment for using Visual Fortran tools
c:\data\tessa>cd..
C:\data>cd..
C:\>cd pdxpop1.1j
C:\pdxpop1.1j>cd paccarbo
C:\pdxpop1.1j\PACCARBO>nmfe5 022.ctl 022.res
WARNINGS AND ERRORS (IF ANY) FOR PROBLEM 1
(WARNING 2) NM-TRAN INFERS THAT THE DATA ARE POPULATION.
(WARNING 48) DES-DEFINED ITEMS ARE COMPUTED ONLY WHEN EVENT TIME
INCREASES. E.G., DISPLAYED VALUES ASSOCIATED WITH THE FIRST EVENT RECORD
OF AN INDIVIDUAL RECORD ARE COMPUTED WITH (THE LAST ADVANCE TO) AN EVENT
TIME OF THE PRIOR INDIVIDUAL RECORD.
Stop - Program terminated.
1 bestand(en) gekopieerd.
fsubs.for(7) : Error: Syntax error, found END-OF-STATEMENT when expecting one of
: ( > ...
READ(77,*) TID(I),TCOV(I,1), TCOV(I,2), TCOV(I,3), TCOV(I,4),
--------------------------------------------------------------------------------
-------^
fsubs.for(8) : Error: Syntax error, found END-OF-STATEMENT when expecting one of
: ( > ...
TCOV(I,5), TCOV(I,6), TCOV(I,7), TCOV(I,8), TCOV(I,9),
--------------------------------------------------------------------------------
-------^
fsubs.for(9) : Error: Syntax error, found END-OF-STATEMENT when expecting one of
: ( > ...
TCOV(I,10), TCOV(I,11), TCOV(I,12), TCOV(I,13), TCOV(I,14),
--------------------------------------------------------------------------------
-------^
fsubs.for(11) : Error: Illegal character in statement label field [C]
11 CONTINUE
---^
fsubs.for(11) : Error: Illegal character in statement label field [O]
11 CONTINUE
----^
fsubs.for(11) : Error: A continuation character is illegal on a line with a stat
ement label
11 CONTINUE
-----^
fsubs.for(11) : Error: Syntax error, found END-OF-STATEMENT when expecting one o
f: ( : % . = =>
11 CONTINUE
--------------------------------------------------------------------------------
^
fsubs.for(11) : Error: A specification statement cannot appear in the executable
section.
11 CONTINUE
^
fsubs.for(6) : Error: This statement label cannot be referenced in this context.
[11]
DO 11,I = 1,108
-----------^
fsubs.for(6) : Error: An unterminated block exists.
DO 11,I = 1,108
^
fsubs.for(14) : Error: Syntax error, found END-OF-STATEMENT when expecting one o
f: *
SUBROUTINE GET_SUB(NEWIND,ID,CURID,COV1,COV2,COV3,COV4,COV5,
--------------------------------------------------------------------------------
^
fsubs.for(15) : Error: Syntax error, found END-OF-STATEMENT when expecting one o
f: *
COV6,COV7,COV8,COV9,COV10,COV11,COV12,COV13,COV14,COV15,
--------------------------------------------------------------------------------
-------^
fsubs.for(18) : Error: Syntax error, found END-OF-STATEMENT when expecting one o
f: %FILL
REAL ID,CURID,COV1,COV2,COV3,COV4,COV5,COV6,COV7,COV8,COV9,
--------------------------------------------------------------------------------
-------^
fsubs.for(24) : Error: Illegal character in statement label field [C]
12 CONTINUE
---^
fsubs.for(24) : Error: Illegal character in statement label field [O]
12 CONTINUE
----^
fsubs.for(24) : Error: A continuation character is illegal on a line with a stat
ement label
12 CONTINUE
-----^
fsubs.for(24) : Error: Syntax error, found END-OF-STATEMENT when expecting one o
f: ( : % . = =>
12 CONTINUE
--------------------------------------------------------------------------------
^
fsubs.for(24) : Error: A specification statement cannot appear in the executable
section.
12 CONTINUE
^
fsubs.for(27) : Error: This statement label cannot be referenced in this context
. [12]
GOTO 12
-------------^
fsubs.for(175) : Error: Syntax error, found END-OF-STATEMENT when expecting one
of: ( * :: , : ) + . - % . ** / // .LT. < .LE. <= .EQ. == .NE. ...
XICOV2,ICOV3,ICOV4,ICOV5,ICOV
----------------------------------^
fsubs.for(167) : Error: A specification statement cannot appear in the executabl
e section.
REAL CURID,ICOV1,ICOV2,ICOV3,ICOV4,ICOV5,ICOV6,ICOV7,ICOV8,ICOV9,ICOV10
------^
fsubs.for(168) : Error: A specification statement cannot appear in the executabl
e section.
LOGICAL READ
------^
fsubs.for(169) : Error: The highest data type rank permitted is INTEGER(KIND=4).
[READ]
IF(.NOT.READ)THEN
--------------^
fsubs.for(169) : Error: The highest data type rank permitted is INTEGER(KIND=4).
IF(.NOT.READ)THEN
---------^
fsubs.for(169) : Error: A logical data type is required in this context.
IF(.NOT.READ)THEN
---------^
No nonmem execution.
There may be error messages in file fl32.txt.
C:\pdxpop1.1j\PACCARBO>
SUBROUTINE READER()
COMMON /READ/ TID,TCOV
REAL TID(139),TCOV(139,15)
OPEN(UNIT = 77,FILE = 'DATA.PRN')
REWIND 77
DO 11,I = 1,139
READ(77,*) TID(I),TCOV(I,1), TCOV(I,2), TCOV(I,3), TCOV(I,4),
TCOV(I,5), TCOV(I,6), TCOV(I,7), TCOV(I,8), TCOV(I,9),
TCOV(I,10), TCOV(I,11), TCOV(I,12), TCOV(I,13), TCOV(I,14),
TCOV(I,15), TCOV(I,16)
11 CONTINUE
END
SUBROUTINE GET_SUB(NEWIND,ID,CURID,COV1,COV2,COV3,COV4,COV5,
COV6,COV7,COV8,COV9,COV10,COV11,COV12,COV13,COV14,COV15,
COV16,MID)
COMMON /READ/ TID,TCOV
REAL ID,CURID,COV1,COV2,COV3,COV4,COV5,COV6,COV7,COV8,COV9,
COV10,COV11,COV12,COV13,COV14,COV15,COV16,MID
INTEGER NEWIND
REAL TID(139),TCOV(139,15)
IF (NEWIND.EQ.1) THEN
CURID=1
12 CONTINUE
IF(ID.GT.TID(CURID)) THEN
CURID =CURID + 1
GOTO 12
END IF
ELSEIF(NEWIND.EQ.0) THEN
CURID = 1
END IF
COV1 = TCOV(CURID,1)
COV2 = TCOV(CURID,2)
COV3 = TCOV(CURID,3)
COV4 = TCOV(CURID,4)
COV5 = TCOV(CURID,5)
COV6 = TCOV(CURID,6)
COV7 = TCOV(CURID,7)
COV8 = TCOV(CURID,8)
COV9 = TCOV(CURID,9)
COV10 = TCOV(CURID,10)
COV11 = TCOV(CURID,11)
COV12 = TCOV(CURID,12)
COV13 = TCOV(CURID,13)
COV14 = TCOV(CURID,14)
COV15 = TCOV(CURID,15)
COV16 = TCOV(CURID,16)
MID = TID(CURID)
END
C
C ID BSL V1 KM1 VM1 KM2 VM2 K21 V3 Q K13 K31 CL V10 QCA V11
2 5.7 9.4 1.38 23.8 1.1 303.9 2.36 6068.6 21.3 2.266 0.004 6.11 14.25 6.83 9.09
3 4.2 8.4 1.38 10.9 0.7 303.7 2.36 3945.1 26.4 3.143 0.007 8.14 14.43 10.55 10.55
4 3.3 8.4 1.38 19 0.8 285.2 2.36 4758.2 19.8 2.357 0.004 8.1 8.3 3.87 5.36
5 5.6 8.7 1.38 20.6 0.7 300 2.36 4502 18.6 2.138 0.004 9.36 18.24 9.68 9.53
6 4.5 10.9 1.38 19.2 0.8 304 2.36 4619.3 21.4 1.963 0.005 9.15 10.07 5.12 7.73
8 7.1 7.6 1.38 17.8 0.8 288 2.36 4627.4 46.7 6.145 0.01 9.51 16.32 8.07 12.28
10 5.8 7.6 1.38 13.3 0.7 282 2.36 3159.8 30.8 4.053 0.01 7.43 13.59 6.71 9.47
11 4.2 8.3 1.38 23.8 0.7 308.1 2.36 4461.5 18.1 2.181 0.004 6.44 15.17 7.88 7.92
12 3.2 9.6 1.38 8.7 0.8 303.3 2.36 3743.4 26.1 2.719 0.007 6.99 15.84 5.51 12.06
13 4.7 7.1 1.38 9.9 0.9 284.4 2.36 4408.8 14.3 2.014 0.003 6.31 14.4 7.99 7.05
14 3.8 7.3 1.38 3.7 1.2 298.7 2.36 3791.2 10 1.37 0.003 4.04 15.67 5.45 6.52
15 4.4 9.1 1.38 6.8 0.7 311.3 2.36 3336.7 20.4 2.242 0.006 4.73 13.7 13.7 8.56
16 2.3 7.7 1.38 17.6 0.9 266 2.36 4378.9 13.8 1.792 0.003 7.46 10.63 6.45 8.39
18 3 7 1.38 10 0.8 273.4 2.36 3893.2 20.1 2.871 0.005 5.91 10.44 6.51 11.97
19 4 8.1 1.38 17.6 0.8 283.1 2.36 4342.1 12.7 1.568 0.003 4.81 11.77 6.6 8.06
20 8.1 5 1.38 17.3 0.5 264.8 2.36 3837.4 38.5 7.7 0.01 8.57 11.67 4.7 7.38
22 2.9 6.3 1.38 18.6 0.7 261.9 2.36 4115.9 16.9 2.683 0.004 4.19 12.33 9.56 9.9
23 3.2 6.2 1.38 20.7 0.7 266.3 2.36 4336.1 13.7 2.21 0.003 6.96 11.34 4.24 10.2
24 6.3 10.3 1.38 14.8 1.1 290.5 2.36 4659.9 8.6 0.835 0.002 5.9 12.06 5.75 7.81
25 5 9.3 1.38 16.3 1 301.8 2.36 5207.2 17.9 1.925 0.003 6.36 11.96 5.79 10.99
26 3.2 8.9 1.38 8 0.7 315 2.36 3753.7 19.2 2.157 0.005 4.71 9.46 4.03 9.43
27 3.7 8.8 1.38 17.3 0.7 300.3 2.36 4540.3 13.6 1.545 0.003 6.56 13.68 5.46 10.02
28 5.9 7.3 1.38 5.5 1 319.9 2.36 3410.9 18 2.466 0.005 8.13 14.16 12.07 10.16
29 4 7.6 1.38 11.3 0.8 289.7 2.36 3772.1 37.2 4.895 0.01 3.94 3.1 1.92 9.15
31 4 8.4 1.38 13.5 0.8 287.3 2.36 4480.6 18.4 2.19 0.004 7.2 12.16 7.23 10.28
33 4.6 4.2 1.38 2.5 1 303.2 2.36 3516.9 6.8 1.619 0.002 4.77 10.9 11.49 5.8
34 3.3 8.9 1.38 9.4 0.8 293 2.36 4005.1 23.8 2.674 0.006 7.27 12.32 4.49 10.29
35 3.8 6.7 1.38 12.1 1 284.4 2.36 4801.4 17.1 2.552 0.004 5.66 11.87 5.11 8.65
36 6.1 12.3 1.38 25.8 1 303.5 2.36 5491.9 22.3 1.813 0.004 9.9 14.62 2.15 10.22
36 6.1 12.3 1.38 25.8 1 303.5 2.36 5491.9 22.3 1.813 0.004 9.9 14.62 2.15 10.22
37 6.7 7.9 1.38 10.6 0.9 290.1 2.36 4396.4 19.4 2.456 0.004 5.18 14.87 8.35 8.19
38 5.4 12.7 1.38 13 0.8 302.6 2.36 4423.4 15.5 1.22 0.004 5.7 10.78 6.43 8.07
39 7.1 7.5 1.38 13.4 0.9 280.5 2.36 4559.6 12.9 1.72 0.003 5.58 13.37 7.13 7.23
43 3.9 6.2 1.38 9.8 0.9 277.6 2.36 4428.5 16.1 2.597 0.004 5.33 8.85 4.97 6.06
44 4.2 7.6 1.38 9.6 1.3 266.1 2.36 4565.2 9.8 1.289 0.002 5.48 8.04 5.65 6.03
45 3.5 10. 1.38 18.9 0.7 312.3 2.36 4600.7 12.8 1.255 0.003 8.36 14.07 5.43 7.24
46 9 5.8 1.38 6.8 0.9 282.8 2.36 3560.2 23.9 4.121 0.007 9.3 17.22 5.41 10.07
48 3.6 8 1.38 11.2 1 276.3 2.36 4545.7 13.3 1.663 0.003 8.24 11.87 5.7 7.41
50 4.5 6.2 1.38 9.9 0.9 287.9 2.36 4535.3 13.6 2.194 0.003 7.3 12.22 6.51 7.22
51 2.2 8.8 1.38 7.6 1.1 290.2 2.36 4355.2 13.2 1.5 0.003 6.07 12.15 5.8 8.19
52 2.7 8.1 1.38 16 0.8 292.9 2.36 4520.1 30.7 3.79 0.007 10.32 15.27 5.21 9.73
53 4.9 11.9 1.38 16.7 1 287 2.36 4601.8 13.4 1.126 0.003 4.89 9.46 6.08 6.86
56 3.7 10.4 1.38 25.3 1 290.3 2.36 4641.8 15.2 1.462 0.003 7.48 9.71 5.81 6
57 4.5 7 1.38 9.6 0.8 290.6 2.36 4041 21.2 3.029 0.005 7.37 11.6 5.96 8.06
59 4.6 8.2 1.38 7.4 0.9 291.3 2.36 4296.6 12.3 1.5 0.003 15.42 3.07 1.07 9.58
60 5.9 8.5 1.38 6.7 0.7 303.3 2.36 3349 20.3 2.388 0.006 7.08 13.89 7.32 10.45
61 7.2 4.6 1.38 8.5 0.8 265.4 2.36 4410.5 13.5 2.935 0.003 5.67 9.34 5.6 6.23
62 5.9 8.5 1.38 11 0.9 294.4 2.36 4432.8 18.7 2.2 0.004 7.92 10.18 5.59 6.92
63 2.7 20.4 1.38 12.7 0.6 315.2 2.36 4126.8 20.6 1.01 0.005 5.44 11.76 7.36 8.96
66 3.6 10.4 1.38 11.9 0.7 309 2.36 4282.7 20.4 1.962 0.005 5.01 11.69 6.81 8.18
68 6.5 7.2 1.38 30.2 0.9 281.1 2.36 4594.3 18.8 2.611 0.004 10.87 11.24 7.62 6.69
70 3.2 3.6 1.38 5.2 0.9 271.2 2.36 4332.9 6.4 1.778 0.001 10.24 9.8 4.63 9.04
71 3.6 11.8 1.38 7.6 0.9 297.1 2.36 4369.8 9.4 0.797 0.002 6.65 22.82 11.91 9.8
73 7.6 9.2 1.38 21.2 0.9 285.8 2.36 4617.8 14.5 1.576 0.003 8.32 9.91 3.42 8.09
75 3.1 6.4 1.38 21.5 0.8 281.9 2.36 4574.9 16.9 2.641 0.004 8.52 10.56 3.05 7.85
75 3.1 6.4 1.38 21.5 0.8 281.9 2.36 4574.9 16.9 2.641 0.004 8.52 10.56 3.05 7.85
76 6.3 7.1 1.38 6.1 0.8 291.4 2.36 4371.4 7.7 1.085 0.002 6.6 11.66 6.6 6.79
78 3.5 7 1.38 12.5 0.9 282.7 2.36 4471 24.9 3.557 0.006 10.64 13.48 8 8.9
79 4.4 6.5 1.38 6.2 1 285.1 2.36 4153.4 14.4 2.215 0.003 6.96 10.21 4.06 7.02
80 3.6 9.1 1.38 10.3 0.8 305.8 2.36 4139.7 25.2 2.769 0.006 6.73 13.07 6.77 8.93
81 2.7 8 1.38 15 0.8 285.8 2.36 4438.5 22.6 2.825 0.005 10.61 15.45 6.97 8.71
83 6.2 6.8 1.38 7 0.9 285.3 2.36 4309.2 10.6 1.559 0.002 6.99 13.57 9.34 8.57
84 6.9 3.4 1.38 9.1 1 251.4 2.36 4530.4 5.9 1.735 0.001 7.37 11.6 5.96 8.06
86 2.7 7.6 1.38 10.4 0.9 285.8 2.36 4396.4 20.2 2.658 0.005 8.45 14.62 7.64 8.45
87 4.8 11.3 1.38 18.1 1 291.1 2.36 4692.6 18.4 1.628 0.004 9.39 11.42 5.98 10.35
88 5.8 5.4 1.38 7.3 1.1 263.3 2.36 4151.4 19.7 3.648 0.005 7.23 13.36 6.42 7.07
91 5.5 9.6 1.38 12 0.8 302.6 2.36 4362.2 17.6 1.833 0.004 6.53 13.37 6.94 9.82
92 5.9 8.4 1.38 16.9 0.9 285 2.36 4603.2 18.6 2.214 0.004 7.76 11.49 5.75 7.39
94 3.1 7.6 1.38 26.1 0.8 285 2.36 4564.3 24.4 3.211 0.005 9 8.93 2.49 8.01
96 7.3 6.2 1.38 15.2 0.9 277.7 2.36 4606.5 16.3 2.629 0.004 8.28 6.76 6.14 8.37
97 8.2 8.6 1.38 12.5 0.9 293.4 2.36 4446.4 25.6 2.977 0.006 10.24 9.29 11.67 12.47
98 5.1 6.4 1.38 10.1 0.8 281.7 2.36 4318 19.2 3 0.004 6.09 8.96 5.73 6.5
99 8.9 8.3 1.38 6.5 0.9 300 2.36 4227.1 12.9 1.554 0.003 6 11.77 6.27 6.48
100 3.8 7 1.38 30.8 0.7 282.9 2.36 4270.3 35.3 5.043 0.008 8.95 12.15 3.51 9
101 4.3 8.9 1.38 5.4 0.9 330.3 2.36 3818.8 15.8 1.775 0.004 6.37 14.31 6.64 7.74
102 2.4 6.2 1.38 22.6 0.8 273.2 2.36 4624.7 24.5 3.952 0.005 8.04 11.49 4.51 9.57
103 3.2 8.1 1.38 16.9 0.9 275.2 2.36 4579.8 15 1.852 0.003 11.85 7.29 6.1 5.65
107 2 7.3 1.38 12.8 0.7 277.1 2.36 4314.3 15.4 2.11 0.004 6.64 10.91 3.69 8.7
108 3.1 7.8 1.38 5 0.8 301.7 2.36 2481.3 24 3.077 0.01 6.68 7.16 8.28 6.48
109 1.9 12.5 1.38 31.6 1.4 278.2 2.36 4952.7 13.6 1.088 0.003 6.71 7.32 9.17 5.95
110 1.5 6.8 1.38 12.7 1 311.7 2.36 5429 19 2.794 0.003 7.38 7.54 8.23 6.36
111 5.3 7.9 1.38 6.4 0.9 315.5 2.36 4148.6 11.6 1.468 0.003 5.11 9.27 4.48 9.79
112 2.2 6.9 1.38 14.5 0.9 296.9 2.36 5779 16.1 2.333 0.003 9 10.15 4.16 7.25
113 4.5 7.7 1.38 7 0.8 303.9 2.36 3117.9 17.6 2.286 0.006 6.07 10.31 3.7 8.27
114 4.2 5.4 1.38 15.4 0.7 283 2.36 4462.4 15.9 2.944 0.004 7.8 12.5 5.87 7.55
115 6 6.9 1.38 23.2 0.8 281.2 2.36 4608.4 21.7 3.145 0.005 8.63 11.7 5.28 7.08
116 2.4 8.3 1.38 14.2 0.8 287.4 2.36 4536.9 19.5 2.349 0.004 8.61 15.79 4.35 10.83
117 2.6 6 1.38 7.9 0.9 292.7 2.36 4026.9 23.2 3.867 0.006 7.2 11.93 6.15 7.47
118 6 7.7 1.38 11.8 1.2 282.2 2.36 5150.5 27.4 3.558 0.005 7.37 11.6 5.96 8.06
119 6.6 8.1 1.38 15.2 0.9 285 2.36 4590.4 18 2.222 0.004 7.12 11.98 4.51 7.17
120 2.2 4.1 1.38 23.3 0.8 224.8 2.36 3085 22.6 5.512 0.007 5.43 9.55 4.17 6.89
121 1.6 7.4 1.38 21.4 0.9 284.4 2.36 4852 26.9 3.635 0.006 6.91 11.7 3.77 8.38
122 4.6 9.8 1.38 15.8 0.7 314.6 2.36 4498.7 24.3 2.48 0.005 7.16 10.19 1.91 8.32
123 3.1 8 1.38 25.5 1.1 302 2.36 5926.3 20.8 2.6 0.004 5.96 9.53 5.35 7.6
124 6.2 7.4 1.38 24 0.5 296.4 2.36 6436.7 21.1 2.851 0.003 9.04 8.95 2.97 5.13
125 5.3 8.5 1.38 14.3 0.7 280.9 2.36 4372.6 20 2.353 0.005 4.84 8.36 5.6 6.45
126 7 8 1.38 8.4 0.8 290.8 2.36 3617.5 18 2.25 0.005 5.19 7.92 6.35 5.92
127 6.6 7.7 1.38 11.9 0.7 285.7 2.36 3709.1 22.8 2.961 0.006 5.31 7.92 8.49 8.22
128 8.4 8.5 1.38 7.9 0.8 299.8 2.36 4283.4 13.4 1.576 0.003 7.87 10.21 5.23 7.89
130 3 6.9 1.38 12.3 0.8 283.2 2.36 4403.7 22.8 3.304 0.005 8.1 10.86 5.89 7.66
133 5.4 8.9 1.38 24.3 0.7 296.6 2.36 4508.8 22.6 2.539 0.005 7.07 8.99 5.24 8.71
135 4 5.7 1.38 13.5 0.7 274.7 2.36 4342.4 21.7 3.807 0.005 6.49 6.98 5.93 8.54
136 3.5 5.3 1.38 5.3 1 286.8 2.36 3960.8 14.1 2.66 0.004 5.24 7.46 3.76 5.99
137 9.4 7.3 1.38 14.4 0.8 289.6 2.36 4442.9 27 3.699 0.006 8.42 12.26 4.17 9.27
138 3.1 5.9 1.38 12.2 0.9 275.8 2.36 4484.9 19.5 3.305 0.004 7.37 11.6 5.96 8.06
139 5.2 7.7 1.38 9.2 0.8 295.7 2.36 3505.1 31.1 4.039 0.009 7.37 11.6 5.96 8.06
exceeding max. number of INPUT-variables
From: mark.e.sale@gsk.com
Subject: Re: [NMusers] exceeding max. number of INPUT-variables
Date: Wed, 24 Aug 2005 10:55:58 -0400
Markus,
You're right, there are a number of error in the Fortran source. It may not be realistic to get
this running without some Fortran expertise. But, the errors include:
1. No continuation line (non-blank character in column 6) for this code (and others)
READ(77,*) TID(I),TCOV(I,1), TCOV(I,2), TCOV(I,3), TCOV(I,4),
TCOV(I,5), TCOV(I,6), TCOV(I,7), TCOV(I,8), TCOV(I,9),
TCOV(I,10), TCOV(I,11), TCOV(I,12), TCOV(I,13), TCOV(I,14),
TCOV(I,15), TCOV(I,16)
Fortran lines must be <=72 columns, if you need continutation lines, you must use the continuation
character in column 6 of the next line.
the line numbers seem to have gotten messed up (likely by the various email servers) in
11 CONTINUE
the 11 needs to go in columns 1 -5, the CONTINUE starting in 7
You also have a number of tab characters in the file, which is not permitted in Fortran.
Finally, the code in the control file
" COMMON/PRCOMG/IDUM1,IDUM2,IMAX,IDUM4,IDUM5
" INTEGER IDUM1,IDUM2,IMAX,IDUM4,IDUM5
" IMAX=1000000
is not related to reading in the data, it is to increase the number of steps permitted by the ODE solver.
Below is an editted Get_subs.for, this compiles, don't know if it is correct. I also suspect that you'll have problems with long
lines in your verbatim code in the control file. Note in the code below the "1"s need to be in column 6 and
the "11" and "12" need to be in columns 1 to 5. Code needs to start in Column 7, the email servers may change this. Tabs
are not permitted.
SUBROUTINE READER()
COMMON /READ/ TID,TCOV
REAL TID(139),TCOV(139,15)
OPEN(UNIT = 77,FILE = 'DATA.PRN')
REWIND 77
DO 11,I = 1,139
READ(77,*) TID(I),TCOV(I,1), TCOV(I,2), TCOV(I,3), TCOV(I,4),
1 TCOV(I,5), TCOV(I,6), TCOV(I,7), TCOV(I,8), TCOV(I,9),
1 TCOV(I,10), TCOV(I,11), TCOV(I,12), TCOV(I,13), TCOV(I,14),
1 TCOV(I,15), TCOV(I,16)
11 CONTINUE
END
SUBROUTINE GET_SUB(NEWIND,ID,CURID,COV1,COV2,COV3,COV4,COV5,
1 COV6,COV7,COV8,COV9,COV10,COV11,COV12,COV13,COV14,COV15,
1 COV16,MID)
COMMON /READ/ TID,TCOV
REAL ID,CURID,COV1,COV2,COV3,COV4,COV5,COV6,COV7,COV8,COV9,
1 COV10,COV11,COV12,COV13,COV14,COV15,COV16,MID
INTEGER NEWIND
REAL TID(139),TCOV(139,15)
IF (NEWIND.EQ.1) THEN
CURID=1
12 CONTINUE
IF(ID.GT.TID(CURID)) THEN
CURID =CURID + 1
GOTO 12
END IF
ELSEIF(NEWIND.EQ.0) THEN
CURID = 1
END IF
COV1 = TCOV(CURID,1)
COV2 = TCOV(CURID,2)
COV3 = TCOV(CURID,3)
COV4 = TCOV(CURID,4)
COV5 = TCOV(CURID,5)
COV6 = TCOV(CURID,6)
COV7 = TCOV(CURID,7)
COV8 = TCOV(CURID,8)
COV9 = TCOV(CURID,9)
COV10 = TCOV(CURID,10)
COV11 = TCOV(CURID,11)
COV12 = TCOV(CURID,12)
COV13 = TCOV(CURID,13)
COV14 = TCOV(CURID,14)
COV15 = TCOV(CURID,15)
COV16 = TCOV(CURID,16)
MID = TID(CURID)
END
Mark Sale M.D.
Global Director, Research Modeling and Simulation
GlaxoSmithKline
919-483-1808
Mobile
919-522-6668
From: "Serge Guzy"
Subject: RE: [NMusers] exceeding max. number of INPUT-variables
Date: Wed, 24 Aug 2005 08:47:41 -0700
I think that your subroutine get subs contained errors.
I debugged it and here the one that can be built.
Hope this helps
Serge Guzy
President POP_PHARM
SUBROUTINE READER()
COMMON /READ/ TID,TCOV
REAL TID(139),TCOV(139,15)
OPEN(UNIT = 77,FILE = 'DATA.PRN')
REWIND 77
DO 11,I = 1,139
READ(77,*) TID(I),TCOV(I,1), TCOV(I,2), TCOV(I,3), TCOV(I,4),
1 TCOV(I,5), TCOV(I,6), TCOV(I,7), TCOV(I,8), TCOV(I,9),
1 TCOV(I,10), TCOV(I,11), TCOV(I,12), TCOV(I,13), TCOV(I,14),
1 TCOV(I,15), TCOV(I,16)
11 CONTINUE
END
SUBROUTINE GET_SUB(NEWIND,ID,CURID,COV1,COV2,COV3,COV4,COV5,
1 COV6,COV7,COV8,COV9,COV10,COV11,COV12,COV13,COV14,COV15,
1 COV16,MID)
COMMON /READ/ TID,TCOV
REAL ID,CURID,COV1,COV2,COV3,COV4,COV5,COV6,COV7,COV8,COV9,
1 COV10,COV11,COV12,COV13,COV14,COV15,COV16,MID
INTEGER NEWIND
REAL TID(139),TCOV(139,15)
IF (NEWIND.EQ.1) THEN
CURID=1
12 CONTINUE
IF(ID.GT.TID(CURID)) THEN
CURID =CURID + 1
GOTO 12
END IF
ELSEIF(NEWIND.EQ.0) THEN
CURID = 1
END IF
COV1 = TCOV(CURID,1)
COV2 = TCOV(CURID,2)
COV3 = TCOV(CURID,3)
COV4 = TCOV(CURID,4)
COV5 = TCOV(CURID,5)
COV6 = TCOV(CURID,6)
COV7 = TCOV(CURID,7)
COV8 = TCOV(CURID,8)
COV9 = TCOV(CURID,9)
COV10 = TCOV(CURID,10)
COV11 = TCOV(CURID,11)
COV12 = TCOV(CURID,12)
COV13 = TCOV(CURID,13)
COV14 = TCOV(CURID,14)
COV15 = TCOV(CURID,15)
COV16 = TCOV(CURID,16)
MID = TID(CURID)
END
From: mark.e.sale@gsk.com
Subject: Re: [NMusers] exceeding max. number of INPUT-variables
Date: Wed, 24 Aug 2005 16:17:31 -0400
Markus,
Here is the NONMEM documentation for using continuation lines in the $PK block verbatim code.
Continuation lines
If blanks occur in positions 1-5 following the initial " and a
non-blank in position 6, this conforms to the usual FORTRAN syn-
tax for continuation lines, and the text is not moved from these
positions. Example:
" X +D/E
I think you'll need this. If I had given this a little more thought, I would have just changed the variable
names so it could use the implicit type for the variables, so the REAL statement wouldn't be needed, although
the call to GET_SUBS is still too long, and will need a continuation line.
Mark Sale M.D.
Global Director, Research Modeling and Simulation
GlaxoSmithKline
919-483-1808
Mobile
919-522-6668
_______________________________________________________