3 |
|
|
4 |
#include "DIC_OPTIONS.h" |
#include "DIC_OPTIONS.h" |
5 |
#include "PTRACERS_OPTIONS.h" |
#include "PTRACERS_OPTIONS.h" |
|
#include "GCHEM_OPTIONS.h" |
|
6 |
|
|
7 |
CBOP |
CBOP |
8 |
C !ROUTINE: DIC_SURFFORCING |
C !ROUTINE: DIC_SURFFORCING |
24 |
#include "PARAMS.h" |
#include "PARAMS.h" |
25 |
#include "GRID.h" |
#include "GRID.h" |
26 |
#include "FFIELDS.h" |
#include "FFIELDS.h" |
27 |
#include "DIC_ABIOTIC.h" |
#include "DIC_VARS.h" |
28 |
|
|
29 |
C !INPUT PARAMETERS: =================================================== |
C !INPUT PARAMETERS: =================================================== |
30 |
C myThid :: thread number |
C myThid :: thread number |
100 |
|
|
101 |
CALL CARBON_COEFFS( |
CALL CARBON_COEFFS( |
102 |
I theta,salt, |
I theta,salt, |
103 |
I bi,bj,iMin,iMax,jMin,jMax) |
I bi,bj,iMin,iMax,jMin,jMax,myThid) |
104 |
C==================================================================== |
C==================================================================== |
105 |
|
|
106 |
DO j=jmin,jmax |
DO j=jmin,jmax |
141 |
I aks(i,j,bi,bj),akb(i,j,bi,bj),akw(i,j,bi,bj), |
I aks(i,j,bi,bj),akb(i,j,bi,bj),akw(i,j,bi,bj), |
142 |
I aksi(i,j,bi,bj),akf(i,j,bi,bj),ff(i,j,bi,bj), |
I aksi(i,j,bi,bj),akf(i,j,bi,bj),ff(i,j,bi,bj), |
143 |
I bt(i,j,bi,bj),st(i,j,bi,bj),ft(i,j,bi,bj), |
I bt(i,j,bi,bj),st(i,j,bi,bj),ft(i,j,bi,bj), |
144 |
U pH(i,j,bi,bj),pCO2(i,j,bi,bj) ) |
U pH(i,j,bi,bj),pCO2(i,j,bi,bj), |
145 |
|
I myThid ) |
146 |
ELSE |
ELSE |
147 |
pCO2(i,j,bi,bj)=0. _d 0 |
pCO2(i,j,bi,bj)=0. _d 0 |
148 |
ENDIF |
ENDIF |