pkg/ebm/ebm_load_climatology.F

d8206d87ee Patr*0001 #include "EBM_OPTIONS.h"
b9b591469d Jean*0002 #ifdef ALLOW_EXCH2
                0003 # include "W2_OPTIONS.h"
                0004 #endif /* ALLOW_EXCH2 */
d8206d87ee Patr*0005
2a68837bc4 Jean*0006 CBOP 0
                0007 C !ROUTINE: EBM_LOAD_CLIMATOLOGY
                0008
                0009 C !INTERFACE:
d8206d87ee Patr*0010       SUBROUTINE EBM_LOAD_CLIMATOLOGY( myThid )
                0011
2a68837bc4 Jean*0012 C     !DESCRIPTION:
                0013 C     *==========================================================*
                0014 C     | S/R EBM_LOAD_CLIMATOLOGY
                0015 C     *==========================================================*
                0016
                0017 C     !USES:
                0018       IMPLICIT NONE
d8206d87ee Patr*0019 C     === Global variables ===
                0020 #include "SIZE.h"
                0021 #include "EEPARAMS.h"
                0022 #include "PARAMS.h"
b9b591469d Jean*0023 #ifdef ALLOW_EXCH2
                0024 # include "W2_EXCH2_SIZE.h"
                0025 # include "W2_EXCH2_TOPOLOGY.h"
                0026 #endif /* ALLOW_EXCH2 */
                0027 #include "GRID.h"
                0028 #include "SET_GRID.h"
d8206d87ee Patr*0029 #include "FFIELDS.h"
                0030 #ifdef ALLOW_EBM
                0031 # include "EBM.h"
                0032 #endif
                0033
2a68837bc4 Jean*0034 C     !INPUT PARAMETERS:
                0035 C     myThid     :: my Thread Id number
                0036       INTEGER myThid
                0037 CEOP
d8206d87ee Patr*0038
                0039 #ifdef ALLOW_EBM
2a68837bc4 Jean*0040 C     !LOCAL VARIABLES:
b9b591469d Jean*0041       INTEGER gridNy
d8206d87ee Patr*0042       INTEGER bi,bj,i,j
                0043       _RL distY,tauX,tauMax,lY
                0044
b9b591469d Jean*0045 #ifdef ALLOW_EXCH2
                0046       gridNy = exch2_mydNy(1)
                0047 #else /* ALLOW_EXCH2 */
                0048       gridNy = Ny
                0049 #endif /* ALLOW_EXCH2 */
                0050
d8206d87ee Patr*0051       IF ( RunoffFile .NE. ' ' ) THEN
                0052        CALL READ_FLD_XY_RS( RunoffFile, ' ', Run, 0, myThid )
                0053       ENDIF
                0054
                0055       IF ( zonalWindFile .EQ. ' ' ) THEN
                0056 C      In cartesian yc, delY and ly are meters.
                0057 C      In spherical polar yc, delY and ly are degrees
                0058        tauMax = 0.1 _d 0
                0059        tauMax = 1.0 * 1./(delR(1)*rhonil)
                0060        lY = 0. _d 0
b9b591469d Jean*0061        DO j=1,gridNy-1
d8206d87ee Patr*0062         lY = lY + delY(j)
                0063        ENDDO
                0064        DO bj = myByLo(myThid), myByHi(myThid)
                0065         DO bi = myBxLo(myThid), myBxHi(myThid)
                0066          DO j=1,sNy
                0067           DO i=1,sNx
74305a5cbb Jean*0068 c          distY = (yC(i,j,bi,bj)-(yC0))/lY
                0069 c jmc:    yC0 has been removed (not used except here,
                0070 c         and since it was not initialised, was generaly zero):
                0071            distY = (yC(i,j,bi,bj)- 0.  )/lY
                0072 c          tauX  = -tauMax*cos(2. _d 0*PI*distY)
d8206d87ee Patr*0073            tauX  = tauMax*sin(PI*distY)
                0074            fu(i,j,bi,bj) = tauX
                0075           ENDDO
                0076          ENDDO
                0077         ENDDO
                0078        ENDDO
                0079       ENDIF
                0080
6637358eea Jean*0081       _EXCH_XY_RS( Run, myThid )
d4bd695091 Jean*0082
a24915ab1a Jean*0083       IF ( plotLevel.GE.debLevC ) THEN
2a68837bc4 Jean*0084         CALL PLOT_FIELD_XYRS( SST, 'Theta Climatology' , 1, myThid )
                0085         CALL PLOT_FIELD_XYRS( SSS, 'Salt  Climatology' , 1, myThid )
                0086         CALL PLOT_FIELD_XYRS( Run, 'Runoff Climatology', 1, myThid )
                0087         CALL PLOT_FIELD_XYRS(fu,'Wind Stress Climatology FU',1,myThid)
                0088         CALL PLOT_FIELD_XYRS(fv,'Wind Stress Climatology FV',1,myThid)
                0089       ENDIF
d8206d87ee Patr*0090
2a68837bc4 Jean*0091 #endif /* ALLOW_EBM */
d8206d87ee Patr*0092
                0093       RETURN
                0094       END