Changeset 13024 for utils/tools_dev_r12970_AGRIF_CMEMS/DOMAINcfg/src/mppini.F90

Timestamp:

2020-06-03T16:26:23+02:00 (4 years ago)

Author:

rblod

Message:

First version of new nesting tools merged with domaincfg, see ticket #2129

File:

• utils/tools_dev_r12970_AGRIF_CMEMS/DOMAINcfg/src/mppini.F90 (modified) (2 diffs)

utils/tools_dev_r12970_AGRIF_CMEMS/DOMAINcfg/src/mppini.F90

@@ -99 +99 @@
          CALL ctl_stop( 'mpp_init: equality  jpni = jpnj = jpnij = 1 is not satisfied',   &
             &           'the domain is lay out for distributed memory computing!' )
+
+#if defined key_agrif
+     IF (.not.agrif_root()) THEN
+        CALL agrif_nemo_init
+     ENDIF
+
+      IF( .NOT. Agrif_Root() ) THEN       ! AGRIF children: specific setting (cf. agrif_user.F90)
+         print *,'nbcellsx = ',nbcellsx,nbghostcells_x
+         print *,'nbcellsy = ',nbcellsy,nbghostcells_y_s,nbghostcells_y_n
+         IF( jpiglo /= nbcellsx + 2 + 2*nbghostcells_x ) THEN
+            IF(lwp) THEN
+               WRITE(numout,*)
+               WRITE(numout,*) 'jpiglo should be: ', nbcellsx + 2 + 2*nbghostcells_x
+            ENDIF        
+            CALL ctl_stop( 'STOP', 'mpp_init: Agrif children requires jpiglo == nbcellsx + 2 + 2*nbghostcells_x' )
+         ENDIF   
+         IF( jpjglo /= nbcellsy + 2 + nbghostcells_y_s + nbghostcells_y_n ) THEN
+            IF(lwp) THEN
+               WRITE(numout,*)
+               WRITE(numout,*) 'jpjglo shoud be: ', nbcellsy + 2 + nbghostcells_y_s + nbghostcells_y_n
+            ENDIF        
+            CALL ctl_stop( 'STOP', &
+                'mpp_init: Agrif children requires jpjglo == nbcellsy + 2 + nbghostcells_y_s + nbghostcells_y_n' )
+         ENDIF   
+         IF( ln_use_jattr )   CALL ctl_stop( 'STOP', 'mpp_init:Agrif children requires ln_use_jattr = .false. ' )
+      ENDIF
+#endif
          !
    END SUBROUTINE mpp_init
@@ -274 +301 @@
             IF(lwp) THEN
                WRITE(numout,*)
-               WRITE(numout,*) 'jpiglo shoud be: ', nbcellsx + 2 + 2*nbghostcells
+               WRITE(numout,*) 'jpiglo should be: ', nbcellsx + 2 + 2*nbghostcells
             ENDIF        
             CALL ctl_stop( 'STOP', 'mpp_init: Agrif children requires jpiglo == nbcellsx + 2 + 2*nbghostcells' )