/[MITgcm]/MITgcm_contrib/ocean_inversion_project/code/ptracers_init.F
ViewVC logotype

Diff of /MITgcm_contrib/ocean_inversion_project/code/ptracers_init.F

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

revision 1.7 by dimitri, Wed Oct 22 16:02:04 2003 UTC revision 1.9 by dimitri, Wed Oct 6 20:26:54 2004 UTC
# Line 18  C !USES: =============================== Line 18  C !USES: ===============================
18        IMPLICIT NONE        IMPLICIT NONE
19  #include "SIZE.h"  #include "SIZE.h"
20  #include "EEPARAMS.h"  #include "EEPARAMS.h"
21    #include "EESUPPORT.h"
22  #include "PARAMS.h"  #include "PARAMS.h"
23  #include "GRID.h"  #include "GRID.h"
24  #include "PTRACERS.h"  #include "PTRACERS.h"
# Line 102  C     Initialize pTracerTakahashi Line 103  C     Initialize pTracerTakahashi
103    
104  C     Normalize pTracerTakahashi so that 1e18 mol/yr is released  C     Normalize pTracerTakahashi so that 1e18 mol/yr is released
105  C     from each model region defined in pTracerMasks.  C     from each model region defined in pTracerMasks.
106  C     It is assumed that each year is 365.25 days (31557600 s)  C     It is assumed that each year is 365.24167 days (31556880 s)
107  C     long and that each month is 2629800 s.  C     long and that each month is 2629740 s.
108    
109        DO iTracer = 1, PTRACERS_num        DO iTracer = 1, PTRACERS_num
110         SumPtracer = 0.         SumPtracer = 0.
# Line 112  C     long and that each month is 262980 Line 113  C     long and that each month is 262980
113           DO j=1,sNy           DO j=1,sNy
114            DO i=1,sNx            DO i=1,sNx
115             DO iMonth=1,12             DO iMonth=1,12
116              SumPtracer = SumPtracer + 2629800. * rA(I,J,bi,bj) *              SumPtracer = SumPtracer + 2629740. * rA(I,J,bi,bj) *
117       &             pTracerMasks    (I,J,iTracer,bi,bj) *       &             pTracerMasks    (I,J,iTracer,bi,bj) *
118       &             pTracerTakahashi(I,J,iMonth ,bi,bj)       &             pTracerTakahashi(I,J,iMonth ,bi,bj)
119             ENDDO             ENDDO
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22