io_nemo_f(filename, slen, MAXBODY, param, ....) character * filename(*);/* input/output file name */ integer slen; /* #characters in ’infile’ */ integer MAXBODY; /* maximun bodies allowed */ character * param(*); /* parameters list */ Descriptionio_nemo_f is a high level C function which allows the user to perform I/O operations on NEMO snapshot files from a FORTRAN program. The use of this function is very simple and doesn’t require any knowledge of how to program in NEMO. Parameters filename=character variable Input/Output file name of the NEMO snapshot that you want to read/save. (this CHARACTER string must be a VARIABLE, not a constant !). slen=integer Lenght of the file name. (We need this information because of the interface bettween FORTRAN and C language) MAXBODY=integer Size of array that you use to get/put the data. Match the size of the array to the maximum number of bodies that you are able to read/save. param=character Param is a string in which you specify what you want to do with the NEMO file. Each choice is defined with a flag separated with ’,’. There are two kind of flags. (1) Information flags allow to specify some actions during the I/O. (2) Variables flags allow to specify what you want to get/put from/into the NEMO file. The FLAGS list is described below. DO NOT FORGET TO ADD A \0 at the end of Param variable as you can see in the example below. (1) Information flagss | save Specify that you want to save data to a NEMO file. r | read Specify that you want to read data from a NEMO file. real4|single Specify that the variables that you use to get/put data from/into NEMO files have been declared in single precision. (real * 4). real8|double Specify that the variables you use to get/put data from/into NEMO files have been declared in double precision. (real * 8). f | n3 Assumes that you declared yours two-dimensionnal FORTRAN arrays (pos, vel, acc) using the number of bodies for the first dimension and the number of space coordinates (always 3D) for the second dimension. (example : real * 4 pos(10000,3) ). c | 3n Assumes that you declared yours two-dimensionnal FORTRAN arrays (pos, vel, acc) using the number of space coordinates (always 3D) for the first dimension and the number of bodies for the second dimension. (example : real * 4 vel(3,10000) ). info | diag Gives some informations during the runtime execution. (2) Variables flags n|nbody Match to the number of bodies. t|time Match to the snapshot time. m|mass Match to the particle masses. x|pos Match to the particle positions. v|vel Match to the particle velocities. p|pot Match to the particle potentials. a|acc Match to the particle accelerations. aux Match to the auxiliary array. k|keys Match to the keys array. d|dens Match to the particle density. e|eps Match to the particle softening. st Select snapshot’s times. You have to put the selected time in a string limited by the character ’#’. Example "100:150#" select the time beween 100 and 150. The string is fully compatible with the expression usually given to nemoinpi (Cf man nemoinpi). sp Select snapshot’s particles. You have to put the selected particles in a string limited by the character ’#’. Example "0:299999#" select the particles beween 0 and 299999. The returned nbody corresponds to the selected nbody. In the example above it would thus be 300000. The string is fully compatible with the expression usually given to nemoinpi (Cf man nemoinpi). ArgumentsThe variable flags specify which data you will get/put from/into the NEMO file, hence after the selection string param, you must insert all the selected variables (arguments) in the same order with which they have been declared in param. Return ValueThe function return 1 if successfull, 0 if the end of snapshot has been reached. ExampleHere is a FORTRAN program to illustrate the use of the function io_nemo_f(). The program reads all the time steps into a NEMO a file, and saves the data in another NEMO file. Source file of the example below is located here: $NEMOSRC/nbody/io_nemo/test_src/nemo_fortran_d_3n.F C----------------------------------------------------------- program nemo_fortran implicit none c maximum # of bodies integer MAXBODY parameter (MAXBODY=100000) real*8 c particle masses. + mass(1:MAXBODY), c particle positions and velocities + pos(3,MAXBODY), vel(3,MAXBODY), c snapshot time + ts c file names character infile*80,outfile*80 integer i,j,k,nbody,io_nemo_f,close_io_nemo_f c get input snapshot filename write(*,*)’Infile name : ’ read(*,’(a80)’)infile c get output snapshot filename write(*,*)’outfile name : ’ read(*,’(a80)’)outfile i = 1 do while (i.gt.0) c read the snapshot up to the end of file c at the end of snapshot, io_nemo_f return 0 i=io_nemo_f(infile,80,MAXBODY,"real8,3n,read,n,m,x,v,t,info\0", $ nbody,mass,pos,vel,ts) c save the snapshot if (i.gt.0) then j=io_nemo_f(outfile,80,MAXBODY,"real8,3n,save,n,m,x,v,t, $ info\0",nbody,mass,pos,vel,ts) endif end do c close the snapshot ’outfile’ k= close_io_nemo_f(outfile,80) end C----------------------------------------------------------- Important Thingsa) Notice in the example above, that in the parameter list, ’n’ matches ’nbody’, ’m’ matches ’mass’, ’x’ matches ’pos’, ’v’ matches ’vel’, ’t’ matches ’ts’. All the variables are in the same order that they have been declared in the param list. Notice also that you must add \0 at the end of Param variable "real8,3n,read,n,m,x,v,t,info\0". b) You must declare all the two-dimensionnal arrays in the same way. That means all the dimensions must be the same for all the arrays, moreover both one-dimensional and two-dimensional array must have the same size for the maximum of bodies. c) All the arrays must be declared in the same floating type. d) During a "read" operation, the function io_nemo_f() return ’0’ if it is the end of the NEMO file. That means that no new values have been read.
# ---------------------------------------- # MAKEFILE to use IO_NEMO_F # # ---------------------------------------- # path for NEMO Library and IO_NEMO_F library LIBS = -L$(NEMOLIB) # - - - - - - - - - - - - - - - - - - - - # compilation with g77 compiler # - - - - - - - - - - - - - - - - - - - - G77FLAGS = -fno-second-underscore -Wno-globals nemo_fortran_g77 : nemo_fortran.F g77 $(G77FLAGS) -o $@ nemo_fortran.F $(LIBS) \ -lnemomaing77 -lnemo -lm # - - - - - - - - - - - - - - - - - - - - # compilation with gfortran compiler # - - - - - - - - - - - - - - - - - - - - GFORTFLAGS = -fno-second-underscore nemo_fortran_g77 : nemo_fortran.F gfortran $(GFORTFLAGS) -o $@ nemo_fortran.F $(LIBS) \ -lnemomaing77 -lnemo -lm # - - - - - - - - - - - - - - - - - - - - # compilation with absoft f77 compiler # - - - - - - - - - - - - - - - - - - - - ABFFLAGS = -B108 -K nemo_fortran_ABSOFT : nemo_fortran.F $(ABSOFT)/bin/f77 $(ABFFLAGS) -o $@ nemo_fortran.F $(LIBS) \ -lnemo -lm # ----------------------------------------
15-Jun-95 V1.0 : created JCL 21-Jun-95 V1.10: bugs fixes JCL 12-Dec-95 V1.11: possibility to close file JCL 11-Mar-96 V1.12: acceleration I/O added JCL 04-Apr-97 V1.13: generic real format JCL 07-Apr-97 V1.14: manual created JCL 19-Jul-02 V1.20: io_nemo/io_nemo_f unified JCL 18-Mar-04 V1.21: bugs fixed, softening added JCL 03-Mar-05 V1.30: code cleaning, valgrind safe JCL 24-Apr-06 V1.31: memory leak fixed JCL 19-Jun-06 V1.32: happy gfortran JCL 29-May-07 V1.42: handle snapshot with different JCL nbodies 29-Feb-08 V1.50: add Aux and Dens array JCL
Table of Contents