Back to home page

MITgcm

 
 

    


File indexing completed on 2020-01-15 06:10:57 UTC

view on githubraw file Latest commit 15f808dc on 2019-09-13 16:02:45 UTC
6451e229ef Jean*0001 #include "FLT_OPTIONS.h"
                0002 
                0003       SUBROUTINE FLT_INIT_FIXED( myThid  )
                0004 
                0005 C     ==================================================================
                0006 C     SUBROUTINE FLT_INIT_FIXED
                0007 C     ==================================================================
a11169c200 Jean*0008 C     o initializes Fixed array
6451e229ef Jean*0009 C     ==================================================================
                0010 
a11169c200 Jean*0011 C     !USES:
                0012       IMPLICIT NONE
                0013 
6451e229ef Jean*0014 #include "SIZE.h"
                0015 #include "EEPARAMS.h"
                0016 #include "PARAMS.h"
                0017 #include "GRID.h"
ad773b031f Oliv*0018 #ifdef ALLOW_EXCH2
135f7f7fe0 Jean*0019 #include "W2_EXCH2_SIZE.h"
ad773b031f Oliv*0020 #include "W2_EXCH2_PARAMS.h"
                0021 #endif
4c099ab430 Oliv*0022 #include "FLT_SIZE.h"
6451e229ef Jean*0023 #include "FLT.h"
                0024 
                0025 C     == routine arguments ==
                0026 C     myThid - thread number for this instance of the routine.
                0027       INTEGER myThid
                0028 
                0029 C     == local variables ==
                0030       INTEGER bi, bj
                0031       _RL xLo, xHi, yLo, yHi
                0032       CHARACTER*(MAX_LEN_MBUF) msgBuf
                0033 
                0034 C     == end of interface ==
                0035 
                0036 C Set domain dependent constants:
                0037       xLo = -xG( 1,  1,1,1)
                0038       xHi = xG(sNx+1,1,1,1)
                0039       yLo = -yG(1, 1,  1,1)
                0040       yHi = yG(1,sNy+1,1,1)
                0041       DO bj=myByLo(myThid),myByHi(myThid)
                0042        DO bi=myBxLo(myThid),myBxHi(myThid)
                0043          xLo = MAX( xLo, -xG(1,1,bi,bj) )
                0044          xHi = MAX( xHi, xG(sNx+1,1,bi,bj) )
                0045          yLo = MAX( yLo, -yG(1,1,bi,bj) )
                0046          yHi = MAX( yHi, yG(1,sNy+1,bi,bj) )
                0047        ENDDO
                0048       ENDDO
6637358eea Jean*0049       _GLOBAL_MAX_RL( xLo, myThid )
                0050       _GLOBAL_MAX_RL( xHi, myThid )
                0051       _GLOBAL_MAX_RL( yLo, myThid )
                0052       _GLOBAL_MAX_RL( yHi, myThid )
6451e229ef Jean*0053 
                0054       _BEGIN_MASTER(myThid)
                0055        flt_xLo = -xLo
                0056        flt_xHi =  xHi
                0057        flt_yLo = -yLo
                0058        flt_yHi =  yHi
                0059        WRITE(msgBuf,'(A,2F18.8)')
                0060      &    '  FLT_xLo,xHi=', xLo, xHi
                0061          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
                0062      &                       SQUEEZE_RIGHT, myThid )
                0063          WRITE(msgBuf,'(A,2F18.8)')
                0064      &    '  FLT_yLo,yHi=', yLo, yHi
                0065          CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
                0066      &                       SQUEEZE_RIGHT, myThid )
                0067       _END_MASTER(myThid)
                0068 
f5995a4aae Gael*0069 #ifndef DEVEL_FLT_EXCH2
ad773b031f Oliv*0070 #ifdef ALLOW_EXCH2
                0071         WRITE(msgBuf,'(2A)') 'FLT_INIT_FIXED:',
15f808dcb0 gael*0072      &    ' exch2 is not yet fully supported in pkg/flt;'
                0073         CALL PRINT_ERROR( msgBuf, myThid )
                0074         WRITE(msgBuf,'(2A)') 'FLT_INIT_FIXED:',
                0075      &    ' if you want to try using it anyway then recompile with'
                0076         CALL PRINT_ERROR( msgBuf, myThid )
                0077         WRITE(msgBuf,'(2A)') 'FLT_INIT_FIXED:',
                0078      &    ' "#define DEVEL_FLT_EXCH2" added in FLT_OPTIONS.h'
ad773b031f Oliv*0079         CALL PRINT_ERROR( msgBuf, myThid )
                0080         STOP 'ABNORMAL END: S/R FLT_INIT_FIXED'
                0081 #endif
f5995a4aae Gael*0082 #endif
ad773b031f Oliv*0083 
6451e229ef Jean*0084       RETURN
                0085       END