     np=Hp->ny;
     nt=Hp->nx;
     dx=Hp->dy;
-    dy=Hp->dx;
     dCh='y';
     setVHalo(bndT,bndB);
     toPrimY<<<BL_TH((np+4)*nt,nTh)>>>(d_q,d_u);

    np=Hp->nx;
    nt=Hp->ny;
    dx=Hp->dx;
    dy=Hp->dy;
    dCh='x';
    setHHalo(Hp->bndL,Hp->bndR);
    toPrimX<<<BL_TH((np+4)*nt,nTh)>>>(d_q,d_u);


void runPass(double dt, int dir){
  int np,nt;
  double dx;
  char dCh;
  char outLab[30];

