Produces a list of derivative matrix d(omaga x zeta)/dzeta, where omega is the rotation speed of the A FoR, ASSUMING constant panel norm.
Each list is such that: - the ii-th element is associated to the ii-th bound surface collocation point, and will contain a sub-list such that:
the j-th element of the sub-list is the dAIC_dzeta matrices w.r.t. the
zeta d.o.f. of the j-th bound surface.
Hence, DAIC*[ii][jj] will have size K_ii x Kzeta_jj
call: ncDOmegaZetavert = nc_domegazetadzeta(Surfs,Surfs_star)