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 <[log in to unmask]> 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);
>
|