1 |
C $Header: /u/gcmpack/MITgcm/pkg/dic/dic_store_fluxco2.F,v 1.1 2009/01/05 15:27:32 dfer Exp $ |
2 |
C $Name: $ |
3 |
|
4 |
#include "DIC_OPTIONS.h" |
5 |
|
6 |
CStartOfInterface |
7 |
SUBROUTINE DIC_STORE_FLUXCO2( myTime, |
8 |
I myIter, |
9 |
I myThid ) |
10 |
C *==========================================================* |
11 |
C | SUBROUTINE DIC_STORE_FLUXCO2 | |
12 |
C | o Routine for controlling storage of air-sea CO2 flux to | |
13 |
C | to the coupling layer. | |
14 |
C |==========================================================* |
15 |
C | This version talks to the MIT Coupler. It uses the MIT | |
16 |
C | Coupler "checkpoint1" library calls. | |
17 |
C *==========================================================* |
18 |
IMPLICIT NONE |
19 |
|
20 |
C == Global variables == |
21 |
#include "SIZE.h" |
22 |
#include "EEPARAMS.h" |
23 |
#include "PARAMS.h" |
24 |
C == DIC variables |
25 |
#include "DIC_VARS.h" |
26 |
#ifdef COMPONENT_MODULE |
27 |
C == Global variables for coupling interface == |
28 |
#include "OCNCPL.h" |
29 |
#include "OCNIDS.h" |
30 |
#include "CPL_PARAMS.h" |
31 |
#endif |
32 |
|
33 |
C == Routine arguments == |
34 |
C myThid :: Thread number for this instance of the routine |
35 |
C myIter :: Current timestep number |
36 |
C myTime :: Current model time |
37 |
_RL myTime |
38 |
INTEGER myIter |
39 |
INTEGER myThid |
40 |
CEndOfInterface |
41 |
|
42 |
#ifdef COMPONENT_MODULE |
43 |
C == Local variables == |
44 |
C i,j :: Loop counters |
45 |
C bi,bj :: Tile index |
46 |
INTEGER i,j,bi,bj |
47 |
|
48 |
#ifdef USE_ATMOSCO2 |
49 |
IF ( ocnCpl_exchange_DIC ) THEN |
50 |
DO bj=myByLo(myThid),myByHi(myThid) |
51 |
DO bi=myBxLo(myThid),myBxHi(myThid) |
52 |
DO j=1,sNy |
53 |
DO i=1,sNx |
54 |
fluxCO2cpl(i,j,bi,bj) = fluxCO2(i,j,bi,bj) |
55 |
ENDDO |
56 |
ENDDO |
57 |
ENDDO |
58 |
ENDDO |
59 |
ENDIF |
60 |
#endif |
61 |
|
62 |
#endif /* COMPONENT_MODULE */ |
63 |
|
64 |
RETURN |
65 |
END |