X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc-M%2Fenecalc1.F;h=fee94cf97cb71bb4a710980dc9fcb35217fbf2d2;hb=d101c97dea752458d76055fdbae49c26fff03c1f;hp=ab0a40869b4b9dc6ac613f8adba9a673982a4fe4;hpb=3a74889c3f6514588d40ec326a78f33b3663f5be;p=unres.git diff --git a/source/wham/src-M/enecalc1.F b/source/wham/src-M/enecalc1.F index ab0a408..fee94cf 100644 --- a/source/wham/src-M/enecalc1.F +++ b/source/wham/src-M/enecalc1.F @@ -60,6 +60,8 @@ enddo enddo enddo + write (iout,*) "indstart(me1),indend(me1)" + &,indstart(me1),indend(me1) do i=indstart(me1),indend(me1) #else do iparm=1,nParmSet @@ -71,6 +73,7 @@ enddo do i=1,ntot #endif + read(ientout,rec=i,err=101) & ((csingle(l,k),l=1,3),k=1,nres), & ((csingle(l,k+nres),l=1,3),k=nnt,nct), @@ -154,7 +157,10 @@ c & " kfac",kfac,"quot",quot," fT",fT & wcorr,wcorr5,wcorr6,wturn4,wturn3,wturn6,wel_loc, & wtor_d,wsccor,wbond #endif +C write (iout,*) "tuz przed energia" call etotal(energia(0),fT) +C write (iout,*) "tuz za energia" + #ifdef DEBUG write (iout,*) "Conformation",i write (iout,'(8f10.5)') ((c(l,k),l=1,3),k=1,nres) @@ -220,6 +226,7 @@ c call intout endif endif endif +C write (iout,*) "Czy tu dochodze" potE(iii+1,iparm)=energia(0) do k=1,21 enetb(k,iii+1,iparm)=energia(k)