Forgot to add a separate case for the 2D cases. Everything works now

Author: danieljvickers

examples/2D_zero_circ_vortex/case.py

@@ -85,7 +85,7 @@
             "patch_icpp(2)%vel(2)": 0.0,
             "patch_icpp(2)%pres": 0.0,
             "patch_icpp(2)%alpha_rho(1)": 0.0,
-            "patch_icpp(2)%hcid": 280,
+            "patch_icpp(2)%hcid": 282,
             "patch_icpp(2)%alpha(1)": 1.0,
             "patch_icpp(2)%alter_patch(1)": "T",
             # CBC Inflow / Outflow

src/pre_process/include/2dHardcodedIC.fpp

@@ -163,7 +163,7 @@
 
     case (280)
         ! This is patch is hard-coded for test suite optimization used in the
-        ! 2D_zero_circ_vortex case:
+        ! 2D_isentropicvortex case:
         ! This analytic patch uses geometry 2
         if (patch_id == 1) then
             q_prim_vf(E_idx)%sf(i, j, 0) = 1d0*(1d0 - (1d0/1d0)*(5d0/(2d0*3.141592653589793))*(5d0/(8d0*1d0*(1.4 + 1d0)*3.141592653589793))*exp(2d0*1d0*(1d0 - (x_cc(i) - patch_icpp(1)%x_centroid)**2d0- (y_cc(j) - patch_icpp(1)%y_centroid)**2d0)))**(1.4 + 1d0)
@@ -181,6 +181,17 @@
             q_prim_vf(contxb + 0)%sf(i, j, 0) = 1*(1 - 0.5*(1.4 - 1)*(0.4)**2*exp(0.5*(1 - sqrt(x_cc(i)**2 + y_cc(j)**2))))**(1 / (1.4 - 1))
         end if
 
+    case (282)
+        ! This is patch is hard-coded for test suite optimization used in the
+        ! 2D_zero_circ_vortex case:
+        ! This analytic patch uses geometry 2
+        if (patch_id == 2) then
+            q_prim_vf(E_idx)%sf(i, j, 0) = 101325*(1 - 0.5*(1.4 - 1)*(0.1 / 0.3)**2*exp(0.5*(1 - sqrt(x_cc(i)**2 + y_cc(j)**2))))**(1.4 / (1.4 - 1))
+            q_prim_vf(contxb + 0)%sf(i, j, 0) = 1*(1 - 0.5*(1.4 - 1)*(0.1 / 0.3)**2*exp(0.5*(1 - sqrt(x_cc(i)**2 + y_cc(j)**2))))**(1 / (1.4 - 1))
+            q_prim_vf(momxb + 0)%sf(i, j, 0) = 112.99092883944267*(1 - (0.1 / 0.3))*y_cc(j)*exp(0.5*(1 - sqrt(x_cc(i)**2 + y_cc(j)**2))) 
+            q_prim_vf(momxb + 1)%sf(i, j, 0) = 112.99092883944267*((0.1 / 0.3))*x_cc(i)*exp(0.5*(1 - sqrt(x_cc(i)**2 + y_cc(j)**2)))
+        end if
+
     case default
         if (proc_rank == 0) then
             call s_int_to_str(patch_id, iStr)