Hello,
Once you get the output file from cpptraj, you could first use awk to
process your output file,
converting number of frames to ns. Then you can plot the processed file
with gnuplot.
awk '{print $1/100.0 , $2}' output > newoutput, 100.0 means you wrote
out every 10 ps.
All the best,
Qinghua
On 01/10/2018 02:44 PM, Rana Rehan Khalid wrote:
> Dear Amber Users
>
> Through CCPTRAJ DSSP analysis we get .gnu file. The plot x-axis range
> [0.000:7002.000], this data separated with xtic values as
> 0,1000,2000,.......7000 (frame numbers) on x axis. While i want to
> change/multiply/transform this frame number into nano second, so that my
> x-axis values showed as 0,1,2......14(nano second) after multiplying with
> .002. here is some lines of my .gnu files
>
>
> enter code here
>
> set pm3d map corners2color c1
> set ytics 1.000, 1.000
> set ytics("HM1_187.O2D-SER_137.OG-HG" 1.000,"HM1_187.O1D-SER_137.OG-HG"
> 2.000)
> set xlabel "Frame"`enter code here`
> set ylabel ""
> set yrange [ 0.000: 4.000]
> set xrange [ 0.000:7002.000]
> splot "-" with pm3d title "hhhh.gnu"
> 1.000 1.000 1
> 1.000 2.000 0
> 1.000 3.000 0``
>
> 2.000 1.000 1
> 2.000 2.000 1
> 2.000 3.000 0
>
> 3.000 1.000 0
> 3.000 2.000 1
> 3.000 3.000 0
>
> 4.000 1.000 0
> 4.000 2.000 0
> 4.000 3.000 0
>
> 5.000 1.000 1
> 5.000 2.000 0
> 5.000 3.000 0
>
> 6.000 1.000 1
> 6.000 2.000 0
> 6.000 3.000 0
> and so on.
> kindly guide how I can change my x axis tics to nano second rather then
> frames numbers. Thanks.
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Jan 10 2018 - 06:00:13 PST