Commit b63132be authored by hb's avatar hb
Browse files

Extensions to compiler option SLICE_MODEL for better representation of zero...

Extensions to compiler option SLICE_MODEL for better representation of zero gradients in y-direction
parent af4c532c
!$Id: momentum.F90,v 1.8 2005-10-06 09:54:00 hb Exp $
!$Id: momentum.F90,v 1.9 2006-01-28 20:07:52 hb Exp $
#include "cppdefs.h"
!-----------------------------------------------------------------------
!BOP
......@@ -28,6 +28,9 @@
! Original author(s): Hans Burchard & Karsten Bolding
!
! $Log: momentum.F90,v $
! Revision 1.9 2006-01-28 20:07:52 hb
! Extensions to compiler option SLICE_MODEL for better representation of zero gradients in y-direction
!
! Revision 1.8 2005-10-06 09:54:00 hb
! added support for vertical slice model - via -DSLICE_MODEL
!
......@@ -168,6 +171,12 @@
U=(U-dtm*(g*DU*zx+dry_u*(-tausu/rho_0-fV+UEx+SlUx+Slr)))/(1+dtm*ru/DU)
end where
#ifdef SLICE_MODEL
do i=imin,imax
U(i,3)=U(i,2)
end do
#endif
! now u is calculated
call update_2d_halo(U,U,au,imin,jmin,imax,jmax,U_TAG)
call wait_halo(U_TAG)
......
!$Id: internal_pressure.F90,v 1.11 2005-11-17 13:50:22 kbk Exp $
!$Id: internal_pressure.F90,v 1.12 2006-01-28 20:07:54 hb Exp $
#include "cppdefs.h"
!-----------------------------------------------------------------------
!BOP
......@@ -52,6 +52,9 @@
! Original author(s): Hans Burchard & Karsten Bolding
!
! $Log: internal_pressure.F90,v $
! Revision 1.12 2006-01-28 20:07:54 hb
! Extensions to compiler option SLICE_MODEL for better representation of zero gradients in y-direction
!
! Revision 1.11 2005-11-17 13:50:22 kbk
! fixes to compile with gfortran
!
......@@ -242,6 +245,14 @@
stop 'do_internal_pressure()'
end select
#ifdef SLICE_MODEL
do i=iimin,iimax
do k=kmin(i,2),kmax
idpdx(i,3,k)=idpdx(i,2,k)
end do
end do
#endif
#ifdef SUBSTR_INI_PRESS
if (first) then
first = .false.
......
!$Id: uu_momentum_3d.F90,v 1.5 2004-07-28 14:58:18 hb Exp $
!$Id: uu_momentum_3d.F90,v 1.6 2006-01-28 20:07:54 hb Exp $
#include "cppdefs.h"
!-----------------------------------------------------------------------
!BOP
......@@ -50,6 +50,9 @@
! Original author(s): Hans Burchard & Karsten Bolding
!
! $Log: uu_momentum_3d.F90,v $
! Revision 1.6 2006-01-28 20:07:54 hb
! Extensions to compiler option SLICE_MODEL for better representation of zero gradients in y-direction
!
! Revision 1.5 2004-07-28 14:58:18 hb
! Changing subroutine calling order via MUDFLAT
!
......@@ -241,7 +244,16 @@
end if
end if
end do
end do
#ifdef SLICE_MODEL
do i=iimin,iimax
do k=kumin(i,2),kmax
uu(i,3,k)=uu(i,2,k)
end do
end do
#endif
! Update the halo zones
call update_3d_halo(uu,uu,au,iimin,jjmin,iimax,jjmax,kmax,U_TAG)
......
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