small updates

Author: wilfonba

docs/documentation/authors.md

@@ -17,3 +17,4 @@ We try to maintain a list of current and past developers:
 + M. Rodriguez 
 + K. Schmidmayer 
 + J. Spratt 
++ B. Wilfong    

docs/documentation/case.md

@@ -205,15 +205,14 @@ end if
 
 #### Hard Coded Patches
 
-Some patch configurations are not adequatley handeled with the above analytic variable definitions. In this case, a hard coded patch can be used. Hard coded patches can be added by adding additional hard coded patch identifers to `src/pre_process/include/1[2,3]dHardcodedIC.fpp`. For example, to add a 2D Hardcoded patch with an id of 201, one would add the following to `src/pre_process/include/2dHardcodedIC.fpp`
+Some patch configurations are not adequatley handeled with the above analytic variable definitions. In this case, a hard coded patch can be used. Hard coded patches can be added by adding additional hard coded patch identifers to `src/pre_process/include/1[2,3]dHardcodedIC.fpp`. For example, to add a 2D Hardcoded patch with an id of 200, one would add the following to `src/pre_process/include/2dHardcodedIC.fpp`
 
 ```f90
-    if (patch_icpp%hcid == 201) then
+    case(200)
         ! Primitive variables assignment
-    end if
 ```
 
-and use `patch_icpp(i)%geometry = y` and `patch_icpp(i)%hcid = 201` in the input file. As a convention, any hard coded patches that are part of the MFC master branch should be identied as 1[2,3]xx where the first digit indiates the number of dimensions.
+and use `patch_icpp(i)%geometry = 7` and `patch_icpp(i)%hcid = 200` in the input file. Additional variables can be declared in `Hardcoded1[2,3]DVariables` and used in `hardcoded1[2,3]D`. As a convention, any hard coded patches that are part of the MFC master branch should be identied as 1[2,3]xx where the first digit indiates the number of dimensions.
 
 #### Parameter Descriptions
 

examples/2D_hardcodied_ic/README.md

@@ -0,0 +1,9 @@
+# 2D Hardcodied IC Example
+
+## Initial Condition
+
+![Initial Condition](initial.png)
+
+## Result
+
+![Result](result.png)

examples/2D_hardcodied_ic/case.py

@@ -26,10 +26,10 @@
     'm'                            : 199,
     'n'                            : 199,
     'p'                            : 0,
-    'dt'                           : 5.E-07,
+    'dt'                           : 5.E-06,
     't_step_start'                 : 0,
-    't_step_stop'                  : 100000,
-    't_step_save'                  : 1000,
+    't_step_stop'                  : 2000,
+    't_step_save'                  : 200,
     # ==========================================================
 
     # Simulation Algorithm Parameters ==========================
@@ -50,9 +50,9 @@
     'wave_speeds'                  : 1,
     'avg_state'                    : 2,
     'bc_x%beg'                     : -2,
-    'bc_x%end'                     : -2,
-    'bc_y%beg'                     : -6,
-    'bc_y%end'                     : -6,
+    'bc_x%end'                     : -7,
+    'bc_y%beg'                     : -2,
+    'bc_y%end'                     : -7,
     # ==========================================================
 
     # Formatted Database Files Structure Parameters ============
@@ -71,7 +71,7 @@
     'patch_icpp(1)%length_y'       : 8.,
     'patch_icpp(1)%vel(1)'         : 0.,
     'patch_icpp(1)%vel(2)'         : 0.,
-    'patch_icpp(1)%pres'           : 100,
+    'patch_icpp(1)%pres'           : 10,
     'patch_icpp(1)%alpha_rho(1)'   : (1-eps)*1000,
     'patch_icpp(1)%alpha_rho(2)'   : eps*1,
     'patch_icpp(1)%alpha(1)'       : 1-eps,

examples/2D_hardcodied_ic/initial.png

@@ -127,7 +127,7 @@ subroutine s_check_input_file() ! --------------------------------------
 
     end subroutine s_check_input_file ! ------------------------------------
 
-        subroutine s_perform_time_step(t_step)
+    subroutine s_perform_time_step(t_step)
 
         integer, intent(INOUT) :: t_step 
         if (proc_rank == 0) then

examples/2D_hardcodied_ic/result.png

@@ -6,8 +6,12 @@
 
 #:def Hardcoded1D()
 
-    if (patch_icpp(patch_id)%hcid == 100) then
-
-    end if
+    select case(patch_icpp(patch_id)%hcid)
+        case(100)
+            ! Put your variable assignments here
+        case default
+            call s_int_to_str(patch_id, iStr)
+            call s_mpi_abort("Invalid hcid specified for patch " // trim(iStr))
+    end select
 
 #:enddef

src/post_process/m_start_up.f90

@@ -7,20 +7,24 @@
 #:enddef
 
 #:def Hardcoded2D()
-
-    if (patch_icpp(patch_id)%hcid == 200) then ! 2D_hardcoded_ic example case
-
-        if (y_cc(j) <= (-x_cc(i)**3 + 1)**(1d0/3d0)) then
-            ! Volume Fractions
-            q_prim_vf(advxb)%sf(i, j, 0) = eps
-            q_prim_vf(advxe)%sf(i, j, 0) = 1d0-eps
-            ! Denssities
-            q_prim_vf(contxb)%sf(i, j, 0) = eps*1000d0
-            q_prim_vf(contxe)%sf(i, j, 0) = (1d0-eps)*1d0
-            ! Pressure
-            q_prim_vf(E_idx)%sf(i, j, 0) = 1d0
-        end if
-
-    end if
+    
+    select case(patch_icpp(patch_id)%hcid) ! 2D_hardcoded_ic example case
+        case(200)
+            if (y_cc(j) <= (-x_cc(i)**3 + 1)**(1d0/3d0)) then
+                ! Volume Fractions
+                q_prim_vf(advxb)%sf(i, j, 0) = eps
+                q_prim_vf(advxe)%sf(i, j, 0) = 1d0-eps
+                ! Denssities
+                q_prim_vf(contxb)%sf(i, j, 0) = eps*1000d0
+                q_prim_vf(contxe)%sf(i, j, 0) = (1d0-eps)*1d0
+                ! Pressure
+                q_prim_vf(E_idx)%sf(i, j, 0) = 1000d0
+            end if
+        case default
+            if (proc_rank == 0) then
+                call s_int_to_str(patch_id, iStr)
+                call s_mpi_abort("Invalid hcid specified for patch " // trim(iStr))
+            end if
+    end select
 
 #:enddef

src/pre_process/include/1dHardcodedIC.fpp

@@ -6,8 +6,12 @@
 
 #:def Hardcoded3D()
 
-    if (patch_icpp(patch_id)%hcid == 300) then
-
-    end if
+    select case(patch_icpp(patch_id)%hcid)
+        case(300)
+            ! Put your variable assignments here
+        case default
+            call s_int_to_str(patch_id, iStr)
+            call s_mpi_abort("Invalid hcid specified for patch " // trim(iStr))
+    end select
 
 #:enddef