User Tools

Site Tools


comp-sa

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision
Previous revision
comp-sa [2021/06/09 08:44] admincomp-sa [2023/10/12 11:04] (current) edit
Line 5: Line 5:
 Instruction to obtain surface area per lipid of **__multicomponent system__** \\ Instruction to obtain surface area per lipid of **__multicomponent system__** \\
 Original Authors: Xiaohong Zhuang and Viviana Monje\\ Original Authors: Xiaohong Zhuang and Viviana Monje\\
-Updated Authors: Yalun Yu and Jeffery Klauda+Updated Authors: Yalun Yu and Jeffery Klauda\\ 
 +Updated rtfpsf folder: Joshua Lucker
  
-The script is in DT2 path: **/lustre/jbklauda/scripts/area_mult_lipids** \\+The script is in ZT1 path: **/afs/shell.umd.edu/project/energybio/shared/jbklauda/scripts/area_mult_lipids** \\
 or using zip file: {{ area_mult_lipids.gz |}} or using zip file: {{ area_mult_lipids.gz |}}
  
Line 112: Line 113:
 **>> area-a**:  area-a is input file for top leaflet, area-b is for bottom leaflet **>> area-a**:  area-a is input file for top leaflet, area-b is for bottom leaflet
  
-**>> 2**: The lipid types of 2 are used instead of 8. (1-atom represented choline/steroland other 3-atoms represented lipids, e.g. glycerol lipids, sphingo lipids). So if you have sterols set this two , if none set to 1.+**>> 2**: This is the dimension of the tessellationwhich is always 2.
  
 **>> 2196**: "image atoms" = 9*(#chol+3*#lip)=9*(28+3*72) **>> 2196**: "image atoms" = 9*(#chol+3*#lip)=9*(28+3*72)
Line 124: Line 125:
  
 **9. calc_avg_vertical.py & calc_avg_final.py:** \\ **9. calc_avg_vertical.py & calc_avg_final.py:** \\
-Update the number of each lipid component per leaflet in the same order as in lipid.dat but remove the sterol flag.+Update the number of each lipid component per leaflet in the same order as in lipid.dat. If you lack sterol, the first number in the Nlist MUST be zero due to the assumption in calc_avg_vertical.py.
 <code> <code>
 # Update the number of each lipid component per leaflet in the same order as in area-a.inp # Update the number of each lipid component per leaflet in the same order as in area-a.inp
 Nlist = [28,12,13,9,5,22,7,4] Nlist = [28,12,13,9,5,22,7,4]
 </code> </code>
 +
 +**10. rtfpsf.str:** \\
 +The file names in the rtfpsf.str file will most likely need to be updated.
 +
  
 **Run the scripts:** \\ **Run the scripts:** \\
comp-sa.1623242690.txt.gz · Last modified: 2021/06/09 08:44 by admin