!:physfluxInitialize:!
        CASE (Initialize)
          ! Allocate local work space
          ALLOCATE(delta(NrVars),coef(mwaves),                                          &
                   sqrtrho(1-mbc:mx+mbc),pres(1-mbc:mx+mbc),rho(1-mbc:mx+mbc),          &
                   u(1-mbc:mx+mbc),v(1-mbc:mx+mbc),h(1-mbc:mx+mbc),c(1-mbc:mx+mbc),     &
                   uR(2-mbc:mx+mbc),vR(2-mbc:mx+mbc),hR(2-mbc:mx+mbc),cR(2-mbc:mx+mbc), &
                   rhosq2(2-mbc:mx+mbc),g1c2(2-mbc:mx+mbc),u2v2(2-mbc:mx+mbc),          &
                   euv(2-mbc:mx+mbc),                                                   &
                   STAT=iError)
          IF (iError/=OK) THEN
            PRINT *,'Cannot allocate work space in problem module, physflux'
            STOP
          END IF