parameter (im=72,jm=59,kseg=24,mm=kseg+2,kb=21)
       real kth,kthD,kp,kaveD
       common/blkwave/
C pom variables
     &    dx(im,jm),dy(im,jm),z(kb),dz(kb),zz(kb),dzz(kb),
     &    d(im,jm),fsm(im,jm),dum(im,jm),dvm(im,jm),
     &    u10x(im,jm),u10y(im,jm),u10(im,jm),u(im,jm,kb),v(im,jm,kb),
     &    iint,iprint,iend,
C wave variables
     &    ent(im,jm),thtav(im,jm),Hs(im,jm),
     &    aren(im,jm),areninv(im,jm),x(im),y(jm), 
     &    cs(mm),sn(mm),theta(mm),thetam(mm),csm(mm),snm(mm),
     &    enw(jm,mm),ens(im,mm),ene(jm,mm),enn(im,mm),
     &    sigthw(jm,mm),sigthe(jm,mm),
     &    sigthf(im,jm,mm),sigth(im,jm,mm),sigthb(im,jm,mm),
     &    cth(im,jm,mm),cg(im,jm,mm),kth(im,jm,mm),kthD(im,jm,mm),
     &    kaveD(im,jm),cgx(im,jm,mm),cgy(im,jm,mm),
     &    age(im,jm),ageinv(im,jm),tht_wnd(im,jm),
     &    sigp(im,jm),ustw(im,jm),cd(im,jm),
     &    cp(im,jm),kp(im,jm), cthth(im,jm,mm),
     &    enb(im,jm,mm),en(im,jm,mm),enf(im,jm,mm),  
     &    sdis(im,jm,mm),swin(im,jm,mm),
     &    tpx0(im,jm),tpy0(im,jm),ttx0(im,jm),tty0(im,jm),
     &    qdis(im,jm),Hm(im,jm),bdis1(im,jm,mm),bdis2(im,jm,mm),
     &    sdistot(im,jm),bdistot(im,jm),bdiss(mm),
     &    tps2d(im,jm)