fork download 
  1.                         open(unit=u_live,file=livename,status='old')
  2.                         i=0
  3.                         do
  4.                                 i=i+1
  5.                                 read(u_live,*,IOSTAT=iostatus) p(1:ndims,i),l(i)
  6.                                 if(iostatus<0) then
  7.                                         i=i-1
  8.                                         if(i<nlive) then
  9.                                                 write(*,*)"ERROR: live points file has less than ",nlive," points."
  10.                                                 write(*,*)"Aborting"
  11. #ifdef MPI
  12.                                                 call MPI_ABORT(MPI_COMM_WORLD,errcode)
  13. #endif
  14.                                                 stop
  15.                                         endif
  16.                                         exit
  17.                                 endif
  18.                                 if(i>nlive) then
  19.                                         write(*,*)"ERROR: live points file has greater than ",nlive," points."
  20.                                         write(*,*)"Aborting"
  21. #ifdef MPI
  22.                                         call MPI_ABORT(MPI_COMM_WORLD,errcode)
  23. #endif
  24.                                         stop
  25.                                 endif
  26.                         enddo
  27.                         close(u_live)
  28.  
  29. ! ERROR: live points file has less than        20000  points. ndims=           8 
  30. !  i=       10000
  31. ! Aborting
  32. !
  33. ! "live points file" is LargeLogMasseslive.points and i give the output of wc for it below
  34. !
  35. !  [akulasuj@qb1 HiggsNS]$ wc LargeLogMasseslive.points 
  36. !    20000  160000 4500000 LargeLogMasseslive.points
  37. !
  38.  
Not running #stdin #stdout 0s 0KB
comments (?)
stdin 
Standard input is empty
stdout 
Standard output is empty
https://ideone.com/PLPIm
language:
Fortran (gfortran 8.3)
created:
14 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!