pkg/bulk_force/bulkf_init_varia.F

6d54cf9ca1 Ed H*0001 #include "BULK_FORCE_OPTIONS.h"
7753507405 Curt*0002
e33db7ef40 Jean*0003 CBOP
                0004 C     !ROUTINE: BULKF_INIT_VARIA
                0005 C     !INTERFACE:
                0006       SUBROUTINE BULKF_INIT_VARIA( myThid )
                0007
                0008 C     *==========================================================*
                0009 C     | SUBROUTINE BULKF_INIT_VARIA
                0010 C     | o Initialize bulk formula arrays
                0011 C     *==========================================================*
                0012 C     \ev
7753507405 Curt*0013
e33db7ef40 Jean*0014 C     !USES:
                0015       IMPLICIT NONE
7753507405 Curt*0016 C     === Global variables ===
                0017 #include "SIZE.h"
                0018 #include "EEPARAMS.h"
679d149d01 Jean*0019 #include "BULKF_PARAMS.h"
7753507405 Curt*0020 #include "BULKF.h"
e33db7ef40 Jean*0021 #include "BULKF_TAVE.h"
7753507405 Curt*0022 #ifdef CONSERV_BULKF
                0023 #include "BULKF_CONSERV.h"
                0024 #endif
                0025
e33db7ef40 Jean*0026 C     !INPUT/OUTPUT PARAMETERS:
                0027 C     myThid ::  Number of this instance of BULKF_INIT_VARIA
7753507405 Curt*0028       INTEGER myThid
e33db7ef40 Jean*0029 CEOP
7753507405 Curt*0030
6a1d3c464b Jean*0031 #ifdef ALLOW_BULK_FORCE
7753507405 Curt*0032 C     == Local variables ==
6bf801d687 Jean*0033 C     bi,bj  :: tile indices
                0034 C     i,j    :: Loop counters
7753507405 Curt*0035       INTEGER bi, bj
6bf801d687 Jean*0036       INTEGER i, j
7753507405 Curt*0037
679d149d01 Jean*0038 C     set up bulk formula arrays to zero
                0039       DO bj = myByLo(myThid), myByHi(myThid)
                0040        DO bi = myBxLo(myThid), myBxHi(myThid)
                0041          DO j=1-Oly,sNy+Oly
                0042            DO i=1-Olx,sNx+Olx
6bf801d687 Jean*0043              Tair  (i,j,bi,bj) = 0. _d 0
                0044              Qair  (i,j,bi,bj) = 0. _d 0
                0045              Solar (i,j,bi,bj) = 0. _d 0
                0046              flwdwn(i,j,bi,bj) = 0. _d 0
                0047              cloud (i,j,bi,bj) = 0. _d 0
                0048              wspeed(i,j,bi,bj) = 0. _d 0
                0049              uwind (i,j,bi,bj) = 0. _d 0
                0050              vwind (i,j,bi,bj) = 0. _d 0
                0051              runoff(i,j,bi,bj) = 0. _d 0
9c2a091b47 Jean*0052 #ifdef ALLOW_FORMULA_AIM
6bf801d687 Jean*0053              thAir (i,j,bi,bj) = 0. _d 0
9c2a091b47 Jean*0054 #endif
6bf801d687 Jean*0055              qnetch(i,j,bi,bj) = 0. _d 0
                0056              empch (i,j,bi,bj) = 0. _d 0
679d149d01 Jean*0057            ENDDO
                0058          ENDDO
adedc74ecc Alis*0059 #ifdef ALLOW_TIMEAVE
6bf801d687 Jean*0060          BULK_timeAve(bi,bj) = 0. _d 0
                0061          DO j=1-Oly,sNy+Oly
679d149d01 Jean*0062            DO i=1-Olx,sNx+Olx
6bf801d687 Jean*0063               BULK_Qnet_Ave    (i,j,bi,bj) = 0. _d 0
                0064               BULK_EmPmR_Ave   (i,j,bi,bj) = 0. _d 0
                0065               BULK_fu_Ave      (i,j,bi,bj) = 0. _d 0
                0066               BULK_fv_Ave      (i,j,bi,bj) = 0. _d 0
                0067               BULK_latent_Ave  (i,j,bi,bj) = 0. _d 0
                0068               BULK_sensible_Ave(i,j,bi,bj) = 0. _d 0
                0069               BULK_evap_Ave    (i,j,bi,bj) = 0. _d 0
                0070               BULK_flwup_Ave   (i,j,bi,bj) = 0. _d 0
                0071               BULK_flwupnet_Ave(i,j,bi,bj) = 0. _d 0
                0072               BULK_solar_Ave   (i,j,bi,bj) = 0. _d 0
                0073               BULK_ssq_Ave     (i,j,bi,bj) = 0. _d 0
adedc74ecc Alis*0074            ENDDO
6bf801d687 Jean*0075          ENDDO
adedc74ecc Alis*0076 #endif /* ALLOW_TIMEAVE */
7753507405 Curt*0077 #ifdef CONSERV_BULKF
679d149d01 Jean*0078          DO j=1-Oly,sNy+Oly
                0079            DO i=1-Olx,sNx+Olx
6bf801d687 Jean*0080               CONS_Qnet (i,j,bi,bj) = 0. _d 0
                0081               CONS_EmPmR(i,j,bi,bj) = 0. _d 0
7753507405 Curt*0082            ENDDO
679d149d01 Jean*0083          ENDDO
                0084 #endif
                0085        ENDDO
                0086       ENDDO
7753507405 Curt*0087
                0088 #ifdef CONSERV_BULKF
6bf801d687 Jean*0089       constim = 0. _d 0
7753507405 Curt*0090 #endif
                0091
6a1d3c464b Jean*0092 #endif /* ALLOW_BULK_FORCE */
7753507405 Curt*0093
                0094       RETURN
                0095       END