42 |
INTEGER iMin,iMax,jMin,jMax, bi, bj |
INTEGER iMin,iMax,jMin,jMax, bi, bj |
43 |
|
|
44 |
C !OUTPUT PARAMETERS: =================================================== |
C !OUTPUT PARAMETERS: =================================================== |
45 |
c GDC :: tendency term due to air-sea exchange |
c GDC :: surface flux due to air-sea exchange |
46 |
_RL GDC(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL GDC(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
47 |
|
|
48 |
#ifdef ALLOW_PTRACERS |
#ifdef ALLOW_PTRACERS |
172 |
C in salinity. Thus, also increase in other surface tracers |
C in salinity. Thus, also increase in other surface tracers |
173 |
C (i.e. positive virtual flux into surface layer) |
C (i.e. positive virtual flux into surface layer) |
174 |
C ...so here, VirtualFLux = dC/dt! |
C ...so here, VirtualFLux = dC/dt! |
175 |
VirtualFlux(i,j)=gsm_DIC*surfaceTendencyS(i,j,bi,bj)/gsm_s |
VirtualFlux(i,j)=gsm_DIC*surfaceForcingS(i,j,bi,bj)/gsm_s |
176 |
c OR |
c OR |
177 |
c let virtual flux be zero |
c let virtual flux be zero |
178 |
c VirtualFlux(i,j)=0.d0 |
c VirtualFlux(i,j)=0.d0 |
187 |
DO j=1-OLy,sNy+OLy |
DO j=1-OLy,sNy+OLy |
188 |
DO i=1-OLx,sNx+OLx |
DO i=1-OLx,sNx+OLx |
189 |
GDC(i,j)= maskC(i,j,kLev,bi,bj)*( |
GDC(i,j)= maskC(i,j,kLev,bi,bj)*( |
190 |
& FluxCO2(i,j,bi,bj)*recip_drF(kLev) |
& FluxCO2(i,j,bi,bj) |
191 |
& + VirtualFlux(i,j) |
& + VirtualFlux(i,j) |
192 |
& ) |
& ) |
193 |
ENDDO |
ENDDO |