pkg/atm2d/subtract_means.F

09a6f3668a Jeff*0001 #include "ctrparam.h"
                0002 #include "ATM2D_OPTIONS.h"
                0003
                0004 C     !INTERFACE:
                0005       SUBROUTINE SUBTRACT_MEANS( wght0, wght1,
                0006      &                           intime0, intime1, myThid )
                0007 C     *==========================================================*
                0008 C     | Subtract out (1D) mean values of wind variables          |
                0009 c     | (anomaly coupling)                                       |
                0010 C     *==========================================================*
                0011       IMPLICIT NONE
                0012
                0013 #include "ATMSIZE.h"
                0014 #include "SIZE.h"
                0015 #include "EEPARAMS.h"
                0016 #include "ATM2D_VARS.h"
                0017
                0018 C     !INPUT/OUTPUT PARAMETERS:
                0019 C     === Routine arguments ===
                0020 C     myThid - Thread no. that called this routine.
                0021       _RL  wght0
                0022       _RL  wght1
                0023       INTEGER intime0
                0024       INTEGER intime1
                0025       INTEGER myThid
                0026
                0027 C     LOCAL VARIABLES:
                0028       INTEGER j_atm
                0029
                0030 C     Interpolate and subtract 1D mean wind files
9274434acc Jean*0031
09a6f3668a Jeff*0032       DO j_atm=2,jm0-1
                0033         atm_tauu(j_atm)= atm_tauu(j_atm) -
                0034      &       (wght0*atau(j_atm,intime0) + wght1*atau(j_atm,intime1))
                0035         atm_tauv(j_atm)= atm_tauv(j_atm) -
                0036      &       (wght0*atav(j_atm,intime0) + wght1*atav(j_atm,intime1))
                0037         atm_windspeed(j_atm)= atm_windspeed(j_atm) -
                0038      &       (wght0*awind(j_atm,intime0) + wght1*awind(j_atm,intime1))
                0039       ENDDO
                0040
                0041       RETURN
                0042       END