Changeset 434 for trunk/NEMO/OPA_SRC/DOM/domhgr.F90

Timestamp:

2006-04-10T17:46:12+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_044 : CT : update the passive tracers TOP component and the standard GYRE configuration

File:

• trunk/NEMO/OPA_SRC/DOM/domhgr.F90 (modified) (4 diffs)

trunk/NEMO/OPA_SRC/DOM/domhgr.F90

@@ -13 +13 @@
    USE phycst          ! physical constants
    USE in_out_manager  ! I/O manager
+   USE lib_mpp
 
    IMPLICIT NONE
@@ -106 +107 @@
          ztj, zuj, zvj, zfj,         &  !
          zphi0, zbeta, znorme,       &  !
-         zarg, zf0
+         zarg, zf0, zminff, zmaxff
       REAL(wp) ::   &
          zlam1, zcos_alpha, zim1 , zjm1 , ze1, ze1deg,   &
@@ -442 +443 @@
          
          IF(lwp) WRITE(numout,*) 
-         IF(lwp) WRITE(numout,*) ' Beta-plane: Beta parameter = constant = ', ff(1,1)
-         IF(lwp) WRITE(numout,*) ' Coriolis parameter varies from ', ff(1,1),' to ', ff(1,jpj)
-
-      CASE ( 5 )                     ! beta-plane and rotated domain
+         IF(lwp) WRITE(numout,*) ' Beta-plane: Beta parameter = constant = ', ff(nldi,nldj)
+         IF(lwp) WRITE(numout,*) ' Coriolis parameter varies from ', ff(nldi,nldj),' to ', ff(nldi,nlej)
+         IF( lk_mpp ) THEN 
+            zminff=ff(nldi,nldj)
+            zmaxff=ff(nldi,nlej)
+            CALL mpp_min( zminff )   ! min over the global domain
+            CALL mpp_max( zmaxff )   ! max over the global domain
+            IF(lwp) WRITE(numout,*) ' Coriolis parameter varies globally from ', zminff,' to ', zmaxff
+         END IF
+
+      CASE ( 5 )                     ! beta-plane and rotated domain (gyre configuration)
 
          zbeta = 2. * omega * COS( rad * ppgphi0 ) / ra                     ! beta at latitude ppgphi0
@@ -453 +461 @@
          ff(:,:) = ( zf0 + zbeta * ABS( gphif(:,:) - zphi0 ) * rad * ra )   ! f = f0 +beta* y ( y=0 at south)
 
-         IF(lwp) WRITE(numout,*) '          Beta-plane: Beta parameter = constant = ', ff(1,1)
-         IF(lwp) WRITE(numout,*) '                      Coriolis parameter varies from ', ff(1,1),' to ', ff(1,jpj)
+         IF(lwp) WRITE(numout,*) ' Beta-plane and rotated domain : '
+         IF(lwp) WRITE(numout,*) ' Coriolis parameter varies in this processor from ', ff(nldi,nldj),' to ', ff(nldi,nlej)
+         IF( lk_mpp ) THEN 
+            zminff=ff(nldi,nldj)
+            zmaxff=ff(nldi,nlej)
+            CALL mpp_min( zminff )   ! min over the global domain
+            CALL mpp_max( zmaxff )   ! max over the global domain
+            IF(lwp) WRITE(numout,*) ' Coriolis parameter varies globally from ', zminff,' to ', zmaxff
+         END IF
 
       END SELECT