program eplf_hf PROVIDE ao_prim_num_max call write_grid_eplf() IRP_IF MPI integer :: ierr call MPI_FINALIZE(ierr) IRP_ENDIF end