Skip to content

Commit

Permalink
added analysis
Browse files Browse the repository at this point in the history
  • Loading branch information
simongravelle committed Jan 4, 2024
1 parent cdee745 commit 3620687
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 3620687

Please sign in to comment.