1 |
C $Header$ |
C $Header$ |
2 |
C $Name$ |
C $Name$ |
3 |
|
|
4 |
|
#include "PACKAGES_CONFIG.h" |
5 |
#include "CPP_OPTIONS.h" |
#include "CPP_OPTIONS.h" |
6 |
|
|
7 |
CBOP |
CBOP |
54 |
C == Local variables == |
C == Local variables == |
55 |
C msgBuf - Used for informational I/O. |
C msgBuf - Used for informational I/O. |
56 |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
57 |
|
#ifdef ALLOW_EXF |
58 |
|
INTEGER i, j, bi, bj |
59 |
|
#endif /* ALLOW_EXF */ |
60 |
CEOP |
CEOP |
61 |
|
|
62 |
C-- Set up vertical grid and coordinate system |
C-- Set up vertical grid and coordinate system |
78 |
_END_MASTER(myThid) |
_END_MASTER(myThid) |
79 |
ENDIF |
ENDIF |
80 |
|
|
81 |
|
#ifdef ALLOW_EXF |
82 |
|
C-- exf_interp assumes that 0 <= xG, xC <= 360 |
83 |
|
C This is a quick fix until this assumption is relaxed |
84 |
|
C and the interpolation weights are pre-computed. |
85 |
|
DO bj = myByLo(myThid), myByHi(myThid) |
86 |
|
DO bi = myBxLo(myThid), myBxHi(myThid) |
87 |
|
DO J=1-Oly,sNy+Oly |
88 |
|
DO I=1-Olx,sNx+Olx |
89 |
|
IF ( xG(I,J,bi,bj) .LT. 0 ) |
90 |
|
& xG(I,J,bi,bj) = xG(I,J,bi,bj) + 360 |
91 |
|
IF ( xC(I,J,bi,bj) .LT. 0 ) |
92 |
|
& xC(I,J,bi,bj) = xC(I,J,bi,bj) + 360 |
93 |
|
ENDDO |
94 |
|
ENDDO |
95 |
|
ENDDO |
96 |
|
ENDDO |
97 |
|
#endif /* ALLOW_EXF */ |
98 |
|
|
99 |
C-- Write certain grid data to files (useful for creating netCDF |
C-- Write certain grid data to files (useful for creating netCDF |
100 |
C and general post-analysis) |
C and general post-analysis) |
101 |
CALL WRITE_FLD_XY_RS( 'XC',' ',XC,0,myThid) |
CALL WRITE_FLD_XY_RS( 'XC',' ',XC,0,myThid) |