Commit 50d80604 authored by hb's avatar hb
Browse files

Changing subroutine calling order via MUDFLAT

parent 3665d881
!$Id: m3d.F90,v 1.19 2004-06-15 08:25:57 kbk Exp $
!$Id: m3d.F90,v 1.20 2004-07-28 14:58:18 hb Exp $
#include "cppdefs.h"
!-----------------------------------------------------------------------
!BOP
......@@ -53,7 +53,10 @@
! Original author(s): Karsten Bolding & Hans Burchard
!
! $Log: m3d.F90,v $
! Revision 1.19 2004-06-15 08:25:57 kbk
! Revision 1.20 2004-07-28 14:58:18 hb
! Changing subroutine calling order via MUDFLAT
!
! Revision 1.19 2004/06/15 08:25:57 kbk
! added supoort for spm - Ruiz
!
! Revision 1.18 2004/05/04 09:23:51 kbk
......@@ -389,6 +392,11 @@
!
! !DESCRIPTION:
! A wrapper to call all 3D related subroutines in one subroutine.
! if #MUDFLAT is defined, then the sequence of velocity equations
! and coordinate construction is made such that drying and flooding
! is stable. If #MUDFLAT is not defined, then adaptive grids with
! Lagrangian component are supported. Both, drying and flooding and
! Lagrangian coordinates does not go together.
!
! !REVISION HISTORY:
! See log for module
......@@ -408,6 +416,9 @@
#ifndef NO_BAROCLINIC
if (bdy3d) call do_bdy_3d(0,T)
#endif
#ifdef MUDFLAT
call coordinates(vert_cord,cord_relax,maxdepth)
#endif
#ifndef NO_BOTTFRIC
if (kmax .gt. 1) then
call bottom_friction_3d()
......@@ -429,7 +440,9 @@
call uu_momentum_3d(bdy3d)
ufirst=.true.
end if
#ifndef MUDFLAT
call coordinates(vert_cord,cord_relax,maxdepth)
#endif
if (kmax .gt. 1) then
call ww_momentum_3d()
end if
......
!$Id: uu_momentum_3d.F90,v 1.4 2004-04-20 16:49:37 hb Exp $
!$Id: uu_momentum_3d.F90,v 1.5 2004-07-28 14:58:18 hb Exp $
#include "cppdefs.h"
!-----------------------------------------------------------------------
!BOP
......@@ -10,6 +10,11 @@
!
! !DESCRIPTION:
!
! Three-dimensional velocity equation in eastern direction.
! If #MUDFLAT is defined, fitting of profiles is made with
! respect to the new surface elevation, otherwise to the
! old surface elevation.
!
! !USES:
use parameters, only: g,avmmol,rho_0
use domain, only: iimin,iimax,jjmin,jjmax,kmax,H,HU,min_depth
......@@ -45,7 +50,10 @@
! Original author(s): Hans Burchard & Karsten Bolding
!
! $Log: uu_momentum_3d.F90,v $
! Revision 1.4 2004-04-20 16:49:37 hb
! Revision 1.5 2004-07-28 14:58:18 hb
! Changing subroutine calling order via MUDFLAT
!
! Revision 1.4 2004/04/20 16:49:37 hb
! call to coordinates moved for better consistency (see JMB)
!
! Revision 1.3 2003/04/23 12:16:34 kbk
......@@ -219,8 +227,11 @@
do k=kumin(i,j),kmax
ResInt=ResInt+Res(k)
end do
#ifdef MUDFLAT
Diff=(Uint(i,j)-ResInt)/(ssun(i,j)+HU(i,j))
#else
Diff=(Uint(i,j)-ResInt)/(ssuo(i,j)+HU(i,j))
#endif
do k=kumin(i,j),kmax
uu(i,j,k)=Res(k) +hun(i,j,k)*Diff
......
!$Id: vv_momentum_3d.F90,v 1.5 2004-04-20 16:49:37 hb Exp $
!$Id: vv_momentum_3d.F90,v 1.6 2004-07-28 14:58:18 hb Exp $
#include "cppdefs.h"
!-----------------------------------------------------------------------
!BOP
......@@ -10,6 +10,11 @@
!
! !DESCRIPTION:
!
! Three-dimensional velocity equation in northern direction.
! If #MUDFLAT is defined, fitting of profiles is made with
! respect to the new surface elevation, otherwise to the
! old surface elevation.
!
! !USES:
use parameters, only: g,avmmol,rho_0
use domain, only: iimin,iimax,jjmin,jjmax,kmax,H,HV,min_depth
......@@ -45,7 +50,10 @@
! Original author(s): Hans Burchard & Karsten Bolding
!
! $Log: vv_momentum_3d.F90,v $
! Revision 1.5 2004-04-20 16:49:37 hb
! Revision 1.6 2004-07-28 14:58:18 hb
! Changing subroutine calling order via MUDFLAT
!
! Revision 1.5 2004/04/20 16:49:37 hb
! call to coordinates moved for better consistency (see JMB)
!
! Revision 1.4 2003/06/29 17:06:23 kbk
......@@ -221,7 +229,11 @@
ResInt=ResInt+Res(k)
end do
#ifdef MUDFLAT
Diff=(Vint(i,j)-ResInt)/(ssvn(i,j)+HV(i,j))
#else
Diff=(Vint(i,j)-ResInt)/(ssvo(i,j)+HV(i,j))
#endif
do k=kvmin(i,j),kmax
......
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