Commit 44d6f609 authored by kbk's avatar kbk
Browse files

molecular diffusion for salt and temp (manuel)

parent 7762093e
!$Id: salinity.F90,v 1.6 2003-09-13 10:52:21 kbk Exp $
!$Id: salinity.F90,v 1.7 2003-12-16 16:00:46 kbk Exp $
#include "cppdefs.h"
!-----------------------------------------------------------------------
!BOP
......@@ -39,7 +39,10 @@
! Original author(s): Karsten Bolding & Hans Burchard
!
! $Log: salinity.F90,v $
! Revision 1.6 2003-09-13 10:52:21 kbk
! Revision 1.7 2003-12-16 16:00:46 kbk
! molecular diffusion for salt and temp (manuel)
!
! Revision 1.6 2003/09/13 10:52:21 kbk
! changed field_no to salt_field_no and temp_field_no
!
! Revision 1.5 2003/08/03 08:13:09 kbk
......@@ -268,6 +271,7 @@ STDERR 'salinity= ',iimin,iimax,i+ioff,iextr/2
#else
use domain, only: dx,dy,ard1
#endif
use parameters, only: avmols
IMPLICIT NONE
!
! !INPUT PARAMETERS:
......@@ -366,8 +370,8 @@ STDERR 'salinity= ',iimin,iimax,i+ioff,iextr/2
if (kmax.gt.1) then
! Auxilury terms, old and new time level,
do k=1,kmax-1
auxo(k)=2.*(1-cnpar)*dt*nuh(i,j,k)/(hn(i,j,k+1)+hn(i,j,k))
auxn(k)=2.* cnpar *dt*nuh(i,j,k)/(hn(i,j,k+1)+hn(i,j,k))
auxo(k)=2.*(1-cnpar)*dt*(nuh(i,j,k)+avmols)/(hn(i,j,k+1)+hn(i,j,k))
auxn(k)=2.* cnpar *dt*(nuh(i,j,k)+avmols)/(hn(i,j,k+1)+hn(i,j,k))
end do
! Matrix elements for surface layer
......
!$Id: temperature.F90,v 1.6 2003-09-13 10:52:21 kbk Exp $
!$Id: temperature.F90,v 1.7 2003-12-16 16:00:46 kbk Exp $
#include "cppdefs.h"
!-----------------------------------------------------------------------
!BOP
......@@ -35,7 +35,10 @@
! Original author(s): Karsten Bolding & Hans Burchard
!
! $Log: temperature.F90,v $
! Revision 1.6 2003-09-13 10:52:21 kbk
! Revision 1.7 2003-12-16 16:00:46 kbk
! molecular diffusion for salt and temp (manuel)
!
! Revision 1.6 2003/09/13 10:52:21 kbk
! changed field_no to salt_field_no and temp_field_no
!
! Revision 1.5 2003/08/03 08:13:09 kbk
......@@ -211,6 +214,7 @@ temp_field_no=1
#else
use domain, only: dx,dy,ard1
#endif
use parameters, only: avmolt
IMPLICIT NONE
!
! !INPUT PARAMETERS:
......@@ -288,8 +292,8 @@ temp_field_no=1
if (kmax.gt.1) then
! Auxilury terms, old and new time level,
do k=1,kmax-1
auxo(k)=2.*(1-cnpar)*dt*nuh(i,j,k)/(hn(i,j,k+1)+hn(i,j,k))
auxn(k)=2.* cnpar *dt*nuh(i,j,k)/(hn(i,j,k+1)+hn(i,j,k))
auxo(k)=2.*(1-cnpar)*dt*(nuh(i,j,k)+avmolt)/(hn(i,j,k+1)+hn(i,j,k))
auxn(k)=2.* cnpar *dt*(nuh(i,j,k)+avmolt)/(hn(i,j,k+1)+hn(i,j,k))
end do
! Matrix elements for surface layer
......
!$Id: parameters.F90,v 1.2 2003-04-23 12:02:43 kbk Exp $
!$Id: parameters.F90,v 1.3 2003-12-16 16:00:46 kbk Exp $
#include "cppdefs.h"
!-----------------------------------------------------------------------
!BOP
......@@ -24,6 +24,8 @@
REALTYPE, parameter :: cp = 3985.
REALTYPE, parameter :: kappa = 0.4
REALTYPE, parameter :: avmmol = 1.8e-6
REALTYPE, parameter :: avmolt= 1.4e-7
REALTYPE, parameter :: avmols = 1.1e-9
!
! Turbulence related constants - see www.gotm.net
!
......@@ -31,7 +33,10 @@
! Original author(s): Karsten Bolding & Hans Burchard
!
! $Log: parameters.F90,v $
! Revision 1.2 2003-04-23 12:02:43 kbk
! Revision 1.3 2003-12-16 16:00:46 kbk
! molecular diffusion for salt and temp (manuel)
!
! Revision 1.2 2003/04/23 12:02:43 kbk
! cleaned code + TABS to spaces
!
! Revision 1.1.1.1 2002/05/02 14:01:19 gotm
......
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