This command (make-and-draw-map "../refined01.mtz" "DELFWT" "PHDELFWT" "" 0 1)
fails to make the difference map. Coot appears to report that the mtz-file (produced by refmac) does not have the DELFWT/PHDELFWT columns, which it does. The 2fo-fc map works just fine. Moreover, I can import the difference map from the gui and the columns in question are found. Coot reports the following valid_labels(../refined01.mtz,DELFWT,PHDELFWT,,0) returns 0 WARNING:: label(s) not found in mtz file ../refined01.mtz DELFWT PHDELFWT I noticed that gui uses make-and-draw-map-with-reso-with-refmac-params and calls the columns "/unknown/unknown030610/DELFWT" and "/unknown/unknown030610/PHDELWT", which must be the CRYSTAL, DATASET records from the mtz file. I want to load several maps from a script knowing nothing but file names and column names, so what am I doing wrong? Ed. PS. I just found the solution for refmac output, which is auto-read-make-and-draw-maps. Works fine and I don't need to specify column names. Still, shouldn't make-and-draw-map work too? -- Coot verendus est
