pkg/bulk_force/bulkf_ave.F

6d54cf9ca1 Ed H*0001 #include "BULK_FORCE_OPTIONS.h"
7753507405 Curt*0002
dd80d278b6 Jean*0003 CBOP
7753507405 Curt*0004 C     !ROUTINE: BULKF_AVE
                0005 C     !INTERFACE:
dd80d278b6 Jean*0006       SUBROUTINE BULKF_AVE(bi,bj,myThid)
                0007
                0008 C     !DESCRIPTION: \bv
7753507405 Curt*0009 C     *==========================================================*
                0010 C     | S/R  BULKF_AVE
                0011 C     | o save values for timeaveraging
                0012 C     |   for sea/ice surface and atmosphere
                0013 C     *==========================================================*
dd80d278b6 Jean*0014 C     \ev
7753507405 Curt*0015
                0016 C     !USES:
                0017       IMPLICIT NONE
                0018 C     == Global variables ==
                0019 #include "SIZE.h"
                0020 #include "EEPARAMS.h"
                0021 #include "PARAMS.h"
                0022 #include "FFIELDS.h"
                0023 #include "BULKF.h"
f4245d1665 Curt*0024 #include "BULKF_INT.h"
dd80d278b6 Jean*0025 #include "BULKF_TAVE.h"
7753507405 Curt*0026
dd80d278b6 Jean*0027 C     !INPUT/OUTPUT PARAMETERS:
                0028       INTEGER bi,bj
                0029       INTEGER myThid
                0030 CEOP
7753507405 Curt*0031
                0032 #ifdef ALLOW_BULK_FORCE
                0033 #ifdef ALLOW_TIMEAVE
                0034
dd80d278b6 Jean*0035       INTEGER i,j
7753507405 Curt*0036
                0037 C--     Time-average
dd80d278b6 Jean*0038       DO j = 1,sNy
                0039         DO i = 1,sNx
6bf801d687 Jean*0040           BULK_Qnet_Ave(i,j,bi,bj)     = BULK_Qnet_Ave(i,j,bi,bj)
                0041      &                   +Qnet(i,j,bi,bj)*deltaTclock
7753507405 Curt*0042 c    &          +lambdaThetaClimRelax*
                0043 c    &           (theta(i,j,1,bi,bj)-SST(i,j,bi,bj))/
                0044 c    &           (recip_Cp*recip_rhoNil*recip_dRf(1))*deltaTclock
6bf801d687 Jean*0045           BULK_EmPmR_Ave(i,j,bi,bj)    = BULK_EmPmR_Ave(i,j,bi,bj)
                0046      &                   +EmPmR(i,j,bi,bj)*deltaTclock
                0047           BULK_fu_Ave(i,j,bi,bj)       = BULK_fu_Ave(i,j,bi,bj)
                0048      &                   +fu(i,j,bi,bj)*deltaTclock
                0049           BULK_fv_Ave(i,j,bi,bj)       = BULK_fv_Ave(i,j,bi,bj)
                0050      &                   +fv(i,j,bi,bj)*deltaTclock
                0051           BULK_latent_Ave(i,j,bi,bj)   = BULK_latent_Ave(i,j,bi,bj)
                0052      &                   +flh(i,j,bi,bj)*deltaTclock
                0053           BULK_sensible_Ave(i,j,bi,bj) = BULK_sensible_Ave(i,j,bi,bj)
                0054      &                   +fsh(i,j,bi,bj)*deltaTclock
                0055           BULK_evap_Ave(i,j,bi,bj)     = BULK_evap_Ave(i,j,bi,bj)
                0056      &                   +evap(i,j,bi,bj)*deltaTclock
                0057           BULK_flwup_Ave(i,j,bi,bj)    = BULK_flwup_Ave(i,j,bi,bj)
                0058      &                   +flwup(i,j,bi,bj)*deltaTclock
                0059           BULK_flwupnet_Ave(i,j,bi,bj) = BULK_flwupnet_Ave(i,j,bi,bj)
                0060      &                   +flwupnet(i,j,bi,bj)*deltaTclock
                0061           BULK_solar_Ave(i,j,bi,bj)    = BULK_solar_Ave(i,j,bi,bj)
                0062      &                   +fswnet(i,j,bi,bj)*deltaTclock
                0063           BULK_ssq_Ave(i,j,bi,bj)      = BULK_ssq_Ave(i,j,bi,bj)
                0064      &                   +savssq(i,j,bi,bj)*deltaTclock
dd80d278b6 Jean*0065         ENDDO
                0066       ENDDO
6bf801d687 Jean*0067       BULK_timeAve(bi,bj) = BULK_timeAve(bi,bj)+deltaTclock
7753507405 Curt*0068
                0069 #endif  /*ALLOW_TIMEAVE*/
                0070 #endif  /*ALLOW_BULK_FORCE*/
dd80d278b6 Jean*0071
6a1d3c464b Jean*0072       RETURN
                0073       END
7753507405 Curt*0074