X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_CSA%2Freadrtns_csa.F;h=9ab45102c12aa482966d11726fe08d40fe94de89;hb=4d5f4674e86e6debc71c19e681c3c57e075f231d;hp=6e222d42a68ea1184a2950521a890f3fb84bbdad;hpb=951f0fe4ce7a66542c398754524d8ae3d9b5dc5a;p=unres.git diff --git a/source/unres/src_CSA/readrtns_csa.F b/source/unres/src_CSA/readrtns_csa.F index 6e222d4..9ab4510 100644 --- a/source/unres/src_CSA/readrtns_csa.F +++ b/source/unres/src_CSA/readrtns_csa.F @@ -455,6 +455,9 @@ c print *,'Finished reading pdb data' call contact(.false.,ncont_ref,icont_ref,co) if (sideadd) then +C Following 2 lines for diagnostics; comment out if not needed + write (iout,*) "Before sideadd" + call intout if(me.eq.king.or..not.out1file) & write(iout,*)'Adding sidechains' maxsi=1000 @@ -472,6 +475,11 @@ c call gen_side(iti,theta(i+1),alph(i),omeg(i),fail) endif enddo endif +C 10/03/12 Adam: Recalculate coordinates with new side chain positions + call chainbuild +C Following 2 lines for diagnostics; comment out if not needed + write (iout,*) "After sideadd" + call intout endif if (indpdb.eq.0) then @@ -577,13 +585,8 @@ C Juyong:READ read_info C READ fragment information!! C both routines should be in dfa.F file!! -<<<<<<< HEAD if (.not. (wdfa_dist.eq.0.0 .and. wdfa_tor.eq.0.0 .and. & wdfa_nei.eq.0.0 .and. wdfa_beta.eq.0.0)) then -======= - if (.not. (wdfa_dist.eq.0.0 .and. wdfa_tor.eq.0.0 - & .and. wdfa_nei.eq.0.0 .and. wdfa_beta.eq.0.0)) then ->>>>>>> adam call init_dfa_vars print*, 'init_dfa_vars finished!' call read_dfa_info