Skip to content

Commit

Permalink
added analysis
Browse files Browse the repository at this point in the history
Former-commit-id: 3620687
  • Loading branch information
simongravelle committed Jan 4, 2024
1 parent f1cd5fd commit 2c60fab
Show file tree
Hide file tree
Showing 11 changed files with 4 additions and 3 deletions.
7 changes: 4 additions & 3 deletions examples/illustrations/bulk-water-tip3p/analysis.ipynb
Original file line number Diff line number Diff line change
Expand Up @@ -100,16 +100,17 @@
"outputs": [],
"source": [
"while 1<2: \n",
" # high res\n",
" for T in [280, 290, 300, 310, 320]:\n",
" datapath = git_path + \"/data/bulk-water-tip3p/bulk-water/raw-data/N4000_\" + str(T) + \"K/\"\n",
" xtc = [datapath+\"prod1.xtc\", datapath+\"prod2.xtc\"]\n",
" u = mda.Universe(datapath+\"prod1.tpr\", xtc)\n",
" dt = np.round(u.trajectory.dt,2)\n",
" hydrogen = u.select_atoms(\"name HW1 HW2\")\n",
" n_hydrogen = hydrogen.n_atoms\n",
" intra = nmrformd.NMR(u, atom_group = hydrogen, neighbor_group = hydrogen, type_analysis = \"intra_molecular\", number_i=5)\n",
" inter = nmrformd.NMR(u, atom_group = hydrogen, neighbor_group = hydrogen, type_analysis = \"inter_molecular\", number_i=1)\n",
" intra = nmrformd.NMR(u, atom_group = hydrogen, neighbor_group = hydrogen,\n",
" type_analysis = \"intra_molecular\", number_i=5)\n",
" inter = nmrformd.NMR(u, atom_group = hydrogen, neighbor_group = hydrogen,\n",
" type_analysis = \"inter_molecular\", number_i=1)\n",
" n_intra = save_result(intra, name = \"N4000_intra_T\" + str(T) + \"K\")\n",
" n_inter = save_result(inter, name = \"N4000_inter_T\" + str(T) + \"K\")\n",
" print(n_hydrogen, T, n_intra)"
Expand Down
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.

0 comments on commit 2c60fab

Please sign in to comment.