openbc: remove call to boundary on init, warn if synturb parameters a…

…re bad
dalesteam · Dec 9, 2024 · 10c9cc2 · 10c9cc2
1 parent c7d95e1
commit 10c9cc2
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/src/modstartup.f90 b/src/modstartup.f90
@@ -844,7 +844,6 @@ subroutine readinitfiles
 
       call thermodynamics
       call surface
-      call boundary
 
       if ( lopenbc ) then
         call openboundary_ghost()

diff --git a/src/modsynturb.f90 b/src/modsynturb.f90
@@ -72,6 +72,9 @@ subroutine initsynturb
         nxturb = int(dx/dxturb*real(itot));
         nyturb = int(dy/dyturb*real(jtot));
         nzturb = kmax
+        if (nxturb < 1 .or. nyturb < 1) then
+           STOP 'nxturb or nyturb < 1, perhaps dxturb or dyturb is too large.'
+        end if
         lambdas = merge(lambda,lambdas,lambdas==-1.)
         lambdasxyz = (/merge(lambdas,lambdas_x,lambdas_x==-1.), &
                     & merge(lambdas,lambdas_y,lambdas_y==-1.), &