#!/usr/bin/env bash

tnl-grid-setup --dimensions 2 \
               --origin-x 0.0 \
               --origin-y 0.0 \
               --proportions-x 1.0 \
               --proportions-y 1.0 \
               --size-x 400 \
               --size-y 400

#tnl-init --test-function sin-wave \
#         --output-file init.tnl
#              --boundary-conditions-type neumann \
#              --boundary-conditions-constant 0 \

tnl-euler-2d-dbg --time-discretisation explicit \
              --time-step 5.0e-4 \
              --boundary-conditions-type myneumann \
              --discrete-solver euler \
              --snapshot-period 0.03 \
              --final-time 0.3 \
              --NW-density 0.5323 \
              --NW-velocityX 1.206 \
              --NW-velocityY 0.0 \
              --NW-pressure 0.3 \
              --NE-density 1.5 \
              --NE-velocityX 0.0 \
              --NE-velocityY 0.0 \
              --NE-pressure 1.5 \
              --SW-density 0.138 \
              --SW-velocityX 1.206 \
              --SW-velocityY 1.206 \
              --SW-pressure 0.029 \
              --SE-density 0.5323 \
              --SE-velocityX 0 \
              --SE-velocityY 1.206 \
              --SE-pressure 0.3 \
              --gamma 1.4 \
              --riemann-border 0.5 \

tnl-view --mesh mesh.tnl --input-files *tnl     
