pkg/atm_compon_interf/atm_store_uvelground.F

C $Header: $
C $Name:   $  

#include "AIM_OPTIONS.h"

CStartOfInterface
      SUBROUTINE ATM_STORE_UVELGROUND( bi,bj,
     I                           myCurrentTime, 
     I                           myCurrentIter, 
     I                           myThid )
C     /==========================================================\
C     | SUBROUTINE ATM_STORE_UVELGROUND                          |
C     | o Routine for saving atmos. lower boundary wind speed    |
C     |   for export to coupling layer.                          |
C     |==========================================================|
C     | This version interfaces to the MITgcm AIMPHYS package.   |
C     \==========================================================/
      IMPLICIT NONE

C     == Global variables ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "CPL_PARAMS.h"
#include "DYNVARS.h"
C     == Global variables for coupling interface ==
#include "ATMCPL.h"

C     == Routine arguments ==
C     bi,bj         - Tile index
C     myThid        - Thread number for this instance of the routine
C     myCurrentIter - Current timestep number
C     myCurrentTime - Current model time
      INTEGER bi, bj
      _RL     myCurrentTime
      INTEGER myCurrentIter
      INTEGER myThid
CEndOfInterface

C     == Local variables ==
C     I,J,K - Loop counters
      INTEGER I,J,K
      _RL cplTimeFraction

C     o Accumulate zonal wind speed that will be exported to the 
C       coupling layer. Speed is in m/s with positive 
C       to the east. Wind speed is in COMMON /DYNVARS_R/,
C       DYNVARS.h variable uVel using a south-west-c-grid staggering.
      cplTimeFraction = 1. _d 0 / DFLOAT(cplSendFrq_iter)
c     DO bj=myByLo(myThid),myByHi(myThid)
c      DO bi=myBxLo(myThid),myBxHi(myThid)
        uVelGroundTime(bi,bj) = uVelGroundTime(bi,bj) + cplTimeFraction
        DO J=1,sNy
         DO I=1,sNx
          uVelGround(I,J,bi,bj) = 
     &    uVelGround(I,J,bi,bj) + uVel(I,J,1,bi,bj)*cplTimeFraction
         ENDDO
        ENDDO
c      ENDDO
c     ENDDO

      RETURN
      END
1	C $Header: $
2	C $Name: $
3
4	#include "AIM_OPTIONS.h"
5
6	CStartOfInterface
7	SUBROUTINE ATM_STORE_UVELGROUND( bi,bj,
8	I myCurrentTime,
9	I myCurrentIter,
10	I myThid )
11	C /==========================================================\
12	C \| SUBROUTINE ATM_STORE_UVELGROUND \|
13	C \| o Routine for saving atmos. lower boundary wind speed \|
14	C \| for export to coupling layer. \|
15	C \|==========================================================\|
16	C \| This version interfaces to the MITgcm AIMPHYS package. \|
17	C \==========================================================/
18	IMPLICIT NONE
19
20	C == Global variables ==
21	#include "SIZE.h"
22	#include "EEPARAMS.h"
23	#include "PARAMS.h"
24	#include "CPL_PARAMS.h"
25	#include "DYNVARS.h"
26	C == Global variables for coupling interface ==
27	#include "ATMCPL.h"
28
29	C == Routine arguments ==
30	C bi,bj - Tile index
31	C myThid - Thread number for this instance of the routine
32	C myCurrentIter - Current timestep number
33	C myCurrentTime - Current model time
34	INTEGER bi, bj
35	_RL myCurrentTime
36	INTEGER myCurrentIter
37	INTEGER myThid
38	CEndOfInterface
39
40	C == Local variables ==
41	C I,J,K - Loop counters
42	INTEGER I,J,K
43	_RL cplTimeFraction
44
45	C o Accumulate zonal wind speed that will be exported to the
46	C coupling layer. Speed is in m/s with positive
47	C to the east. Wind speed is in COMMON /DYNVARS_R/,
48	C DYNVARS.h variable uVel using a south-west-c-grid staggering.
49	cplTimeFraction = 1. _d 0 / DFLOAT(cplSendFrq_iter)
50	c DO bj=myByLo(myThid),myByHi(myThid)
51	c DO bi=myBxLo(myThid),myBxHi(myThid)
52	uVelGroundTime(bi,bj) = uVelGroundTime(bi,bj) + cplTimeFraction
53	DO J=1,sNy
54	DO I=1,sNx
55	uVelGround(I,J,bi,bj) =
56	& uVelGround(I,J,bi,bj) + uVel(I,J,1,bi,bj)*cplTimeFraction
57	ENDDO
58	ENDDO
59	c ENDDO
60	c ENDDO
61
62	RETURN
63	END