New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
domcfg.F90 in trunk/NEMO/OPA_SRC/DOM – NEMO

source: trunk/NEMO/OPA_SRC/DOM/domcfg.F90 @ 65

Last change on this file since 65 was 65, checked in by opalod, 20 years ago

CT : BUGFIX039 : Indices and comments correction for zoom functionalities
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 10.3 KB
RevLine 
[3]1MODULE domcfg
2   !!==============================================================================
3   !!                       ***  MODULE domcfg   ***
4   !! Ocean initialization : domain configuration initialization
5   !!==============================================================================
6
7   !!----------------------------------------------------------------------
8   !!   dom_cfg        : initialize the domain configuration
9   !!----------------------------------------------------------------------
10   !! * Modules used
11   USE dom_oce         ! ocean space and time domain
12   USE phycst          ! physical constants
13   USE in_out_manager  ! I/O manager
14   USE lib_mpp         ! distributed memory computing library
15   USE solisl          ! ???
16
17   IMPLICIT NONE
18   PRIVATE
19
20   !! * Routine accessibility
21   PUBLIC dom_cfg        ! called by opa.F90
22   !!----------------------------------------------------------------------
23   !!   OPA 9.0 , LODYC-IPSL  (2003)
24   !!----------------------------------------------------------------------
25
26CONTAINS
27
28   SUBROUTINE dom_cfg
29      !!----------------------------------------------------------------------
30      !!                  ***  ROUTINE dom_cfg  ***
31      !!                   
32      !! ** Purpose :   set the domain configuration
33      !!
34      !! ** Method  :
35      !!
36      !! History :
37      !!   9.0  !  03-09  (G. Madec)  Original code
38      !!----------------------------------------------------------------------
39      !! * Local declarations
40      INTEGER ::   iconf = 0         ! temporary integers
41      !!----------------------------------------------------------------------
42
43      IF(lwp) THEN
44         WRITE(numout,*)
45         WRITE(numout,*) 'dom_cfg : set the ocean configuration'
46         WRITE(numout,*) '~~~~~~~      ocean model configuration used :',   &
47            &                             ' cp_cfg = ', cp_cfg, ' jp_cfg = ', jp_cfg
48      ENDIF
49
50      ! Global domain boundary conditions
51      ! ---------------------------------
52      IF(lwp) THEN
53         WRITE(numout,*) '          global domain lateral boundaries'
54
55         IF( jperio == 0 ) WRITE(numout,*) '             jperio= 0, closed'
56         IF( jperio == 1 ) WRITE(numout,*) '             jperio= 1, cyclic east-west'
57         IF( jperio == 2 ) WRITE(numout,*) '             jperio= 2, equatorial symmetric'
58         IF( jperio == 3 ) WRITE(numout,*) '             jperio= 3, north fold with T-point pivot'
59         IF( jperio == 4 ) WRITE(numout,*) '             jperio= 4, cyclic east-west and',   &
60                                                                  ' north fold with T-point pivot'
61         IF( jperio == 5 ) WRITE(numout,*) '             jperio= 5, north fold with F-point pivot'
62         IF( jperio == 6 ) WRITE(numout,*) '             jperio= 6, cyclic east-west and',   &
63                                                                  ' north fold with F-point pivot'
64      ENDIF
65      IF( jperio <  0 .OR. jperio > 6 ) THEN
66          IF(lwp) WRITE(numout,cform_err)
67          IF(lwp) WRITE(numout,*) 'jperio is out of range'
68          nstop = nstop + 1
69      ENDIF
70
71
72      ! global domain versus zoom and/or local domain
73      ! ---------------------------------------------
74
75      CALL dom_glo 
76
77   END SUBROUTINE dom_cfg
78
79
80   SUBROUTINE dom_glo
81      !!----------------------------------------------------------------------
82      !!                     ***  ROUTINE dom_glo  ***
83      !!
84      !! ** Purpose :   initialization for global domain, zoom and local domain
85      !!
86      !! ** Method  :   
87      !!
88      !! ** Action  : - mig  , mjg :
89      !!              - mi0  , mi1   :
90      !!              - mj0, , mj1   :
91      !!
92      !! History :
93      !!   8.5  !  02-08  (G. Madec)    Original code
94      !!----------------------------------------------------------------------
95      !! * Local declarations
96      INTEGER ::   ji, jj            ! dummy loop argument
97      !!----------------------------------------------------------------------
98
99      ! Local domain
100      ! ============
101
[65]102      ! local domain indices ==> data domain indices
[3]103      DO ji = 1, jpi
104        mig(ji) = ji + jpizoom - 1 + nimpp - 1
105      END DO
106      DO jj = 1, jpj
107        mjg(jj) = jj + jpjzoom - 1 + njmpp - 1
108      END DO
109
[65]110      ! data domain indices ==> local domain indices
111      ! (return (m.0,m.1)=(1,0) if data domain gridpoint is to the west/south of the
112      ! local domain, or (m.0,m.1)=(jp.+1,jp.) to the east/north of local domain.
113      DO ji = 1, jpidta
114        mi0(ji) = MAX( 1, MIN( ji - jpizoom + 1 - nimpp + 1, jpi+1 ) )
115        mi1(ji) = MAX( 0, MIN( ji - jpizoom + 1 - nimpp + 1, jpi   ) )
[3]116      END DO
[65]117      DO jj = 1, jpjdta
118        mj0(jj) = MAX( 1, MIN( jj - jpjzoom + 1 - njmpp + 1, jpj+1 ) )
119        mj1(jj) = MAX( 0, MIN( jj - jpjzoom + 1 - njmpp + 1, jpj   ) )
[3]120      END DO
121
122      IF(lwp) THEN
123         WRITE(numout,*)
[65]124         WRITE(numout,*) 'dom_glo : domain: data / local '
[3]125         WRITE(numout,*) '~~~~~~~ '
126         WRITE(numout,*) '          data input domain    : jpidta = ', jpidta,   &
127            &                                            ' jpjdta = ', jpjdta, ' jpkdta = ', jpkdta
128         WRITE(numout,*) '          global or zoom domain: jpiglo = ', jpiglo,   &
129            &                                            ' jpjglo = ', jpjglo, ' jpk    = ', jpk
130         WRITE(numout,*) '          local domain         : jpi    = ', jpi   ,   &
131            &                                            ' jpj    = ', jpj   , ' jpk    = ', jpk
132         WRITE(numout,*)
133         WRITE(numout,*) '          south-west indices    jpizoom = ', jpizoom,   &
134            &                                           ' jpjzoom = ', jpjzoom
135         WRITE(numout,*)
[65]136         WRITE(numout,*) '          conversion local  ==> data i-index domain'
[3]137         WRITE(numout,25)              (mig(ji),ji = 1,jpi)
138         WRITE(numout,*)
[65]139         WRITE(numout,*) '          conversion data   ==> local  i-index domain'
[3]140         WRITE(numout,*) '             starting index'
[65]141         WRITE(numout,25)              (mi0(ji),ji = 1,jpidta)
[3]142         WRITE(numout,*) '             ending index'
[65]143         WRITE(numout,25)              (mi1(ji),ji = 1,jpidta)
[3]144         WRITE(numout,*)
[65]145         WRITE(numout,*) '          conversion local  ==> data i-index domain'
[3]146         WRITE(numout,25)              (mjg(jj),jj = 1,jpj)
147         WRITE(numout,*)
[65]148         WRITE(numout,*) '          conversion data  ==> local i-index domain'
[3]149         WRITE(numout,*) '             starting index'
[65]150         WRITE(numout,25)              (mj0(jj),jj = 1,jpjdta)
[3]151         WRITE(numout,*) '             ending index'
[65]152         WRITE(numout,25)              (mj1(jj),jj = 1,jpjdta)
[3]153      ENDIF
154 25   FORMAT( 100(10x,19i4,/) )
155
156      ! Zoom domain
157      ! ===========
158
159      ! zoom control
[65]160      IF( jpiglo + jpizoom - 1  >  jpidta .OR.   &
161          jpjglo + jpjzoom - 1  >  jpjdta      ) THEN
[3]162         IF(lwp)WRITE(numout,cform_err)
163         IF(lwp)WRITE(numout,*)' global or zoom domain exceed the data domain ! '
164         nstop = nstop + 1
165      ENDIF
166
167      ! set zoom flag
168      IF ( jpiglo < jpidta .OR. jpjglo < jpjdta )   lzoom = .TRUE.
169
170      ! set zoom type flags
[65]171      IF( lzoom .AND. jpizoom /= 1 )   lzoom_w = .TRUE.                     !
[3]172      IF( lzoom .AND. jpjzoom /= 1 )   lzoom_s = .TRUE.
[65]173      IF( lzoom .AND. jpiglo + jpizoom -1 /= jpidta )   lzoom_e = .TRUE.
[3]174      IF( lzoom .AND. jpjglo + jpjzoom -1 /= jpjdta )   lzoom_n = .TRUE.
175
176      IF(lwp) THEN
177         WRITE(numout,*)
178         WRITE(numout,*) '          zoom flags : '
179         WRITE(numout,*) '             lzoom   = ', lzoom  , ' (T = zoom, F = global )'
180         WRITE(numout,*) '             lzoom_e = ', lzoom_e, ' (T = forced closed east  boundary)'
181         WRITE(numout,*) '             lzoom_w = ', lzoom_w, ' (T = forced closed west  boundary)'
182         WRITE(numout,*) '             lzoom_s = ', lzoom_s, ' (T = forced closed South boundary)'
183         WRITE(numout,*) '             lzoom_n = ', lzoom_n, ' (T = forced closed North boundary)'
184      ENDIF
185      IF(  ( lzoom_e .OR. lzoom_w )  .AND.  ( jperio == 1 .OR. jperio == 4 .OR. jperio == 6 )  ) THEN
186         IF(lwp)WRITE(numout,cform_err)
187         IF(lwp)WRITE(numout,*)' Your zoom choice is inconsistent with east-west cyclic boundary condition'
188         nstop = nstop + 1
189      ENDIF
190      IF(  lzoom_n  .AND.  ( 3 <= jperio .AND. jperio <= 6 )  ) THEN
191         IF(lwp)WRITE(numout,cform_err)
192         IF(lwp)WRITE(numout,*)' Your zoom choice is inconsistent with North fold boundary condition'
193         nstop = nstop + 1
194      ENDIF
[65]195      IF(  lzoom  .AND.  lk_isl  ) THEN
[3]196         IF(lwp)WRITE(numout,cform_err)
197         IF(lwp)WRITE(numout,*)' key_islands and zoom are not allowed'
198         nstop = nstop + 1
199      ENDIF
200
201      ! Pre-defined arctic/antarctic zoom of ORCA configuration flag
202      IF( cp_cfg == "orca" ) THEN
203         SELECT CASE ( jp_cfg )
204         !                                        ! =======================
205         CASE ( 2 )                               !  ORCA_R2 configuration
206            !                                     ! =======================
207            IF(  jpiglo  == 142    .AND. jpjglo  ==  53 .AND.   &
208               & jpizoom ==  21    .AND. jpjzoom ==  97         )   lzoom_arct = .TRUE.
209            IF(  jpiglo  == jpidta .AND. jpjglo  ==  50 .AND.   &
210               & jpizoom ==   1    .AND. jpjzoom ==   1         )   lzoom_anta = .TRUE.
211            !                                     ! =======================
212         CASE ( 05 )                              !  ORCA_R05 configuration
213            !                                     ! =======================
214            IF(  jpiglo  == 562    .AND. jpjglo  == 202 .AND.   &
215               & jpizoom ==  81    .AND. jpjzoom == 301         )   lzoom_arct = .TRUE.
216            IF(  jpiglo  == jpidta .AND. jpjglo  == 187 .AND.   &
217               & jpizoom ==   1    .AND. jpjzoom ==   1         )   lzoom_anta = .TRUE.
218         END SELECT
219         !
220         IF(lwp) WRITE(numout,*) '          ORCA configuration: antarctic/arctic zoom flags : '
221         IF(lwp) WRITE(numout,*) '             lzoom_arct = ', lzoom_arct, ' (T=   arctic zoom, F=global)'
222         IF(lwp) WRITE(numout,*) '             lzoom_anta = ', lzoom_anta, ' (T=antarctic zoom, F=global)'
223         !
224      ENDIF
225         
226   END SUBROUTINE dom_glo
227
228   !!======================================================================
229END MODULE domcfg
Note: See TracBrowser for help on using the repository browser.