diff --git a/hamocc/mo_read_shelfmask.F90 b/hamocc/mo_read_shelfmask.F90
index 2d237231..6bc66ae0 100644
--- a/hamocc/mo_read_shelfmask.F90
+++ b/hamocc/mo_read_shelfmask.F90
@@ -105,7 +105,7 @@ subroutine ini_read_shelfmask(kpie,kpje,kbnd,pbath,omask)
       ncstat=nf90_close(ncid)
     else
       ! reconstruct shelf sea mask from internal bathymetry
-      !$OMP DO PARALLEL PRIVATE (i,j)
+      !$OMP PARALLEL DO PRIVATE (i,j)
       do j=1,kpje
         do i=1,kpie
           if((omask(i,j) > 0.5) .and. (pbath(i,j) <= shelfbreak_depth)) then
@@ -116,7 +116,7 @@ subroutine ini_read_shelfmask(kpie,kpje,kbnd,pbath,omask)
       !$OMP END PARALLEL DO
     endif
 
-    !$OMP DO PARALLEL PRIVATE (i,j)
+    !$OMP PARALLEL DO PRIVATE (i,j)
     ! Eventually fill the logical shelfsea mask field
     do j = 1,kpje
       do i = 1,kpie