CHARMM Development Project Forums CHARMM Community Script Archive Sample scripts for basic tasks in CHARMM

Hi,
It appears the problem was I was running a parrallel compiled vesion of charmm, instead of the serial version.

When I run the serial version the scripts seem to work fine.

Thanks for your time and help,
Ed

hi, ED

The inputfile is your input file name, not just the "****.pdb". You misunderstand.

The templete solvent-box.str was used. Anybody has comments on the following?

size 35x35x95 is Ok, but 45x45x95 gives error:

CHARMM> read sequence tip3 @NSLAB
Parameter: NSLAB -> "1944"

CHARMM> gene wat3 noangle nodihe
NO PATCHING WILL BE DONE ON THE FIRST RESIDUE
NO PATCHING WILL BE DONE ON THE LAST RESIDUE
GENPSF> Segment 5 has been generated. Its identifier is ^@^@^@^@^@^@^@^@.


***** LEVEL -4 WARNING FROM *****
***** DATA STRUCTURE REALLOCATED IN ALLDT
******************************************
BOMLEV ( -2) IS REACHED - TERMINATING. WRNLEV IS 5

The size of the system that can be created using these scripts depends on the size parameter used when CHARMM was installed at your site. To resolve this you have to re-install CHARMM using eg LARGE or XXLARGE. Note that even if the final number of atoms/residues could be handled by your present CHARMM version, the scripts may still fail at an intermediate stage when water molecules are being added. This issue has been discussed before.

what fortran version I must use for make scrips?

FORTRAN is usually used to make "programs", not scripts. If you are modifying CHARMM source code FORTRAN77 is the standard for versions before c33b1. For c33b1 and later versions use FORTRAN90.
Tip: Use the POST function instead of REPLY when creating a post on a new topic.

Hi lennart,
Thanks for your sharing scripts. I am a beginner for charmm. Your scripts really benifit me much. But I have some questions for solvent-shell.str. I think I need your helps. In the DO_X part in solvent-shell.str, you write:

! do the same slab again, but translated -1A and rotated 30 degrees around the x-axis
read sequence tip3 @NSLAB
gene wat3 noangle nodihe
coor dupl sele segid wat0 end sele segid wat3 end
coor rota xdir 1.0 phi 30.0 sele segid wat3 end
coor transl xdir -1.0 sele segid wat3 end

I don't clearly understand why make the slab translated -1A and rotated 30 degrees around the x-axis. Why not just let the slab slide along the x-axisdo without "translated -1A and rotated 30 degrees around the x-axis" ?

Thanks.

Ardie

Just to increase the chance that we fill all possible cavities and crevices with water; there may be low density regions in the water box we started from.

Thanks. But I am still confused. Why do you rotate it by just 30 rather than by other value?

If I had picked any other value you could still have asked the same question, so perhaps the number 30 is not so special.
You can change it to any other value that better suits your mood, or remove that part of the script if you want.