source/cluster/wham/src/COMMON.MCM

   1 C... Following COMMON block contains general variables controlling the MC/MCM
   2 C... procedure
   3 c-----------------------------------------------------------------------------
   4       double precision Tcur,Tmin,Tmax,TstepH,TstepC,RanFract,
   5      &        overlap_cut,e_up,delte
   6       integer nstepH,nstepC,maxacc,maxgen,maxtrial,maxtrial_iter,
   7      &        maxrepm,ngen,ntrial,ntherm,nrepm,neneval,nsave,maxoverlap,
   8      &        nsave_part,max_mcm_it,nsweep,print_mc
   9       logical print_stat,print_int
  10       common /mcm/ Tcur,Tmin,Tmax,TstepH,TstepC,Rbol,betbol,RanFract,
  11      & overlap_cut,e_up,delte,
  12      & nstepH,nstepC,maxacc,maxgen,maxtrial,maxtrial_iter,maxrepm,
  13      & maxoverlap,ntrial,max_mcm_it,
  14      & ngen,ntherm,nrepm,neneval,nsave,nsave_part(max_cg_procs),nsweep,
  15      & print_mc,print_stat,print_int
  16 c-----------------------------------------------------------------------------
  17 C... The meaning of the above variables is as follows:
  18 C... Tcur,Tmin,Tmax - Current,minimum and maximum temperature, respectively;
  19 C... NstepC,NStepH - Number of cooling and heating steps, respectively;
  20 C... TstepH,TstepC - factors by which T is multiplied in order to be
  21 C...                 increased or decreased.
  22 C... betbol - Boltzmann's inverse temperature (1/(Rbol*Tcur));
  23 C... Rbol - the gas constant;
  24 C... RanFract - the chance that a new conformation will be random-generated;
  25 C... maxacc - maximum number of accepted conformations;
  26 C... maxgen,ngen - Maximum and current number of generated conformations;
  27 C... maxtrial,ntrial - maximum number of trials before temperature is increased
  28 C...                   and current number of trials, respectively;
  29 C... maxrepm,nrepm - maximum number of allowed minima repetition and current
  30 C...                 number of minima repetitions, respectively;
  31 C... maxoverlap - max. # of overlapping confs generated in a single iteration;
  32 C... neneval - number of energy evaluations;
  33 C... nsave - number of confs. in the backup array;
  34 C... nsweep - the number of macroiterations in generating the distributions.
  35 c------------------------------------------------------------------------------
  36 C... Following COMMON block contains variables controlling motion.
  37 c------------------------------------------------------------------------------
  38       double precision sumpro_type,sumpro_bond
  39       integer koniecl, Nbm,MaxSideMove,nmove,moves(-1:MaxMoveType+1),
  40      &   moves_acc(-1:MaxMoveType+1),nacc_tot,nacc_part(0:MaxProcs)
  41       common /move/ sumpro_type(0:MaxMoveType),sumpro_bond(0:maxres),
  42      & koniecl,Nbm,MaxSideMove,nmove,nbond_move(maxres),
  43      & nbond_acc(maxres),moves,moves_acc
  44       common /accept_stats/ nacc_tot,nacc_part
  45       integer nwindow,winstart,winend,winlen
  46       common /windows/ nwindow,winstart(maxres),winend(maxres),
  47      &        winlen(maxres)
  48       character*16 MovTypID
  49       common /moveID/ MovTypID(-1:MaxMoveType+1)
  50 c------------------------------------------------------------------------------
  51 C... koniecl - the number of bonds to be considered "end bonds" subjected to
  52 C...          end moves;
  53 C... Nbm - The maximum length of N-bond segment to be moved;
  54 C... MaxSideMove - maximum number of side chains subjected to local moves
  55 C...               simultaneously;
  56 C... nmove - the current number of attempted moves;
  57 C... nbond_move(*) array that stores the total numbers of 2-bond,3-bond,...
  58 C...            moves;
  59 C... nendmove - number of endmoves;
  60 C... nbackmove - number of backbone moves;
  61 C... nsidemove - number of local side chain moves;
  62 C... sumpro_type(*) - array that stores the lower and upper boundary of the
  63 C...                  random-number range that determines the type of move
  64 C...                  (N-bond, backbone or side chain);
  65 C... sumpro_bond(*) - array that stores the probabilities to perform bond
  66 C...                  moves of consecutive segment length.
  67 C... winstart(*) - the starting position of the perturbation window;
  68 C... winend(*) -  the end position of the perturbation window;
  69 C... winlen(*) - length of the perturbation window;
  70 C... nwindow - the number of perturbation windows (0 - entire chain).