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	jscott	1.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	jscott	1.2	C \| Subtract out (1D) mean values of wind variables \|
9			c \| (anomaly coupling) \|
10	jscott	1.1	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