X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_CSA%2Freadrtns_csa.F;fp=source%2Funres%2Fsrc_CSA%2Freadrtns_csa.F;h=6e222d42a68ea1184a2950521a890f3fb84bbdad;hb=eb0ee5aa9234e4a56580f967363a14e0accfba01;hp=df144f4e7531dfe7012acc09c8b51f93bd4b67f7;hpb=e6be6f3432696603440b91ad5a5732214f003264;p=unres.git diff --git a/source/unres/src_CSA/readrtns_csa.F b/source/unres/src_CSA/readrtns_csa.F index df144f4..6e222d4 100644 --- a/source/unres/src_CSA/readrtns_csa.F +++ b/source/unres/src_CSA/readrtns_csa.F @@ -138,7 +138,7 @@ C Set up the time limit (caution! The time must be input in minutes!) call readi(controlcard,'MAXGEN',maxgen,10000) call readi(controlcard,'MAXOVERLAP',maxoverlap,1000) call readi(controlcard,"KDIAG",kdiag,0) - call readi(controlcard,"RESCALE_MODE",rescale_mode,1) + call readi(controlcard,"RESCALE_MODE",rescale_mode,0) if(me.eq.king .or. .not. out1file .and. fg_rank.eq.0) & write (iout,*) "RESCALE_MODE",rescale_mode split_ene=index(controlcard,'SPLIT_ENE').gt.0 @@ -380,7 +380,7 @@ C bad(i,1)=scalscp*bad(i,1) bad(i,2)=scalscp*bad(i,2) enddo - call rescale_weights(t_bath) +c call rescale_weights(t_bath) if(me.eq.king.or..not.out1file) & write (iout,22) wsc,wscp,welec,wvdwpp,wbond,wang,wscloc,wtor, & wtor_d,wstrain,wel_loc,wcorr,wcorr5,wcorr6,wsccor,wturn3, @@ -577,8 +577,13 @@ 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 @@ -1154,6 +1159,9 @@ c!bankt call readi(minimcard,'MINFUN',minfun,maxmin) call reada(minimcard,'TOLF',tolf,1.0D-2) call reada(minimcard,'RTOLF',rtolf,1.0D-4) + print_min_stat=min0(index(minimcard,'PRINT_MIN_STAT'),1) + print_min_res=min0(index(minimcard,'PRINT_MIN_RES'),1) + print_min_ini=min0(index(minimcard,'PRINT_MIN_INI'),1) write (iout,'(/80(1h*)/20x,a/80(1h*))') & 'Options in energy minimization:' write (iout,'(4(a,i5),a,1pe14.5,a,1pe14.5)')