
 type (iosystem_desc_t)     :: iosystem
 integer (i4)               :: dims(1)
 integer (i4)               :: compdof
 type (io_desc_t)           :: iodesc
 ...
 !---------------------------------------
 ! Initializing the decomposition on PE 0
 !---------------------------------------
 dims(1) = 8
 compdof = (/2,4,5/)
 call PIO_initdecomp(iosystem,PIO_double,dims,compdof,iodesc)
