/[MITgcm]/MITgcm/pkg/atm_compon_interf/atm_store_airco2.F
ViewVC logotype

Diff of /MITgcm/pkg/atm_compon_interf/atm_store_airco2.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by jmc, Mon Jun 18 21:28:52 2007 UTC revision 1.2 by dfer, Mon Jan 5 15:19:17 2009 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "PACKAGES_CONFIG.h"  #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6    #ifdef ALLOW_AIM
7    # include "AIM_OPTIONS.h"
8    #endif
9    
10  CBOP  CBOP
11  C     !ROUTINE: ATM_STORE_AIRCO2  C     !ROUTINE: ATM_STORE_AIRCO2
# Line 24  C     !USES: Line 27  C     !USES:
27        IMPLICIT NONE        IMPLICIT NONE
28    
29  C     == Global variables ==  C     == Global variables ==
30  #include "SIZE.h"  #ifdef ALLOW_AIM
31    # include "AIM_SIZE.h"
32    #else
33    # include "SIZE.h"
34    #endif
35    
36  #include "EEPARAMS.h"  #include "EEPARAMS.h"
37  #include "PARAMS.h"  #include "PARAMS.h"
38  #include "CPL_PARAMS.h"  #include "CPL_PARAMS.h"
39  c#ifdef ALLOW_AIM  #ifdef ALLOW_AIM
40  c#include "AIM_SIZE.h"  #include "AIM_CO2.h"
41  c#include "com_physvar.h"  #endif
 c#endif  
42  C     == Global variables for coupling interface ==  C     == Global variables for coupling interface ==
43  #include "ATMCPL.h"  #include "ATMCPL.h"
44    
# Line 49  C     myThid :: my Thread Id number Line 55  C     myThid :: my Thread Id number
55  CEOP  CEOP
56    
57  #ifdef ALLOW_AIM  #ifdef ALLOW_AIM
58    #ifdef ALLOW_AIM_CO2
59  C     == Local variables ==  C     == Local variables ==
60  C     i,j :: Loop counters  C     i,j :: Loop counters
61        INTEGER i,j, ij        INTEGER i,j, ij
62        _RL cplTimeFraction        _RL cplTimeFraction
63    
64  C     o Accumulate Sea-Ice Mass from thSIce pkg that will be exported  C     o Accumulate atmospheric CO2 from Aim pkg that will be exported
65  C       to the coupling layer. seaIceMass is per surface unit, in kg/m2.  C       to the coupling layer.
66         cplTimeFraction = 1. _d 0 / DFLOAT(cplSendFrq_iter)         cplTimeFraction = 1. _d 0 / DFLOAT(cplSendFrq_iter)
67         airCO2Time(bi,bj) = airCO2Time(bi,bj) + cplTimeFraction         airCO2Time(bi,bj) = airCO2Time(bi,bj) + cplTimeFraction
68         DO j=1,sNy         DO j=1,sNy
69           DO i=1,sNx           DO i=1,sNx
70            ij = i + (j-1)*sNx            ij = i + (j-1)*sNx
71            airCO2(i,j,bi,bj) = airCO2(i,j,bi,bj)            airCO2(i,j,bi,bj) = airCO2(i,j,bi,bj)
72  c    &                      + aim_CO2(ij,myThid)*cplTimeFraction       &                      + aim_CO2(ij,myThid)*cplTimeFraction
73           ENDDO           ENDDO
74         ENDDO         ENDDO
75    
76    #endif
77  #endif /* ALLOW_AIM */  #endif /* ALLOW_AIM */
78    
79        RETURN        RETURN
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22