+#ifdef DEBUG
+ write (iout,*) "dudconst, duscdiff, dugamma,dutheta"
+ write (iout,*) "wumb",wumb
+ do i=1,nct
+ write (iout,'(i5,3f10.5,5x,3f10.5,5x,2f10.5)')
+ & i,(dudconst(j,i),j=1,3),(duscdiff(j,i),j=1,3),
+ & dugamma(i),dutheta(i)
+ enddo
+#endif