
 type (iosystem_desc_t)     :: iosystem
 integer (i4)               :: dims(1)
 integer (kind=PIO_OFFSET)  :: start(1), count(1)
 type (io_desc_t)           :: iodesc
 ...
 !---------------------------------------
 ! Initializing the decomposition on PE 0
 !---------------------------------------
 dims(1) = 8
 start(1) = 3
 count(1) = 3
 call PIO_initdecomp(iosystem,PIO_double,dims,start,count,iodesc)
