c1: CENTER ATOMS=1-200 

DENSITY SPECIES=2-201 LABEL=d1

d2: INCYLINDER ATOM=c1 DATA=d1 DIRECTION=Z RADIUS={TANH R_0=1.5} SIGMA=0.1 LOWER=-0.1 UPPER=0.1

# DENSITY SPECIES=2-201 LABEL=d1n
# d2n: INCYLINDER ATOM=c1 DATA=d1n DIRECTION=Z RADIUS={TANH R_0=1.5} SIGMA=0.1 LOWER=-0.1 UPPER=0.1 NUMERICAL_DERIVATIVES

d2s: INSPHERE ATOM=c1 DATA=d1 RADIUS={TANH R_0=1.5}

PRINT ARG=d2,d2s FILE=colvar FMT=%8.4f 
DUMPDERIVATIVES ARG=d2,d2s FILE=derivatives FMT=%8.4f


