Commit 7762093e authored by kbk's avatar kbk
Browse files

back ground viscosity and diffusivity (manuel)

parent 55079a24
!$Id: gotm.F90,v 1.5 2003-05-05 15:51:57 kbk Exp $
!$Id: gotm.F90,v 1.6 2003-12-16 15:58:54 kbk Exp $
#include "cppdefs.h"
!-----------------------------------------------------------------------
!BOP
......@@ -15,6 +15,7 @@
use domain, only: iimin,iimax,jjmin,jjmax,kmax,az,min_depth,crit_depth
use variables_2d, only: D,zub,zvb,z
use variables_3d, only: dt,kmin,ho,hn,tke,eps,SS,NN,num,nuh,taus,taub
use variables_3d, only: avmback,avhback
use turbulence, only: do_turbulence,cde
use turbulence, only: tke1d => tke, eps1d => eps, L1d => L
use turbulence, only: num1d => num, nuh1d => nuh
......@@ -30,7 +31,10 @@
! Original author(s): Karsten Bolding & Hans Burchard
!
! $Log: gotm.F90,v $
! Revision 1.5 2003-05-05 15:51:57 kbk
! Revision 1.6 2003-12-16 15:58:54 kbk
! back ground viscosity and diffusivity (manuel)
!
! Revision 1.5 2003/05/05 15:51:57 kbk
! proper update of halo zones for num and nuh
!
! Revision 1.4 2003/04/23 12:16:34 kbk
......@@ -96,6 +100,11 @@
write(debug,*) 'gotm() # ',Ncall
#endif
#ifdef CONSTANT_VISCOSITY
num(i,j,:) = avmback
nuh(i,j,:) = avhback
#else
do j=jjmin,jjmax
do i=iimin,iimax
......@@ -138,23 +147,16 @@
tke(i,j,:) = tke1d
eps(i,j,:) = eps1d
num(i,j,:) = num1d
nuh(i,j,:) = nuh1d
#ifdef HAIDVOGEL_TEST
num(i,j,:) = 0.e-10
nuh(i,j,:) = 0.e-10
#endif
#ifdef NOMADS_TEST
num(i,j,:) = _ZERO_
nuh(i,j,:) = _ZERO_
#endif
num(i,j,:) = num1d + avmback
nuh(i,j,:) = nuh1d + avhback
end if
end do
end do
call update_3d_halo(num,num,az,iimin,jjmin,iimax,jjmax,kmax,H_TAG)
call wait_halo(H_TAG)
#endif
call update_3d_halo(nuh,nuh,az,iimin,jjmin,iimax,jjmax,kmax,H_TAG)
call wait_halo(H_TAG)
......
!$Id: m3d.F90,v 1.7 2003-09-12 16:23:38 kbk Exp $
!$Id: m3d.F90,v 1.8 2003-12-16 15:58:54 kbk Exp $
#include "cppdefs.h"
!-----------------------------------------------------------------------
!BOP
......@@ -33,15 +33,12 @@
use variables_3d
use advection_3d, only: init_advection_3d
use bdy_3d, only: init_bdy_3d, do_bdy_3d
#ifdef PARALLEL
use halo_mpi
#endif
IMPLICIT NONE
!
! !PUBLIC DATA MEMBERS:
integer :: M=1
REALTYPE :: cord_relax=_ZERO_
REALTYPE :: cord_relax=_ZERO_
logical :: calc_temp=.true.
logical :: calc_salt=.true.
logical :: calc_spm=.false.
......@@ -53,7 +50,10 @@ use halo_mpi
! Original author(s): Karsten Bolding & Hans Burchard
!
! $Log: m3d.F90,v $
! Revision 1.7 2003-09-12 16:23:38 kbk
! Revision 1.8 2003-12-16 15:58:54 kbk
! back ground viscosity and diffusivity (manuel)
!
! Revision 1.7 2003/09/12 16:23:38 kbk
! fixed order of use statement - now compiles on with IFC 7.1
!
! Revision 1.6 2003/08/28 15:13:57 kbk
......@@ -200,7 +200,8 @@ use halo_mpi
M,cnpar,cord_relax, &
bdy3d,bdyfmt_3d,bdyramp_3d,bdyfile_3d, &
vel_hor_adv,vel_ver_adv,vel_strang, &
calc_temp,calc_salt,calc_spm
calc_temp,calc_salt,calc_spm, &
avmback,avhback
!
!EOP
!-------------------------------------------------------------------------
......
!$Id: variables_3d.F90,v 1.4 2003-04-23 12:16:34 kbk Exp $
!$Id: variables_3d.F90,v 1.5 2003-12-16 15:58:54 kbk Exp $
#include "cppdefs.h"
!-----------------------------------------------------------------------
!BOP
......@@ -20,6 +20,7 @@
!
! !PUBLIC DATA MEMBERS:
REALTYPE :: dt,cnpar=0.9
REALTYPE :: avmback=_ZERO_,avhback=_ZERO_
!
#ifdef STATIC
#include "static_3d.h"
......@@ -33,7 +34,10 @@
! Original author(s): Karsten Bolding & Hans Burchard
!
! $Log: variables_3d.F90,v $
! Revision 1.4 2003-04-23 12:16:34 kbk
! Revision 1.5 2003-12-16 15:58:54 kbk
! back ground viscosity and diffusivity (manuel)
!
! Revision 1.4 2003/04/23 12:16:34 kbk
! cleaned code + TABS to spaces
!
! Revision 1.3 2003/04/07 16:22:31 kbk
......
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