Subroutine get_pesname USE FileUnits_Module Use CommonInfo_Module Use PES_Module Implicit None !==================================================================== ! This routine reads the type of potential energy surface that will ! be used !===================================================================== ! I N P U T !==================================================================== ! I N T E R N A L S INTEGER :: istat !===================================================================== ! N A M E L I S T S NAMELIST / pes / pesname OPEN(UNIT=Input_Unit,FILE=Trim(inputfile1),IOSTAT=istat,STATUS='old',ACTION='read') IF (istat.ne.0) THEN !WRITE(Msg_Unit,*)"InputFile1=",inputfile1 STOP 'get_pesname - OPEN failed' ENDIF READ(Input_Unit,pes) CLOSE(UNIT=Input_Unit,IOSTAT=istat,STATUS='keep') IF (istat.ne.0) STOP 'get_pesname - CLOSE failed' !===================================================================== ! Wirte pesname to file and screen Write(Out_Unit,21) 'Potential Energy Surface Info:','pesname', TRIM(pesname) Write(*,22) 'pesname', TRIM(pesname) 21 FORMAT(/1X,a,/1X,50('-'),/1X,a,10('.'),a) 22 FORMAT(1X,a,8('.'),a) Return End Subroutine get_pesname