Merge pull request #1 from smullangi3/sickle-cell-pr

comp-physics · Aug 1, 2023 · ea82c01 · ea82c01
2 parents 6de9935 + 793c37d
commit ea82c01
Show file tree

Hide file tree

Showing 14 changed files with 357,288 additions and 4 deletions.
diff --git a/case/initcond.F90 b/case/initcond.F90
@@ -154,7 +154,7 @@ subroutine Recenter_Cells_and_Walls
         rbc%x(:,:,1) = rbc%x(:,:,1) + 0.5*Lb(1)
         rbc%x(:,:,2) = rbc%x(:,:,2) + 0.5*Lb(2)
         print*,'Zc', sum(rbc%x(:,:,3))/real(rbc%nlat*rbc%nlon)
-   end do
+    end do
 
     ! walls
     ! find the bounding box
@@ -179,4 +179,4 @@ subroutine Recenter_Cells_and_Walls
 
 end subroutine Recenter_Cells_and_Walls
 
-end program InitCond
+end program InitCond
diff --git a/case_different_celltypes/D/.gittouch b/case_different_celltypes/D/.gittouch
diff --git a/case_different_celltypes/Input/SickleCell.dat b/case_different_celltypes/Input/SickleCell.dat
diff --git a/case_different_celltypes/Input/cylinder.jou b/case_different_celltypes/Input/cylinder.jou
@@ -0,0 +1,10 @@
+create cylinder radius 3.0 height 13.33
+
+curve 1 interval size 0.15
+curve 2 interval size 0.15
+surface 1 scheme trimesh
+
+mesh surface 1
+
+block 1 surface 1
+export mesh "cylinder_D_6.0_L_13.33.g" dimension 3 block 1 overwrite
diff --git a/case_different_celltypes/Input/cylinder_D_3.5_L_13.06.g b/case_different_celltypes/Input/cylinder_D_3.5_L_13.06.g
diff --git a/case_different_celltypes/Input/cylinder_D_6.0_L_13.33.g b/case_different_celltypes/Input/cylinder_D_6.0_L_13.33.g