Changeset 3875 for branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/trcnam_pisces.F90

Timestamp:

2013-04-18T16:38:06+02:00 (11 years ago)

Author:

clevy

Message:

Configuration Setting/Step? 1, see ticket:#1074

File:

• branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/trcnam_pisces.F90 (modified) (4 diffs)

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/trcnam_pisces.F90

@@ -49 +49 @@
       !!
       INTEGER :: jl, jn
+      INTEGER :: ios                 ! Local integer output status for namelist read
       TYPE(DIAG), DIMENSION(jp_pisces_2d)  :: pisdia2d
       TYPE(DIAG), DIMENSION(jp_pisces_3d)  :: pisdia3d
@@ -69 +70 @@
 #endif
       IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~'
-      CALL ctl_opn( numnatp, TRIM( clname ), 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
-
+      CALL ctl_opn( numnatp_ref, TRIM( clname )//'_ref', 'OLD'    , 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
+      CALL ctl_opn( numnatp_cfg, TRIM( clname )//'_cfg', 'OLD'    , 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
+      CALL ctl_opn( numonp     , 'output.namelist.pis' , 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
       !
       IF( .NOT.lk_iomput .AND. ln_diatrc ) THEN
@@ -88 +90 @@
          END DO
 
-         REWIND( numnatp )               ! 
-         READ  ( numnatp, nampisdia )
+         REWIND( numnatp_ref )              ! Namelist nampisdia in reference namelist : Pisces diagnostics
+         READ  ( numnatp_ref, nampisdia, IOSTAT = ios, ERR = 901)
+901      IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisdia in reference namelist', lwp )
+
+         REWIND( numnatp_cfg )              ! Namelist nampisdia in configuration namelist : Pisces diagnostics
+         READ  ( numnatp_cfg, nampisdia, IOSTAT = ios, ERR = 902 )
+902      IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisdia in configuration namelist', lwp )
+         WRITE ( numonp, nampisdia )
 
          DO jl = 1, jp_pisces_2d
@@ -140 +148 @@
          END DO
 
-         REWIND( numnatp )
-         READ  ( numnatp, nampisdbi )
+         REWIND( numnatp_ref )              ! Namelist nampisdbi in reference namelist : Pisces add. diagnostics
+         READ  ( numnatp_ref, nampisdbi, IOSTAT = ios, ERR = 903)
+903      IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisdbi in reference namelist', lwp )
+
+         REWIND( numnatp_cfg )              ! Namelist nampisdbi in configuration namelist : Pisces add. diagnostics
+         READ  ( numnatp_cfg, nampisdbi, IOSTAT = ios, ERR = 904 )
+904      IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisdbi in configuration namelist', lwp )
+         WRITE ( numonp, nampisdbi )
 
          DO jl = 1, jp_pisces_trd