Commit 3c06cb28 authored by kbk's avatar kbk
Browse files

using HALO in loop boundaries

parent ccce9bf4
!$Id: depth_update.F90,v 1.4 2003-05-12 09:22:39 kbk Exp $ !$Id: depth_update.F90,v 1.5 2003-06-18 08:27:41 kbk Exp $
#include "cppdefs.h" #include "cppdefs.h"
!----------------------------------------------------------------------- !-----------------------------------------------------------------------
!BOP !BOP
...@@ -26,7 +26,10 @@ ...@@ -26,7 +26,10 @@
! Original author(s): Hans Burchard & Karsten Bolding ! Original author(s): Hans Burchard & Karsten Bolding
! !
! $Log: depth_update.F90,v $ ! $Log: depth_update.F90,v $
! Revision 1.4 2003-05-12 09:22:39 kbk ! Revision 1.5 2003-06-18 08:27:41 kbk
! using HALO in loop boundaries
!
! Revision 1.4 2003/05/12 09:22:39 kbk
! no use of update_2d_halo, expand loop boundaries instead ! no use of update_2d_halo, expand loop boundaries instead
! !
! Revision 1.3 2003/04/23 12:09:43 kbk ! Revision 1.3 2003/04/23 12:09:43 kbk
...@@ -75,8 +78,8 @@ ...@@ -75,8 +78,8 @@
D = z+H D = z+H
! U-points ! U-points
do j=jmin-1,jmax+1 do j=jmin-HALO,jmax+HALO
do i=imin-1,imax+1 do i=imin-HALO,imax+HALO-1
#ifdef USE_MASK #ifdef USE_MASK
if(au(i,j) .gt. 0) then if(au(i,j) .gt. 0) then
#endif #endif
...@@ -90,8 +93,8 @@ ...@@ -90,8 +93,8 @@
end do end do
! V-points ! V-points
do j=jmin-1,jmax+1 do j=jmin-HALO,jmax+HALO-1
do i=imin-1,imax+1 do i=imin-HALO,imax+HALO
#ifdef USE_MASK #ifdef USE_MASK
if(av(i,j) .gt. 0) then if(av(i,j) .gt. 0) then
#endif #endif
......
!$Id: uv_depths.F90,v 1.6 2003-05-12 09:22:39 kbk Exp $ !$Id: uv_depths.F90,v 1.7 2003-06-18 08:27:41 kbk Exp $
#include "cppdefs.h" #include "cppdefs.h"
!----------------------------------------------------------------------- !-----------------------------------------------------------------------
!BOP !BOP
...@@ -25,7 +25,10 @@ ...@@ -25,7 +25,10 @@
! Original author(s): Hans Burchard & Karsten Bolding ! Original author(s): Hans Burchard & Karsten Bolding
! !
! $Log: uv_depths.F90,v $ ! $Log: uv_depths.F90,v $
! Revision 1.6 2003-05-12 09:22:39 kbk ! Revision 1.7 2003-06-18 08:27:41 kbk
! using HALO in loop boundaries
!
! Revision 1.6 2003/05/12 09:22:39 kbk
! no use of update_2d_halo, expand loop boundaries instead ! no use of update_2d_halo, expand loop boundaries instead
! !
! Revision 1.5 2003/04/23 12:09:44 kbk ! Revision 1.5 2003/04/23 12:09:44 kbk
...@@ -61,15 +64,14 @@ ...@@ -61,15 +64,14 @@
write(0,*) 'uv_depths() # ',Ncall write(0,*) 'uv_depths() # ',Ncall
#endif #endif
do j=jmin-1,jmax+1 do j=jmin-HALO,jmax+HALO
do i=imin-1,imax+1 do i=imin-HALO,imax+HALO-1
#ifdef MIN_VEL_DEPTH
HU(i,j)=min(H(i,j),H(i+1,j))
HV(i,j)=min(H(i,j),H(i,j+1))
#else
HU(i,j)=0.5*(H(i,j)+H(i+1,j)) HU(i,j)=0.5*(H(i,j)+H(i+1,j))
end do
end do
do j=jmin-HALO,jmax+HALO-1
do i=imin-HALO,imax+HALO
HV(i,j)=0.5*(H(i,j)+H(i,j+1)) HV(i,j)=0.5*(H(i,j)+H(i,j+1))
#endif
end do end do
end do end do
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment