
            imax=nsoa*nsob*nsoc
            jmax=nsod
            kmax=ncod
            i = 0
