pkg/atm_compon_interf/atm_store_dynvars.F

C $Header: /u/gcmpack/MITgcm/pkg/atm_compon_interf/atm_store_dynvars.F,v 1.2 2015/11/12 01:01:57 jmc Exp $
C $Name:  $

#include "ATM_CPL_OPTIONS.h"

CBOP
C     !ROUTINE: ATM_STORE_DYNVARS
C     !INTERFACE:
      SUBROUTINE ATM_STORE_DYNVARS(
     I                     bi, bj,
     I                     myTime, myIter, myThid )

C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE ATM_STORE_DYNVARS
C     | o Routine for saving dynamical fields (e.g.,  atmos.
C     |   sea-level pressure ) for export to coupling layer.
C     *==========================================================*
C     | This version interfaces to the main model
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE

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

C     !INPUT/OUTPUT PARAMETERS:
C     bi, bj    :: Tile indices
C     myTime    :: Current time in simulation (s)
C     myIter    :: Current iteration number
C     myThid    :: My Thread Id. number
      INTEGER bi, bj
      _RL     myTime
      INTEGER myIter
      INTEGER myThid

C     !LOCAL VARIABLES:
C     i, j      :: Loop counters
      INTEGER i,j
c     _RL cplTimeFraction
CEOP

C     o Accumulate Sea-level pressure anomaly from dynamical core
C       to the coupling layer. Eta is in Pa
c      cplTimeFraction = 1. _d 0 / DFLOAT(cplSendFrq_iter)
c      atmSLPrTime(bi,bj) = atmSLPrTime(bi,bj) + cplTimeFraction
       atmSLPrTime(bi,bj) = 1. _d 0
       DO j=1,sNy
        DO i=1,sNx
c         atmSLPr(i,j,bi,bj) = atmSLPr(i,j,bi,bj)
c    &                 + etaN(i,j,bi,bj)*cplTimeFraction
          atmSLPr(i,j,bi,bj) = etaN(i,j,bi,bj)
        ENDDO
       ENDDO

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/atm_compon_interf/atm_store_dynvars.F,v 1.2 2015/11/12 01:01:57 jmc Exp $
2	C $Name: $
3
4	#include "ATM_CPL_OPTIONS.h"
5
6	CBOP
7	C !ROUTINE: ATM_STORE_DYNVARS
8	C !INTERFACE:
9	SUBROUTINE ATM_STORE_DYNVARS(
10	I bi, bj,
11	I myTime, myIter, myThid )
12
13	C !DESCRIPTION: \bv
14	C ==========================================================
15	C \| SUBROUTINE ATM_STORE_DYNVARS
16	C \| o Routine for saving dynamical fields (e.g., atmos.
17	C \| sea-level pressure ) for export to coupling layer.
18	C ==========================================================
19	C \| This version interfaces to the main model
20	C ==========================================================
21	C \ev
22
23	C !USES:
24	IMPLICIT NONE
25
26	C == Global variables ==
27	#include "SIZE.h"
28	#include "EEPARAMS.h"
29	#include "PARAMS.h"
30	#include "CPL_PARAMS.h"
31	C == Global variables (main model)
32	#include "DYNVARS.h"
33	C == Global variables for coupling interface ==
34	#include "ATMCPL.h"
35
36	C !INPUT/OUTPUT PARAMETERS:
37	C bi, bj :: Tile indices
38	C myTime :: Current time in simulation (s)
39	C myIter :: Current iteration number
40	C myThid :: My Thread Id. number
41	INTEGER bi, bj
42	_RL myTime
43	INTEGER myIter
44	INTEGER myThid
45
46	C !LOCAL VARIABLES:
47	C i, j :: Loop counters
48	INTEGER i,j
49	c _RL cplTimeFraction
50	CEOP
51
52	C o Accumulate Sea-level pressure anomaly from dynamical core
53	C to the coupling layer. Eta is in Pa
54	c cplTimeFraction = 1. _d 0 / DFLOAT(cplSendFrq_iter)
55	c atmSLPrTime(bi,bj) = atmSLPrTime(bi,bj) + cplTimeFraction
56	atmSLPrTime(bi,bj) = 1. _d 0
57	DO j=1,sNy
58	DO i=1,sNx
59	c atmSLPr(i,j,bi,bj) = atmSLPr(i,j,bi,bj)
60	c & + etaN(i,j,bi,bj)*cplTimeFraction
61	atmSLPr(i,j,bi,bj) = etaN(i,j,bi,bj)
62	ENDDO
63	ENDDO
64
65	RETURN
66	END