24 |
#include "OBCS.h" |
#include "OBCS.h" |
25 |
#include "ORLANSKI.h" |
#include "ORLANSKI.h" |
26 |
|
|
27 |
C SPK 6/2/00: Added radiative OBC's for salinity. |
C SPK 6/2/00: Added radiative OBCs for salinity. |
28 |
C SPK 6/6/00: Changed calculation of OB*w. When K=1, the |
C SPK 6/6/00: Changed calculation of OB*w. When K=1, the |
29 |
C upstream value is used. For example on the eastern OB: |
C upstream value is used. For example on the eastern OB: |
30 |
C IF (K.EQ.1) THEN |
C IF (K.EQ.1) THEN |
38 |
C 3) (Sonya Legg) Changed application of uVel and vVel. |
C 3) (Sonya Legg) Changed application of uVel and vVel. |
39 |
C uVel on the western OB is actually applied at I_obc+1 |
C uVel on the western OB is actually applied at I_obc+1 |
40 |
C while vVel on the southern OB is applied at J_obc+1. |
C while vVel on the southern OB is applied at J_obc+1. |
41 |
C 4) (Sonya Legg) Added templates for forced OB's. |
C 4) (Sonya Legg) Added templates for forced OBs. |
42 |
C |
C |
43 |
C SPK 7/17/00: Non-uniform resolution is now taken into account in diagnosing |
C SPK 7/17/00: Non-uniform resolution is now taken into account in diagnosing |
44 |
C phase speeds and time-stepping OB values. CL is still the |
C phase speeds and time-stepping OB values. CL is still the |
87 |
/* cvelTimeScale is averaging period for phase speed in sec. */ |
/* cvelTimeScale is averaging period for phase speed in sec. */ |
88 |
|
|
89 |
fracCVEL = deltaT/cvelTimeScale /* fraction of new phase speed used*/ |
fracCVEL = deltaT/cvelTimeScale /* fraction of new phase speed used*/ |
90 |
f1 = fracCVEL /* don't change this. Set cvelTimeScale */ |
f1 = fracCVEL /* dont change this. Set cvelTimeScale */ |
91 |
f2 = 1.0-fracCVEL /* don't change this. set cvelTimeScale */ |
f2 = 1.0-fracCVEL /* dont change this. set cvelTimeScale */ |
92 |
|
|
93 |
C Northern OB (Orlanski Radiation Condition) |
C Northern OB (Orlanski Radiation Condition) |
94 |
DO K=1,Nr |
DO K=1,Nr |