1 |
jmc |
1.2 |
C $Header: $ |
2 |
|
|
C $Name: $ |
3 |
|
|
|
4 |
jscott |
1.1 |
#include "ctrparam.h" |
5 |
|
|
#include "ATM2D_OPTIONS.h" |
6 |
|
|
|
7 |
|
|
C !INTERFACE: |
8 |
|
|
SUBROUTINE SUBTRACT_MEANS( wght0, wght1, |
9 |
|
|
& intime0, intime1, myThid ) |
10 |
|
|
C *==========================================================* |
11 |
|
|
C | Subtract out (1D) mean values of wind variables | |
12 |
|
|
c | (anomaly coupling) | |
13 |
|
|
C *==========================================================* |
14 |
|
|
IMPLICIT NONE |
15 |
|
|
|
16 |
|
|
#include "ATMSIZE.h" |
17 |
|
|
#include "SIZE.h" |
18 |
|
|
#include "EEPARAMS.h" |
19 |
|
|
#include "ATM2D_VARS.h" |
20 |
|
|
|
21 |
|
|
C !INPUT/OUTPUT PARAMETERS: |
22 |
|
|
C === Routine arguments === |
23 |
|
|
C myThid - Thread no. that called this routine. |
24 |
|
|
_RL wght0 |
25 |
|
|
_RL wght1 |
26 |
|
|
INTEGER intime0 |
27 |
|
|
INTEGER intime1 |
28 |
|
|
INTEGER myThid |
29 |
|
|
|
30 |
|
|
C LOCAL VARIABLES: |
31 |
|
|
INTEGER j_atm |
32 |
|
|
|
33 |
|
|
C Interpolate and subtract 1D mean wind files |
34 |
jmc |
1.2 |
|
35 |
jscott |
1.1 |
DO j_atm=2,jm0-1 |
36 |
|
|
atm_tauu(j_atm)= atm_tauu(j_atm) - |
37 |
|
|
& (wght0*atau(j_atm,intime0) + wght1*atau(j_atm,intime1)) |
38 |
|
|
atm_tauv(j_atm)= atm_tauv(j_atm) - |
39 |
|
|
& (wght0*atav(j_atm,intime0) + wght1*atav(j_atm,intime1)) |
40 |
|
|
atm_windspeed(j_atm)= atm_windspeed(j_atm) - |
41 |
|
|
& (wght0*awind(j_atm,intime0) + wght1*awind(j_atm,intime1)) |
42 |
|
|
ENDDO |
43 |
|
|
|
44 |
|
|
RETURN |
45 |
|
|
END |