/[MITgcm]/MITgcm/pkg/exch2/w2_set_single_facet.F
ViewVC logotype

Diff of /MITgcm/pkg/exch2/w2_set_single_facet.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by jmc, Tue May 12 19:40:33 2009 UTC revision 1.2 by jmc, Mon Jun 29 03:14:27 2009 UTC
# Line 17  C     Set-up simple single facet (domain Line 17  C     Set-up simple single facet (domain
17  C     !USES:  C     !USES:
18        IMPLICIT NONE        IMPLICIT NONE
19    
20  C      Tile toplogy settings data structures  C      Tile topology settings data structures
21  #include "SIZE.h"  #include "SIZE.h"
22  #include "EEPARAMS.h"  #include "EEPARAMS.h"
23  #include "W2_EXCH2_SIZE.h"  #include "W2_EXCH2_SIZE.h"
# Line 32  C               (Note: not relevant sinc Line 32  C               (Note: not relevant sinc
32  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
33  C     === Local variables ===  C     === Local variables ===
34  C     msgBuf     :: Informational/error meesage buffer  C     msgBuf     :: Informational/error meesage buffer
35  c     CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
36          INTEGER j
37  CEOP  CEOP
38    
39          WRITE(msgBuf,'(2A,I3,A)') 'W2_SET_SINGLE_FACET:',
40         &              ' preDefTopol=', preDefTopol, ' selected'
41          CALL PRINT_MESSAGE( msgBuf, W2_oUnit,SQUEEZE_RIGHT,myThid )
42    
43    C--   Number of facets:
44        nFacets = 1        nFacets = 1
 C-    set connection graph (topology)  
 C-    need to set/check facet dims  
45    
46        STOP 'ABNORMAL END: S/R W2_SET_SINGLE_FACET empty'  C--   Set Facet Edge connections (topology graph) ignoring any previous
47    C     setting from data.exch2 (Edges order: N,S,E,W <==> 1,2,3,4 )
48    C     face 1 N(=1) edge connects to face 1 S(=2) edge:
49          facet_link(1,1) = 1.2
50    C     face 1 S(=2) edge connects to face 1 N(=1) edge:
51          facet_link(2,1) = 1.1
52    C     face 1 E(=3) edge connects to face 1 W(=4) edge:
53          facet_link(3,1) = 1.4
54    C     face 1 W(=4) edge connects to face 1 E(=3) edge:
55          facet_link(4,1) = 1.3
56    
57    C--   Facet dimension: take the 1rst 2 numbers from facet_dims (if correct)
58          IF ( facet_dims(1).EQ.0 .AND. facet_dims(2).EQ.0 ) THEN
59    C-    Default: take global dimension from SIZE.h (will fail with blank tiles)
60            facet_dims(1) = Nx
61            facet_dims(2) = Ny
62          ENDIF
63          IF ( facet_dims(1).LE.0 .OR. facet_dims(2).LE.0 ) THEN
64             WRITE(msgBuf,'(2A,2I5)') 'W2_SET_SINGLE_FACET:',
65         &     ' unvalid 1rst 2 dimensions:', facet_dims(1), facet_dims(2)
66             CALL PRINT_MESSAGE( msgBuf, W2_oUnit,SQUEEZE_RIGHT,myThid )
67             CALL PRINT_ERROR( msgBuf, myThid )
68             CALL ALL_PROC_DIE( myThid )
69             STOP 'ABNORMAL END: S/R W2_SET_SINGLE_FACET: unvalid dims'
70          ENDIF
71          DO j=3,W2_maxNbFacets*2
72           IF ( facet_dims(j).NE.0 ) THEN
73             WRITE(msgBuf,'(2A,I5)') 'W2_SET_SINGLE_FACET:',
74         &     ' no more than 2 dims (X,Y) expected for single facet'
75             CALL PRINT_MESSAGE( msgBuf, W2_oUnit,SQUEEZE_RIGHT,myThid )
76             CALL PRINT_ERROR( msgBuf, myThid )
77             CALL ALL_PROC_DIE( myThid )
78             STOP 'ABNORMAL END: S/R W2_SET_SINGLE_FACET: unexpected dims'
79           ENDIF
80          ENDDO
81    
82        RETURN        RETURN
83        END        END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22