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

