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

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

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

revision 1.2 by jmc, Tue Jul 13 18:03:31 2004 UTC revision 1.10 by jmc, Sat Jan 2 23:07:39 2010 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4  cswdcost -- add sunroutine ---  #include "DIC_OPTIONS.h"
 #include "CPP_OPTIONS.h"  
 #include "GCHEM_OPTIONS.h"  
   
5    
6  CStartOfInterFace  CStartOfInterFace
7        SUBROUTINE DIC_AVER(        SUBROUTINE DIC_AVER(
8       I           myTime,myIter,myThid)       I           myTime,myIter,myThid)
9    
10  C     /==========================================================\  C     *==========================================================*
11  C     | SUBROUTINE DIC_AVER  i                            |  C     | SUBROUTINE DIC_AVER
12  C     |==========================================================|  C     *==========================================================*
13        IMPLICIT NONE        IMPLICIT NONE
14    
15  C     == GLobal variables ==  C     == GLobal variables ==
# Line 22  C     == GLobal variables == Line 19  C     == GLobal variables ==
19  #include "PARAMS.h"  #include "PARAMS.h"
20  #include "GRID.h"  #include "GRID.h"
21  #include "PTRACERS_SIZE.h"  #include "PTRACERS_SIZE.h"
22  #include "PTRACERS.h"  #include "PTRACERS_FIELDS.h"
23  #include "GCHEM.h"  #include "DIC_VARS.h"
 #include "DIC_ABIOTIC.h"  
24  #ifdef DIC_BIOTIC  #ifdef DIC_BIOTIC
 #include "DIC_BIOTIC.h"  
25  #include "DIC_DIAGS.h"  #include "DIC_DIAGS.h"
26  #include "DIC_COST.h"  #include "DIC_COST.h"
27  #endif  #endif
 #ifdef ALLOW_SEAICE  
 #include "ICE.h"  
 #endif  
28    
29  C     == Routine arguments ==  C     == Routine arguments ==
       INTEGER myIter  
30        _RL myTime        _RL myTime
31          INTEGER myIter
32        INTEGER myThid        INTEGER myThid
33    
34  #ifdef ALLOW_DIC_COST  #ifdef ALLOW_DIC_COST
35    
36  C     == Local variables ==  C     == Local variables ==
       LOGICAL  DIFFERENT_MULTIPLE  
       EXTERNAL DIFFERENT_MULTIPLE  
37        INTEGER i, j, bi, bj, k, it        INTEGER i, j, bi, bj, k, it
38        _RL rdt        _RL rdt
39        INTEGER nForcingPeriods,Imytm,Ifprd,Ifcyc,Iftm        INTEGER nForcingPeriods,Imytm,Ifprd,Ifcyc,Iftm
40    
41  cswddmonth -- end-  Cswddmonth -- end-
42  c  
43  c averages  C averages
44        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
45         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
46            OBS_timetave(bi,bj) = OBS_timetave(bi,bj)+deltaTclock
47          DO k=1,Nr          DO k=1,Nr
          OBS_timetave(bi,bj,k)=OBS_timetave(bi,bj,k)+  
      &                         deltaTclock  
48           DO j=1-OLy,sNy+OLy           DO j=1-OLy,sNy+OLy
49            DO i=1-OLx,sNx+OLx            DO i=1-OLx,sNx+OLx
50              po4ann(i,j,k,bi,bj)=po4ann(i,j,k,bi,bj)+              po4ann(i,j,k,bi,bj)=po4ann(i,j,k,bi,bj)+
# Line 67  c averages Line 56  c averages
56          ENDDO          ENDDO
57         ENDDO         ENDDO
58        ENDDO        ENDDO
59  cswdmonth-add--  Cswdmonth-add--
60        rdt=1. _d 0 / deltaTclock        rdt=1. _d 0 / deltaTclock
61        nForcingPeriods=int(externForcingCycle/externForcingPeriod+0.5)        nForcingPeriods=int(externForcingCycle/externForcingPeriod+0.5)
62        Imytm=int(myTime*rdt+0.5)        Imytm=int(myTime*rdt+0.5)
# Line 78  cswdmonth-add-- Line 67  cswdmonth-add--
67  c     print*,'QQ timing check', mytime, myIter, it  c     print*,'QQ timing check', mytime, myIter, it
68        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
69         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
70           OBSM_timetave(bi,bj,it)=OBSM_timetave(bi,bj,it)+           OBSM_timetave(it,bi,bj) = OBSM_timetave(it,bi,bj)+deltaTclock
      &                         deltaTclock  
71           DO j=1-OLy,sNy+OLy           DO j=1-OLy,sNy+OLy
72            DO i=1-OLx,sNx+OLx            DO i=1-OLx,sNx+OLx
73              po4lev1(i,j,it,bi,bj)=po4lev1(i,j,it,bi,bj)+              po4lev1(i,j,it,bi,bj)=po4lev1(i,j,it,bi,bj)+
# Line 103  cQQ  &               PTRACER(i,j,4,bi,bj Line 91  cQQ  &               PTRACER(i,j,4,bi,bj
91         ENDDO         ENDDO
92        ENDDO        ENDDO
93    
94  cswdmonth-end--  Cswdmonth-end--
95    
96    #endif /* ALLOW_DIC_COST */
97    
 #endif  
 c  
98        RETURN        RETURN
99        END        END
 cswd -- end added subroutine --  
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22