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
0006
0007
a11169c200 Jean*0008
6451e229ef Jean*0009
0010
a11169c200 Jean*0011
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
0026
0027 INTEGER myThid
0028
0029
0030 INTEGER bi, bj
0031 _RL xLo, xHi, yLo, yHi
0032 CHARACTER*(MAX_LEN_MBUF) msgBuf
0033
0034
0035
0036
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