Changeset 72 for trunk/NEMO/OPA_SRC/DOM/domain.F90

Timestamp:

2004-04-22T14:23:01+02:00 (20 years ago)

Author:

opalod

Message:

CT : BUGFIX046 : Bug correction in mpp case, use mpp_max(ze1max ... instead of mpp_min(ze1max)

File:

• trunk/NEMO/OPA_SRC/DOM/domain.F90 (modified) (7 diffs)

trunk/NEMO/OPA_SRC/DOM/domain.F90

@@ -100 +100 @@
          hv(:,:) = hv(:,:) + fse3v(:,:,jk) * vmask(:,:,jk)
       END DO
-#else
+#elif defined key_dynspg_rl
       ! Inverse of the local depth
       hur(:,:) = fse3u(:,:,1)             ! Lower bound : thickness of the first model level
@@ -161 +161 @@
 
       IF(lwp) THEN
-         WRITE(numout,*) 'namrun'
-         WRITE(numout,*) ' job number                      no        = ', no
-         WRITE(numout,*) ' experiment name for output      cexper    = ', cexper
-         WRITE(numout,*) ' restart logical                 ln_rstart = ', ln_rstart
-         WRITE(numout,*) ' control of time step            nrstdt    = ', nrstdt
-         WRITE(numout,*) ' number of the first time step   nit000    = ', nit000
-         WRITE(numout,*) ' number of the last time step    nitend    = ', nitend
-         WRITE(numout,*) ' initial calendar date aammjj    ndate0    = ', ndate0
-         WRITE(numout,*) ' leap year calendar (0/1)        nleapy    = ', nleapy
-         WRITE(numout,*) ' initial state output            ninist    = ', ninist
-         WRITE(numout,*) ' level of print                  nprint    = ', nprint
-         WRITE(numout,*) ' frequency of restart file       nstock    = ', nstock
-         WRITE(numout,*) ' frequency of output file        nwrite    = ', nwrite
-         WRITE(numout,*) ' runoff option                   nrunoff   = ', nrunoff
-         WRITE(numout,*) ' run control (for debugging)     ln_ctl    = ', ln_ctl
-      ENDIF
-
-      l_ctl = ln_ctl .AND. lwp             ! trend control print on the 1st processor only
-
-      ndastp = ndate0                     ! Assign initial date to current date
+         WRITE(numout,*) '        Namelist namrun'
+         WRITE(numout,*) '           job number                      no        = ', no
+         WRITE(numout,*) '           experiment name for output      cexper    = ', cexper
+         WRITE(numout,*) '           restart logical                 ln_rstart = ', ln_rstart
+         WRITE(numout,*) '           control of time step            nrstdt    = ', nrstdt
+         WRITE(numout,*) '           number of the first time step   nit000    = ', nit000
+         WRITE(numout,*) '           number of the last time step    nitend    = ', nitend
+         WRITE(numout,*) '           initial calendar date aammjj    ndate0    = ', ndate0
+         WRITE(numout,*) '           leap year calendar (0/1)        nleapy    = ', nleapy
+         WRITE(numout,*) '           initial state output            ninist    = ', ninist
+         WRITE(numout,*) '           level of print                  nprint    = ', nprint
+         WRITE(numout,*) '           frequency of restart file       nstock    = ', nstock
+         WRITE(numout,*) '           frequency of output file        nwrite    = ', nwrite
+         WRITE(numout,*) '           runoff option                   nrunoff   = ', nrunoff
+         WRITE(numout,*) '           run control (for debugging)     ln_ctl    = ', ln_ctl
+      ENDIF
+
+      l_ctl = ln_ctl .AND. lwp       ! trend control print on the 1st processor only
+
+      ndastp = ndate0                ! Assign initial date to current date
 
 
@@ -202 +202 @@
       CASE (  1 ) 
          CALL ioconf_calendar('gregorian')
-         IF(lwp) WRITE(numout,*)'The IOIPSL calendar is "gregorian", i.e. leap year'
+         IF(lwp) WRITE(numout,*) '           The IOIPSL calendar is "gregorian", i.e. leap year'
       CASE (  0 )
          CALL ioconf_calendar('noleap')
-         IF(lwp) WRITE(numout,*)'The IOIPSL calendar is "noleap", i.e. no leap year'
+         IF(lwp) WRITE(numout,*) '           The IOIPSL calendar is "noleap", i.e. no leap year'
       CASE ( 30 )
          CALL ioconf_calendar('360d')
-         IF(lwp) WRITE(numout,*)'The IOIPSL calendar is "360d", i.e. 360 days in a year'
+         IF(lwp) WRITE(numout,*) '           The IOIPSL calendar is "360d", i.e. 360 days in a year'
       END SELECT
 #endif
@@ -227 +227 @@
       END SELECT
       IF(lwp) THEN
-         WRITE(numout,*) ' '
-         WRITE(numout,*) ' nb of days per year      raajj = ', raajj,' days'
-         WRITE(numout,*) ' nb of seconds per year   raass = ', raass, ' s'
-         WRITE(numout,*) ' nb of seconds per month  rmoss = ', rmoss, ' s'
+         WRITE(numout,*)
+         WRITE(numout,*) '           nb of days per year      raajj = ', raajj,' days'
+         WRITE(numout,*) '           nb of seconds per year   raass = ', raass, ' s'
+         WRITE(numout,*) '           nb of seconds per month  rmoss = ', rmoss, ' s'
       ENDIF
 
@@ -238 +238 @@
 
       IF(lwp) THEN
-         WRITE(numout,*) 'namdom'
-         WRITE(numout,*) ' '
-         WRITE(numout,*) ' flag read/compute bathymetry   ntopo     = ', ntopo
-         WRITE(numout,*) ' minimum thickness of partial   e3zps_min = ', e3zps_min, ' (m)'
-         WRITE(numout,*) '    step level                  e3zps_rat = ', e3zps_rat
-         WRITE(numout,*) ' flag read/compute coordinates  ngrid     = ', ngrid
-         WRITE(numout,*) ' flag write mesh/mask file(s)   nmsh      = ', nmsh
-         WRITE(numout,*) '      = 0   no file created                 '
-         WRITE(numout,*) '      = 1   mesh_mask                       '
-         WRITE(numout,*) '      = 2   mesh and mask                   '
-         WRITE(numout,*) '      = 3   mesh_hgr, msh_zgr and mask      '
-         WRITE(numout,*) ' acceleration of converge       nacc      = ', nacc
-         WRITE(numout,*) ' asselin time filter parameter  atfp      = ', atfp
-         WRITE(numout,*) ' time step                      rdt       = ', rdt
-         WRITE(numout,*) ' minimum time step on tracers   rdtmin    = ', rdtmin
-         WRITE(numout,*) ' maximum time step on tracers   rdtmax    = ', rdtmax
-         WRITE(numout,*) ' depth variation tracer step    rdth      = ', rdth
+         WRITE(numout,*)
+         WRITE(numout,*) '        Namelist namdom'
+         WRITE(numout,*) '           flag read/compute bathymetry   ntopo     = ', ntopo
+         WRITE(numout,*) '           minimum thickness of partial   e3zps_min = ', e3zps_min, ' (m)'
+         WRITE(numout,*) '              step level                  e3zps_rat = ', e3zps_rat
+         WRITE(numout,*) '           flag read/compute coordinates  ngrid     = ', ngrid
+         WRITE(numout,*) '           flag write mesh/mask file(s)   nmsh      = ', nmsh
+         WRITE(numout,*) '                = 0   no file created                 '
+         WRITE(numout,*) '                = 1   mesh_mask                       '
+         WRITE(numout,*) '                = 2   mesh and mask                   '
+         WRITE(numout,*) '                = 3   mesh_hgr, msh_zgr and mask      '
+         WRITE(numout,*) '           acceleration of converge       nacc      = ', nacc
+         WRITE(numout,*) '           asselin time filter parameter  atfp      = ', atfp
+         WRITE(numout,*) '           time step                      rdt       = ', rdt
+         WRITE(numout,*) '           minimum time step on tracers   rdtmin    = ', rdtmin
+         WRITE(numout,*) '           maximum time step on tracers   rdtmax    = ', rdtmax
+         WRITE(numout,*) '           depth variation tracer step    rdth      = ', rdth
          IF( lk_ice_lim ) THEN
-            WRITE(numout,*) ' ice model coupling frequency   nfice  = ', nfice
+            WRITE(numout,*) '           ice model coupling frequency      nfice  = ', nfice
             nfbulk = nfice
          ENDIF
-         WRITE(numout,*) ' bulk computation frequency     nfbulk = ', nfbulk, ' = nfice if ice model used'
-         WRITE(numout,*) ' flag closed sea or not         nclosea= ', nclosea
-
+         WRITE(numout,*) '           bulk computation frequency       nfbulk  = ', nfbulk, ' = nfice if ice model used'
+         WRITE(numout,*) '           flag closed sea or not           nclosea = ', nclosea
       ENDIF
 
@@ -271 +270 @@
       READ  ( numnam, namcla )
       IF(lwp) THEN
-         WRITE(numout,*) ' namcla'
-         WRITE(numout,*) ' cross land advection        n_cla           = ',n_cla
+         WRITE(numout,*)
+         WRITE(numout,*) '        Namelist namcla'
+         WRITE(numout,*) '           cross land advection        n_cla        = ',n_cla
       ENDIF
 
@@ -307 +307 @@
       IF( lk_mpp )   CALL mpp_min( ze1min )   ! min over the global domain
       IF( lk_mpp )   CALL mpp_min( ze2min )
-      IF( lk_mpp )   CALL mpp_min( ze1max )
-      IF( lk_mpp )   CALL mpp_min( ze2max )
+      IF( lk_mpp )   CALL mpp_max( ze1max )   ! max over the global domain
+      IF( lk_mpp )   CALL mpp_max( ze2max )
 
       iloc  = MINLOC( e1t(:,:), mask = tmask(:,:,1) == 1.e0 )