$! awa 970408 ver 1.1
$! converted from unix example
$!
$! mapmask - map/mask extend program 
$!
$!
$if f$search(" CCP4_SCR:toxd_fc.map") .eqs. "" then goto term
$!
$ass/user CCP4_SCR:toxd_fc.map mapin
$ass/user CCP4_SCR:toxd_ins1.ext mapout
$ass/user ccp4_master:[ccp4.examples.toxd]toxd.pdb xyzin
$mapmask   
BORDER 5.0
END
$!
$! Or:
$!
$ass/user CCP4_SCR:toxd_fc.map mapin
$ass/user CCP4_SCR:toxd_ins2.ext mapout
$mapmask    
XYZLIM 0 79 0 79 0 19
END
$!
$! To Resection A Map:
$!
$ass/user CCP4_SCR:toxd_fc.map mapin
$ass/user CCP4_SCR:toxd_inszyx.map mapout
$mapmask  
AXIS Z Y X
END
$!
$! To Rescale A Map: (r'=2r+1)
$!
$ass/user CCP4_SCR:toxd_fc.map mapin
$ass/user CCP4_SCR:toxd_insscaled.map mapout
mapmask   
SCALE FACTOR 2.0 1.0
END
$!
$! To Make A Mask From A Map (Mask Covers Region Of Map > 0.5):
$!
$ass/user CCP4_SCR:toxd_fc.map mapin
$ass/user CCP4_SCR:toxd_ins1.msk mskout
$mapmask                  
END
$!
$! To Make A Mask From A Map To Cover 70% Of The Unit Cell:
$!
$ass/user CCP4_SCR:toxd_fc.map mapin
$ass/user CCP4_SCR:toxd_ins2.msk mskout
$mapmask                   
MASK FRAC 0.7
END
$!
$! To Make A Map From A Mask:
$!
$ass/user CCP4_SCR:toxd_ins1.msk  mapin
$ass/user CCP4_SCR:toxd_from_msk.map mapout
$mapmask   
END
$! 
$! To print out a mask section z=5 of all points above 0.9:
$! 
$ass/user CCP4_SCR:toxd_fc.map mapin
$mapmask   
MASK CUT 0.9
PRINT MASK Z 5 5
END
$!
$! To print out alternate map x-sections, scaled x100:
$!
$ass/user CCP4_SCR:toxd_fc.map mapin
$mapmask  
SCALE FACTOR 100
PRINT MAP X 0 24 2
END
$!
$!To Make A Solvent Envelope, 2.5 Angstroms About Atoms.
$! 
$! Use sfall to construct density for all points up to but
$! not exceeding 2.5 A from atoms. Then use mapmask to select
$! this region.
$! 
$! Step 1: sfall
$!
$ass/user ccp4_master:[ccp4.examples.toxd]toxd.pdb xyzin
$ass/user CCP4_SCR:junk.map mapout
$sfall    
MODE ATMMAP		
VDWR 2.5	
SYMM P212121		 
END
$! 
$! Step 2: mapmask
$!
$ass/user CCP4_SCR:junk.map mapin
$ass/user CCP4_SCR:bar.msk mapout
$mapmask                
MASK CUT 0.01		
END
$!
$exit
$term:
$write sys$output "! run sf_calc.com first"