/[MITgcm]/MITgcm/pkg/fizhi/fizhi_rayleigh.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/fizhi_rayleigh.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by molod, Fri Oct 22 14:52:14 2004 UTC revision 1.3 by molod, Fri Oct 22 16:33:39 2004 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3                                                                                                  
4  #include "FIZHI_OPTIONS.h"  #include "FIZHI_OPTIONS.h"
5        subroutine rayleigh(myid,pres,pkap,pekap,zsurf,u,v,t,s,im,jm,lm,        subroutine rayleigh(myid,pres,pkap,pekap,zsurf,u,v,t,s,im,jm,lm,
6       .                                                bi,bj,rfu,rfv,rft)       .                                                bi,bj,rfu,rfv,rft)
7  C **********************************************************************  C **********************************************************************
8  C  C
9  C  PURPOSE  C  PURPOSE
10  C     To implement Rayleigh Friction  C     Rayleigh Friction -- Linear Drag (Strong Damping) Above 70 Km
11  C  C
12  C  ARGUMENTS   DESCRIPTION  C  ARGUMENTS   DESCRIPTION
13  C  C
# Line 50  C ************************************** Line 50  C **************************************
50        _RL rf(im,jm,lm)        _RL rf(im,jm,lm)
51        _RL z(im,jm,lm)        _RL z(im,jm,lm)
52        _RL dz(im,jm,lm)        _RL dz(im,jm,lm)
53        _RL cpog, cpinv, virtcon, getcon        _RL cpog, cpinv, virtcon, getcon, dampcoef
54    
55  C **********************************************************************  C **********************************************************************
56  C ****   APPLY RAYLEIGH FRICTION TO WIND (INCLUDE HEATING)           ***  C ****   APPLY RAYLEIGH FRICTION TO WIND (INCLUDE HEATING)           ***
# Line 59  C ************************************** Line 59  C **************************************
59        cpog = getcon('CP')/getcon('GRAVITY')        cpog = getcon('CP')/getcon('GRAVITY')
60        cpinv = 1.0/getcon('CP')        cpinv = 1.0/getcon('CP')
61        virtcon = getcon('VIRTCON')        virtcon = getcon('VIRTCON')
62          dampcoef = 2./3.
63    
64        do L=1,lm        do L=1,lm
65         do j=1,jm         do j=1,jm
# Line 95  C ************************************** Line 96  C **************************************
96        do L=1,lm        do L=1,lm
97         do j=1,jm         do j=1,jm
98         do i=1,im         do i=1,im
99          rf(i,j,L) = (2./3.)*(1+tanh((z(i,j,L)-80000.)/5000.))/86400.          rf(i,j,L) = dampcoef*(1+tanh((z(i,j,L)-80000.)/5000.))/86400.
100          rfu(i,j,L) = - rf(i,j,L) * u(i,j,L)          rfu(i,j,L) = - rf(i,j,L) * u(i,j,L)
101          rfv(i,j,L) = - rf(i,j,L) * v(i,j,L)          rfv(i,j,L) = - rf(i,j,L) * v(i,j,L)
102          rft(i,j,L) = -(u(i,j,L)*rfu(i,j,L) + v(i,j,L)*rfv(i,j,L) )*cpinv          rft(i,j,L) = -(u(i,j,L)*rfu(i,j,L) + v(i,j,L)*rfv(i,j,L) )*cpinv
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22