# Color a map by local resolution computed by ResMap
#   Initial version of script courtesy of Tom Goddard, UCSF.

# Open both volumes and hide the ResMap volume.
set bg_color white
open #0 run2_half1_class001_unfil.mrc
open #1 run2_half1_class001_unfil_resmap.mrc
volume #1 hide

# Color the original map with the values from the ResMap output.
scolor #0 volume #1 cmap 11.00,blue:12.80,cyan:14.60,green:16.40,yellow:18.20,orange:20.00,red:20.01,gray



# OPTIONAL: ResMap Slice Animation.

# Show midway slice of the original map with contour level below the minimum map value.
volume #0 planes z,100 step 1 level -1 style surface

# Show a smooth transparent contour surface indicating the structure boundaries.
vop gaussian #0 sDev 5 model #2
volume #2 level 0.02 step 1 color .9,.7,.7,.5

# Zoom out a bit and tilt to a nice viewing angle.
turn x -45
turn y -30
turn z -30
scale 0.5

# Cycle through planes from N/2 to 4N/5 up to N/5 and back to N/2.
volume #0 planes z,100,160,0.25
wait 240
volume #0 planes z,160,40,0.25
wait 480
volume #0 planes z,40,100,0.25