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; calls routine to convert/get |
C     | the atmos fluxes, do ther derivative term calculation,   |
C     | then add restoring and fixed fluxes, if any.             |
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     iloop - loop counter for main loop (coupled periods)
C     aloop - loop counter for atm time steps (within a coupled per.)
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 ! comes back true if time to do a file re-load
      _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( 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; calls routine to convert/get \|
8	C \| the atmos fluxes, do ther derivative term calculation, \|
9	C \| then add restoring and fixed fluxes, if any. \|
10	C ==========================================================
11	IMPLICIT NONE
12
13	C === Global Atmos/Ocean/Seaice Interface Variables ===
14	#include "ATMSIZE.h"
15	#include "SIZE.h"
16	#include "EEPARAMS.h"
17	#include "ATM2D_VARS.h"
18
19
20	C !INPUT/OUTPUT PARAMETERS:
21	C === Routine arguments ===
22	C iloop - loop counter for main loop (coupled periods)
23	C aloop - loop counter for atm time steps (within a coupled per.)
24	C inMonth - current month
25	C myThid - Thread no. that called this routine.
26	INTEGER iloop
27	INTEGER aloop
28	INTEGER inMonth
29	INTEGER myIter
30	INTEGER myThid
31
32	C LOCAL VARIABLES:
33	_RL wght0,wght1
34	INTEGER intime0,intime1
35	LOGICAL iftime ! comes back true if time to do a file re-load
36	_RL curTime ! current time (s) of simulation starting from 0.0
37
38	curTime = float(iloop-1)dtcouplo + (float(aloop)-0.5)dtatmo
39
40	PRINT *,'Executing atm2ocn main',curTime
41	CALL CALC_FILELOAD( curTime, iloop, aloop, wght0, wght1,
42	& intime0, intime1, ifTime, myThid)
43
44	CALL READ_ATMOS( inMonth,myThid)
45
46	IF ( (atmosTauuFile.NE.' ').OR.
47	& (atmosTauvFile.NE.' ').OR.(atmosWindFile.NE.' ')) THEN
48	CALL SUBTRACT_MEANS( wght0, wght1,
49	& intime0, intime1, myThid)
50	ENDIF
51
52	CALL CALC_1DTO2D( myThid)
53
54	CALL FIXED_FLUX_ADD( wght0,wght1,
55	& intime0,intime1,ifTime,myIter,myThid)
56
57	CALL RELAX_ADD( wght0,wght1,
58	& intime0,intime1,ifTime,myIter,myThid)
59
60	RETURN
61	END