|
39 | 39 | " dataLabel = \"ion current\",\n", |
40 | 40 | " mode = 'append'\n", |
41 | 41 | ")\n", |
42 | | - "# Have fixed axis for better visualization\n", |
43 | | - "#plot[\"ax\"].set_xlim(0, 100)\n", |
| 42 | + "# fixed y-axis for visualization\n", |
44 | 43 | "plot[\"ax\"].set_ylim(-25,25)\n", |
45 | 44 | "\n", |
46 | 45 | "u = mda.Universe(NAMD_TOPOL, \"imd://10.206.72.141:8888\", buffersize = 1024*1024*1024)\n", |
|
61 | 60 | " # correct for PBC jumps\n", |
62 | 61 | " deltas[deltas>zdist/2.] -= zdist\n", |
63 | 62 | " deltas[deltas<-zdist/2.] += zdist\n", |
64 | | - " cd = deltas*charges\n", |
65 | | - " current = np.sum(cd)/zdist/dt*unitConv\n", |
| 63 | + " Qd = deltas*charges\n", |
| 64 | + " current = np.sum(Qd)/zdist/dt*unitConv\n", |
66 | 65 | " ionRef = ionPos\n", |
67 | 66 | " plot['update'](time, current)\n", |
68 | 67 | " i = i + 1" |
|
84 | 83 | " dataLabels = [\"cation current\", \"anion current\"],\n", |
85 | 84 | " mode = 'append'\n", |
86 | 85 | ")\n", |
87 | | - "# Have fixed axis for better visualization\n", |
88 | | - "#plot[\"ax\"].set_xlim(0, 100)\n", |
| 86 | + "# fixed y-axis for visualization\n", |
89 | 87 | "plot[\"ax\"].set_ylim(-25,25)\n", |
90 | 88 | "\n", |
91 | 89 | "u = mda.Universe(NAMD_TOPOL, \"imd://10.206.72.141:8888\", buffersize = 1024*1024*1024)\n", |
|
0 commit comments