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"
                0021 #ifdef CONSERV_BULKF
                0022 #include "BULKF_CONSERV.h"
                0023 #endif
                0024
e33db7ef40 Jean*0025 C     !INPUT/OUTPUT PARAMETERS:
                0026 C     myThid ::  Number of this instance of BULKF_INIT_VARIA
7753507405 Curt*0027       INTEGER myThid
e33db7ef40 Jean*0028 CEOP
7753507405 Curt*0029
6a1d3c464b Jean*0030 #ifdef ALLOW_BULK_FORCE
7753507405 Curt*0031 C     == Local variables ==
6bf801d687 Jean*0032 C     bi,bj  :: tile indices
                0033 C     i,j    :: Loop counters
7753507405 Curt*0034       INTEGER bi, bj
6bf801d687 Jean*0035       INTEGER i, j
7753507405 Curt*0036
679d149d01 Jean*0037 C     set up bulk formula arrays to zero
                0038       DO bj = myByLo(myThid), myByHi(myThid)
                0039        DO bi = myBxLo(myThid), myBxHi(myThid)
b7411f1a84 Jean*0040          DO j=1-OLy,sNy+OLy
                0041            DO i=1-OLx,sNx+OLx
6bf801d687 Jean*0042              Tair  (i,j,bi,bj) = 0. _d 0
                0043              Qair  (i,j,bi,bj) = 0. _d 0
                0044              Solar (i,j,bi,bj) = 0. _d 0
                0045              flwdwn(i,j,bi,bj) = 0. _d 0
                0046              cloud (i,j,bi,bj) = 0. _d 0
                0047              wspeed(i,j,bi,bj) = 0. _d 0
                0048              uwind (i,j,bi,bj) = 0. _d 0
                0049              vwind (i,j,bi,bj) = 0. _d 0
                0050              runoff(i,j,bi,bj) = 0. _d 0
9c2a091b47 Jean*0051 #ifdef ALLOW_FORMULA_AIM
6bf801d687 Jean*0052              thAir (i,j,bi,bj) = 0. _d 0
9c2a091b47 Jean*0053 #endif
6bf801d687 Jean*0054              qnetch(i,j,bi,bj) = 0. _d 0
                0055              empch (i,j,bi,bj) = 0. _d 0
679d149d01 Jean*0056            ENDDO
                0057          ENDDO
7753507405 Curt*0058 #ifdef CONSERV_BULKF
b7411f1a84 Jean*0059          DO j=1-OLy,sNy+OLy
                0060            DO i=1-OLx,sNx+OLx
6bf801d687 Jean*0061               CONS_Qnet (i,j,bi,bj) = 0. _d 0
                0062               CONS_EmPmR(i,j,bi,bj) = 0. _d 0
7753507405 Curt*0063            ENDDO
679d149d01 Jean*0064          ENDDO
                0065 #endif
                0066        ENDDO
                0067       ENDDO
7753507405 Curt*0068
                0069 #ifdef CONSERV_BULKF
6bf801d687 Jean*0070       constim = 0. _d 0
7753507405 Curt*0071 #endif
                0072
6a1d3c464b Jean*0073 #endif /* ALLOW_BULK_FORCE */
7753507405 Curt*0074
                0075       RETURN
                0076       END