PRO test1_p3,totmom,q=q,h=h
; IDL Version 4.0.1 (Win32 x86)
; Journal File for user@Microsoft Windows Host
; Working directory: C:\RSI\IDL40
; Date: Fri Feb 02 10:01:27 1996

gam=9810.D0
IF NOT KEYWORD_SET(q) THEN q=0.55D0 ELSE q=DOUBLE(q)
l=1.0D0
IF NOT KEYWORD_SET(h) THEN h=(7.D0/6.D0)*l ELSE h=DOUBLE(h)

abot=(1.D0-q)*l*l
atop=q*l*l

ybartop=h+q*l/2.D0
ybarbot=h+q*l+(1.D0-q)*l/2.D0

ftop=gam*ybartop*atop
fbot=gam*ybarbot*abot

ibartop=l*((q*l)^3)/12.D0
ibarbot=l*(((1.D0-q)*l)^3)/12.D0

ycptop=ybartop+ibartop/(ybartop*atop)
ycpbot=ybarbot+ibarbot/(ybarbot*abot)

toparm=(h+q*l)-ycptop
botarm=ycpbot-(h+q*l)

cwmom=ftop*toparm
ccwmom=fbot*botarm
totmom=ccwmom-cwmom

relmom=totmom/(cwmom+ccwmom)

RETURN

END