Hi Alex Wouldn't you be better off using the API (http://www.ccp4.ac.uk/dist/html/C_library/cmaplib_8h.html): that way your program will have a good chance of still working in another 10 years' time!
Cheers -- Ian On 2 July 2012 14:41, Soares, Alexei <[email protected]> wrote: > Hello, Alex here. > > Ten years ago, I read CCP4 maps this way, but it no longer works ... can > anyone tell me why? Note that map1 is an integer so that multiple headers > can be read. Also, supposedly the map header is text, but a simple "head > -1" dumps lots of digital garbage from the map header. > > Thanks, and sorry for the beginner question! > > // ******************** READ MAP HEADER **************************** > > fread( &MAP_H[map1].NC , sizeof(long) , 1, read1); > fread( &MAP_H[map1].NR , sizeof(long) , 1, read1); > fread( &MAP_H[map1].NS , sizeof(long) , 1, read1); > > fread( &MAP_H[map1].MODE , sizeof(long) , 1, read1); > > fread( &MAP_H[map1].NCSTART, sizeof(long) , 1, read1); > fread( &MAP_H[map1].NRSTART, sizeof(long) , 1, read1); > fread( &MAP_H[map1].NSSTART, sizeof(long) , 1, read1); > > fread( &MAP_H[map1].NX , sizeof(long) , 1, read1); > fread( &MAP_H[map1].NY , sizeof(long) , 1, read1); > fread( &MAP_H[map1].NZ , sizeof(long) , 1, read1); > > fread( MAP_H[map1].CELL , sizeof(float), 6, read1); > > fread( &MAP_H[map1].MAPC , sizeof(long) , 1, read1); > fread( &MAP_H[map1].MAPR , sizeof(long) , 1, read1); > fread( &MAP_H[map1].MAPS , sizeof(long) , 1, read1); > > fread( &MAP_H[map1].AMIN , sizeof(long) , 1, read1); > fread( &MAP_H[map1].AMAX , sizeof(long) , 1, read1); > fread( &MAP_H[map1].AMEAN , sizeof(long) , 1, read1); > > fread( &MAP_H[map1].ISPG , sizeof(long) , 1, read1); > > fread( &MAP_H[map1].NSY , sizeof(long) , 1, read1); > > fread( MAP_H[map1].REST , sizeof(float), 32, read1); > > fread( MAP_H[map1].LAB , sizeof(char) , LAB_LEN, read1); > fread( MAP_H[map1].SYM , sizeof(char) , MAP_H[map1].NSY, read1); >
