1 |
C $Header: $ |
2 |
C $Name: $ |
3 |
|
4 |
#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 |
|
35 |
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 |