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);

Reply via email to