X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc%2Fparmread.F;h=36730f5a7bfa12512c7c19cb17ff8d5f5cae3004;hb=b3332a0737da32abddd105126f2c43a69da0419d;hp=aae7c015c36683853f83a43e5bec7ff11e2cc43e;hpb=34d9af87d8a30caeb12b841c0310ced100c44970;p=unres.git diff --git a/source/wham/src/parmread.F b/source/wham/src/parmread.F index aae7c01..36730f5 100644 --- a/source/wham/src/parmread.F +++ b/source/wham/src/parmread.F @@ -961,7 +961,7 @@ C C C Define the constants of the disulfide bridge C - ebr=-5.50D0 +c ebr=-5.50D0 c c Old arbitrary potential - commented out. c @@ -972,21 +972,21 @@ c Constants of the disulfide-bond potential determined based on the RHF/6-31G** c energy surface of diethyl disulfide. c A. Liwo and U. Kozlowska, 11/24/03 c - D0CM = 3.78d0 - AKCM = 15.1d0 - AKTH = 11.0d0 - AKCT = 12.0d0 - V1SS =-1.08d0 - V2SS = 7.61d0 - V3SS = 13.7d0 +c D0CM = 3.78d0 +c AKCM = 15.1d0 +c AKTH = 11.0d0 +c AKCT = 12.0d0 +c V1SS =-1.08d0 +c V2SS = 7.61d0 +c V3SS = 13.7d0 - if (lprint) then - write (iout,'(/a)') "Disulfide bridge parameters:" - write (iout,'(a,f10.2)') 'S-S bridge energy: ',ebr - write (iout,'(2(a,f10.2))') 'd0cm:',d0cm,' akcm:',akcm - write (iout,'(2(a,f10.2))') 'akth:',akth,' akct:',akct - write (iout,'(3(a,f10.2))') 'v1ss:',v1ss,' v2ss:',v2ss, - & ' v3ss:',v3ss - endif +c if (lprint) then +c write (iout,'(/a)') "Disulfide bridge parameters:" +c write (iout,'(a,f10.2)') 'S-S bridge energy: ',ebr +c write (iout,'(2(a,f10.2))') 'd0cm:',d0cm,' akcm:',akcm +c write (iout,'(2(a,f10.2))') 'akth:',akth,' akct:',akct +c write (iout,'(3(a,f10.2))') 'v1ss:',v1ss,' v2ss:',v2ss, +c & ' v3ss:',v3ss +c endif return end