pkg/aim_compon_interf/atm_store_qshortwave.F

C $Header: $
C $Name:   $  

#include "AIM_OPTIONS.h"

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

C     == Global variables ==
#include "AIM_SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "CPL_PARAMS.h"
C     == Global variables for coupling interface ==
#include "ATMCPL.h"
C     == AIMPHYS specific global data ==
#include "com_physvar.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,II - Loop counters
      INTEGER I,J,K,II
      _RL cplTimeFraction

C     o Accumulate heat flux that will be exported to the coupling layer.
C     +ve into atmos., W/m^2.
      cplTimeFraction = 1. _d 0 / DFLOAT(cplSendFrq_iter)
c     DO bj=myByLo(myThid),myByHi(myThid)
c      DO bi=myBxLo(myThid),myBxHi(myThid)
        QshortwaveTime(bi,bj) = QshortwaveTime(bi,bj) + cplTimeFraction
        DO J=1,sNy
         DO I=1,sNx
          II = (sNx)*(J-1)+I
          Qshortwave(I,J,bi,bj) = 
     &    Qshortwave(I,J,bi,bj) - SSR(II,myThid)*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_QSHORTWAVE(bi,bj,
8	I myCurrentTime,
9	I myCurrentIter,
10	I myThid )
11	C /==========================================================\
12	C \| SUBROUTINE ATM_STORE_QSHORTWAVE \|
13	C \| o Routine for saving atmos. lower bdy shortwave heat flux\|
14	C \| 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 "AIM_SIZE.h"
22	#include "EEPARAMS.h"
23	#include "PARAMS.h"
24	#include "CPL_PARAMS.h"
25	C == Global variables for coupling interface ==
26	#include "ATMCPL.h"
27	C == AIMPHYS specific global data ==
28	#include "com_physvar.h"
29
30
31	C == Routine arguments ==
32	C bi,bj - Tile index
33	C myThid - Thread number for this instance of the routine
34	C myCurrentIter - Current timestep number
35	C myCurrentTime - Current model time
36	INTEGER bi, bj
37	_RL myCurrentTime
38	INTEGER myCurrentIter
39	INTEGER myThid
40	CEndOfInterface
41
42	C == Local variables ==
43	C I,J,K,II - Loop counters
44	INTEGER I,J,K,II
45	_RL cplTimeFraction
46
47	C o Accumulate heat flux that will be exported to the coupling layer.
48	C +ve into atmos., W/m^2.
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	QshortwaveTime(bi,bj) = QshortwaveTime(bi,bj) + cplTimeFraction
53	DO J=1,sNy
54	DO I=1,sNx
55	II = (sNx)*(J-1)+I
56	Qshortwave(I,J,bi,bj) =
57	& Qshortwave(I,J,bi,bj) - SSR(II,myThid)*cplTimeFraction
58	ENDDO
59	ENDDO
60	c ENDDO
61	c ENDDO
62
63	RETURN
64	END