/[MITgcm]/MITgcm/pkg/dic/alk_surfforcing.F
ViewVC logotype

Diff of /MITgcm/pkg/dic/alk_surfforcing.F

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

revision 1.4 by jmc, Sun Jul 18 01:13:50 2004 UTC revision 1.5 by stephd, Fri Aug 6 19:50:52 2004 UTC
# Line 29  C  myIter               :: current times Line 29  C  myIter               :: current times
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)
# Line 57  C Positive EminusPforV => loss of water Line 56  C Positive EminusPforV => loss of water
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
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22