!--------------------------------------------------------------------------
   !---------------------------- coriolis ----------------------------------
   !
   DO iq=1,nqdyn
      !$OMP DO SCHEDULE(STATIC)
      DO edge = 1, edge_num
         ij_left = left(edge)
         ij_right = right(edge)
         DO l = 1, llm
            Ftheta(l,edge) = .5*(theta(l,ij_left,iq)+theta(l,ij_right,iq))*hflux(l,edge)
         END DO
      END DO
      !$OMP END DO
      !$OMP DO SCHEDULE(STATIC)
      DO ij = 1, primal_num
         DO l = 1, llm
            divF=0.
            DO iedge = 1, primal_deg(ij)
               edge = primal_edge(iedge,ij)
               divF = divF + Ftheta(l,edge)*primal_ne(iedge,ij)
            END DO
            dtheta_rhodz(l,ij,iq) = -divF / Ai(ij)
         END DO
      END DO
      !$OMP END DO
   END DO ! iq
   !
   !$OMP DO SCHEDULE(STATIC)
   DO ij = 1, primal_num
      DO l = 1, llm
         divF=0.
         DO iedge = 1, primal_deg(ij)
            edge = primal_edge(iedge,ij)
            divF = divF + hflux(l,edge)*primal_ne(iedge,ij)
         END DO
         convm(l,ij) = -divF / Ai(ij)
      END DO
   END DO
   !$OMP END DO
   !
   !$OMP DO SCHEDULE(STATIC)
   DO edge = 1, edge_num
      DO l = 1, llm
         du_trisk=0.
         DO itrisk = 1, trisk_deg(edge)
            edge_trisk = trisk(itrisk,edge)
            du_trisk = du_trisk + wee(itrisk,edge)*hflux(l,edge_trisk)*(qu(l,edge)+qu(l,edge_trisk))
         END DO
         du(l,edge) = du(l,edge) + .5*du_trisk
      END DO
   END DO
   !$OMP END DO
   !---------------------------- coriolis ----------------------------------
   !--------------------------------------------------------------------------