

      der(1) = exp( thetalog(t,1) * b(2,je) + thetalog(t,2)*b(3,je)) 

      der(2) = b(1,je)*der(1) * thetalog(t,1) 

      der(3) = b(1,je)*der(1) * thetalog(t,2)

      