jscott/pkg_atm2d/subtract_means.F

#include "ctrparam.h"
#include "ATM2D_OPTIONS.h"

C     !INTERFACE:
      SUBROUTINE SUBTRACT_MEANS( wght0, wght1,
     &                           intime0, intime1, myThid )
C     *==========================================================*
C     | Subtract out (1D) mean values of wind variables          |
c     | (anomaly coupling)                                       |
C     *==========================================================*
      IMPLICIT NONE

#include "ATMSIZE.h"
#include "SIZE.h"
#include "EEPARAMS.h"
#include "ATM2D_VARS.h"

C     !INPUT/OUTPUT PARAMETERS:
C     === Routine arguments ===
C     myThid - Thread no. that called this routine.
      _RL  wght0
      _RL  wght1
      INTEGER intime0
      INTEGER intime1
      INTEGER myThid

C     LOCAL VARIABLES:
      INTEGER j_atm

C     Interpolate and subtract 1D mean wind files
      
      DO j_atm=2,jm0-1
        atm_tauu(j_atm)= atm_tauu(j_atm) -
     &       (wght0*atau(j_atm,intime0) + wght1*atau(j_atm,intime1))
        atm_tauv(j_atm)= atm_tauv(j_atm) -
     &       (wght0*atav(j_atm,intime0) + wght1*atav(j_atm,intime1))
        atm_windspeed(j_atm)= atm_windspeed(j_atm) -
     &       (wght0*awind(j_atm,intime0) + wght1*awind(j_atm,intime1))
      ENDDO

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