5258       subroutine ggrad2omgrad(ggrad,omgrad,omop,urop,uiop,  argument
5281      &     urop(ndim,ndim,2,2),  uiop(ndim,ndim,2,2)
5418                             uiui_t = uiop(ip,it,imp,imt)*
5419      &                               uiop(iq,it,imq,imt)
5420                             uiur_t = uiop(ip,it,imp,imt)*
5423      &                               uiop(iq,it,imq,imt)
5430                               uiui_u = uiop(ir,iu,imr,imu)*
5431      &                                 uiop(is,iu,ims,imu)
5432                               uiur_u = uiop(ir,iu,imr,imu)*
5435      &                                 uiop(is,iu,ims,imu)
5485       subroutine ggrad2ugrad(ggrad,urgrad,omop,urop,uiop,  argument
5508      &     urop(ndim,ndim,2,2),  uiop(ndim,ndim,2,2)
5655                             uiui_rs = uiop(ir,iu,imr,imu)*
5656      &                                  uiop(is,iu,ims,imu)
5657                             uiur_rs = uiop(ir,iu,imr,imu)*
5660      &                                  uiop(is,iu,ims,imu)
5664                             uiui_pq = uiop(ip,iu,imp,imu)*
5665      &                                  uiop(iq,iu,imq,imu)
5666                             uiur_pq = uiop(ip,iu,imp,imu)*
5669      &                                  uiop(iq,iu,imq,imu)
5701                               ui = uiop(ir,iw,imr,imw)
5707                               ui = - uiop(is,iw,ims,imw)
5714                               ui = uiop(ip,iw,imp,imw)
5720                               ui = - uiop(iq,iw,imq,imw)
5766       subroutine ggrad2ugrad_old(ggrad,urgrad,omop,urop,uiop,  argument
5789      &     urop(ndim,ndim,2,2),  uiop(ndim,ndim,2,2)
5937                               uiui_rs = uiop(ir,iu,imr,imu)*
5938      &                                  uiop(is,iu,ims,imu)
5939                               uiur_rs = uiop(ir,iu,imr,imu)*
5942      &                                  uiop(is,iu,ims,imu)
5946                               uiui_pq = uiop(ip,iu,imp,imu)*
5947      &                                  uiop(iq,iu,imq,imu)
5948                               uiur_pq = uiop(ip,iu,imp,imu)*
5951      &                                  uiop(iq,iu,imq,imu)
5983                                     ui = uiop(ir,iw,imr,imw)
5986                                     ui = ui - uiop(is,iw,ims,imw)
5995                                     ui = uiop(ip,iw,imp,imw)
5998                                     ui = ui - uiop(iq,iw,imq,imw)
6047       subroutine uou2g(omop,urop,uiop,gop,  argument
6073      &     uiop(ndim,ndim,2,2)
6118             call wrtmat2(uiop(1,1,imp,imq),ndim,ndim,ndim,ndim)
6243                             uiui_t = uiop(ip,it,imp,imt)*
6244      &                               uiop(iq,it,imq,imt)
6245                             uiur_t = uiop(ip,it,imp,imt)*
6248      &                               uiop(iq,it,imq,imt)
6254                               uiui_u = uiop(ir,iu,imr,imu)*
6255      &                                 uiop(is,iu,ims,imu)
6256                               uiur_u = uiop(ir,iu,imr,imu)*
6259      &                                 uiop(is,iu,ims,imu)