Reorder initial conditions for consistency. (#134)

CaNS-World · Jan 15, 2025 · ceae6b5 · ceae6b5
1 parent a39b687
commit ceae6b5
Showing 1 changed file with 15 additions and 15 deletions.
diff --git a/src/initflow.f90 b/src/initflow.f90
@@ -114,6 +114,21 @@ subroutine initflow(inivel,bcvel,ng,lo,l,dl,zc,zf,dzc,dzf,visc, &
           end do
         end do
       end do
+    case('tgw')
+      do k=1,n(3)
+        do j=1,n(2)
+          yc = (j+lo(2)-1-.5)*dl(2)
+          yf = (j+lo(2)-1-.0)*dl(2)
+          do i=1,n(1)
+            xc = (i+lo(1)-1-.5)*dl(1)
+            xf = (i+lo(1)-1-.0)*dl(1)
+            u(i,j,k) =  cos(xf)*sin(yc)*uref
+            v(i,j,k) = -sin(xc)*cos(yf)*uref
+            w(i,j,k) = 0.
+            p(i,j,k) = -(cos(2.*xc)+cos(2.*yc))/4.*uref**2
+          end do
+        end do
+      end do
     case('ant')
       !
       ! see M. Antuono, JFM 890, A23 (2020)
@@ -137,21 +152,6 @@ subroutine initflow(inivel,bcvel,ng,lo,l,dl,zc,zf,dzc,dzf,visc, &
           end do
         end do
       end do
-    case('tgw')
-      do k=1,n(3)
-        do j=1,n(2)
-          yc = (j+lo(2)-1-.5)*dl(2)
-          yf = (j+lo(2)-1-.0)*dl(2)
-          do i=1,n(1)
-            xc = (i+lo(1)-1-.5)*dl(1)
-            xf = (i+lo(1)-1-.0)*dl(1)
-            u(i,j,k) =  cos(xf)*sin(yc)*uref
-            v(i,j,k) = -sin(xc)*cos(yf)*uref
-            w(i,j,k) = 0.
-            p(i,j,k) = -(cos(2.*xc)+cos(2.*yc))/4.*uref**2
-          end do
-        end do
-      end do
     case('pdc','hdc')
       lref  = l(3)/2.
       if(trim(inivel) /= 'pdc') lref = 2.*lref