29 |
C myTime :: current time |
C myTime :: current time |
30 |
C bi,bj :: tile indices |
C bi,bj :: tile indices |
31 |
C PTR_ALK :: alkalinity field |
C PTR_ALK :: alkalinity field |
|
C GALK :: surface flux of alkalinity |
|
32 |
INTEGER myIter, myThid |
INTEGER myIter, myThid |
33 |
_RL myTime |
_RL myTime |
34 |
_RL PTR_ALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr) |
_RL PTR_ALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr) |
56 |
C in salinity. Thus, also increase in other surface tracers |
C in salinity. Thus, also increase in other surface tracers |
57 |
C (i.e. positive virtual flux into surface layer) |
C (i.e. positive virtual flux into surface layer) |
58 |
cQQ GALK(i,j)=gsm_ALK*ftotflux(i,j,bi,bj)/gsm_S |
cQQ GALK(i,j)=gsm_ALK*ftotflux(i,j,bi,bj)/gsm_S |
59 |
GALK(i,j)=gsm_ALK*surfaceForcingS(i,j,bi,bj)/gsm_S |
GALK(i,j)=gsm_ALK*surfaceForcingS(i,j,bi,bj)* |
60 |
|
& recip_drF(1)/gsm_S |
61 |
c OR |
c OR |
62 |
c let virtual flux be zero |
c let virtual flux be zero |
63 |
c GALK(i,j)=0.d0 |
c GALK(i,j)=0.d0 |