AMBER: Calculating cumulative average

From: Francesco Pietra <chiendarret.yahoo.com>
Date: Thu, 2 Aug 2007 10:00:10 -0700 (PDT)

Who would be so kind to adapt the awk script below (from thermodynamic
integration tutorial A6)

BEGIN{ isd=0; epsum=0; npts=0 }

/Molecular dynamics:/{
    while( $1 != "nstlim") getline;
    split($3,a,",");
    nstlim = a[1];
}

/DV\/DL,/{ isd=1}

/NSTEP/{ if(isd == 1 ){nstep = $3 ; tim = $6; }}

/EPtot/ {
    if(isd == 1 ) {
        if( nstep >= nstart ) {
            npts += 1;
            eptot = $9;
            epsum = epsum + $9;
            print tim,eptot,epsum/npts ;
        }
        isd = 0;
    }
}


to get cumulative average for both PM3ESCF and EPtot for QM/MM (tutorial A2). I
should certainly learn the language for such scripts, though I am immediately
interested in the reliability of PM3ESCF in relation to nstlim#.

Thanks a lot

francesco pietra


      ____________________________________________________________________________________
Shape Yahoo! in your own image. Join our Network Research Panel today! http://surveylink.yahoo.com/gmrs/yahoo_panel_invite.asp?a=7


-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Aug 05 2007 - 06:07:36 PDT
Custom Search