Changeset 496 for trunk/NEMO/OFF_SRC/opa.F90

Timestamp:

2006-09-12T12:59:38+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_063:CE:integration of the control print option for debugging

File:

• trunk/NEMO/OFF_SRC/opa.F90 (modified) (3 diffs)

trunk/NEMO/OFF_SRC/opa.F90

@@ -79 +79 @@
       !! * Local declarations
       INTEGER ::   istp       ! time step index
+      CHARACTER (len=20) ::   namelistname
+      CHARACTER (len=28) ::   file_out
       CHARACTER (len=64) ::        &
          cform_aaa="( /, 'AAAAAAAA', / ) "     ! flag for output listing
+
       !!----------------------------------------------------------------------
-      
-      
+
       ! Initializations
       ! ===============
-      
+
+      file_out = 'ocean.output'
+
       ! open listing and namelist units
-      IF ( numout /= 0 .AND. numout /= 6 ) THEN 
-         OPEN( UNIT=numout, FILE='ocean.output', FORM='FORMATTED' )
-      ENDIF
-
-      OPEN( UNIT=numnam, FILE='namelist', FORM='FORMATTED', STATUS='OLD' )
-
-      IF(lwp) THEN
+      IF ( numout /= 0 .AND. numout /= 6 ) THEN
+         CALL ctlopn( numout, file_out, 'UNKNOWN', 'FORMATTED',   &
+            &         'SEQUENTIAL', 1, numout, .FALSE., 1 )
+      ENDIF
+
+      namelistname = 'namelist'
+      CALL ctlopn( numnam, namelistname, 'OLD', 'FORMATTED', 'SEQUENTIAL',   &
+         &           1, numout, .FALSE., 1 )
+
+      WRITE(numout,*)
+      WRITE(numout,*) '                 L O D Y C - I P S L'
+      WRITE(numout,*) '                     O P A model'
+      WRITE(numout,*) '            Ocean General Circulation Model'
+      WRITE(numout,*) '               version OPA 9.0  (2005) '
+      WRITE(numout,*)
+      WRITE(numout,*)
+
+      ! Nodes selection
+      narea = mynode()
+      narea = narea + 1    ! mynode return the rank of proc (0 --> jpnij -1 )
+      lwp   = narea == 1
+
+      IF( lk_mpp )   THEN
+         CLOSE( numout )       ! standard model output file
+         WRITE(file_out,FMT="('ocean.output_',I4.4)") narea-1
+         IF ( numout /= 0 .AND. numout /= 6 ) THEN
+            CALL ctlopn( numout, file_out, 'UNKNOWN', 'FORMATTED',   &
+                 &         'SEQUENTIAL', 1, numout, .FALSE., 1 )
+         ENDIF
+         !
          WRITE(numout,*)
          WRITE(numout,*) '                 L O D Y C - I P S L'
          WRITE(numout,*) '                     O P A model'
          WRITE(numout,*) '            Ocean General Circulation Model'
-         WRITE(numout,*) '               version OPA 9.0  (2003)'
+         WRITE(numout,*) '               version OPA 9.0  (2005) '
+         WRITE(numout,*) '                   MPI Ocean output '
          WRITE(numout,*)
-      ENDIF
-
-      ! Nodes selection
-      narea = mynode()
-      narea = narea + 1    ! mynode return the rank of proc (0 --> jpnij -1 )
-      lwp   = narea == 1
+         WRITE(numout,*)
+      ENDIF
+
 
       !                                     ! ============================== !
@@ -135 +160 @@
                        CALL bn2( tn, sn, rn2 )              ! before Brunt-Vaisala frequency
 
-      IF( lk_zps    )   CALL zps_hde( nit000, tn, sn, rhd,  &  ! Partial steps: before Horizontal DErivative
+      IF( ln_zps    )   CALL zps_hde( nit000, tn, sn, rhd,  &  ! Partial steps: before Horizontal DErivative
                                           gtu, gsu, gru, &  ! of t, s, rd at the bottom ocean level
                                           gtv, gsv, grv )
@@ -152 +177 @@
       CALL tra_qsr_init                         ! Solar radiation penetration
 
+#if ! defined key_off_degrad
       CALL ldf_tra_init                         ! Lateral ocean tracer physics
+#endif 
 
       CALL zdf_init                             ! Vertical ocean physics