integer maxtwin_l parameter (maxtwin_l=24) integer nmaxpart_l parameter (nmaxpart_l = 3) real scale_ov0,b_ov0,scale_bulk0,b_bulk0,baniso_ov0(6) real scale_part0(nmaxpart_l),b_part0(nmaxpart_l) real twin_frac0(maxtwin_l) common /save_ls_p/scale_ov0,b_ov0,scale_bulk0,b_bulk0, & baniso_ov0,scale_part0,b_part0, & twin_frac0