On 04/04/2016 05:18 PM, Keller, Jacob wrote:
> Dear Crystallographers,
>
> I have CADed a number of files together in very parallel ways, but one of 13 gives me the error:
>
> CAD: *** Datum less than VRSET_MAGIC, it MUST be changed ***
>
> Does anyone know what this means?
>
> JPK
>
There may be a clue in here-
grep VRSET_MAGIC $CCP4/src/*.f
/raid/bin/ccp4-6.1.3/src/cad.f: + VRSET_MAGIC
/raid/bin/ccp4-6.1.3/src/cad.f: VRSET_MAGIC = -9.0E+12
/raid/bin/ccp4-6.1.3/src/cad.f: IF (RNGMTZ(1,I) .LE. VRSET_MAGIC) CALL CCPERR(1,
/raid/bin/ccp4-6.1.3/src/cad.f: + ' *** Datum less than VRSET_MAGIC, it MUST be changed ***')
/raid/bin/ccp4-6.1.3/src/cad.f: RECOUT(JDO170) = VRSET_MAGIC
/raid/bin/ccp4-6.1.3/src/cad.f: RECOUT(JDO370) = VRSET_MAGIC
/raid/bin/ccp4-6.1.3/src/cad.f: IF (RECWORK(JDO390,JDO380) .GT. VRSET_MAGIC)
/raid/bin/ccp4-6.1.3/src/cad.f: IF( RECOUT(JDO420).GT.VRSET_MAGIC) RECSUM = RECSUM + 1
/raid/bin/ccp4-6.1.3/src/cad.f:C---- Reset MNF from VRSET_MAGIC to VAL_MISSING to write out
/raid/bin/ccp4-6.1.3/src/cad.f: CALL RESET_MAGIC(0,RECOUT,RECIN,NUM_COLS_OUT,VRSET_MAGIC,
/raid/bin/ccp4-6.1.3/src/cad.f: CALL RESET_MAGIC(0,RECOUT,BDATA,NUM_COLS_OUT,VRSET_MAGIC,VALM)
/
It looks like you have a large negative datum that needs to be set to the missin-number flag value in order to proceed.
>
>
> *******************************************
> Jacob Pearson Keller, PhD
> Looger Lab/HHMI Janelia Research Campus
> 19700 Helix Dr, Ashburn, VA 20147
> email: [log in to unmask]
> *******************************************
>
|