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

Timestamp:

2006-10-05T17:53:31+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_072 : CT : - lights modifications to ensure good control prints for the restartability & reproductibility tests

• build all ocean.output_XXXX files only if logical namelist ln_ctl is TRUE, in all case keep the initial ocean.output file (require a new namelist namctl)
• conditionning some prints using nprint parameter

File:

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

trunk/NEMO/OPA_SRC/DOM/domhgr.F90

@@ -365 +365 @@
          gphi0 = zphi1 + zsin_alpha * ze1deg * FLOAT( jpjglo-2 )
 
-         IF(lwp) WRITE(numout,*) 'ze1', ze1, 'cosalpha', zcos_alpha, 'sinalpha', zsin_alpha
-         IF(lwp) WRITE(numout,*) 'ze1deg', ze1deg, 'glam0', glam0, 'gphi0', gphi0
+         IF( nprint==1 .AND. lwp )   THEN
+            WRITE(numout,*) '          ze1', ze1, 'cosalpha', zcos_alpha, 'sinalpha', zsin_alpha
+            WRITE(numout,*) '          ze1deg', ze1deg, 'glam0', glam0, 'gphi0', gphi0
+         ENDIF
 
          DO jj = 1, jpj
@@ -404 +406 @@
       ! ----------------
 
-      IF(lwp .AND. .NOT.ln_rstart ) THEN
+      IF( lwp .AND. .NOT.ln_rstart ) THEN
          WRITE(numout,*)
          WRITE(numout,*) '          longitude and e1 scale factors'
@@ -468 +470 @@
          ff(:,:) = ( zf0  + zbeta * gphif(:,:) * 1.e+3 )                        ! f = f0 +beta* y ( y=0 at south)
          
-         IF(lwp) WRITE(numout,*) 
-         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(lwp) THEN
+            WRITE(numout,*) 
+            WRITE(numout,*) '          Beta-plane: Beta parameter = constant = ', ff(nldi,nldj)
+            WRITE(numout,*) '          Coriolis parameter varies from ', ff(nldi,nldj),' to ', ff(nldi,nlej)
+         ENDIF
          IF( lk_mpp ) THEN 
             zminff=ff(nldi,nldj)
@@ -476 +480 @@
             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
+            IF(lwp) WRITE(numout,*) '          Coriolis parameter varies globally from ', zminff,' to ', zmaxff
          END IF
 
@@ -487 +491 @@
          ff(:,:) = ( zf0 + zbeta * ABS( gphif(:,:) - zphi0 ) * rad * ra )   ! f = f0 +beta* y ( y=0 at south)
 
-         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(lwp) THEN
+            WRITE(numout,*) 
+            WRITE(numout,*) '          Beta-plane and rotated domain : '
+            WRITE(numout,*) '          Coriolis parameter varies in this processor from ', ff(nldi,nldj),' to ', ff(nldi,nlej)
+         ENDIF
+
          IF( lk_mpp ) THEN 
             zminff=ff(nldi,nldj)
@@ -494 +502 @@
             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
+            IF(lwp) WRITE(numout,*) '          Coriolis parameter varies globally from ', zminff,' to ', zmaxff
          END IF