/[MITgcm]/MITgcm/pkg/bulk_force/bulkf_ave.F
ViewVC logotype

Diff of /MITgcm/pkg/bulk_force/bulkf_ave.F

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

revision 1.7 by jmc, Fri Nov 4 01:28:38 2005 UTC revision 1.8 by jmc, Tue May 9 18:59:48 2006 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "BULK_FORCE_OPTIONS.h"  #include "BULK_FORCE_OPTIONS.h"
5    
6    CBOP
7  C     !ROUTINE: BULKF_AVE  C     !ROUTINE: BULKF_AVE
8  C     !INTERFACE:  C     !INTERFACE:
9        SUBROUTINE BULKF_AVE(bi,bj,mythid)        SUBROUTINE BULKF_AVE(bi,bj,myThid)
10    
11    C     !DESCRIPTION: \bv
12  C     *==========================================================*  C     *==========================================================*
13  C     | S/R  BULKF_AVE  C     | S/R  BULKF_AVE
14  C     | o save values for timeaveraging  C     | o save values for timeaveraging
15  C     |   for sea/ice surface and atmosphere  C     |   for sea/ice surface and atmosphere
16  C     *==========================================================*  C     *==========================================================*
17    C     \ev
18    
19  C     !USES:  C     !USES:
20        IMPLICIT NONE        IMPLICIT NONE
# Line 21  C     == Global variables == Line 25  C     == Global variables ==
25  #include "FFIELDS.h"  #include "FFIELDS.h"
26  #include "BULKF.h"  #include "BULKF.h"
27  #include "BULKF_INT.h"  #include "BULKF_INT.h"
28  #include "BULKF_DIAG.h"  #include "BULKF_TAVE.h"
 c     == routine arguments ==  
29    
30        integer mythid  C     !INPUT/OUTPUT PARAMETERS:
31        integer bi,bj        INTEGER bi,bj
32          INTEGER myThid
33    CEOP
34    
35  #ifdef ALLOW_BULK_FORCE  #ifdef ALLOW_BULK_FORCE
36  #ifdef ALLOW_TIMEAVE  #ifdef ALLOW_TIMEAVE
37    
38        integer i,j        INTEGER i,j
39    
40  C--     Time-average  C--     Time-average
41            do j = 1,sny        DO j = 1,sNy
42              do i = 1,snx          DO i = 1,sNx
43               BULK_Qnet_AVE(i,j,bi,bj)=BULK_Qnet_AVE(i,j,bi,bj)               BULK_Qnet_AVE(i,j,bi,bj)=BULK_Qnet_AVE(i,j,bi,bj)
44       &                       +Qnet(i,j,bi,bj)*deltaTclock       &                       +Qnet(i,j,bi,bj)*deltaTclock
45  c    &          +lambdaThetaClimRelax*  c    &          +lambdaThetaClimRelax*
# Line 60  c    &           (recip_Cp*recip_rhoNil* Line 65  c    &           (recip_Cp*recip_rhoNil*
65       &                       +fswnet(i,j,bi,bj)*deltaTclock       &                       +fswnet(i,j,bi,bj)*deltaTclock
66               BULK_ssq_AVE(i,j,bi,bj)=BULK_ssq_AVE(i,j,bi,bj)               BULK_ssq_AVE(i,j,bi,bj)=BULK_ssq_AVE(i,j,bi,bj)
67       &                       +savssq(i,j,bi,bj)*deltaTclock       &                       +savssq(i,j,bi,bj)*deltaTclock
68              enddo          ENDDO
69             enddo        ENDDO
70             BULKF_timeave(1,bi,bj)=BULKF_timeave(1,bi,bj)+deltaTclock        BULKF_timeave(1,bi,bj)=BULKF_timeave(1,bi,bj)+deltaTclock
71    
72    
73  #endif  /*ALLOW_TIMEAVE*/  #endif  /*ALLOW_TIMEAVE*/
74  #endif  /*ALLOW_BULK_FORCE*/  #endif  /*ALLOW_BULK_FORCE*/
75    
76        RETURN        RETURN
77        END        END
78    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22