get_pcyl
cur_dtype->dtype_pcyl = get_pcyl(ncyl, cur_dtype->dtype_acyl);
d->dtype_pcyl = get_pcyl(d->dtype_ncyl, d->dtype_acyl);
int get_pcyl(int n_cyls, int p_cyls);