/[MITgcm]/MITgcm/pkg/ptracers/ptracers_forcing.F
ViewVC logotype

Diff of /MITgcm/pkg/ptracers/ptracers_forcing.F

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

revision 1.5 by heimbach, Sun Oct 26 01:10:34 2003 UTC revision 1.6 by jmc, Sun Jul 18 01:12:04 2004 UTC
# Line 20  C !USES: =============================== Line 20  C !USES: ===============================
20  #include "SIZE.h"  #include "SIZE.h"
21  #include "EEPARAMS.h"  #include "EEPARAMS.h"
22  #include "PARAMS.h"  #include "PARAMS.h"
23  #include "FFIELDS.h"  c #include "FFIELDS.h"
24  #include "DYNVARS.h"  c #include "DYNVARS.h"
25  #include "GRID.h"  #include "GRID.h"
26    
27  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
# Line 57  CEOP Line 57  CEOP
57         kSurface = 1         kSurface = 1
58        endif        endif
59    
60  C--   Surface forcing term is surfaceTendencyPtr, previously computed  C--   Surface forcing term is surfaceForcingPtr, previously computed
61  C     by PTRACERS_FORCING because it is needed by KPP_TRANSPORT_PTR.  C     by PTRACERS_FORCING because it is needed by KPP_TRANSPORT_PTR.
62    
63         IF ( k .EQ. kSurface ) THEN         IF ( k .EQ. kSurface ) THEN
64          DO j=jMin,jMax          DO j=jMin,jMax
65           DO i=iMin,iMax           DO i=iMin,iMax
66            gPtracer(i,j,k,bi,bj) = gPtracer(i,j,k,bi,bj) +            gPtracer(i,j,k,bi,bj) = gPtracer(i,j,k,bi,bj)
67       &           maskC(i,j,k,bi,bj) *       &          + surfPtracer(i,j,bi,bj)
68       &           surfPtracer(i,j,bi,bj)       &           *recip_drF(k)*recip_hFacC(i,j,k,bi,bj)
69           ENDDO           ENDDO
70          ENDDO          ENDDO
71         ENDIF         ENDIF
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22