This function generate an incomplete split-plot design for given number of main plot treatments, number of subplot treatments, number of blocks and number of main plot treatments in each block
ispd.imcs(v1,b,k1,v2)
number of main plot treatments
number of blocks
number of main plot treatments in each block
number of subplot treatments
A list containing parameters, design layout and column layout of design