simCRSM function

simulate data according to CRSM