/[MITgcm]/MITgcm/pkg/aim/aim_calc_diags.F
ViewVC logotype

Diff of /MITgcm/pkg/aim/aim_calc_diags.F

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

revision 1.6 by jmc, Tue Sep 25 19:53:57 2001 UTC revision 1.7 by jmc, Fri Sep 27 20:05:11 2002 UTC
# Line 9  C     /================================= Line 9  C     /=================================
9  C     | SUBROUTINE AIM_CALC_DIAGS                                |  C     | SUBROUTINE AIM_CALC_DIAGS                                |
10  C     | o Calculate AIM diagnostics                              |  C     | o Calculate AIM diagnostics                              |
11  C     \==========================================================/  C     \==========================================================/
12        IMPLICIT rEAL*8 (A-H,O-Z)        IMPLICIT NONE
13    
14  C     === Global variables ===  C     === Global variables ===
15  #include "SIZE.h"  C-- size for MITgcm & Physics package :
16    #include "AIM_SIZE.h"
17    
18  #include "EEPARAMS.h"  #include "EEPARAMS.h"
19  #include "PARAMS.h"  #include "PARAMS.h"
 #include "AIM_DIAGS.h"  
20    
21  #ifdef ALLOW_AIM  #include "AIM_DIAGS.h"
 #include "atparam0.h"  
 #include "atparam1.h"  
       INTEGER NGP  
       INTEGER NLON  
       INTEGER NLAT  
       INTEGER NLEV  
       PARAMETER ( NLON=IX, NLAT=IL, NLEV=KX, NGP=NLON*NLAT )  
22  #include "com_physvar.h"  #include "com_physvar.h"
 #endif /* ALLOW_AIM */  
23    
24        LOGICAL  DIFFERENT_MULTIPLE        LOGICAL  DIFFERENT_MULTIPLE
25        EXTERNAL DIFFERENT_MULTIPLE        EXTERNAL DIFFERENT_MULTIPLE
# Line 53  C     == Local variables == Line 46  C     == Local variables ==
46        DO J=1,sNy        DO J=1,sNy
47         DO I=1,sNx         DO I=1,sNx
48          I2 = sNx*(J-1)+I          I2 = sNx*(J-1)+I
49  c       USTRtave(i,j,bi,bj)   = USTRtave(i,j,bi,bj)   + USTR(I2,3,myThid)*DDTT  c       USTRtave(i,j,bi,bj)   = USTRtave(i,j,bi,bj)
50  c       VSTRtave(i,j,bi,bj)   = VSTRtave(i,j,bi,bj)   + VSTR(I2,3,myThid)*DDTT  c    &                        + USTR(I2,3,myThid)*DDTT
51          TSRtave(i,j,bi,bj)    = TSRtave(i,j,bi,bj)    + TSR(I2,myThid)*DDTT  c       VSTRtave(i,j,bi,bj)   = VSTRtave(i,j,bi,bj)
52          OLRtave(i,j,bi,bj)    = OLRtave(i,j,bi,bj)    + OLR(I2,myThid)*DDTT  c    &                        + VSTR(I2,3,myThid)*DDTT
53          SSRtave(i,j,bi,bj)    = SSRtave(i,j,bi,bj)    + SSR(I2,myThid)*DDTT          TSRtave(i,j,bi,bj)    = TSRtave(i,j,bi,bj)
54          SLRtave(i,j,bi,bj)    = SLRtave(i,j,bi,bj)    + SLR(I2,myThid)*DDTT       &                        + TSR(I2,myThid)*DDTT
55          SHFtave(i,j,bi,bj)    = SHFtave(i,j,bi,bj)    + SHF(I2,3,myThid)*DDTT          OLRtave(i,j,bi,bj)    = OLRtave(i,j,bi,bj)
56          EVAPtave(i,j,bi,bj)   = EVAPtave(i,j,bi,bj)   + EVAP(I2,3,myThid)*DDTT       &                        + OLR(I2,myThid)*DDTT
57          PRECNVtave(i,j,bi,bj) = PRECNVtave(i,j,bi,bj) + PRECNV(I2,myThid)*DDTT          SSRtave(i,j,bi,bj)    = SSRtave(i,j,bi,bj)
58          PRECLStave(i,j,bi,bj) = PRECLStave(i,j,bi,bj) + PRECLS(I2,myThid)*DDTT       &                        + SSR(I2,myThid)*DDTT
59          CLOUDCtave(i,j,bi,bj) = CLOUDCtave(i,j,bi,bj) + CLOUDC(I2,myThid)*DDTT          SLRtave(i,j,bi,bj)    = SLRtave(i,j,bi,bj)
60         &                        + SLR(I2,myThid)*DDTT
61            SHFtave(i,j,bi,bj)    = SHFtave(i,j,bi,bj)
62         &                        + SHF(I2,3,myThid)*DDTT
63            EVAPtave(i,j,bi,bj)   = EVAPtave(i,j,bi,bj)
64         &                        + EVAP(I2,3,myThid)*DDTT
65            PRECNVtave(i,j,bi,bj) = PRECNVtave(i,j,bi,bj)
66         &                        + PRECNV(I2,myThid)*DDTT
67            PRECLStave(i,j,bi,bj) = PRECLStave(i,j,bi,bj)
68         &                        + PRECLS(I2,myThid)*DDTT
69            CLOUDCtave(i,j,bi,bj) = CLOUDCtave(i,j,bi,bj)
70         &                        + CLOUDC(I2,myThid)*DDTT
71         ENDDO         ENDDO
72        ENDDO        ENDDO
73    
74  C-    Keep record of how much time has been integrated over  C-    Keep record of how much time has been integrated over
75        DO K=1,Nr        DO K=1,Nr
76          AIM_TimeAve(k,bi,bj)=AIM_TimeAve(k,bi,bj)+DDTT          AIM_TimeAve(k,bi,bj)=AIM_TimeAve(k,bi,bj)+DDTT
77        ENDDO                                                                                            ENDDO
78    
79    
80  #endif /* ALLOW_TIMEAVE */  #endif /* ALLOW_TIMEAVE */
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22