      MODULE NMPRD4P
      USE SIZES, ONLY: DPSIZE
      USE NMPRD4,ONLY: VRBL
      IMPLICIT NONE
      SAVE
      REAL(KIND=DPSIZE), DIMENSION (:),POINTER ::COM
      REAL(KIND=DPSIZE), POINTER ::MU_1,MU_2,MU_3,MU_4,MU_5,MU_6,MU_7
      REAL(KIND=DPSIZE), POINTER ::ka,ke,Vd,EC50,n,delta,c,q,d,e,s,Y
      REAL(KIND=DPSIZE), POINTER ::A000043,A000045,A000167,A000168
      REAL(KIND=DPSIZE), POINTER ::A000169,A000170,A000171,A000172
      REAL(KIND=DPSIZE), POINTER ::A000173,A000174,A000175,A000176
      REAL(KIND=DPSIZE), POINTER ::A000177,A000178,A000179,A000180
      REAL(KIND=DPSIZE), POINTER ::E000007,F000081,E000009,E000010
      REAL(KIND=DPSIZE), POINTER ::F000082,F000083,E000014,E000015
      REAL(KIND=DPSIZE), POINTER ::F000088,F000087,F000086,E000023
      REAL(KIND=DPSIZE), POINTER ::E000024,E000022,F000089,F000090
      REAL(KIND=DPSIZE), POINTER ::E000046,E000047,E000048,F000121
      REAL(KIND=DPSIZE), POINTER ::F000122,F000120,F000119,F000123
      REAL(KIND=DPSIZE), POINTER ::D000001,D000002,D000003,D000004
      REAL(KIND=DPSIZE), POINTER ::D000005,D000006,D000007,D000008
      REAL(KIND=DPSIZE), POINTER ::D000009,D000010,D000011,D000012
      REAL(KIND=DPSIZE), POINTER ::D000013,D000014,D000015,D000016
      REAL(KIND=DPSIZE), POINTER ::D000017,D000018,D000019,D000020
      REAL(KIND=DPSIZE), POINTER ::D000021,D000022,D000023,D000024
      REAL(KIND=DPSIZE), POINTER ::D000025,D000026,D000027,D000028
      REAL(KIND=DPSIZE), POINTER ::D000029,D000030,D000031,D000032
      REAL(KIND=DPSIZE), POINTER ::D000033,D000034,D000035,D000404
      REAL(KIND=DPSIZE), POINTER ::D000411,D000412,D000413,D000414
      REAL(KIND=DPSIZE), POINTER ::D000415,D000416,D000417,D000405
      REAL(KIND=DPSIZE), POINTER ::D000418,D000419,D000420,D000421
      REAL(KIND=DPSIZE), POINTER ::D000422,D000423,D000406,D000428
      REAL(KIND=DPSIZE), POINTER ::D000427,D000426,D000425,D000424
      REAL(KIND=DPSIZE), POINTER ::D000407,D000429,D000430,D000431
      REAL(KIND=DPSIZE), POINTER ::D000432,D000408,D000433,D000434
      REAL(KIND=DPSIZE), POINTER ::D000435,D000409,D000436,D000437
      REAL(KIND=DPSIZE), POINTER ::D000410,D000438,C000036,C000035
      CONTAINS
      SUBROUTINE ASSOCNMPRD4
      COM=>VRBL
      MU_1=>COM(000001);MU_2=>COM(000002);MU_3=>COM(000003)
      MU_4=>COM(000004);MU_5=>COM(000005);MU_6=>COM(000006)
      MU_7=>COM(000007);ka=>COM(000008);ke=>COM(000009)
      Vd=>COM(000010);EC50=>COM(000011);n=>COM(000012)
      delta=>COM(000013);c=>COM(000014);q=>COM(000015);d=>COM(000016)
      e=>COM(000017);s=>COM(000018);Y=>COM(000019)
      A000043=>COM(000020);A000045=>COM(000021);A000167=>COM(000022)
      A000168=>COM(000023);A000169=>COM(000024);A000170=>COM(000025)
      A000171=>COM(000026);A000172=>COM(000027);A000173=>COM(000028)
      A000174=>COM(000029);A000175=>COM(000030);A000176=>COM(000031)
      A000177=>COM(000032);A000178=>COM(000033);A000179=>COM(000034)
      A000180=>COM(000035);E000007=>COM(000036);F000081=>COM(000037)
      E000009=>COM(000038);E000010=>COM(000039);F000082=>COM(000040)
      F000083=>COM(000041);E000014=>COM(000042);E000015=>COM(000043)
      F000088=>COM(000044);F000087=>COM(000045);F000086=>COM(000046)
      E000023=>COM(000047);E000024=>COM(000048);E000022=>COM(000049)
      F000089=>COM(000050);F000090=>COM(000051);E000046=>COM(000052)
      E000047=>COM(000053);E000048=>COM(000054);F000121=>COM(000055)
      F000122=>COM(000056);F000120=>COM(000057);F000119=>COM(000058)
      F000123=>COM(000059);D000001=>COM(000060);D000002=>COM(000061)
      D000003=>COM(000062);D000004=>COM(000063);D000005=>COM(000064)
      D000006=>COM(000065);D000007=>COM(000066);D000008=>COM(000067)
      D000009=>COM(000068);D000010=>COM(000069);D000011=>COM(000070)
      D000012=>COM(000071);D000013=>COM(000072);D000014=>COM(000073)
      D000015=>COM(000074);D000016=>COM(000075);D000017=>COM(000076)
      D000018=>COM(000077);D000019=>COM(000078);D000020=>COM(000079)
      D000021=>COM(000080);D000022=>COM(000081);D000023=>COM(000082)
      D000024=>COM(000083);D000025=>COM(000084);D000026=>COM(000085)
      D000027=>COM(000086);D000028=>COM(000087);D000029=>COM(000088)
      D000030=>COM(000089);D000031=>COM(000090);D000032=>COM(000091)
      D000033=>COM(000092);D000034=>COM(000093);D000035=>COM(000094)
      D000404=>COM(000095);D000411=>COM(000096);D000412=>COM(000097)
      D000413=>COM(000098);D000414=>COM(000099);D000415=>COM(000100)
      D000416=>COM(000101);D000417=>COM(000102);D000405=>COM(000103)
      D000418=>COM(000104);D000419=>COM(000105);D000420=>COM(000106)
      D000421=>COM(000107);D000422=>COM(000108);D000423=>COM(000109)
      D000406=>COM(000110);D000428=>COM(000111);D000427=>COM(000112)
      D000426=>COM(000113);D000425=>COM(000114);D000424=>COM(000115)
      D000407=>COM(000116);D000429=>COM(000117);D000430=>COM(000118)
      D000431=>COM(000119);D000432=>COM(000120);D000408=>COM(000121)
      D000433=>COM(000122);D000434=>COM(000123);D000435=>COM(000124)
      D000409=>COM(000125);D000436=>COM(000126);D000437=>COM(000127)
      D000410=>COM(000128);D000438=>COM(000129);C000036=>COM(000130)
      C000035=>COM(000131)
      END SUBROUTINE ASSOCNMPRD4
      END MODULE NMPRD4P
      SUBROUTINE MODEL (IDNO,NCM,NPAR,IR,IATT,LINK)                           
      USE PRMOD_CHAR, ONLY: NAME                                              
      USE SIZES,     ONLY: DPSIZE,ISIZE,SD
      USE PRDIMS,    ONLY: GPRD,HPRD,GERD,HERD,GPKD
      INTEGER(KIND=ISIZE) :: IDNO,NCM,NPAR,IR,IATT,LINK,I,J                   
      DIMENSION :: IATT(IR,*),LINK(IR,*)                                      
      SAVE
      INTEGER(KIND=ISIZE), DIMENSION (5,7) :: MOD
      CHARACTER(LEN=SD), DIMENSION(5) :: CMOD
      DATA (MOD(I,  1),I=  1,  5)/&
      1,1,1,1,1 /
      DATA (MOD(I,  2),I=  1,  5)/&
      1,1,1,1,1 /
      DATA (MOD(I,  3),I=  1,  5)/&
      1,1,1,1,1 /
      DATA (MOD(I,  4),I=  1,  5)/&
      1,0,0,0,0 /
      DATA (MOD(I,  5),I=  1,  5)/&
      1,0,0,0,0 /
      DATA (MOD(I,  6),I=  1,  5)/&
      0,0,0,0,0 /
      DATA (MOD(I,  7),I=  1,  5)/&
      0,0,0,0,0 /
      DATA (CMOD(I),I=  1,  5) &
      /'COMP 1','COMP 2','COMP 3','COMP 4','COMP 5'/
      FORALL (I=1:5) NAME(I)=CMOD(I)
      FORALL (I=1:5,J=1:7) IATT(I,J)=MOD(I,J)
      IDNO=9999                                                               
      NCM=  5
      NPAR=011
      RETURN
      END
      SUBROUTINE PK(ICALL,IDEF,THETA,IREV,EVTREC,NVNT,INDXS,IRGG,GG,NETAS)      
      USE NMPRD4P
      USE SIZES,     ONLY: DPSIZE,ISIZE
      USE PRDIMS,    ONLY: GPRD,HPRD,GERD,HERD,GPKD
      USE NMBAYES_REAL,    ONLY: PRIOR_INFO
      USE PRDATA, ONLY: MXSTEP=>MXSTP01
      USE PRLS01_INT, ONLY: ITASK_=>ITASK
      USE PRLS01_REAL, ONLY: STOP_TIME
      USE NMBAYES_REAL, ONLY: ADDL_ACTUAL,ADDL_TIMEDIFF,ADDL_TIME
      USE NMPRD_INT, ONLY: PRED_IGNORE_DATA,PRED_IGNORE_DATA_TEST
      USE PRCOM_INT, ONLY: DES_DER=>MITER
      USE NMPRD_REAL,ONLY: ETA,EPS                                            
      USE NMPRD_INT, ONLY: MSEC=>ISECDER,MFIRST=>IFRSTDER,COMACT,COMSAV,IFIRSTEM,FIRSTEM
      USE NMPRD_INT, ONLY: MDVRES,ETASXI,NPDE_MODE,NOFIRSTDERCODE
      USE NMPRD_REAL, ONLY: DV_LOQ,CDF_L,DV_LAQ,CDF_LA
      USE NMPRD_INT, ONLY: IQUIT
      USE PROCM_INT, ONLY: NEWIND=>PNEWIF                                       
      USE PROCM_INT, ONLY: A_0FLG,A_UFLG
      USE PRMOD_REAL,ONLY: A_0,DA_0,D2A_0
      USE NMBAYES_REAL, ONLY: LDF                                             
      USE PRDDESLVU, ONLY: PASTZERO
      IMPLICIT REAL(KIND=DPSIZE) (A-Z)                                          
      REAL(KIND=DPSIZE) :: EVTREC                                               
      SAVE
      INTEGER(KIND=ISIZE) :: ICALL,IDEF,IREV,NVNT,INDXS,IRGG,NETAS              
      DIMENSION :: IDEF(7,*),THETA(*),EVTREC(IREV,*),INDXS(*),GG(IRGG,GPKD+1,*) 
      FIRSTEM=IFIRSTEM
      IF (ICALL <= 1) THEN                                                      
      CALL ASSOCNMPRD4
      IDEF(   1,0001)=  -9
      IDEF(   1,0002)=  -1
      IDEF(   1,0003)=   1
      IDEF(   1,0004)=   0
      IDEF(   2,0003)=   0
      IDEF(   2,0004)=   0
      CALL GETETA(ETA)                                                          
      IF (IQUIT == 1) RETURN                                                    
      RETURN                                                                    
      ENDIF                                                                     
      IF (NEWIND /= 2) THEN
      IF (ICALL == 4) THEN
      CALL SIMETA(ETA)
      ELSE
      CALL GETETA(ETA)
      ENDIF
      IF (IQUIT == 1) RETURN
      ENDIF
 !  level            0
      MU_1=DLOG(THETA(001)) 
      MU_2=DLOG(THETA(002)) 
      MU_3=DLOG(THETA(003)) 
      MU_4=DLOG(THETA(004)) 
      MU_5=DLOG(THETA(005)) 
      MU_6=DLOG(THETA(006)) 
      MU_7=DLOG(THETA(007)) 
      B000001=MU_1+ETA(001) 
      B000002=DEXP(B000001) 
      ka=B000002 
      B000003=MU_2+ETA(002) 
      B000004=DEXP(B000003) 
      ke=B000004 
      B000005=MU_3+ETA(003) 
      B000006=DEXP(B000005) 
      Vd=B000006 
      B000007=MU_4+ETA(004) 
      B000008=DEXP(B000007) 
      EC50=B000008 
      B000009=MU_5+ETA(005) 
      B000010=DEXP(B000009) 
      n=B000010 
      B000011=MU_6+ETA(006) 
      B000012=DEXP(B000011) 
      delta=B000012 
      B000013=MU_7+ETA(007) 
      B000014=DEXP(B000013) 
      c=B000014 
      q=100.D0 
      d=0.001D0 
      e=0.0000001D0 
      s=20000.D0 
      IF(A_0FLG == 1)THEN 
      A_0(1)=0.D0 
      ENDIF 
      IF(A_0FLG == 1)THEN 
      A_0(2)=0.D0 
      ENDIF 
      IF(A_0FLG == 1)THEN 
      B000015=c*delta 
      B000016=q*e 
      A_0(3)=B000015/B000016 
      ENDIF 
      IF(A_0FLG == 1)THEN 
      B000020=s*e*q-d*c*delta 
      B000021=q*delta*e 
      A_0(4)=B000020/B000021 
      ENDIF 
      IF(A_0FLG == 1)THEN 
      B000036=s*e*q-d*c*delta 
      B000037=c*delta*e 
      A_0(5)=B000036/B000037 
      B000045=delta 
      ENDIF 
      P000001=ka 
      P000002=ke 
      P000003=e 
      P000004=d 
      P000005=s 
      P000006=delta 
      P000007=EC50 
      P000008=n 
      P000009=Vd 
      P000010=q 
      P000011=c 
      IF (FIRSTEM == 1) THEN
!                      A000043 = DERIVATIVE OF Vd W.R.T. ETA(003)
      A000043=B000006 
      IF(A_0FLG == 1)THEN 
!                      A000066 = DERIVATIVE OF B000015 W.R.T. ETA(007)
      A000066=delta*B000014 
!                      A000067 = DERIVATIVE OF B000015 W.R.T. ETA(006)
      A000067=c*B000012 
      B000019=1.D0/B000016 
!                      A000070 = DERIVATIVE OF A_0(3) W.R.T. ETA(006)
      A000070=B000019*A000067 
!                      A000071 = DERIVATIVE OF A_0(3) W.R.T. ETA(007)
      A000071=B000019*A000066 
      DA_0(003,006)=A000070 
      DA_0(003,007)=A000071 
      ENDIF 
      IF(A_0FLG == 1)THEN 
      B000022=-d*delta 
!                      A000082 = DERIVATIVE OF B000020 W.R.T. ETA(007)
      A000082=B000022*B000014 
      B000023=-d*c 
!                      A000083 = DERIVATIVE OF B000020 W.R.T. ETA(006)
      A000083=B000023*B000012 
      B000027=q*e 
!                      A000086 = DERIVATIVE OF B000021 W.R.T. ETA(006)
      A000086=B000027*B000012 
      B000030=1.D0/B000021 
!                      A000087 = DERIVATIVE OF A_0(4) W.R.T. ETA(006)
      A000087=B000030*A000083 
!                      A000088 = DERIVATIVE OF A_0(4) W.R.T. ETA(007)
      A000088=B000030*A000082 
      B000031=-B000020/B000021/B000021 
!                      A000089 = DERIVATIVE OF A_0(4) W.R.T. ETA(006)
      A000089=B000031*A000086+A000087 
      DA_0(004,006)=A000089 
      DA_0(004,007)=A000088 
      ENDIF 
      IF(A_0FLG == 1)THEN 
      B000038=-d*delta 
!                      A000114 = DERIVATIVE OF B000036 W.R.T. ETA(007)
      A000114=B000038*B000014 
      B000039=-d*c 
!                      A000115 = DERIVATIVE OF B000036 W.R.T. ETA(006)
      A000115=B000039*B000012 
      B000043=B000045*e 
!                      A000118 = DERIVATIVE OF B000037 W.R.T. ETA(007)
      A000118=B000043*B000014 
      B000044=c*e 
!                      A000119 = DERIVATIVE OF B000037 W.R.T. ETA(006)
      A000119=B000044*B000012 
      B000048=1.D0/B000037 
!                      A000122 = DERIVATIVE OF A_0(5) W.R.T. ETA(006)
      A000122=B000048*A000115 
!                      A000123 = DERIVATIVE OF A_0(5) W.R.T. ETA(007)
      A000123=B000048*A000114 
      B000049=-B000036/B000037/B000037 
!                      A000124 = DERIVATIVE OF A_0(5) W.R.T. ETA(006)
      A000124=B000049*A000119+A000122 
!                      A000125 = DERIVATIVE OF A_0(5) W.R.T. ETA(007)
      A000125=B000049*A000118+A000123 
      DA_0(005,006)=A000124 
      DA_0(005,007)=A000125 
      ENDIF 
!                      A000167 = DERIVATIVE OF P000001 W.R.T. ETA(001)
      A000167=B000002 
!                      A000169 = DERIVATIVE OF P000002 W.R.T. ETA(002)
      A000169=B000004 
!                      A000171 = DERIVATIVE OF P000006 W.R.T. ETA(006)
      A000171=B000012 
!                      A000173 = DERIVATIVE OF P000007 W.R.T. ETA(004)
      A000173=B000008 
!                      A000175 = DERIVATIVE OF P000008 W.R.T. ETA(005)
      A000175=B000010 
!                      A000177 = DERIVATIVE OF P000009 W.R.T. ETA(003)
      A000177=A000043 
!                      A000179 = DERIVATIVE OF P000011 W.R.T. ETA(007)
      A000179=B000014 
      GG(0001,1,1)=P000001
      GG(0001,0002,1)=A000167
      GG(0002,1,1)=P000002
      GG(0002,0003,1)=A000169
      GG(0003,1,1)=P000003
      GG(0004,1,1)=P000004
      GG(0005,1,1)=P000005
      GG(0006,1,1)=P000006
      GG(0006,0007,1)=A000171
      GG(0007,1,1)=P000007
      GG(0007,0005,1)=A000173
      GG(0008,1,1)=P000008
      GG(0008,0006,1)=A000175
      GG(0009,1,1)=P000009
      GG(0009,0004,1)=A000177
      GG(0010,1,1)=P000010
      GG(0011,1,1)=P000011
      GG(0011,0008,1)=A000179
      ELSE
      GG(0001,1,1)=P000001
      GG(0002,1,1)=P000002
      GG(0003,1,1)=P000003
      GG(0004,1,1)=P000004
      GG(0005,1,1)=P000005
      GG(0006,1,1)=P000006
      GG(0007,1,1)=P000007
      GG(0008,1,1)=P000008
      GG(0009,1,1)=P000009
      GG(0010,1,1)=P000010
      GG(0011,1,1)=P000011
      ENDIF
      IF (MSEC == 1) THEN
!                      A000045 = DERIVATIVE OF A000043 W.R.T. ETA(003)
      A000045=B000006 
      IF(A_0FLG == 1)THEN 
!                      A000072 = DERIVATIVE OF A000066 W.R.T. ETA(007)
      A000072=delta*B000014 
!                      A000073 = DERIVATIVE OF A000067 W.R.T. ETA(007)
      A000073=B000012*B000014 
!                      A000074 = DERIVATIVE OF A000067 W.R.T. ETA(006)
      A000074=c*B000012 
!                      A000079 = DERIVATIVE OF A000070 W.R.T. ETA(006)
      A000079=B000019*A000074 
!                      A000080 = DERIVATIVE OF A000070 W.R.T. ETA(007)
      A000080=B000019*A000073 
!                      A000081 = DERIVATIVE OF A000071 W.R.T. ETA(007)
      A000081=B000019*A000072 
      D2A_0(003,006,006)=A000079 
      D2A_0(003,007,006)=A000080 
      D2A_0(003,007,007)=A000081 
      ENDIF 
      IF(A_0FLG == 1)THEN 
!                      A000090 = DERIVATIVE OF B000022 W.R.T. ETA(006)
      A000090=-d*B000012 
!                      A000091 = DERIVATIVE OF A000082 W.R.T. ETA(007)
      A000091=B000022*B000014 
!                      A000092 = DERIVATIVE OF B000023 W.R.T. ETA(007)
      A000092=-d*B000014 
!                      A000093 = DERIVATIVE OF A000083 W.R.T. ETA(007)
      A000093=B000012*A000092 
!                      A000094 = DERIVATIVE OF A000083 W.R.T. ETA(006)
      A000094=B000023*B000012 
!                      A000099 = DERIVATIVE OF A000086 W.R.T. ETA(006)
      A000099=B000027*B000012 
      B000032=-1.D0/B000021/B000021 
!                      A000100 = DERIVATIVE OF B000030 W.R.T. ETA(006)
      A000100=B000032*A000086 
!                      A000101 = DERIVATIVE OF A000087 W.R.T. ETA(006)
      A000101=A000083*A000100 
!                      A000102 = DERIVATIVE OF A000087 W.R.T. ETA(006)
      A000102=B000030*A000094+A000101 
!                      A000103 = DERIVATIVE OF A000087 W.R.T. ETA(007)
      A000103=B000030*A000093 
!                      A000104 = DERIVATIVE OF A000088 W.R.T. ETA(007)
      A000104=B000030*A000091 
      B000033=-1.D0/B000021/B000021 
!                      A000105 = DERIVATIVE OF B000031 W.R.T. ETA(006)
      A000105=B000033*A000083 
!                      A000106 = DERIVATIVE OF B000031 W.R.T. ETA(007)
      A000106=B000033*A000082 
      B000034=B000020/B000021/B000021/B000021 
!                      A000107 = DERIVATIVE OF B000031 W.R.T. ETA(006)
      A000107=B000034*A000086+A000105 
      B000035=B000020/B000021/B000021/B000021 
!                      A000108 = DERIVATIVE OF B000031 W.R.T. ETA(006)
      A000108=B000035*A000086+A000107 
!                      A000109 = DERIVATIVE OF A000089 W.R.T. ETA(006)
      A000109=A000086*A000108 
!                      A000110 = DERIVATIVE OF A000089 W.R.T. ETA(007)
      A000110=A000086*A000106 
!                      A000111 = DERIVATIVE OF A000089 W.R.T. ETA(006)
      A000111=B000031*A000099+A000109 
!                      A000112 = DERIVATIVE OF A000089 W.R.T. ETA(007)
      A000112=A000103+A000110 
!                      A000113 = DERIVATIVE OF A000089 W.R.T. ETA(006)
      A000113=A000102+A000111 
      D2A_0(004,006,006)=A000113 
      D2A_0(004,007,006)=A000112 
      D2A_0(004,007,007)=A000104 
      ENDIF 
      IF(A_0FLG == 1)THEN 
!                      A000126 = DERIVATIVE OF B000038 W.R.T. ETA(006)
      A000126=-d*B000012 
!                      A000127 = DERIVATIVE OF A000114 W.R.T. ETA(007)
      A000127=B000038*B000014 
!                      A000128 = DERIVATIVE OF B000039 W.R.T. ETA(007)
      A000128=-d*B000014 
!                      A000129 = DERIVATIVE OF A000115 W.R.T. ETA(007)
      A000129=B000012*A000128 
!                      A000130 = DERIVATIVE OF A000115 W.R.T. ETA(006)
      A000130=B000039*B000012 
!                      A000135 = DERIVATIVE OF B000043 W.R.T. ETA(006)
      A000135=e*B000012 
!                      A000136 = DERIVATIVE OF A000118 W.R.T. ETA(007)
      A000136=B000043*B000014 
!                      A000137 = DERIVATIVE OF B000044 W.R.T. ETA(007)
      A000137=e*B000014 
!                      A000138 = DERIVATIVE OF A000119 W.R.T. ETA(007)
      A000138=B000012*A000137 
!                      A000139 = DERIVATIVE OF A000119 W.R.T. ETA(006)
      A000139=B000044*B000012 
      B000050=-1.D0/B000037/B000037 
!                      A000144 = DERIVATIVE OF B000048 W.R.T. ETA(006)
      A000144=B000050*A000119 
!                      A000145 = DERIVATIVE OF B000048 W.R.T. ETA(007)
      A000145=B000050*A000118 
!                      A000146 = DERIVATIVE OF A000122 W.R.T. ETA(007)
      A000146=A000115*A000145 
!                      A000147 = DERIVATIVE OF A000122 W.R.T. ETA(006)
      A000147=A000115*A000144 
!                      A000148 = DERIVATIVE OF A000122 W.R.T. ETA(006)
      A000148=B000048*A000130+A000147 
!                      A000149 = DERIVATIVE OF A000122 W.R.T. ETA(007)
      A000149=B000048*A000129+A000146 
!                      A000150 = DERIVATIVE OF A000123 W.R.T. ETA(007)
      A000150=A000114*A000145 
!                      A000151 = DERIVATIVE OF A000123 W.R.T. ETA(007)
      A000151=B000048*A000127+A000150 
      B000051=-1.D0/B000037/B000037 
!                      A000152 = DERIVATIVE OF B000049 W.R.T. ETA(006)
      A000152=B000051*A000115 
!                      A000153 = DERIVATIVE OF B000049 W.R.T. ETA(007)
      A000153=B000051*A000114 
      B000052=B000036/B000037/B000037/B000037 
!                      A000154 = DERIVATIVE OF B000049 W.R.T. ETA(006)
      A000154=B000052*A000119+A000152 
!                      A000155 = DERIVATIVE OF B000049 W.R.T. ETA(007)
      A000155=B000052*A000118+A000153 
      B000053=B000036/B000037/B000037/B000037 
!                      A000156 = DERIVATIVE OF B000049 W.R.T. ETA(006)
      A000156=B000053*A000119+A000154 
!                      A000157 = DERIVATIVE OF B000049 W.R.T. ETA(007)
      A000157=B000053*A000118+A000155 
!                      A000158 = DERIVATIVE OF A000124 W.R.T. ETA(007)
      A000158=A000119*A000157 
!                      A000159 = DERIVATIVE OF A000124 W.R.T. ETA(006)
      A000159=A000119*A000156 
!                      A000160 = DERIVATIVE OF A000124 W.R.T. ETA(006)
      A000160=B000049*A000139+A000159 
!                      A000161 = DERIVATIVE OF A000124 W.R.T. ETA(007)
      A000161=B000049*A000138+A000158 
!                      A000162 = DERIVATIVE OF A000124 W.R.T. ETA(007)
      A000162=A000149+A000161 
!                      A000163 = DERIVATIVE OF A000124 W.R.T. ETA(006)
      A000163=A000148+A000160 
!                      A000164 = DERIVATIVE OF A000125 W.R.T. ETA(007)
      A000164=A000118*A000157 
!                      A000165 = DERIVATIVE OF A000125 W.R.T. ETA(007)
      A000165=B000049*A000136+A000164 
!                      A000166 = DERIVATIVE OF A000125 W.R.T. ETA(007)
      A000166=A000151+A000165 
      D2A_0(005,006,006)=A000163 
      D2A_0(005,007,006)=A000162 
      D2A_0(005,007,007)=A000166 
      ENDIF 
!                      A000168 = DERIVATIVE OF A000167 W.R.T. ETA(001)
      A000168=B000002 
!                      A000170 = DERIVATIVE OF A000169 W.R.T. ETA(002)
      A000170=B000004 
!                      A000172 = DERIVATIVE OF A000171 W.R.T. ETA(006)
      A000172=B000012 
!                      A000174 = DERIVATIVE OF A000173 W.R.T. ETA(004)
      A000174=B000008 
!                      A000176 = DERIVATIVE OF A000175 W.R.T. ETA(005)
      A000176=B000010 
!                      A000178 = DERIVATIVE OF A000177 W.R.T. ETA(003)
      A000178=A000045 
!                      A000180 = DERIVATIVE OF A000179 W.R.T. ETA(007)
      A000180=B000014 
      GG(0001,0002,0002)=A000168
      GG(0002,0003,0003)=A000170
      GG(0006,0007,0007)=A000172
      GG(0007,0005,0005)=A000174
      GG(0008,0006,0006)=A000176
      GG(0009,0004,0004)=A000178
      GG(0011,0008,0008)=A000180
      ENDIF
      RETURN
      END
      SUBROUTINE ERROR (ICALL,IDEF,THETA,IREV,EVTREC,NVNT,INDXS,F,G,HH)       
      USE NMPRD4P
      USE SIZES,     ONLY: DPSIZE,ISIZE
      USE PRDIMS,    ONLY: GPRD,HPRD,GERD,HERD,GPKD
      USE NMPRD_INT, ONLY: PRED_IGNORE_DATA,PRED_IGNORE_DATA_TEST
      USE NMPRD_REAL,ONLY: ETA,EPS                                            
      USE NMPRD_INT, ONLY: MSEC=>ISECDER,MFIRST=>IFRSTDER,IQUIT,IFIRSTEM,FIRSTEM
      USE NMPRD_INT, ONLY: MDVRES,ETASXI,NPDE_MODE,NOFIRSTDERCODE
      USE NMPRD_REAL, ONLY: DV_LOQ,CDF_L,DV_LAQ,CDF_LA
      USE NMBAYES_REAL, ONLY: ADDL_ACTUAL,ADDL_TIMEDIFF,ADDL_TIME
      USE PRLS01_INT, ONLY: ITASK_=>ITASK
      USE PRLS01_REAL, ONLY: STOP_TIME
      USE NMPRD_INT, ONLY: NEWL2
      USE PROCM_INT, ONLY: NEWIND=>PNEWIF                                       
      USE PROCM_REAL,ONLY: TSTATE
      USE PROCM_REAL,ONLY: A=>AMNT,DAETA,D2AETA
      IMPLICIT REAL(KIND=DPSIZE) (A-Z)                                        
      REAL(KIND=DPSIZE) :: EVTREC                                             
      SAVE
      INTEGER(KIND=ISIZE) :: ICALL,IDEF,IREV,NVNT,INDXS                       
      DIMENSION :: IDEF(*),THETA(*),EVTREC(IREV,*),INDXS(*)                   
      REAL(KIND=DPSIZE) :: G(GERD,*),HH(HERD,*)                               
      FIRSTEM=IFIRSTEM
      IF (ICALL <= 1) THEN                                                    
      CALL ASSOCNMPRD4
      IDEF(2)=-1
      IDEF(3)=001
      RETURN
      ENDIF
      IF (ICALL == 4) THEN
      IF (NEWL2 == 1) THEN
      CALL SIMEPS(EPS)
      IF (IQUIT == 1) RETURN
      ENDIF
      ENDIF
 !  level            0
      CMT=EVTREC(NVNT,010)
      Q000000=0.D0 
      Q000001=1.D0 
      IF(CMT == 2.D0)THEN 
      Q000000=1.D0 
      Q000001=0.D0 
      B000003=A(2)/Vd+EPS(001) 
      ELSE 
      B000011=DLOG10(A(5)) 
      B000012=B000011+EPS(002) 
      ENDIF 
      B000015=Q000000*B000003+Q000001*B000012 
!                      C000033 = DERIVATIVE OF B000015 W.R.T. EPS(001)
      C000033=Q000000 
!                      C000034 = DERIVATIVE OF B000015 W.R.T. EPS(002)
      C000034=Q000001 
      Y=B000015 
!                      C000035 = DERIVATIVE OF Y W.R.T. EPS(002)
      C000035=C000034 
!                      C000036 = DERIVATIVE OF Y W.R.T. EPS(001)
      C000036=C000033 
      IF (FIRSTEM == 1) THEN !1
      IF(CMT == 2.D0)THEN 
      B000004=1.D0/Vd 
!                      D000103 = DERIVATIVE OF B000003 W.R.T. ETA(007)
      D000103=B000004*DAETA(002,007) 
!                      D000104 = DERIVATIVE OF B000003 W.R.T. ETA(006)
      D000104=B000004*DAETA(002,006) 
!                      D000105 = DERIVATIVE OF B000003 W.R.T. ETA(005)
      D000105=B000004*DAETA(002,005) 
!                      D000106 = DERIVATIVE OF B000003 W.R.T. ETA(004)
      D000106=B000004*DAETA(002,004) 
!                      D000107 = DERIVATIVE OF B000003 W.R.T. ETA(003)
      D000107=B000004*DAETA(002,003) 
!                      D000108 = DERIVATIVE OF B000003 W.R.T. ETA(002)
      D000108=B000004*DAETA(002,002) 
!                      D000109 = DERIVATIVE OF B000003 W.R.T. ETA(001)
      D000109=B000004*DAETA(002,001) 
      B000005=-A(2)/Vd/Vd 
!                      D000110 = DERIVATIVE OF B000003 W.R.T. ETA(003)
      D000110=B000005*A000043+D000107 
      ELSE 
      B000013=.434294481903D0/A(5) 
!                      D000201 = DERIVATIVE OF B000011 W.R.T. ETA(007)
      D000201=B000013*DAETA(005,007) 
!                      D000202 = DERIVATIVE OF B000011 W.R.T. ETA(006)
      D000202=B000013*DAETA(005,006) 
!                      D000203 = DERIVATIVE OF B000011 W.R.T. ETA(005)
      D000203=B000013*DAETA(005,005) 
!                      D000204 = DERIVATIVE OF B000011 W.R.T. ETA(004)
      D000204=B000013*DAETA(005,004) 
!                      D000205 = DERIVATIVE OF B000011 W.R.T. ETA(003)
      D000205=B000013*DAETA(005,003) 
!                      D000206 = DERIVATIVE OF B000011 W.R.T. ETA(002)
      D000206=B000013*DAETA(005,002) 
!                      D000207 = DERIVATIVE OF B000011 W.R.T. ETA(001)
      D000207=B000013*DAETA(005,001) 
      ENDIF 
!                      D000306 = DERIVATIVE OF B000015 W.R.T. ETA(003)
      D000306=Q000000*D000110 
!                      D000307 = DERIVATIVE OF B000015 W.R.T. ETA(001)
      D000307=Q000000*D000109 
!                      D000308 = DERIVATIVE OF B000015 W.R.T. ETA(002)
      D000308=Q000000*D000108 
!                      D000309 = DERIVATIVE OF B000015 W.R.T. ETA(004)
      D000309=Q000000*D000106 
!                      D000310 = DERIVATIVE OF B000015 W.R.T. ETA(005)
      D000310=Q000000*D000105 
!                      D000311 = DERIVATIVE OF B000015 W.R.T. ETA(006)
      D000311=Q000000*D000104 
!                      D000312 = DERIVATIVE OF B000015 W.R.T. ETA(007)
      D000312=Q000000*D000103 
!                      D000313 = DERIVATIVE OF B000015 W.R.T. ETA(007)
      D000313=Q000001*D000201+D000312 
!                      D000314 = DERIVATIVE OF B000015 W.R.T. ETA(006)
      D000314=Q000001*D000202+D000311 
!                      D000315 = DERIVATIVE OF B000015 W.R.T. ETA(005)
      D000315=Q000001*D000203+D000310 
!                      D000316 = DERIVATIVE OF B000015 W.R.T. ETA(004)
      D000316=Q000001*D000204+D000309 
!                      D000317 = DERIVATIVE OF B000015 W.R.T. ETA(003)
      D000317=Q000001*D000205+D000306 
!                      D000318 = DERIVATIVE OF B000015 W.R.T. ETA(002)
      D000318=Q000001*D000206+D000308 
!                      D000319 = DERIVATIVE OF B000015 W.R.T. ETA(001)
      D000319=Q000001*D000207+D000307 
!                      D000404 = DERIVATIVE OF Y W.R.T. ETA(001)
      D000404=D000319 
!                      D000405 = DERIVATIVE OF Y W.R.T. ETA(002)
      D000405=D000318 
!                      D000406 = DERIVATIVE OF Y W.R.T. ETA(003)
      D000406=D000317 
!                      D000407 = DERIVATIVE OF Y W.R.T. ETA(004)
      D000407=D000316 
!                      D000408 = DERIVATIVE OF Y W.R.T. ETA(005)
      D000408=D000315 
!                      D000409 = DERIVATIVE OF Y W.R.T. ETA(006)
      D000409=D000314 
!                      D000410 = DERIVATIVE OF Y W.R.T. ETA(007)
      D000410=D000313 
      G(001,1)=D000404
      G(002,1)=D000405
      G(003,1)=D000406
      G(004,1)=D000407
      G(005,1)=D000408
      G(006,1)=D000409
      G(007,1)=D000410
      ENDIF !1
      HH(001,1)=C000036
      HH(002,1)=C000035
      IF (MSEC == 1) THEN
      IF(CMT == 2.D0)THEN 
      B000007=-1.D0/Vd/Vd 
!                      D000112 = DERIVATIVE OF B000004 W.R.T. ETA(003)
      D000112=B000007*A000043 
!                      D000113 = DERIVATIVE OF D000103 W.R.T. ETA(007)
      D000113=B000004*D2AETA(002,007,007) 
!                      D000114 = DERIVATIVE OF D000104 W.R.T. ETA(007)
      D000114=B000004*D2AETA(002,007,006) 
!                      D000115 = DERIVATIVE OF D000104 W.R.T. ETA(006)
      D000115=B000004*D2AETA(002,006,006) 
!                      D000116 = DERIVATIVE OF D000105 W.R.T. ETA(007)
      D000116=B000004*D2AETA(002,007,005) 
!                      D000117 = DERIVATIVE OF D000105 W.R.T. ETA(006)
      D000117=B000004*D2AETA(002,006,005) 
!                      D000118 = DERIVATIVE OF D000105 W.R.T. ETA(005)
      D000118=B000004*D2AETA(002,005,005) 
!                      D000119 = DERIVATIVE OF D000106 W.R.T. ETA(007)
      D000119=B000004*D2AETA(002,007,004) 
!                      D000120 = DERIVATIVE OF D000106 W.R.T. ETA(006)
      D000120=B000004*D2AETA(002,006,004) 
!                      D000121 = DERIVATIVE OF D000106 W.R.T. ETA(005)
      D000121=B000004*D2AETA(002,005,004) 
!                      D000122 = DERIVATIVE OF D000106 W.R.T. ETA(004)
      D000122=B000004*D2AETA(002,004,004) 
!                      D000123 = DERIVATIVE OF D000107 W.R.T. ETA(003)
      D000123=DAETA(002,003)*D000112 
!                      D000124 = DERIVATIVE OF D000107 W.R.T. ETA(007)
      D000124=B000004*D2AETA(002,007,003) 
!                      D000125 = DERIVATIVE OF D000107 W.R.T. ETA(006)
      D000125=B000004*D2AETA(002,006,003) 
!                      D000126 = DERIVATIVE OF D000107 W.R.T. ETA(005)
      D000126=B000004*D2AETA(002,005,003) 
!                      D000127 = DERIVATIVE OF D000107 W.R.T. ETA(004)
      D000127=B000004*D2AETA(002,004,003) 
!                      D000128 = DERIVATIVE OF D000107 W.R.T. ETA(003)
      D000128=B000004*D2AETA(002,003,003)+D000123 
!                      D000129 = DERIVATIVE OF D000108 W.R.T. ETA(003)
      D000129=DAETA(002,002)*D000112 
!                      D000130 = DERIVATIVE OF D000108 W.R.T. ETA(007)
      D000130=B000004*D2AETA(002,007,002) 
!                      D000131 = DERIVATIVE OF D000108 W.R.T. ETA(006)
      D000131=B000004*D2AETA(002,006,002) 
!                      D000132 = DERIVATIVE OF D000108 W.R.T. ETA(005)
      D000132=B000004*D2AETA(002,005,002) 
!                      D000133 = DERIVATIVE OF D000108 W.R.T. ETA(004)
      D000133=B000004*D2AETA(002,004,002) 
!                      D000134 = DERIVATIVE OF D000108 W.R.T. ETA(003)
      D000134=B000004*D2AETA(002,003,002)+D000129 
!                      D000135 = DERIVATIVE OF D000108 W.R.T. ETA(002)
      D000135=B000004*D2AETA(002,002,002) 
!                      D000136 = DERIVATIVE OF D000109 W.R.T. ETA(003)
      D000136=DAETA(002,001)*D000112 
!                      D000137 = DERIVATIVE OF D000109 W.R.T. ETA(007)
      D000137=B000004*D2AETA(002,007,001) 
!                      D000138 = DERIVATIVE OF D000109 W.R.T. ETA(006)
      D000138=B000004*D2AETA(002,006,001) 
!                      D000139 = DERIVATIVE OF D000109 W.R.T. ETA(005)
      D000139=B000004*D2AETA(002,005,001) 
!                      D000140 = DERIVATIVE OF D000109 W.R.T. ETA(004)
      D000140=B000004*D2AETA(002,004,001) 
!                      D000141 = DERIVATIVE OF D000109 W.R.T. ETA(003)
      D000141=B000004*D2AETA(002,003,001)+D000136 
!                      D000142 = DERIVATIVE OF D000109 W.R.T. ETA(002)
      D000142=B000004*D2AETA(002,002,001) 
!                      D000143 = DERIVATIVE OF D000109 W.R.T. ETA(001)
      D000143=B000004*D2AETA(002,001,001) 
      B000008=-1.D0/Vd/Vd 
!                      D000144 = DERIVATIVE OF B000005 W.R.T. ETA(007)
      D000144=B000008*DAETA(002,007) 
!                      D000145 = DERIVATIVE OF B000005 W.R.T. ETA(006)
      D000145=B000008*DAETA(002,006) 
!                      D000146 = DERIVATIVE OF B000005 W.R.T. ETA(005)
      D000146=B000008*DAETA(002,005) 
!                      D000147 = DERIVATIVE OF B000005 W.R.T. ETA(004)
      D000147=B000008*DAETA(002,004) 
!                      D000148 = DERIVATIVE OF B000005 W.R.T. ETA(003)
      D000148=B000008*DAETA(002,003) 
!                      D000149 = DERIVATIVE OF B000005 W.R.T. ETA(002)
      D000149=B000008*DAETA(002,002) 
!                      D000150 = DERIVATIVE OF B000005 W.R.T. ETA(001)
      D000150=B000008*DAETA(002,001) 
      B000009=A(2)/Vd/Vd/Vd 
!                      D000151 = DERIVATIVE OF B000005 W.R.T. ETA(003)
      D000151=B000009*A000043+D000148 
      B000010=A(2)/Vd/Vd/Vd 
!                      D000152 = DERIVATIVE OF B000005 W.R.T. ETA(003)
      D000152=B000010*A000043+D000151 
!                      D000153 = DERIVATIVE OF D000110 W.R.T. ETA(003)
      D000153=A000043*D000152 
!                      D000154 = DERIVATIVE OF D000110 W.R.T. ETA(004)
      D000154=A000043*D000147 
!                      D000155 = DERIVATIVE OF D000110 W.R.T. ETA(005)
      D000155=A000043*D000146 
!                      D000156 = DERIVATIVE OF D000110 W.R.T. ETA(006)
      D000156=A000043*D000145 
!                      D000157 = DERIVATIVE OF D000110 W.R.T. ETA(007)
      D000157=A000043*D000144 
!                      D000158 = DERIVATIVE OF D000110 W.R.T. ETA(003)
      D000158=B000005*A000045+D000153 
!                      D000159 = DERIVATIVE OF D000110 W.R.T. ETA(003)
      D000159=D000128+D000158 
!                      D000160 = DERIVATIVE OF D000110 W.R.T. ETA(004)
      D000160=D000127+D000154 
!                      D000161 = DERIVATIVE OF D000110 W.R.T. ETA(005)
      D000161=D000126+D000155 
!                      D000162 = DERIVATIVE OF D000110 W.R.T. ETA(006)
      D000162=D000125+D000156 
!                      D000163 = DERIVATIVE OF D000110 W.R.T. ETA(007)
      D000163=D000124+D000157 
      ELSE 
      B000014=-.434294481903D0/A(5)/A(5) 
!                      D000215 = DERIVATIVE OF B000013 W.R.T. ETA(007)
      D000215=B000014*DAETA(005,007) 
!                      D000216 = DERIVATIVE OF B000013 W.R.T. ETA(006)
      D000216=B000014*DAETA(005,006) 
!                      D000217 = DERIVATIVE OF B000013 W.R.T. ETA(005)
      D000217=B000014*DAETA(005,005) 
!                      D000218 = DERIVATIVE OF B000013 W.R.T. ETA(004)
      D000218=B000014*DAETA(005,004) 
!                      D000219 = DERIVATIVE OF B000013 W.R.T. ETA(003)
      D000219=B000014*DAETA(005,003) 
!                      D000220 = DERIVATIVE OF B000013 W.R.T. ETA(002)
      D000220=B000014*DAETA(005,002) 
!                      D000221 = DERIVATIVE OF B000013 W.R.T. ETA(001)
      D000221=B000014*DAETA(005,001) 
!                      D000222 = DERIVATIVE OF D000201 W.R.T. ETA(007)
      D000222=DAETA(005,007)*D000215 
!                      D000223 = DERIVATIVE OF D000201 W.R.T. ETA(007)
      D000223=B000013*D2AETA(005,007,007)+D000222 
!                      D000224 = DERIVATIVE OF D000202 W.R.T. ETA(006)
      D000224=DAETA(005,006)*D000216 
!                      D000225 = DERIVATIVE OF D000202 W.R.T. ETA(007)
      D000225=DAETA(005,006)*D000215 
!                      D000226 = DERIVATIVE OF D000202 W.R.T. ETA(007)
      D000226=B000013*D2AETA(005,007,006)+D000225 
!                      D000227 = DERIVATIVE OF D000202 W.R.T. ETA(006)
      D000227=B000013*D2AETA(005,006,006)+D000224 
!                      D000228 = DERIVATIVE OF D000203 W.R.T. ETA(005)
      D000228=DAETA(005,005)*D000217 
!                      D000229 = DERIVATIVE OF D000203 W.R.T. ETA(006)
      D000229=DAETA(005,005)*D000216 
!                      D000230 = DERIVATIVE OF D000203 W.R.T. ETA(007)
      D000230=DAETA(005,005)*D000215 
!                      D000231 = DERIVATIVE OF D000203 W.R.T. ETA(007)
      D000231=B000013*D2AETA(005,007,005)+D000230 
!                      D000232 = DERIVATIVE OF D000203 W.R.T. ETA(006)
      D000232=B000013*D2AETA(005,006,005)+D000229 
!                      D000233 = DERIVATIVE OF D000203 W.R.T. ETA(005)
      D000233=B000013*D2AETA(005,005,005)+D000228 
!                      D000234 = DERIVATIVE OF D000204 W.R.T. ETA(004)
      D000234=DAETA(005,004)*D000218 
!                      D000235 = DERIVATIVE OF D000204 W.R.T. ETA(005)
      D000235=DAETA(005,004)*D000217 
!                      D000236 = DERIVATIVE OF D000204 W.R.T. ETA(006)
      D000236=DAETA(005,004)*D000216 
!                      D000237 = DERIVATIVE OF D000204 W.R.T. ETA(007)
      D000237=DAETA(005,004)*D000215 
!                      D000238 = DERIVATIVE OF D000204 W.R.T. ETA(007)
      D000238=B000013*D2AETA(005,007,004)+D000237 
!                      D000239 = DERIVATIVE OF D000204 W.R.T. ETA(006)
      D000239=B000013*D2AETA(005,006,004)+D000236 
!                      D000240 = DERIVATIVE OF D000204 W.R.T. ETA(005)
      D000240=B000013*D2AETA(005,005,004)+D000235 
!                      D000241 = DERIVATIVE OF D000204 W.R.T. ETA(004)
      D000241=B000013*D2AETA(005,004,004)+D000234 
!                      D000242 = DERIVATIVE OF D000205 W.R.T. ETA(003)
      D000242=DAETA(005,003)*D000219 
!                      D000243 = DERIVATIVE OF D000205 W.R.T. ETA(004)
      D000243=DAETA(005,003)*D000218 
!                      D000244 = DERIVATIVE OF D000205 W.R.T. ETA(005)
      D000244=DAETA(005,003)*D000217 
!                      D000245 = DERIVATIVE OF D000205 W.R.T. ETA(006)
      D000245=DAETA(005,003)*D000216 
!                      D000246 = DERIVATIVE OF D000205 W.R.T. ETA(007)
      D000246=DAETA(005,003)*D000215 
!                      D000247 = DERIVATIVE OF D000205 W.R.T. ETA(007)
      D000247=B000013*D2AETA(005,007,003)+D000246 
!                      D000248 = DERIVATIVE OF D000205 W.R.T. ETA(006)
      D000248=B000013*D2AETA(005,006,003)+D000245 
!                      D000249 = DERIVATIVE OF D000205 W.R.T. ETA(005)
      D000249=B000013*D2AETA(005,005,003)+D000244 
!                      D000250 = DERIVATIVE OF D000205 W.R.T. ETA(004)
      D000250=B000013*D2AETA(005,004,003)+D000243 
!                      D000251 = DERIVATIVE OF D000205 W.R.T. ETA(003)
      D000251=B000013*D2AETA(005,003,003)+D000242 
!                      D000252 = DERIVATIVE OF D000206 W.R.T. ETA(002)
      D000252=DAETA(005,002)*D000220 
!                      D000253 = DERIVATIVE OF D000206 W.R.T. ETA(003)
      D000253=DAETA(005,002)*D000219 
!                      D000254 = DERIVATIVE OF D000206 W.R.T. ETA(004)
      D000254=DAETA(005,002)*D000218 
!                      D000255 = DERIVATIVE OF D000206 W.R.T. ETA(005)
      D000255=DAETA(005,002)*D000217 
!                      D000256 = DERIVATIVE OF D000206 W.R.T. ETA(006)
      D000256=DAETA(005,002)*D000216 
!                      D000257 = DERIVATIVE OF D000206 W.R.T. ETA(007)
      D000257=DAETA(005,002)*D000215 
!                      D000258 = DERIVATIVE OF D000206 W.R.T. ETA(007)
      D000258=B000013*D2AETA(005,007,002)+D000257 
!                      D000259 = DERIVATIVE OF D000206 W.R.T. ETA(006)
      D000259=B000013*D2AETA(005,006,002)+D000256 
!                      D000260 = DERIVATIVE OF D000206 W.R.T. ETA(005)
      D000260=B000013*D2AETA(005,005,002)+D000255 
!                      D000261 = DERIVATIVE OF D000206 W.R.T. ETA(004)
      D000261=B000013*D2AETA(005,004,002)+D000254 
!                      D000262 = DERIVATIVE OF D000206 W.R.T. ETA(003)
      D000262=B000013*D2AETA(005,003,002)+D000253 
!                      D000263 = DERIVATIVE OF D000206 W.R.T. ETA(002)
      D000263=B000013*D2AETA(005,002,002)+D000252 
!                      D000264 = DERIVATIVE OF D000207 W.R.T. ETA(001)
      D000264=DAETA(005,001)*D000221 
!                      D000265 = DERIVATIVE OF D000207 W.R.T. ETA(002)
      D000265=DAETA(005,001)*D000220 
!                      D000266 = DERIVATIVE OF D000207 W.R.T. ETA(003)
      D000266=DAETA(005,001)*D000219 
!                      D000267 = DERIVATIVE OF D000207 W.R.T. ETA(004)
      D000267=DAETA(005,001)*D000218 
!                      D000268 = DERIVATIVE OF D000207 W.R.T. ETA(005)
      D000268=DAETA(005,001)*D000217 
!                      D000269 = DERIVATIVE OF D000207 W.R.T. ETA(006)
      D000269=DAETA(005,001)*D000216 
!                      D000270 = DERIVATIVE OF D000207 W.R.T. ETA(007)
      D000270=DAETA(005,001)*D000215 
!                      D000271 = DERIVATIVE OF D000207 W.R.T. ETA(007)
      D000271=B000013*D2AETA(005,007,001)+D000270 
!                      D000272 = DERIVATIVE OF D000207 W.R.T. ETA(006)
      D000272=B000013*D2AETA(005,006,001)+D000269 
!                      D000273 = DERIVATIVE OF D000207 W.R.T. ETA(005)
      D000273=B000013*D2AETA(005,005,001)+D000268 
!                      D000274 = DERIVATIVE OF D000207 W.R.T. ETA(004)
      D000274=B000013*D2AETA(005,004,001)+D000267 
!                      D000275 = DERIVATIVE OF D000207 W.R.T. ETA(003)
      D000275=B000013*D2AETA(005,003,001)+D000266 
!                      D000276 = DERIVATIVE OF D000207 W.R.T. ETA(002)
      D000276=B000013*D2AETA(005,002,001)+D000265 
!                      D000277 = DERIVATIVE OF D000207 W.R.T. ETA(001)
      D000277=B000013*D2AETA(005,001,001)+D000264 
      ENDIF 
!                      D000320 = DERIVATIVE OF D000306 W.R.T. ETA(007)
      D000320=Q000000*D000163 
!                      D000321 = DERIVATIVE OF D000306 W.R.T. ETA(006)
      D000321=Q000000*D000162 
!                      D000322 = DERIVATIVE OF D000306 W.R.T. ETA(005)
      D000322=Q000000*D000161 
!                      D000323 = DERIVATIVE OF D000306 W.R.T. ETA(004)
      D000323=Q000000*D000160 
!                      D000324 = DERIVATIVE OF D000306 W.R.T. ETA(003)
      D000324=Q000000*D000159 
!                      D000325 = DERIVATIVE OF D000307 W.R.T. ETA(001)
      D000325=Q000000*D000143 
!                      D000326 = DERIVATIVE OF D000307 W.R.T. ETA(002)
      D000326=Q000000*D000142 
!                      D000327 = DERIVATIVE OF D000307 W.R.T. ETA(003)
      D000327=Q000000*D000141 
!                      D000328 = DERIVATIVE OF D000307 W.R.T. ETA(004)
      D000328=Q000000*D000140 
!                      D000329 = DERIVATIVE OF D000307 W.R.T. ETA(005)
      D000329=Q000000*D000139 
!                      D000330 = DERIVATIVE OF D000307 W.R.T. ETA(006)
      D000330=Q000000*D000138 
!                      D000331 = DERIVATIVE OF D000307 W.R.T. ETA(007)
      D000331=Q000000*D000137 
!                      D000332 = DERIVATIVE OF D000308 W.R.T. ETA(002)
      D000332=Q000000*D000135 
!                      D000333 = DERIVATIVE OF D000308 W.R.T. ETA(003)
      D000333=Q000000*D000134 
!                      D000334 = DERIVATIVE OF D000308 W.R.T. ETA(004)
      D000334=Q000000*D000133 
!                      D000335 = DERIVATIVE OF D000308 W.R.T. ETA(005)
      D000335=Q000000*D000132 
!                      D000336 = DERIVATIVE OF D000308 W.R.T. ETA(006)
      D000336=Q000000*D000131 
!                      D000337 = DERIVATIVE OF D000308 W.R.T. ETA(007)
      D000337=Q000000*D000130 
!                      D000338 = DERIVATIVE OF D000309 W.R.T. ETA(004)
      D000338=Q000000*D000122 
!                      D000339 = DERIVATIVE OF D000309 W.R.T. ETA(005)
      D000339=Q000000*D000121 
!                      D000340 = DERIVATIVE OF D000309 W.R.T. ETA(006)
      D000340=Q000000*D000120 
!                      D000341 = DERIVATIVE OF D000309 W.R.T. ETA(007)
      D000341=Q000000*D000119 
!                      D000342 = DERIVATIVE OF D000310 W.R.T. ETA(005)
      D000342=Q000000*D000118 
!                      D000343 = DERIVATIVE OF D000310 W.R.T. ETA(006)
      D000343=Q000000*D000117 
!                      D000344 = DERIVATIVE OF D000310 W.R.T. ETA(007)
      D000344=Q000000*D000116 
!                      D000345 = DERIVATIVE OF D000311 W.R.T. ETA(006)
      D000345=Q000000*D000115 
!                      D000346 = DERIVATIVE OF D000311 W.R.T. ETA(007)
      D000346=Q000000*D000114 
!                      D000347 = DERIVATIVE OF D000312 W.R.T. ETA(007)
      D000347=Q000000*D000113 
!                      D000348 = DERIVATIVE OF D000313 W.R.T. ETA(007)
      D000348=Q000001*D000223 
!                      D000349 = DERIVATIVE OF D000313 W.R.T. ETA(007)
      D000349=D000347+D000348 
!                      D000350 = DERIVATIVE OF D000314 W.R.T. ETA(007)
      D000350=Q000001*D000226 
!                      D000351 = DERIVATIVE OF D000314 W.R.T. ETA(006)
      D000351=Q000001*D000227 
!                      D000352 = DERIVATIVE OF D000314 W.R.T. ETA(007)
      D000352=D000346+D000350 
!                      D000353 = DERIVATIVE OF D000314 W.R.T. ETA(006)
      D000353=D000345+D000351 
!                      D000354 = DERIVATIVE OF D000315 W.R.T. ETA(007)
      D000354=Q000001*D000231 
!                      D000355 = DERIVATIVE OF D000315 W.R.T. ETA(006)
      D000355=Q000001*D000232 
!                      D000356 = DERIVATIVE OF D000315 W.R.T. ETA(005)
      D000356=Q000001*D000233 
!                      D000357 = DERIVATIVE OF D000315 W.R.T. ETA(007)
      D000357=D000344+D000354 
!                      D000358 = DERIVATIVE OF D000315 W.R.T. ETA(006)
      D000358=D000343+D000355 
!                      D000359 = DERIVATIVE OF D000315 W.R.T. ETA(005)
      D000359=D000342+D000356 
!                      D000360 = DERIVATIVE OF D000316 W.R.T. ETA(007)
      D000360=Q000001*D000238 
!                      D000361 = DERIVATIVE OF D000316 W.R.T. ETA(006)
      D000361=Q000001*D000239 
!                      D000362 = DERIVATIVE OF D000316 W.R.T. ETA(005)
      D000362=Q000001*D000240 
!                      D000363 = DERIVATIVE OF D000316 W.R.T. ETA(004)
      D000363=Q000001*D000241 
!                      D000364 = DERIVATIVE OF D000316 W.R.T. ETA(007)
      D000364=D000341+D000360 
!                      D000365 = DERIVATIVE OF D000316 W.R.T. ETA(006)
      D000365=D000340+D000361 
!                      D000366 = DERIVATIVE OF D000316 W.R.T. ETA(005)
      D000366=D000339+D000362 
!                      D000367 = DERIVATIVE OF D000316 W.R.T. ETA(004)
      D000367=D000338+D000363 
!                      D000368 = DERIVATIVE OF D000317 W.R.T. ETA(007)
      D000368=Q000001*D000247 
!                      D000369 = DERIVATIVE OF D000317 W.R.T. ETA(006)
      D000369=Q000001*D000248 
!                      D000370 = DERIVATIVE OF D000317 W.R.T. ETA(005)
      D000370=Q000001*D000249 
!                      D000371 = DERIVATIVE OF D000317 W.R.T. ETA(004)
      D000371=Q000001*D000250 
!                      D000372 = DERIVATIVE OF D000317 W.R.T. ETA(003)
      D000372=Q000001*D000251 
!                      D000373 = DERIVATIVE OF D000317 W.R.T. ETA(003)
      D000373=D000324+D000372 
!                      D000374 = DERIVATIVE OF D000317 W.R.T. ETA(004)
      D000374=D000323+D000371 
!                      D000375 = DERIVATIVE OF D000317 W.R.T. ETA(005)
      D000375=D000322+D000370 
!                      D000376 = DERIVATIVE OF D000317 W.R.T. ETA(006)
      D000376=D000321+D000369 
!                      D000377 = DERIVATIVE OF D000317 W.R.T. ETA(007)
      D000377=D000320+D000368 
!                      D000378 = DERIVATIVE OF D000318 W.R.T. ETA(007)
      D000378=Q000001*D000258 
!                      D000379 = DERIVATIVE OF D000318 W.R.T. ETA(006)
      D000379=Q000001*D000259 
!                      D000380 = DERIVATIVE OF D000318 W.R.T. ETA(005)
      D000380=Q000001*D000260 
!                      D000381 = DERIVATIVE OF D000318 W.R.T. ETA(004)
      D000381=Q000001*D000261 
!                      D000382 = DERIVATIVE OF D000318 W.R.T. ETA(003)
      D000382=Q000001*D000262 
!                      D000383 = DERIVATIVE OF D000318 W.R.T. ETA(002)
      D000383=Q000001*D000263 
!                      D000384 = DERIVATIVE OF D000318 W.R.T. ETA(007)
      D000384=D000337+D000378 
!                      D000385 = DERIVATIVE OF D000318 W.R.T. ETA(006)
      D000385=D000336+D000379 
!                      D000386 = DERIVATIVE OF D000318 W.R.T. ETA(005)
      D000386=D000335+D000380 
!                      D000387 = DERIVATIVE OF D000318 W.R.T. ETA(004)
      D000387=D000334+D000381 
!                      D000388 = DERIVATIVE OF D000318 W.R.T. ETA(003)
      D000388=D000333+D000382 
!                      D000389 = DERIVATIVE OF D000318 W.R.T. ETA(002)
      D000389=D000332+D000383 
!                      D000390 = DERIVATIVE OF D000319 W.R.T. ETA(007)
      D000390=Q000001*D000271 
!                      D000391 = DERIVATIVE OF D000319 W.R.T. ETA(006)
      D000391=Q000001*D000272 
!                      D000392 = DERIVATIVE OF D000319 W.R.T. ETA(005)
      D000392=Q000001*D000273 
!                      D000393 = DERIVATIVE OF D000319 W.R.T. ETA(004)
      D000393=Q000001*D000274 
!                      D000394 = DERIVATIVE OF D000319 W.R.T. ETA(003)
      D000394=Q000001*D000275 
!                      D000395 = DERIVATIVE OF D000319 W.R.T. ETA(002)
      D000395=Q000001*D000276 
!                      D000396 = DERIVATIVE OF D000319 W.R.T. ETA(001)
      D000396=Q000001*D000277 
!                      D000397 = DERIVATIVE OF D000319 W.R.T. ETA(007)
      D000397=D000331+D000390 
!                      D000398 = DERIVATIVE OF D000319 W.R.T. ETA(006)
      D000398=D000330+D000391 
!                      D000399 = DERIVATIVE OF D000319 W.R.T. ETA(005)
      D000399=D000329+D000392 
!                      D000400 = DERIVATIVE OF D000319 W.R.T. ETA(004)
      D000400=D000328+D000393 
!                      D000401 = DERIVATIVE OF D000319 W.R.T. ETA(003)
      D000401=D000327+D000394 
!                      D000402 = DERIVATIVE OF D000319 W.R.T. ETA(002)
      D000402=D000326+D000395 
!                      D000403 = DERIVATIVE OF D000319 W.R.T. ETA(001)
      D000403=D000325+D000396 
!                      D000411 = DERIVATIVE OF D000404 W.R.T. ETA(001)
      D000411=D000403 
!                      D000412 = DERIVATIVE OF D000404 W.R.T. ETA(002)
      D000412=D000402 
!                      D000413 = DERIVATIVE OF D000404 W.R.T. ETA(003)
      D000413=D000401 
!                      D000414 = DERIVATIVE OF D000404 W.R.T. ETA(004)
      D000414=D000400 
!                      D000415 = DERIVATIVE OF D000404 W.R.T. ETA(005)
      D000415=D000399 
!                      D000416 = DERIVATIVE OF D000404 W.R.T. ETA(006)
      D000416=D000398 
!                      D000417 = DERIVATIVE OF D000404 W.R.T. ETA(007)
      D000417=D000397 
!                      D000418 = DERIVATIVE OF D000405 W.R.T. ETA(002)
      D000418=D000389 
!                      D000419 = DERIVATIVE OF D000405 W.R.T. ETA(003)
      D000419=D000388 
!                      D000420 = DERIVATIVE OF D000405 W.R.T. ETA(004)
      D000420=D000387 
!                      D000421 = DERIVATIVE OF D000405 W.R.T. ETA(005)
      D000421=D000386 
!                      D000422 = DERIVATIVE OF D000405 W.R.T. ETA(006)
      D000422=D000385 
!                      D000423 = DERIVATIVE OF D000405 W.R.T. ETA(007)
      D000423=D000384 
!                      D000424 = DERIVATIVE OF D000406 W.R.T. ETA(007)
      D000424=D000377 
!                      D000425 = DERIVATIVE OF D000406 W.R.T. ETA(006)
      D000425=D000376 
!                      D000426 = DERIVATIVE OF D000406 W.R.T. ETA(005)
      D000426=D000375 
!                      D000427 = DERIVATIVE OF D000406 W.R.T. ETA(004)
      D000427=D000374 
!                      D000428 = DERIVATIVE OF D000406 W.R.T. ETA(003)
      D000428=D000373 
!                      D000429 = DERIVATIVE OF D000407 W.R.T. ETA(004)
      D000429=D000367 
!                      D000430 = DERIVATIVE OF D000407 W.R.T. ETA(005)
      D000430=D000366 
!                      D000431 = DERIVATIVE OF D000407 W.R.T. ETA(006)
      D000431=D000365 
!                      D000432 = DERIVATIVE OF D000407 W.R.T. ETA(007)
      D000432=D000364 
!                      D000433 = DERIVATIVE OF D000408 W.R.T. ETA(005)
      D000433=D000359 
!                      D000434 = DERIVATIVE OF D000408 W.R.T. ETA(006)
      D000434=D000358 
!                      D000435 = DERIVATIVE OF D000408 W.R.T. ETA(007)
      D000435=D000357 
!                      D000436 = DERIVATIVE OF D000409 W.R.T. ETA(006)
      D000436=D000353 
!                      D000437 = DERIVATIVE OF D000409 W.R.T. ETA(007)
      D000437=D000352 
!                      D000438 = DERIVATIVE OF D000410 W.R.T. ETA(007)
      D000438=D000349 
      G(001,002)=D000411
      G(002,002)=D000412
      G(002,003)=D000418
      G(003,002)=D000413
      G(003,003)=D000419
      G(003,004)=D000428
      G(004,002)=D000414
      G(004,003)=D000420
      G(004,004)=D000427
      G(004,005)=D000429
      G(005,002)=D000415
      G(005,003)=D000421
      G(005,004)=D000426
      G(005,005)=D000430
      G(005,006)=D000433
      G(006,002)=D000416
      G(006,003)=D000422
      G(006,004)=D000425
      G(006,005)=D000431
      G(006,006)=D000434
      G(006,007)=D000436
      G(007,002)=D000417
      G(007,003)=D000423
      G(007,004)=D000424
      G(007,005)=D000432
      G(007,006)=D000435
      G(007,007)=D000437
      G(007,008)=D000438
      ENDIF
      F=Y
      RETURN
      END
      SUBROUTINE TOL(NRD,ANRD,NRDC,ANRDC)
      USE SIZES,     ONLY: ISIZE
      INTEGER(KIND=ISIZE) :: NRD(0:*), ANRD(0:*), NRDC(0:*), ANRDC(0:*)
      NRD(1)=12
      RETURN
      END
      SUBROUTINE DES (A,P,T,DADT,IR,DA,DP,DT)                                 
      USE NMPRD4P
      USE SIZES,     ONLY: DPSIZE,ISIZE
      USE PRDIMS,    ONLY: GPRD,HPRD,GERD,HERD,GPKD
      USE NMPRD_INT, ONLY: IERPRD,IERPRDU,NETEXT,IQUIT                        
      USE NMPRD_CHAR,ONLY: ETEXT                                              
      USE NMPRD_INT, ONLY: MSEC=>ISECDER,MFIRST=>IFRSTDER,IFIRSTEM,IFIRSTEMJAC,FIRSTEM
      USE PRCOM_INT, ONLY: DES_DER=>MITER
      USE NMPRD_INT, ONLY: MDVRES,ETASXI,NPDE_MODE,NOFIRSTDERCODE
      USE NMPRD_REAL, ONLY: DV_LOQ,CDF_L,DV_LAQ,CDF_LA
      USE PRMOD_INT, ONLY: ICALL=>ICALLD,IDEFD,IDEFA
      USE PRDDESLVU, ONLY: PASTZERO
      USE NMBAYES_REAL, ONLY: ADDL_ACTUAL,ADDL_TIMEDIFF,ADDL_TIME
      IMPLICIT REAL(KIND=DPSIZE) (A-Z)                                        
      SAVE
      INTEGER(KIND=ISIZE) :: IR                                               
      DIMENSION :: A(*),P(*),DADT(*),DA(IR,*),DP(IR,*),DT(*)                  
      INTEGER(KIND=ISIZE) :: IFIRSTEMJACIN
      FIRSTEM=IFIRSTEM
      IF(DES_DER==1.OR.DES_DER==4) FIRSTEM=1
      IFIRSTEMJACIN=IFIRSTEMJAC
      IF(NOFIRSTDERCODE/=1) THEN
      IFIRSTEMJAC=FIRSTEM
      ELSE
      IFIRSTEMJAC=0
      ENDIF
      IF(IFIRSTEMJACIN==-2) RETURN
      IF (ICALL == 1) THEN
      CALL ASSOCNMPRD4
      IDEFD(1)=  0
      IDEFD(2)=0
      DA(   1,1)=0000014280
      DA(   2,1)=0000028441
      DA(   3,1)=0000028560
      DA(   4,1)=0000042840
      DA(   5,1)=0000043078
      DA(   6,1)=0000057001
      DA(   7,1)=0000057120
      DA(   8,1)=0000057239
      DA(   9,1)=0000071043
      DA(  10,1)=0000071281
      DA(  11,1)=0000071400
      DA(  12,1)=0000071045
      DA(  13,1)=0000071283
      DA(  14,1)=0000043081
      DA(  15,1)=0000057242
      DA(  16,1)=0000071047
      DA(  17,1)=0000042845
      DA(  18,1)=0000057006
      DA(  19,1)=0000014311
      DA(  20,1)=0000028472
      DA(  21,1)=0000028592
      DA(  22,1)=0000071080
      DA(  23,1)=0000071081
      DA(  24,1)=0000071082
      DA(  25,1)=0000071083
      DA(  26,1)=0000042873
      DA(  27,1)=0000042874
      DA(  28,1)=0000057034
      DA(  29,1)=0000057156
      DA(  30,1)=0000071318
      DA(  31,1)=0000071319
      DA(  32,1)=0000071320
      DA(  33,1)=0000071321
      DA(  34,1)=0000043111
      DA(  35,1)=0000057272
      DA(  36,1)=0000071441
      DA(  37,1)=0000000000
      DP(   1,1)=0000014280
      DP(   2,1)=0000028441
      DP(   3,1)=0000028560
      DP(   4,1)=0000042840
      DP(   5,1)=0000042959
      DP(   6,1)=0000043078
      DP(   7,1)=0000057001
      DP(   8,1)=0000057358
      DP(   9,1)=0000071638
      DP(  10,1)=0000071757
      DP(  11,1)=0000071876
      DP(  12,1)=0000071995
      DP(  13,1)=0000072114
      DP(  14,1)=0000071645
      DP(  15,1)=0000071764
      DP(  16,1)=0000071883
      DP(  17,1)=0000072002
      DP(  18,1)=0000071646
      DP(  19,1)=0000071765
      DP(  20,1)=0000071884
      DP(  21,1)=0000072003
      DP(  22,1)=0000071647
      DP(  23,1)=0000071766
      DP(  24,1)=0000071885
      DP(  25,1)=0000072004
      DP(  26,1)=0000071648
      DP(  27,1)=0000071767
      DP(  28,1)=0000071886
      DP(  29,1)=0000000000
      DT(   1)=0000000000
      RETURN
      ENDIF
 !  level            0
 !  level            0
      DADT(1)=-P(001)*A(1) 
      DADT(2)=P(001)*A(1)-P(002)*A(2) 
      B000001=P(003)*A(5)+P(004) 
      DADT(3)=P(005)-A(3)*B000001 
      DADT(4)=P(003)*A(5)*A(3)-P(006)*A(4) 
      IF(P(007) == 0.D0.AND.P(008) <= 0.D0)THEN 
      IERPRD=1
      ETEXT(1)='DES SUBROUTINE: ERROR IN COMPUTATION'
      ETEXT(2)='ATTEMPT TO COMPUTE 0**POWER WITH POWER<=0.'
      RETURN
      ENDIF 
      IF(P(007) <  0.D0)THEN 
      IERPRD=1
      ETEXT(1)='DES SUBROUTINE: ERROR IN COMPUTATION'
      ETEXT(2)='ATTEMPT TO COMPUTE BASE**POWER WITH BASE<0.'
      RETURN
      ENDIF 
      B000005=0.D0 
      IF(P(007) == 0.D0)THEN 
      B000005=1.D0 
      ENDIF 
      B000006=1.D0-B000005 
      B000007=P(007)+B000005 
      B000008=B000007**P(008) 
      B000009=B000006*B000008 
      B000010=A(2)/P(009) 
      B000011=A(2)/P(009) 
      IF(B000010 == 0.D0.AND.P(008) <= 0.D0)THEN 
      IERPRD=1
      ETEXT(1)='DES SUBROUTINE: ERROR IN COMPUTATION'
      ETEXT(2)='ATTEMPT TO COMPUTE 0**POWER WITH POWER<=0.'
      RETURN
      ENDIF 
      IF(B000010 <  0.D0)THEN 
      IERPRD=1
      ETEXT(1)='DES SUBROUTINE: ERROR IN COMPUTATION'
      ETEXT(2)='ATTEMPT TO COMPUTE BASE**POWER WITH BASE<0.'
      RETURN
      ENDIF 
      B000013=0.D0 
      IF(B000010 == 0.D0)THEN 
      B000013=1.D0 
      ENDIF 
      B000014=1.D0-B000013 
      B000015=B000010+B000013 
      B000016=B000015**P(008) 
      B000017=B000014*B000016 
      IF(B000011 == 0.D0.AND.P(008) <= 0.D0)THEN 
      IERPRD=1
      ETEXT(1)='DES SUBROUTINE: ERROR IN COMPUTATION'
      ETEXT(2)='ATTEMPT TO COMPUTE 0**POWER WITH POWER<=0.'
      RETURN
      ENDIF 
      IF(B000011 <  0.D0)THEN 
      IERPRD=1
      ETEXT(1)='DES SUBROUTINE: ERROR IN COMPUTATION'
      ETEXT(2)='ATTEMPT TO COMPUTE BASE**POWER WITH BASE<0.'
      RETURN
      ENDIF 
      B000019=0.D0 
      IF(B000011 == 0.D0)THEN 
      B000019=1.D0 
      ENDIF 
      B000020=1.D0-B000019 
      B000021=B000011+B000019 
      B000022=B000021**P(008) 
      B000023=B000020*B000022 
      B000018=B000023+B000009 
      B000024=1.D0-B000017/B000018 
      DADT(5)=P(010)*B000024*A(4)-P(011)*A(5) 
      IF (FIRSTEM == 1) THEN ! 1
!                      E000007 = DERIVATIVE OF DADT(1) W.R.T. A(001)
      E000007=-P(001) 
!                      F000081 = DERIVATIVE OF DADT(1) W.R.T. P(001)
      F000081=-A(1) 
!                      E000009 = DERIVATIVE OF DADT(2) W.R.T. A(001)
      E000009=P(001) 
!                      E000010 = DERIVATIVE OF DADT(2) W.R.T. A(002)
      E000010=-P(002) 
!                      F000082 = DERIVATIVE OF DADT(2) W.R.T. P(001)
      F000082=A(1) 
!                      F000083 = DERIVATIVE OF DADT(2) W.R.T. P(002)
      F000083=-A(2) 
!                      E000014 = DERIVATIVE OF DADT(3) W.R.T. A(003)
      E000014=-B000001 
!                      E000015 = DERIVATIVE OF DADT(3) W.R.T. A(005)
      E000015=-A(3)*P(003) 
!                      F000086 = DERIVATIVE OF DADT(3) W.R.T. P(005)
      F000086=1.D0 
!                      F000087 = DERIVATIVE OF DADT(3) W.R.T. P(004)
      F000087=-A(3) 
!                      F000088 = DERIVATIVE OF DADT(3) W.R.T. P(003)
      F000088=-A(3)*A(5) 
      B000002=P(003)*A(3) 
!                      E000022 = DERIVATIVE OF DADT(4) W.R.T. A(005)
      E000022=B000002 
      B000003=P(003)*A(5) 
!                      E000023 = DERIVATIVE OF DADT(4) W.R.T. A(003)
      E000023=B000003 
!                      E000024 = DERIVATIVE OF DADT(4) W.R.T. A(004)
      E000024=-P(006) 
      B000004=A(5)*A(3) 
!                      F000089 = DERIVATIVE OF DADT(4) W.R.T. P(003)
      F000089=B000004 
!                      F000090 = DERIVATIVE OF DADT(4) W.R.T. P(006)
      F000090=-A(4) 
      B000025=1.D0/P(009) 
      B000026=1.D0/P(009) 
      B000027=P(008)-1.D0 
      B000028=B000015**B000027 
      B000029=B000028*P(008) 
!                      E000037 = DERIVATIVE OF B000016 W.R.T. A(002)
      E000037=B000029*B000025 
!                      E000038 = DERIVATIVE OF B000017 W.R.T. A(002)
      E000038=B000014*E000037 
      B000030=P(008)-1.D0 
      B000031=B000021**B000030 
      B000032=B000031*P(008) 
!                      E000041 = DERIVATIVE OF B000022 W.R.T. A(002)
      E000041=B000032*B000026 
!                      E000042 = DERIVATIVE OF B000023 W.R.T. A(002)
      E000042=B000020*E000041 
      B000033=-1.D0/B000018 
!                      E000044 = DERIVATIVE OF B000024 W.R.T. A(002)
      E000044=B000033*E000038 
      B000034=B000017/B000018/B000018 
!                      E000045 = DERIVATIVE OF B000024 W.R.T. A(002)
      E000045=B000034*E000042+E000044 
      B000035=P(010)*A(4) 
!                      E000046 = DERIVATIVE OF DADT(5) W.R.T. A(002)
      E000046=B000035*E000045 
      B000036=P(010)*B000024 
!                      E000047 = DERIVATIVE OF DADT(5) W.R.T. A(004)
      E000047=B000036 
!                      E000048 = DERIVATIVE OF DADT(5) W.R.T. A(005)
      E000048=-P(011) 
      B000047=P(008)-1.D0 
      B000048=B000007**B000047 
      B000049=B000048*P(008) 
      B000050=DLOG(B000007) 
      B000051=B000050*B000008 
!                      F000094 = DERIVATIVE OF B000009 W.R.T. P(008)
      F000094=B000006*B000051 
!                      F000095 = DERIVATIVE OF B000009 W.R.T. P(007)
      F000095=B000006*B000049 
      B000052=-A(2)/P(009)/P(009) 
      B000053=-A(2)/P(009)/P(009) 
      B000054=P(008)-1.D0 
      B000055=B000015**B000054 
      B000056=B000055*P(008) 
!                      F000099 = DERIVATIVE OF B000016 W.R.T. P(009)
      F000099=B000056*B000052 
      B000057=DLOG(B000015) 
      B000058=B000057*B000016 
!                      F000101 = DERIVATIVE OF B000017 W.R.T. P(008)
      F000101=B000014*B000058 
!                      F000102 = DERIVATIVE OF B000017 W.R.T. P(009)
      F000102=B000014*F000099 
      B000059=P(008)-1.D0 
      B000060=B000021**B000059 
      B000061=B000060*P(008) 
!                      F000106 = DERIVATIVE OF B000022 W.R.T. P(009)
      F000106=B000061*B000053 
      B000062=DLOG(B000021) 
      B000063=B000062*B000022 
!                      F000108 = DERIVATIVE OF B000023 W.R.T. P(008)
      F000108=B000020*B000063 
!                      F000109 = DERIVATIVE OF B000023 W.R.T. P(009)
      F000109=B000020*F000106 
!                      F000113 = DERIVATIVE OF B000018 W.R.T. P(008)
      F000113=F000094+F000108 
      B000064=-1.D0/B000018 
!                      F000114 = DERIVATIVE OF B000024 W.R.T. P(008)
      F000114=B000064*F000101 
!                      F000115 = DERIVATIVE OF B000024 W.R.T. P(009)
      F000115=B000064*F000102 
      B000065=B000017/B000018/B000018 
!                      F000116 = DERIVATIVE OF B000024 W.R.T. P(008)
      F000116=B000065*F000113+F000114 
!                      F000117 = DERIVATIVE OF B000024 W.R.T. P(007)
      F000117=B000065*F000095 
!                      F000118 = DERIVATIVE OF B000024 W.R.T. P(009)
      F000118=B000065*F000109+F000115 
      B000066=B000024*A(4) 
!                      F000119 = DERIVATIVE OF DADT(5) W.R.T. P(010)
      F000119=B000066 
      B000067=P(010)*A(4) 
!                      F000120 = DERIVATIVE OF DADT(5) W.R.T. P(009)
      F000120=B000067*F000118 
!                      F000121 = DERIVATIVE OF DADT(5) W.R.T. P(007)
      F000121=B000067*F000117 
!                      F000122 = DERIVATIVE OF DADT(5) W.R.T. P(008)
      F000122=B000067*F000116 
!                      F000123 = DERIVATIVE OF DADT(5) W.R.T. P(011)
      F000123=-A(5) 
      ENDIF !1
      IF (MSEC == 1) THEN 
!                      E000008 = DERIVATIVE OF F000081 W.R.T. A(001)
      E000008=-1.D0 
!                      E000011 = DERIVATIVE OF F000082 W.R.T. A(001)
      E000011=1.D0 
!                      E000012 = DERIVATIVE OF F000083 W.R.T. A(002)
      E000012=-1.D0 
!                      E000016 = DERIVATIVE OF E000014 W.R.T. A(005)
      E000016=-P(003) 
!                      E000017 = DERIVATIVE OF E000015 W.R.T. A(003)
      E000017=-P(003) 
!                      E000018 = DERIVATIVE OF F000084 W.R.T. A(005)
      E000018=1.D0 
!                      E000019 = DERIVATIVE OF F000087 W.R.T. A(003)
      E000019=-1.D0 
!                      E000020 = DERIVATIVE OF F000088 W.R.T. A(003)
      E000020=-A(5) 
!                      E000021 = DERIVATIVE OF F000088 W.R.T. A(005)
      E000021=-A(3)*E000018 
!                      E000026 = DERIVATIVE OF E000022 W.R.T. A(003)
      E000026=P(003) 
!                      E000028 = DERIVATIVE OF E000023 W.R.T. A(005)
      E000028=P(003) 
!                      E000031 = DERIVATIVE OF F000089 W.R.T. A(003)
      E000031=A(5) 
!                      E000032 = DERIVATIVE OF F000089 W.R.T. A(005)
      E000032=A(3) 
!                      E000033 = DERIVATIVE OF F000090 W.R.T. A(004)
      E000033=-1.D0 
      B000037=B000027-1.D0 
      B000038=B000015**B000037 
      B000039=B000038*B000027 
!                      E000049 = DERIVATIVE OF B000028 W.R.T. A(002)
      E000049=B000039*B000025 
!                      E000050 = DERIVATIVE OF B000029 W.R.T. A(002)
      E000050=P(008)*E000049 
!                      E000051 = DERIVATIVE OF E000037 W.R.T. A(002)
      E000051=B000025*E000050 
!                      E000052 = DERIVATIVE OF E000038 W.R.T. A(002)
      E000052=B000014*E000051 
!                      E000053 = DERIVATIVE OF E000039 W.R.T. A(002)
      E000053=E000052 
      B000040=B000030-1.D0 
      B000041=B000021**B000040 
      B000042=B000041*B000030 
!                      E000054 = DERIVATIVE OF B000031 W.R.T. A(002)
      E000054=B000042*B000026 
!                      E000055 = DERIVATIVE OF B000032 W.R.T. A(002)
      E000055=P(008)*E000054 
!                      E000056 = DERIVATIVE OF E000041 W.R.T. A(002)
      E000056=B000026*E000055 
!                      E000057 = DERIVATIVE OF E000042 W.R.T. A(002)
      E000057=B000020*E000056 
!                      E000058 = DERIVATIVE OF E000043 W.R.T. A(002)
      E000058=E000057 
      B000043=1.D0/B000018/B000018 
!                      E000059 = DERIVATIVE OF B000033 W.R.T. A(002)
      E000059=B000043*E000042 
!                      E000060 = DERIVATIVE OF E000044 W.R.T. A(002)
      E000060=E000038*E000059 
!                      E000061 = DERIVATIVE OF E000044 W.R.T. A(002)
      E000061=B000033*E000053+E000060 
      B000044=1.D0/B000018/B000018 
!                      E000062 = DERIVATIVE OF B000034 W.R.T. A(002)
      E000062=B000044*E000038 
      B000045=-B000017/B000018/B000018/B000018 
!                      E000063 = DERIVATIVE OF B000034 W.R.T. A(002)
      E000063=B000045*E000042+E000062 
      B000046=-B000017/B000018/B000018/B000018 
!                      E000064 = DERIVATIVE OF B000034 W.R.T. A(002)
      E000064=B000046*E000042+E000063 
!                      E000065 = DERIVATIVE OF E000045 W.R.T. A(002)
      E000065=E000042*E000064 
!                      E000066 = DERIVATIVE OF E000045 W.R.T. A(002)
      E000066=B000034*E000058+E000065 
!                      E000067 = DERIVATIVE OF E000045 W.R.T. A(002)
      E000067=E000061+E000066 
!                      E000069 = DERIVATIVE OF E000046 W.R.T. A(004)
      E000069=E000045*P(010) 
!                      E000070 = DERIVATIVE OF E000046 W.R.T. A(002)
      E000070=B000035*E000067 
!                      E000071 = DERIVATIVE OF B000036 W.R.T. A(002)
      E000071=P(010)*E000045 
!                      E000072 = DERIVATIVE OF E000047 W.R.T. A(002)
      E000072=E000071 
      B000068=B000047-1.D0 
      B000069=B000007**B000068 
      B000070=B000069*B000047 
      B000071=DLOG(B000007) 
      B000072=B000071*B000048 
!                      F000127 = DERIVATIVE OF B000049 W.R.T. P(008)
      F000127=P(008)*B000072 
!                      F000128 = DERIVATIVE OF B000049 W.R.T. P(007)
      F000128=P(008)*B000070 
!                      F000129 = DERIVATIVE OF B000049 W.R.T. P(008)
      F000129=B000048+F000127 
!                      F000130 = DERIVATIVE OF F000092 W.R.T. P(008)
      F000130=F000129 
!                      F000131 = DERIVATIVE OF F000092 W.R.T. P(007)
      F000131=F000128 
      B000073=1.D0/B000007 
!                      F000133 = DERIVATIVE OF B000051 W.R.T. P(007)
      F000133=B000008*B000073 
!                      F000134 = DERIVATIVE OF B000051 W.R.T. P(008)
      F000134=B000050*B000051 
!                      F000135 = DERIVATIVE OF B000051 W.R.T. P(007)
      F000135=B000050*B000049+F000133 
!                      F000136 = DERIVATIVE OF F000093 W.R.T. P(007)
      F000136=F000135 
!                      F000137 = DERIVATIVE OF F000093 W.R.T. P(008)
      F000137=F000134 
!                      F000138 = DERIVATIVE OF F000094 W.R.T. P(008)
      F000138=B000006*F000137 
!                      F000139 = DERIVATIVE OF F000094 W.R.T. P(007)
      F000139=B000006*F000136 
!                      F000140 = DERIVATIVE OF F000095 W.R.T. P(007)
      F000140=B000006*F000131 
!                      F000141 = DERIVATIVE OF F000095 W.R.T. P(008)
      F000141=B000006*F000130 
      B000074=A(2)/P(009)/P(009)/P(009) 
      B000075=A(2)/P(009)/P(009)/P(009) 
!                      F000143 = DERIVATIVE OF B000052 W.R.T. P(009)
      F000143=B000075+B000074 
!                      F000144 = DERIVATIVE OF F000096 W.R.T. P(009)
      F000144=F000143 
      B000076=A(2)/P(009)/P(009)/P(009) 
      B000077=A(2)/P(009)/P(009)/P(009) 
!                      F000146 = DERIVATIVE OF B000053 W.R.T. P(009)
      F000146=B000077+B000076 
!                      F000147 = DERIVATIVE OF F000097 W.R.T. P(009)
      F000147=F000146 
!                      F000148 = DERIVATIVE OF F000098 W.R.T. P(009)
      F000148=F000144 
      B000078=B000054-1.D0 
      B000079=B000015**B000078 
      B000080=B000079*B000054 
!                      F000150 = DERIVATIVE OF B000055 W.R.T. P(009)
      F000150=B000080*B000052 
      B000081=DLOG(B000015) 
      B000082=B000081*B000055 
!                      F000152 = DERIVATIVE OF B000056 W.R.T. P(008)
      F000152=P(008)*B000082 
!                      F000153 = DERIVATIVE OF B000056 W.R.T. P(009)
      F000153=P(008)*F000150 
!                      F000154 = DERIVATIVE OF B000056 W.R.T. P(008)
      F000154=B000055+F000152 
!                      F000155 = DERIVATIVE OF F000099 W.R.T. P(008)
      F000155=B000052*F000154 
!                      F000156 = DERIVATIVE OF F000099 W.R.T. P(009)
      F000156=B000052*F000153 
!                      F000157 = DERIVATIVE OF F000099 W.R.T. P(009)
      F000157=B000056*F000148+F000156 
      B000083=1.D0/B000015 
!                      F000158 = DERIVATIVE OF B000057 W.R.T. P(009)
      F000158=B000083*B000052 
!                      F000159 = DERIVATIVE OF B000058 W.R.T. P(009)
      F000159=B000016*F000158 
!                      F000160 = DERIVATIVE OF B000058 W.R.T. P(008)
      F000160=B000057*B000058 
!                      F000161 = DERIVATIVE OF B000058 W.R.T. P(009)
      F000161=B000057*F000099+F000159 
!                      F000162 = DERIVATIVE OF F000100 W.R.T. P(009)
      F000162=F000161 
!                      F000163 = DERIVATIVE OF F000100 W.R.T. P(008)
      F000163=F000160 
!                      F000164 = DERIVATIVE OF F000101 W.R.T. P(008)
      F000164=B000014*F000163 
!                      F000165 = DERIVATIVE OF F000101 W.R.T. P(009)
      F000165=B000014*F000162 
!                      F000166 = DERIVATIVE OF F000102 W.R.T. P(009)
      F000166=B000014*F000157 
!                      F000167 = DERIVATIVE OF F000102 W.R.T. P(008)
      F000167=B000014*F000155 
!                      F000168 = DERIVATIVE OF F000103 W.R.T. P(008)
      F000168=F000167 
!                      F000169 = DERIVATIVE OF F000103 W.R.T. P(009)
      F000169=F000166 
!                      F000170 = DERIVATIVE OF F000104 W.R.T. P(009)
      F000170=F000165 
!                      F000171 = DERIVATIVE OF F000104 W.R.T. P(008)
      F000171=F000164 
!                      F000172 = DERIVATIVE OF F000105 W.R.T. P(009)
      F000172=F000147 
      B000084=B000059-1.D0 
      B000085=B000021**B000084 
      B000086=B000085*B000059 
!                      F000174 = DERIVATIVE OF B000060 W.R.T. P(009)
      F000174=B000086*B000053 
      B000087=DLOG(B000021) 
      B000088=B000087*B000060 
!                      F000176 = DERIVATIVE OF B000061 W.R.T. P(008)
      F000176=P(008)*B000088 
!                      F000177 = DERIVATIVE OF B000061 W.R.T. P(009)
      F000177=P(008)*F000174 
!                      F000178 = DERIVATIVE OF B000061 W.R.T. P(008)
      F000178=B000060+F000176 
!                      F000179 = DERIVATIVE OF F000106 W.R.T. P(008)
      F000179=B000053*F000178 
!                      F000180 = DERIVATIVE OF F000106 W.R.T. P(009)
      F000180=B000053*F000177 
!                      F000181 = DERIVATIVE OF F000106 W.R.T. P(009)
      F000181=B000061*F000172+F000180 
      B000089=1.D0/B000021 
!                      F000182 = DERIVATIVE OF B000062 W.R.T. P(009)
      F000182=B000089*B000053 
!                      F000183 = DERIVATIVE OF B000063 W.R.T. P(009)
      F000183=B000022*F000182 
!                      F000184 = DERIVATIVE OF B000063 W.R.T. P(008)
      F000184=B000062*B000063 
!                      F000185 = DERIVATIVE OF B000063 W.R.T. P(009)
      F000185=B000062*F000106+F000183 
!                      F000186 = DERIVATIVE OF F000107 W.R.T. P(009)
      F000186=F000185 
!                      F000187 = DERIVATIVE OF F000107 W.R.T. P(008)
      F000187=F000184 
!                      F000188 = DERIVATIVE OF F000108 W.R.T. P(008)
      F000188=B000020*F000187 
!                      F000189 = DERIVATIVE OF F000108 W.R.T. P(009)
      F000189=B000020*F000186 
!                      F000190 = DERIVATIVE OF F000109 W.R.T. P(009)
      F000190=B000020*F000181 
!                      F000191 = DERIVATIVE OF F000109 W.R.T. P(008)
      F000191=B000020*F000179 
!                      F000192 = DERIVATIVE OF F000110 W.R.T. P(008)
      F000192=F000191 
!                      F000193 = DERIVATIVE OF F000110 W.R.T. P(009)
      F000193=F000190 
!                      F000194 = DERIVATIVE OF F000111 W.R.T. P(009)
      F000194=F000189 
!                      F000195 = DERIVATIVE OF F000111 W.R.T. P(008)
      F000195=F000188 
!                      F000196 = DERIVATIVE OF F000112 W.R.T. P(008)
      F000196=F000141 
!                      F000197 = DERIVATIVE OF F000112 W.R.T. P(007)
      F000197=F000140 
!                      F000198 = DERIVATIVE OF F000113 W.R.T. P(007)
      F000198=F000139 
!                      F000199 = DERIVATIVE OF F000113 W.R.T. P(008)
      F000199=F000138 
!                      F000200 = DERIVATIVE OF F000113 W.R.T. P(008)
      F000200=F000195+F000199 
!                      F000201 = DERIVATIVE OF F000113 W.R.T. P(009)
      F000201=F000194 
      B000090=1.D0/B000018/B000018 
!                      F000202 = DERIVATIVE OF B000064 W.R.T. P(008)
      F000202=B000090*F000113 
!                      F000203 = DERIVATIVE OF B000064 W.R.T. P(007)
      F000203=B000090*F000095 
!                      F000204 = DERIVATIVE OF B000064 W.R.T. P(009)
      F000204=B000090*F000109 
!                      F000205 = DERIVATIVE OF F000114 W.R.T. P(009)
      F000205=F000101*F000204 
!                      F000206 = DERIVATIVE OF F000114 W.R.T. P(007)
      F000206=F000101*F000203 
!                      F000207 = DERIVATIVE OF F000114 W.R.T. P(008)
      F000207=F000101*F000202 
!                      F000208 = DERIVATIVE OF F000114 W.R.T. P(008)
      F000208=B000064*F000171+F000207 
!                      F000209 = DERIVATIVE OF F000114 W.R.T. P(009)
      F000209=B000064*F000170+F000205 
!                      F000210 = DERIVATIVE OF F000115 W.R.T. P(009)
      F000210=F000102*F000204 
!                      F000211 = DERIVATIVE OF F000115 W.R.T. P(007)
      F000211=F000102*F000203 
!                      F000212 = DERIVATIVE OF F000115 W.R.T. P(008)
      F000212=F000102*F000202 
!                      F000213 = DERIVATIVE OF F000115 W.R.T. P(009)
      F000213=B000064*F000169+F000210 
!                      F000214 = DERIVATIVE OF F000115 W.R.T. P(008)
      F000214=B000064*F000168+F000212 
      B000091=1.D0/B000018/B000018 
!                      F000215 = DERIVATIVE OF B000065 W.R.T. P(008)
      F000215=B000091*F000101 
!                      F000216 = DERIVATIVE OF B000065 W.R.T. P(009)
      F000216=B000091*F000102 
      B000092=-B000017/B000018/B000018/B000018 
!                      F000217 = DERIVATIVE OF B000065 W.R.T. P(008)
      F000217=B000092*F000113+F000215 
!                      F000218 = DERIVATIVE OF B000065 W.R.T. P(007)
      F000218=B000092*F000095 
!                      F000219 = DERIVATIVE OF B000065 W.R.T. P(009)
      F000219=B000092*F000109+F000216 
      B000093=-B000017/B000018/B000018/B000018 
!                      F000220 = DERIVATIVE OF B000065 W.R.T. P(008)
      F000220=B000093*F000113+F000217 
!                      F000221 = DERIVATIVE OF B000065 W.R.T. P(007)
      F000221=B000093*F000095+F000218 
!                      F000222 = DERIVATIVE OF B000065 W.R.T. P(009)
      F000222=B000093*F000109+F000219 
!                      F000223 = DERIVATIVE OF F000116 W.R.T. P(009)
      F000223=F000113*F000222 
!                      F000224 = DERIVATIVE OF F000116 W.R.T. P(007)
      F000224=F000113*F000221 
!                      F000225 = DERIVATIVE OF F000116 W.R.T. P(008)
      F000225=F000113*F000220 
!                      F000226 = DERIVATIVE OF F000116 W.R.T. P(009)
      F000226=B000065*F000201+F000223 
!                      F000227 = DERIVATIVE OF F000116 W.R.T. P(008)
      F000227=B000065*F000200+F000225 
!                      F000228 = DERIVATIVE OF F000116 W.R.T. P(007)
      F000228=B000065*F000198+F000224 
!                      F000229 = DERIVATIVE OF F000116 W.R.T. P(009)
      F000229=F000209+F000226 
!                      F000230 = DERIVATIVE OF F000116 W.R.T. P(008)
      F000230=F000208+F000227 
!                      F000231 = DERIVATIVE OF F000116 W.R.T. P(007)
      F000231=F000206+F000228 
!                      F000232 = DERIVATIVE OF F000117 W.R.T. P(009)
      F000232=F000095*F000222 
!                      F000233 = DERIVATIVE OF F000117 W.R.T. P(007)
      F000233=F000095*F000221 
!                      F000234 = DERIVATIVE OF F000117 W.R.T. P(008)
      F000234=F000095*F000220 
!                      F000235 = DERIVATIVE OF F000117 W.R.T. P(007)
      F000235=B000065*F000197+F000233 
!                      F000236 = DERIVATIVE OF F000117 W.R.T. P(008)
      F000236=B000065*F000196+F000234 
!                      F000237 = DERIVATIVE OF F000118 W.R.T. P(009)
      F000237=F000109*F000222 
!                      F000238 = DERIVATIVE OF F000118 W.R.T. P(007)
      F000238=F000109*F000221 
!                      F000239 = DERIVATIVE OF F000118 W.R.T. P(008)
      F000239=F000109*F000220 
!                      F000240 = DERIVATIVE OF F000118 W.R.T. P(009)
      F000240=B000065*F000193+F000237 
!                      F000241 = DERIVATIVE OF F000118 W.R.T. P(008)
      F000241=B000065*F000192+F000239 
!                      F000242 = DERIVATIVE OF F000118 W.R.T. P(008)
      F000242=F000214+F000241 
!                      F000243 = DERIVATIVE OF F000118 W.R.T. P(009)
      F000243=F000213+F000240 
!                      F000244 = DERIVATIVE OF F000118 W.R.T. P(007)
      F000244=F000211+F000238 
!                      F000245 = DERIVATIVE OF B000066 W.R.T. P(009)
      F000245=A(4)*F000118 
!                      F000246 = DERIVATIVE OF B000066 W.R.T. P(007)
      F000246=A(4)*F000117 
!                      F000247 = DERIVATIVE OF B000066 W.R.T. P(008)
      F000247=A(4)*F000116 
!                      F000248 = DERIVATIVE OF F000119 W.R.T. P(008)
      F000248=F000247 
!                      F000249 = DERIVATIVE OF F000119 W.R.T. P(007)
      F000249=F000246 
!                      F000250 = DERIVATIVE OF F000119 W.R.T. P(009)
      F000250=F000245 
!                      F000252 = DERIVATIVE OF F000120 W.R.T. P(010)
      F000252=F000118*A(4) 
!                      F000253 = DERIVATIVE OF F000120 W.R.T. P(007)
      F000253=B000067*F000244 
!                      F000254 = DERIVATIVE OF F000120 W.R.T. P(009)
      F000254=B000067*F000243 
!                      F000255 = DERIVATIVE OF F000120 W.R.T. P(008)
      F000255=B000067*F000242 
!                      F000256 = DERIVATIVE OF F000121 W.R.T. P(010)
      F000256=F000117*A(4) 
!                      F000257 = DERIVATIVE OF F000121 W.R.T. P(008)
      F000257=B000067*F000236 
!                      F000258 = DERIVATIVE OF F000121 W.R.T. P(007)
      F000258=B000067*F000235 
!                      F000259 = DERIVATIVE OF F000121 W.R.T. P(009)
      F000259=B000067*F000232 
!                      F000260 = DERIVATIVE OF F000122 W.R.T. P(010)
      F000260=F000116*A(4) 
!                      F000261 = DERIVATIVE OF F000122 W.R.T. P(007)
      F000261=B000067*F000231 
!                      F000262 = DERIVATIVE OF F000122 W.R.T. P(008)
      F000262=B000067*F000230 
!                      F000263 = DERIVATIVE OF F000122 W.R.T. P(009)
      F000263=B000067*F000229 
      B000094=-1.D0/P(009)/P(009) 
!                      E000074 = DERIVATIVE OF F000096 W.R.T. A(002)
      E000074=B000094 
      B000095=-1.D0/P(009)/P(009) 
!                      E000076 = DERIVATIVE OF F000097 W.R.T. A(002)
      E000076=B000095 
!                      E000077 = DERIVATIVE OF F000098 W.R.T. A(002)
      E000077=E000074 
      B000096=B000054-1.D0 
      B000097=B000015**B000096 
      B000098=B000097*B000054 
!                      E000078 = DERIVATIVE OF B000055 W.R.T. A(002)
      E000078=B000098*B000025 
!                      E000079 = DERIVATIVE OF B000056 W.R.T. A(002)
      E000079=P(008)*E000078 
!                      E000080 = DERIVATIVE OF F000099 W.R.T. A(002)
      E000080=B000052*E000079 
!                      E000081 = DERIVATIVE OF F000099 W.R.T. A(002)
      E000081=B000056*E000077+E000080 
      B000099=1.D0/B000015 
!                      E000082 = DERIVATIVE OF B000057 W.R.T. A(002)
      E000082=B000099*B000025 
!                      E000083 = DERIVATIVE OF B000058 W.R.T. A(002)
      E000083=B000016*E000082 
!                      E000084 = DERIVATIVE OF B000058 W.R.T. A(002)
      E000084=B000057*E000037+E000083 
!                      E000085 = DERIVATIVE OF F000100 W.R.T. A(002)
      E000085=E000084 
!                      E000086 = DERIVATIVE OF F000101 W.R.T. A(002)
      E000086=B000014*E000085 
!                      E000087 = DERIVATIVE OF F000102 W.R.T. A(002)
      E000087=B000014*E000081 
!                      E000088 = DERIVATIVE OF F000103 W.R.T. A(002)
      E000088=E000087 
!                      E000089 = DERIVATIVE OF F000104 W.R.T. A(002)
      E000089=E000086 
!                      E000090 = DERIVATIVE OF F000105 W.R.T. A(002)
      E000090=E000076 
      B000100=B000059-1.D0 
      B000101=B000021**B000100 
      B000102=B000101*B000059 
!                      E000091 = DERIVATIVE OF B000060 W.R.T. A(002)
      E000091=B000102*B000026 
!                      E000092 = DERIVATIVE OF B000061 W.R.T. A(002)
      E000092=P(008)*E000091 
!                      E000093 = DERIVATIVE OF F000106 W.R.T. A(002)
      E000093=B000053*E000092 
!                      E000094 = DERIVATIVE OF F000106 W.R.T. A(002)
      E000094=B000061*E000090+E000093 
      B000103=1.D0/B000021 
!                      E000095 = DERIVATIVE OF B000062 W.R.T. A(002)
      E000095=B000103*B000026 
!                      E000096 = DERIVATIVE OF B000063 W.R.T. A(002)
      E000096=B000022*E000095 
!                      E000097 = DERIVATIVE OF B000063 W.R.T. A(002)
      E000097=B000062*E000041+E000096 
!                      E000098 = DERIVATIVE OF F000107 W.R.T. A(002)
      E000098=E000097 
!                      E000099 = DERIVATIVE OF F000108 W.R.T. A(002)
      E000099=B000020*E000098 
!                      E000100 = DERIVATIVE OF F000109 W.R.T. A(002)
      E000100=B000020*E000094 
!                      E000101 = DERIVATIVE OF F000110 W.R.T. A(002)
      E000101=E000100 
!                      E000102 = DERIVATIVE OF F000111 W.R.T. A(002)
      E000102=E000099 
!                      E000103 = DERIVATIVE OF F000113 W.R.T. A(002)
      E000103=E000102 
      B000104=1.D0/B000018/B000018 
!                      E000104 = DERIVATIVE OF B000064 W.R.T. A(002)
      E000104=B000104*E000042 
!                      E000105 = DERIVATIVE OF F000114 W.R.T. A(002)
      E000105=F000101*E000104 
!                      E000106 = DERIVATIVE OF F000114 W.R.T. A(002)
      E000106=B000064*E000089+E000105 
!                      E000107 = DERIVATIVE OF F000115 W.R.T. A(002)
      E000107=F000102*E000104 
!                      E000108 = DERIVATIVE OF F000115 W.R.T. A(002)
      E000108=B000064*E000088+E000107 
      B000105=1.D0/B000018/B000018 
!                      E000109 = DERIVATIVE OF B000065 W.R.T. A(002)
      E000109=B000105*E000038 
      B000106=-B000017/B000018/B000018/B000018 
!                      E000110 = DERIVATIVE OF B000065 W.R.T. A(002)
      E000110=B000106*E000042+E000109 
      B000107=-B000017/B000018/B000018/B000018 
!                      E000111 = DERIVATIVE OF B000065 W.R.T. A(002)
      E000111=B000107*E000042+E000110 
!                      E000112 = DERIVATIVE OF F000116 W.R.T. A(002)
      E000112=F000113*E000111 
!                      E000113 = DERIVATIVE OF F000116 W.R.T. A(002)
      E000113=B000065*E000103+E000112 
!                      E000114 = DERIVATIVE OF F000116 W.R.T. A(002)
      E000114=E000106+E000113 
!                      E000115 = DERIVATIVE OF F000117 W.R.T. A(002)
      E000115=F000095*E000111 
!                      E000116 = DERIVATIVE OF F000118 W.R.T. A(002)
      E000116=F000109*E000111 
!                      E000117 = DERIVATIVE OF F000118 W.R.T. A(002)
      E000117=B000065*E000101+E000116 
!                      E000118 = DERIVATIVE OF F000118 W.R.T. A(002)
      E000118=E000108+E000117 
!                      E000119 = DERIVATIVE OF B000066 W.R.T. A(002)
      E000119=A(4)*E000045 
!                      E000121 = DERIVATIVE OF F000119 W.R.T. A(004)
      E000121=B000024 
!                      E000122 = DERIVATIVE OF F000119 W.R.T. A(002)
      E000122=E000119 
!                      E000124 = DERIVATIVE OF F000120 W.R.T. A(004)
      E000124=F000118*P(010) 
!                      E000125 = DERIVATIVE OF F000120 W.R.T. A(002)
      E000125=B000067*E000118 
!                      E000126 = DERIVATIVE OF F000121 W.R.T. A(004)
      E000126=F000117*P(010) 
!                      E000127 = DERIVATIVE OF F000121 W.R.T. A(002)
      E000127=B000067*E000115 
!                      E000128 = DERIVATIVE OF F000122 W.R.T. A(004)
      E000128=F000116*P(010) 
!                      E000129 = DERIVATIVE OF F000122 W.R.T. A(002)
      E000129=B000067*E000114 
!                      E000130 = DERIVATIVE OF F000123 W.R.T. A(005)
      E000130=-1.D0 
      ENDIF !msec
      IF (FIRSTEM == 1) THEN !2
      DA(   1,1)=E000007
      DA(   2,1)=E000009
      DA(   3,1)=E000010
      DA(   4,1)=E000014
      DA(   5,1)=E000015
      DA(   6,1)=E000023
      DA(   7,1)=E000024
      DA(   8,1)=E000022
      DA(   9,1)=E000046
      DA(  10,1)=E000047
      DA(  11,1)=E000048
      DP(   1,1)=F000081
      DP(   2,1)=F000082
      DP(   3,1)=F000083
      DP(   4,1)=F000088
      DP(   5,1)=F000087
      DP(   6,1)=F000086
      DP(   7,1)=F000089
      DP(   8,1)=F000090
      DP(   9,1)=F000121
      DP(  10,1)=F000122
      DP(  11,1)=F000120
      DP(  12,1)=F000119
      DP(  13,1)=F000123
      ENDIF !2
      IF (MSEC == 1) THEN
      DA(  12,1)=E000070
      DA(  13,1)=E000072
      DA(  14,1)=E000017
      DA(  15,1)=E000026
      DA(  16,1)=E000069
      DA(  17,1)=E000016
      DA(  18,1)=E000028
      DA(  19,1)=E000008
      DA(  20,1)=E000011
      DA(  21,1)=E000012
      DA(  22,1)=E000127
      DA(  23,1)=E000129
      DA(  24,1)=E000125
      DA(  25,1)=E000122
      DA(  26,1)=E000020
      DA(  27,1)=E000019
      DA(  28,1)=E000031
      DA(  29,1)=E000033
      DA(  30,1)=E000126
      DA(  31,1)=E000128
      DA(  32,1)=E000124
      DA(  33,1)=E000121
      DA(  34,1)=E000021
      DA(  35,1)=E000032
      DA(  36,1)=E000130
      DP(  14,1)=F000258
      DP(  15,1)=F000261
      DP(  16,1)=F000253
      DP(  17,1)=F000249
      DP(  18,1)=F000257
      DP(  19,1)=F000262
      DP(  20,1)=F000255
      DP(  21,1)=F000248
      DP(  22,1)=F000259
      DP(  23,1)=F000263
      DP(  24,1)=F000254
      DP(  25,1)=F000250
      DP(  26,1)=F000256
      DP(  27,1)=F000260
      DP(  28,1)=F000252
      ENDIF
      RETURN
      END
      SUBROUTINE FSIZESR(NAME_FSIZES,F_SIZES)
      USE SIZES, ONLY: ISIZE
      INTEGER(KIND=ISIZE), DIMENSION(*) :: F_SIZES
      CHARACTER(LEN=*),    DIMENSION(*) :: NAME_FSIZES
      NAME_FSIZES(01)='LTH'; F_SIZES(01)=7
      NAME_FSIZES(02)='LVR'; F_SIZES(02)=9
      NAME_FSIZES(03)='LVR2'; F_SIZES(03)=0
      NAME_FSIZES(04)='LPAR'; F_SIZES(04)=38
      NAME_FSIZES(05)='LPAR3'; F_SIZES(05)=0
      NAME_FSIZES(06)='NO'; F_SIZES(06)=0
      NAME_FSIZES(07)='MMX'; F_SIZES(07)=1
      NAME_FSIZES(08)='LNP4'; F_SIZES(08)=0
      NAME_FSIZES(09)='LSUPP'; F_SIZES(09)=1
      NAME_FSIZES(10)='LIM7'; F_SIZES(10)=0
      NAME_FSIZES(11)='LWS3'; F_SIZES(11)=0
      NAME_FSIZES(12)='MAXIDS'; F_SIZES(12)=1
      NAME_FSIZES(13)='LIM1'; F_SIZES(13)=0
      NAME_FSIZES(14)='LIM2'; F_SIZES(14)=0
      NAME_FSIZES(15)='LIM3'; F_SIZES(15)=0
      NAME_FSIZES(16)='LIM4'; F_SIZES(16)=0
      NAME_FSIZES(17)='LIM5'; F_SIZES(17)=0
      NAME_FSIZES(18)='LIM6'; F_SIZES(18)=0
      NAME_FSIZES(19)='LIM8'; F_SIZES(19)=0
      NAME_FSIZES(20)='LIM10'; F_SIZES(20)=0
      NAME_FSIZES(21)='LIM11'; F_SIZES(21)=0
      NAME_FSIZES(22)='LIM13'; F_SIZES(22)=0
      NAME_FSIZES(23)='LIM15'; F_SIZES(23)=0
      NAME_FSIZES(24)='LIM16'; F_SIZES(24)=0
      NAME_FSIZES(25)='MAXRECID'; F_SIZES(25)=0
      NAME_FSIZES(26)='PC'; F_SIZES(26)=0
      NAME_FSIZES(27)='PCT'; F_SIZES(27)=1
      NAME_FSIZES(28)='PIR'; F_SIZES(28)=37
      NAME_FSIZES(29)='PD'; F_SIZES(29)=10
      NAME_FSIZES(30)='PAL'; F_SIZES(30)=0
      NAME_FSIZES(31)='MAXFCN'; F_SIZES(31)=0
      NAME_FSIZES(32)='MAXIC'; F_SIZES(32)=0
      NAME_FSIZES(33)='PG'; F_SIZES(33)=0
      NAME_FSIZES(34)='NPOPMIXMAX'; F_SIZES(34)=0
      NAME_FSIZES(35)='MAXOMEG'; F_SIZES(35)=7
      NAME_FSIZES(36)='MAXPTHETA'; F_SIZES(36)=10
      NAME_FSIZES(37)='MAXITER'; F_SIZES(37)=20
      NAME_FSIZES(38)='ISAMPLEMAX'; F_SIZES(38)=0
      NAME_FSIZES(39)='DIMTMP'; F_SIZES(39)=0
      NAME_FSIZES(40)='DIMCNS'; F_SIZES(40)=0
      NAME_FSIZES(41)='DIMNEW'; F_SIZES(41)=0
      NAME_FSIZES(42)='PDT'; F_SIZES(42)=1
      NAME_FSIZES(43)='LADD_MAX'; F_SIZES(43)=0
      NAME_FSIZES(44)='MAXSIDL'; F_SIZES(44)=0
      NAME_FSIZES(45)='NTT'; F_SIZES(45)=7
      NAME_FSIZES(46)='NOMEG'; F_SIZES(46)=7
      NAME_FSIZES(47)='NSIGM'; F_SIZES(47)=2
      NAME_FSIZES(48)='PPDT'; F_SIZES(48)=1
      NAME_FSIZES(49)='PAST_SIZE'; F_SIZES(49)=0
      NAME_FSIZES(50)='MAXNRDS'; F_SIZES(50)=0
      RETURN
      END SUBROUTINE FSIZESR
      SUBROUTINE MUMODEL2(THETA,MU_,ICALL,IDEF,NEWIND,&
      EVTREC,DATREC,IREV,NVNT,INDXS,F,G,H,IRGG,GG,NETAS)
      USE NMPRD4P
      USE SIZES,     ONLY: DPSIZE,ISIZE
      USE PRDIMS,    ONLY: GPRD,HPRD,GERD,HERD,GPKD
      USE NMBAYES_REAL,    ONLY: PRIOR_INFO
      USE PRDATA, ONLY: MXSTEP=>MXSTP01
      USE PRLS01_INT, ONLY: ITASK_=>ITASK
      USE PRLS01_REAL, ONLY: STOP_TIME
      USE NMBAYES_REAL, ONLY: ADDL_ACTUAL,ADDL_TIMEDIFF,ADDL_TIME
      USE NMPRD_INT, ONLY: PRED_IGNORE_DATA,PRED_IGNORE_DATA_TEST
      USE PRCOM_INT, ONLY: DES_DER=>MITER
      USE NMPRD_REAL,ONLY: ETA,EPS
      USE NMPRD_INT, ONLY: MSEC=>ISECDER,MFIRST=>IFRSTDER,COMACT,COMSAV,IFIRSTEM,FIRSTEM
      USE NMPRD_INT, ONLY: MDVRES,ETASXI,NPDE_MODE,NOFIRSTDERCODE
      USE NMPRD_REAL, ONLY: DV_LOQ,CDF_L,DV_LAQ,CDF_LA
      USE NMPRD_INT, ONLY: IQUIT
      USE PROCM_INT, ONLY: A_0FLG,A_UFLG
      USE PRMOD_REAL,ONLY: A_0,DA_0,D2A_0
      USE NMBAYES_REAL, ONLY: LDF
      USE PRDDESLVU, ONLY: PASTZERO
      IMPLICIT REAL(KIND=DPSIZE) (A-Z)
      REAL(KIND=DPSIZE)   :: MU_(*)
      INTEGER NEWIND
      REAL(KIND=DPSIZE) :: EVTREC
      SAVE
      INTEGER(KIND=ISIZE) :: ICALL,IDEF,IREV,NVNT,INDXS,IRGG,NETAS
      DIMENSION :: IDEF(7,*),THETA(*),EVTREC(IREV,*),INDXS(*),GG(IRGG,GPKD+1,*)
      FIRSTEM=IFIRSTEM
      IF (ICALL <= 1) THEN
      CALL ASSOCNMPRD4
      IDEF(   1,0001)=  -9
      IDEF(   1,0002)=  -1
      IDEF(   1,0003)=   1
      IDEF(   1,0004)=   0
      IDEF(   2,0003)=   0
      IDEF(   2,0004)=   0
      CALL GETETA(ETA)
      IF (IQUIT == 1) RETURN
      RETURN
      ENDIF
      IF (NEWIND /= 2) THEN
      IF (ICALL == 4) THEN
      CALL SIMETA(ETA)
      ELSE
      CALL GETETA(ETA)
      ENDIF
      IF (IQUIT == 1) RETURN
      ENDIF
 !  level            0
      MU_1=DLOG(THETA(001))
      MU_(001)=MU_1
      MU_2=DLOG(THETA(002))
      MU_(002)=MU_2
      MU_3=DLOG(THETA(003))
      MU_(003)=MU_3
      MU_4=DLOG(THETA(004))
      MU_(004)=MU_4
      MU_5=DLOG(THETA(005))
      MU_(005)=MU_5
      MU_6=DLOG(THETA(006))
      MU_(006)=MU_6
      MU_7=DLOG(THETA(007))
      MU_(007)=MU_7
       RETURN
      B000001=MU_1+ETA(001)
      B000002=DEXP(B000001)
      ka=B000002
      B000003=MU_2+ETA(002)
      B000004=DEXP(B000003)
      ke=B000004
      B000005=MU_3+ETA(003)
      B000006=DEXP(B000005)
      Vd=B000006
      B000007=MU_4+ETA(004)
      B000008=DEXP(B000007)
      EC50=B000008
      B000009=MU_5+ETA(005)
      B000010=DEXP(B000009)
      n=B000010
      B000011=MU_6+ETA(006)
      B000012=DEXP(B000011)
      delta=B000012
      B000013=MU_7+ETA(007)
      B000014=DEXP(B000013)
      c=B000014
      q=100.D0
      d=0.001D0
      e=0.0000001D0
      s=20000.D0
      IF(A_0FLG == 1)THEN
      A_0(1)=0.D0
      ENDIF
      IF(A_0FLG == 1)THEN
      A_0(2)=0.D0
      ENDIF
      IF(A_0FLG == 1)THEN
      B000015=c*delta
      B000016=q*e
      A_0(3)=B000015/B000016
      ENDIF
      IF(A_0FLG == 1)THEN
      B000020=s*e*q-d*c*delta
      B000021=q*delta*e
      A_0(4)=B000020/B000021
      ENDIF
      IF(A_0FLG == 1)THEN
      B000036=s*e*q-d*c*delta
      B000037=c*delta*e
      A_0(5)=B000036/B000037
      B000045=delta
      ENDIF
      P000001=ka
      P000002=ke
      P000003=e
      P000004=d
      P000005=s
      P000006=delta
      P000007=EC50
      P000008=n
      P000009=Vd
      P000010=q
      P000011=c
      IF (FIRSTEM == 1) THEN
!                      A000043 = DERIVATIVE OF Vd W.R.T. ETA(003)
      A000043=B000006
      IF(A_0FLG == 1)THEN
!                      A000066 = DERIVATIVE OF B000015 W.R.T. ETA(007)
      A000066=delta*B000014
!                      A000067 = DERIVATIVE OF B000015 W.R.T. ETA(006)
      A000067=c*B000012
      B000019=1.D0/B000016
!                      A000070 = DERIVATIVE OF A_0(3) W.R.T. ETA(006)
      A000070=B000019*A000067
!                      A000071 = DERIVATIVE OF A_0(3) W.R.T. ETA(007)
      A000071=B000019*A000066
      DA_0(003,006)=A000070
      DA_0(003,007)=A000071
      ENDIF
      IF(A_0FLG == 1)THEN
      B000022=-d*delta
!                      A000082 = DERIVATIVE OF B000020 W.R.T. ETA(007)
      A000082=B000022*B000014
      B000023=-d*c
!                      A000083 = DERIVATIVE OF B000020 W.R.T. ETA(006)
      A000083=B000023*B000012
      B000027=q*e
!                      A000086 = DERIVATIVE OF B000021 W.R.T. ETA(006)
      A000086=B000027*B000012
      B000030=1.D0/B000021
!                      A000087 = DERIVATIVE OF A_0(4) W.R.T. ETA(006)
      A000087=B000030*A000083
!                      A000088 = DERIVATIVE OF A_0(4) W.R.T. ETA(007)
      A000088=B000030*A000082
      B000031=-B000020/B000021/B000021
!                      A000089 = DERIVATIVE OF A_0(4) W.R.T. ETA(006)
      A000089=B000031*A000086+A000087
      DA_0(004,006)=A000089
      DA_0(004,007)=A000088
      ENDIF
      IF(A_0FLG == 1)THEN
      B000038=-d*delta
!                      A000114 = DERIVATIVE OF B000036 W.R.T. ETA(007)
      A000114=B000038*B000014
      B000039=-d*c
!                      A000115 = DERIVATIVE OF B000036 W.R.T. ETA(006)
      A000115=B000039*B000012
      B000043=B000045*e
!                      A000118 = DERIVATIVE OF B000037 W.R.T. ETA(007)
      A000118=B000043*B000014
      B000044=c*e
!                      A000119 = DERIVATIVE OF B000037 W.R.T. ETA(006)
      A000119=B000044*B000012
      B000048=1.D0/B000037
!                      A000122 = DERIVATIVE OF A_0(5) W.R.T. ETA(006)
      A000122=B000048*A000115
!                      A000123 = DERIVATIVE OF A_0(5) W.R.T. ETA(007)
      A000123=B000048*A000114
      B000049=-B000036/B000037/B000037
!                      A000124 = DERIVATIVE OF A_0(5) W.R.T. ETA(006)
      A000124=B000049*A000119+A000122
!                      A000125 = DERIVATIVE OF A_0(5) W.R.T. ETA(007)
      A000125=B000049*A000118+A000123
      DA_0(005,006)=A000124
      DA_0(005,007)=A000125
      ENDIF
!                      A000167 = DERIVATIVE OF P000001 W.R.T. ETA(001)
      A000167=B000002
!                      A000169 = DERIVATIVE OF P000002 W.R.T. ETA(002)
      A000169=B000004
!                      A000171 = DERIVATIVE OF P000006 W.R.T. ETA(006)
      A000171=B000012
!                      A000173 = DERIVATIVE OF P000007 W.R.T. ETA(004)
      A000173=B000008
!                      A000175 = DERIVATIVE OF P000008 W.R.T. ETA(005)
      A000175=B000010
!                      A000177 = DERIVATIVE OF P000009 W.R.T. ETA(003)
      A000177=A000043
!                      A000179 = DERIVATIVE OF P000011 W.R.T. ETA(007)
      A000179=B000014
      GG(0001,1,1)=P000001
      GG(0001,0002,1)=A000167
      GG(0002,1,1)=P000002
      GG(0002,0003,1)=A000169
      GG(0003,1,1)=P000003
      GG(0004,1,1)=P000004
      GG(0005,1,1)=P000005
      GG(0006,1,1)=P000006
      GG(0006,0007,1)=A000171
      GG(0007,1,1)=P000007
      GG(0007,0005,1)=A000173
      GG(0008,1,1)=P000008
      GG(0008,0006,1)=A000175
      GG(0009,1,1)=P000009
      GG(0009,0004,1)=A000177
      GG(0010,1,1)=P000010
      GG(0011,1,1)=P000011
      GG(0011,0008,1)=A000179
      ELSE
      GG(0001,1,1)=P000001
      GG(0002,1,1)=P000002
      GG(0003,1,1)=P000003
      GG(0004,1,1)=P000004
      GG(0005,1,1)=P000005
      GG(0006,1,1)=P000006
      GG(0007,1,1)=P000007
      GG(0008,1,1)=P000008
      GG(0009,1,1)=P000009
      GG(0010,1,1)=P000010
      GG(0011,1,1)=P000011
      ENDIF
      IF (MSEC == 1) THEN
!                      A000045 = DERIVATIVE OF A000043 W.R.T. ETA(003)
      A000045=B000006
      IF(A_0FLG == 1)THEN
!                      A000072 = DERIVATIVE OF A000066 W.R.T. ETA(007)
      A000072=delta*B000014
!                      A000073 = DERIVATIVE OF A000067 W.R.T. ETA(007)
      A000073=B000012*B000014
!                      A000074 = DERIVATIVE OF A000067 W.R.T. ETA(006)
      A000074=c*B000012
!                      A000079 = DERIVATIVE OF A000070 W.R.T. ETA(006)
      A000079=B000019*A000074
!                      A000080 = DERIVATIVE OF A000070 W.R.T. ETA(007)
      A000080=B000019*A000073
!                      A000081 = DERIVATIVE OF A000071 W.R.T. ETA(007)
      A000081=B000019*A000072
      D2A_0(003,006,006)=A000079
      D2A_0(003,007,006)=A000080
      D2A_0(003,007,007)=A000081
      ENDIF
      IF(A_0FLG == 1)THEN
!                      A000090 = DERIVATIVE OF B000022 W.R.T. ETA(006)
      A000090=-d*B000012
!                      A000091 = DERIVATIVE OF A000082 W.R.T. ETA(007)
      A000091=B000022*B000014
!                      A000092 = DERIVATIVE OF B000023 W.R.T. ETA(007)
      A000092=-d*B000014
!                      A000093 = DERIVATIVE OF A000083 W.R.T. ETA(007)
      A000093=B000012*A000092
!                      A000094 = DERIVATIVE OF A000083 W.R.T. ETA(006)
      A000094=B000023*B000012
!                      A000099 = DERIVATIVE OF A000086 W.R.T. ETA(006)
      A000099=B000027*B000012
      B000032=-1.D0/B000021/B000021
!                      A000100 = DERIVATIVE OF B000030 W.R.T. ETA(006)
      A000100=B000032*A000086
!                      A000101 = DERIVATIVE OF A000087 W.R.T. ETA(006)
      A000101=A000083*A000100
!                      A000102 = DERIVATIVE OF A000087 W.R.T. ETA(006)
      A000102=B000030*A000094+A000101
!                      A000103 = DERIVATIVE OF A000087 W.R.T. ETA(007)
      A000103=B000030*A000093
!                      A000104 = DERIVATIVE OF A000088 W.R.T. ETA(007)
      A000104=B000030*A000091
      B000033=-1.D0/B000021/B000021
!                      A000105 = DERIVATIVE OF B000031 W.R.T. ETA(006)
      A000105=B000033*A000083
!                      A000106 = DERIVATIVE OF B000031 W.R.T. ETA(007)
      A000106=B000033*A000082
      B000034=B000020/B000021/B000021/B000021
!                      A000107 = DERIVATIVE OF B000031 W.R.T. ETA(006)
      A000107=B000034*A000086+A000105
      B000035=B000020/B000021/B000021/B000021
!                      A000108 = DERIVATIVE OF B000031 W.R.T. ETA(006)
      A000108=B000035*A000086+A000107
!                      A000109 = DERIVATIVE OF A000089 W.R.T. ETA(006)
      A000109=A000086*A000108
!                      A000110 = DERIVATIVE OF A000089 W.R.T. ETA(007)
      A000110=A000086*A000106
!                      A000111 = DERIVATIVE OF A000089 W.R.T. ETA(006)
      A000111=B000031*A000099+A000109
!                      A000112 = DERIVATIVE OF A000089 W.R.T. ETA(007)
      A000112=A000103+A000110
!                      A000113 = DERIVATIVE OF A000089 W.R.T. ETA(006)
      A000113=A000102+A000111
      D2A_0(004,006,006)=A000113
      D2A_0(004,007,006)=A000112
      D2A_0(004,007,007)=A000104
      ENDIF
      IF(A_0FLG == 1)THEN
!                      A000126 = DERIVATIVE OF B000038 W.R.T. ETA(006)
      A000126=-d*B000012
!                      A000127 = DERIVATIVE OF A000114 W.R.T. ETA(007)
      A000127=B000038*B000014
!                      A000128 = DERIVATIVE OF B000039 W.R.T. ETA(007)
      A000128=-d*B000014
!                      A000129 = DERIVATIVE OF A000115 W.R.T. ETA(007)
      A000129=B000012*A000128
!                      A000130 = DERIVATIVE OF A000115 W.R.T. ETA(006)
      A000130=B000039*B000012
!                      A000135 = DERIVATIVE OF B000043 W.R.T. ETA(006)
      A000135=e*B000012
!                      A000136 = DERIVATIVE OF A000118 W.R.T. ETA(007)
      A000136=B000043*B000014
!                      A000137 = DERIVATIVE OF B000044 W.R.T. ETA(007)
      A000137=e*B000014
!                      A000138 = DERIVATIVE OF A000119 W.R.T. ETA(007)
      A000138=B000012*A000137
!                      A000139 = DERIVATIVE OF A000119 W.R.T. ETA(006)
      A000139=B000044*B000012
      B000050=-1.D0/B000037/B000037
!                      A000144 = DERIVATIVE OF B000048 W.R.T. ETA(006)
      A000144=B000050*A000119
!                      A000145 = DERIVATIVE OF B000048 W.R.T. ETA(007)
      A000145=B000050*A000118
!                      A000146 = DERIVATIVE OF A000122 W.R.T. ETA(007)
      A000146=A000115*A000145
!                      A000147 = DERIVATIVE OF A000122 W.R.T. ETA(006)
      A000147=A000115*A000144
!                      A000148 = DERIVATIVE OF A000122 W.R.T. ETA(006)
      A000148=B000048*A000130+A000147
!                      A000149 = DERIVATIVE OF A000122 W.R.T. ETA(007)
      A000149=B000048*A000129+A000146
!                      A000150 = DERIVATIVE OF A000123 W.R.T. ETA(007)
      A000150=A000114*A000145
!                      A000151 = DERIVATIVE OF A000123 W.R.T. ETA(007)
      A000151=B000048*A000127+A000150
      B000051=-1.D0/B000037/B000037
!                      A000152 = DERIVATIVE OF B000049 W.R.T. ETA(006)
      A000152=B000051*A000115
!                      A000153 = DERIVATIVE OF B000049 W.R.T. ETA(007)
      A000153=B000051*A000114
      B000052=B000036/B000037/B000037/B000037
!                      A000154 = DERIVATIVE OF B000049 W.R.T. ETA(006)
      A000154=B000052*A000119+A000152
!                      A000155 = DERIVATIVE OF B000049 W.R.T. ETA(007)
      A000155=B000052*A000118+A000153
      B000053=B000036/B000037/B000037/B000037
!                      A000156 = DERIVATIVE OF B000049 W.R.T. ETA(006)
      A000156=B000053*A000119+A000154
!                      A000157 = DERIVATIVE OF B000049 W.R.T. ETA(007)
      A000157=B000053*A000118+A000155
!                      A000158 = DERIVATIVE OF A000124 W.R.T. ETA(007)
      A000158=A000119*A000157
!                      A000159 = DERIVATIVE OF A000124 W.R.T. ETA(006)
      A000159=A000119*A000156
!                      A000160 = DERIVATIVE OF A000124 W.R.T. ETA(006)
      A000160=B000049*A000139+A000159
!                      A000161 = DERIVATIVE OF A000124 W.R.T. ETA(007)
      A000161=B000049*A000138+A000158
!                      A000162 = DERIVATIVE OF A000124 W.R.T. ETA(007)
      A000162=A000149+A000161
!                      A000163 = DERIVATIVE OF A000124 W.R.T. ETA(006)
      A000163=A000148+A000160
!                      A000164 = DERIVATIVE OF A000125 W.R.T. ETA(007)
      A000164=A000118*A000157
!                      A000165 = DERIVATIVE OF A000125 W.R.T. ETA(007)
      A000165=B000049*A000136+A000164
!                      A000166 = DERIVATIVE OF A000125 W.R.T. ETA(007)
      A000166=A000151+A000165
      D2A_0(005,006,006)=A000163
      D2A_0(005,007,006)=A000162
      D2A_0(005,007,007)=A000166
      ENDIF
!                      A000168 = DERIVATIVE OF A000167 W.R.T. ETA(001)
      A000168=B000002
!                      A000170 = DERIVATIVE OF A000169 W.R.T. ETA(002)
      A000170=B000004
!                      A000172 = DERIVATIVE OF A000171 W.R.T. ETA(006)
      A000172=B000012
!                      A000174 = DERIVATIVE OF A000173 W.R.T. ETA(004)
      A000174=B000008
!                      A000176 = DERIVATIVE OF A000175 W.R.T. ETA(005)
      A000176=B000010
!                      A000178 = DERIVATIVE OF A000177 W.R.T. ETA(003)
      A000178=A000045
!                      A000180 = DERIVATIVE OF A000179 W.R.T. ETA(007)
      A000180=B000014
      GG(0001,0002,0002)=A000168
      GG(0002,0003,0003)=A000170
      GG(0006,0007,0007)=A000172
      GG(0007,0005,0005)=A000174
      GG(0008,0006,0006)=A000176
      GG(0009,0004,0004)=A000178
      GG(0011,0008,0008)=A000180
      ENDIF
      RETURN
      END
