jscott/pkg_atm2d/atm2ocn_main.F

#include "ctrparam.h"
#include "ATM2D_OPTIONS.h"
C     !INTERFACE:
      SUBROUTINE ATM2OCN_MAIN( iloop, aloop, inMonth,
     &                                 myIter, myThid )
C     *==========================================================*
C     | Main call after atmos call, calculate ocean/seaice       |
c     | fluxes                                                   |
C     *==========================================================*
        IMPLICIT NONE

C     === Global Atmos/Ocean/Seaice Interface Variables ===
#include "ATMSIZE.h"
#include "SIZE.h"
#include "EEPARAMS.h"
#include "ATM2D_VARS.h"


C     !INPUT/OUTPUT PARAMETERS:
C     === Routine arguments ===
C     inMonth - current month
C     myThid - Thread no. that called this routine.
      INTEGER iloop
      INTEGER aloop
      INTEGER inMonth
      INTEGER myIter
      INTEGER myThid

C     LOCAL VARIABLES:
      _RL wght0,wght1
      INTEGER intime0,intime1
      LOGICAL iftime
      _RL  curTime   ! current time (s) of simulation starting from 0.0

      curTime = float(iloop-1)*dtcouplo + (float(aloop)-0.5)*dtatmo

      PRINT *,'Executing atm2ocn main',curTime
      CALL CALC_FILELOAD( curTime, iloop, aloop, wght0, wght1,
     &                    intime0, intime1, ifTime, myThid)

      CALL READ_ATMOS(inMonth,myThid)

      IF ( (atmosTauuFile.NE.' ').OR.
     &   (atmosTauvFile.NE.' ').OR.(atmosWindFile.NE.' ')) THEN     
        CALL SUBTRACT_MEANS( wght0, wght1,
     &                       intime0, intime1, myThid)
      ENDIF

      CALL CALC_1DTO2D(inMonth,myThid)

      CALL FIXED_FLUX_ADD(wght0,wght1,
     &              intime0,intime1,ifTime,myIter,myThid)

      CALL RELAX_ADD(wght0,wght1,
     &              intime0,intime1,ifTime,myIter,myThid)
 
      RETURN
      END
1	#include "ctrparam.h"
2	#include "ATM2D_OPTIONS.h"
3	C !INTERFACE:
4	SUBROUTINE ATM2OCN_MAIN( iloop, aloop, inMonth,
5	& myIter, myThid )
6	C ==========================================================
7	C \| Main call after atmos call, calculate ocean/seaice \|
8	c \| fluxes \|
9	C ==========================================================
10	IMPLICIT NONE
11
12	C === Global Atmos/Ocean/Seaice Interface Variables ===
13	#include "ATMSIZE.h"
14	#include "SIZE.h"
15	#include "EEPARAMS.h"
16	#include "ATM2D_VARS.h"
17
18
19	C !INPUT/OUTPUT PARAMETERS:
20	C === Routine arguments ===
21	C inMonth - current month
22	C myThid - Thread no. that called this routine.
23	INTEGER iloop
24	INTEGER aloop
25	INTEGER inMonth
26	INTEGER myIter
27	INTEGER myThid
28
29	C LOCAL VARIABLES:
30	_RL wght0,wght1
31	INTEGER intime0,intime1
32	LOGICAL iftime
33	_RL curTime ! current time (s) of simulation starting from 0.0
34
35	curTime = float(iloop-1)dtcouplo + (float(aloop)-0.5)dtatmo
36
37	PRINT *,'Executing atm2ocn main',curTime
38	CALL CALC_FILELOAD( curTime, iloop, aloop, wght0, wght1,
39	& intime0, intime1, ifTime, myThid)
40
41	CALL READ_ATMOS(inMonth,myThid)
42
43	IF ( (atmosTauuFile.NE.' ').OR.
44	& (atmosTauvFile.NE.' ').OR.(atmosWindFile.NE.' ')) THEN
45	CALL SUBTRACT_MEANS( wght0, wght1,
46	& intime0, intime1, myThid)
47	ENDIF
48
49	CALL CALC_1DTO2D(inMonth,myThid)
50
51	CALL FIXED_FLUX_ADD(wght0,wght1,
52	& intime0,intime1,ifTime,myIter,myThid)
53
54	CALL RELAX_ADD(wght0,wght1,
55	& intime0,intime1,ifTime,myIter,myThid)
56
57	RETURN
58	END