Changeset 627 for trunk/NEMO/TOP_SRC/prtctl_trc.F90

Timestamp:

2007-02-27T20:02:43+01:00 (17 years ago)

Author:

opalod

Message:

nemo_v2_bugfix_025:RB: suppress hard-coded logical units except for dimg and coupled

File:

• trunk/NEMO/TOP_SRC/prtctl_trc.F90 (modified) (6 diffs)

trunk/NEMO/TOP_SRC/prtctl_trc.F90

@@ -121 +121 @@
 
          ! Set logical unit
-         j_id = js - narea + 1
+         j_id = numid_trc(js - narea + 1)
          ! Set indices for the SUM control
          IF( .NOT. lsp_area ) THEN
@@ -153 +153 @@
                zsum          = SUM( ztab3d(nictls:nictle,njctls:njctle,1:kdir) &
                   &                 *zmask(nictls:nictle,njctls:njctle,1:kdir) )
-               WRITE(numid_trc(j_id),FMT="(3x,a,' : ',D23.16)") cl(jn), zsum-zvctl
+               WRITE(j_id,FMT="(3x,a,' : ',D23.16)") cl(jn), zsum-zvctl
                tra_ctl(jn,js) = zsum
             ENDDO
@@ -161 +161 @@
                zsum          = SUM( ztab3d(nictls:nictle,njctls:njctle,1:kdir) &
                   &               * zmask(nictls:nictle,njctls:njctle,1:kdir) )
-               WRITE(numid_trc(j_id),FMT="(3x,a,' : ',D23.16)") cl(jn), zsum
+               WRITE(j_id,FMT="(3x,a,' : ',D23.16)") cl(jn), zsum
             END DO
          ENDIF
@@ -201 +201 @@
       ! Loop over each sub-domain, i.e. number of processors ijsplt
       DO js = sind, eind
-         j_id = js - narea + 1
-         WRITE(numid_trc(j_id),*)clinfo
+         j_id = numid_trc(js - narea + 1)
+         WRITE(j_id,*)clinfo
       ENDDO
 
@@ -273 +273 @@
       DO js = sind, eind
          WRITE(clfile_out,FMT=clb_name) js-1
-         j_id = js -narea + 1
-         CALL ctlopn( numid_trc(j_id), clfile_out, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL',   &
+         j_id = numid_trc(js -narea + 1)
+         CALL ctlopn( j_id, clfile_out, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL',   &
             &         1, numout, .FALSE., 1 )
-         WRITE(numid_trc(j_id),*)
-         WRITE(numid_trc(j_id),*) '                 L O D Y C - I P S L'
-         WRITE(numid_trc(j_id),*) '                     O P A model'
-         WRITE(numid_trc(j_id),*) '            Ocean General Circulation Model'
-         WRITE(numid_trc(j_id),*) '               version OPA 9.0  (2005) '
-         WRITE(numid_trc(j_id),*)
-         WRITE(numid_trc(j_id),*) '                   PROC number: ', js
-         WRITE(numid_trc(j_id),*)
-         WRITE(numid_trc(j_id),FMT="(19x,a20)")cl_run
+         WRITE(j_id,*)
+         WRITE(j_id,*) '                 L O D Y C - I P S L'
+         WRITE(j_id,*) '                     O P A model'
+         WRITE(j_id,*) '            Ocean General Circulation Model'
+         WRITE(j_id,*) '               version OPA 9.0  (2005) '
+         WRITE(j_id,*)
+         WRITE(j_id,*) '                   PROC number: ', js
+         WRITE(j_id,*)
+         WRITE(j_id,FMT="(19x,a20)")cl_run
 
          ! Print the SUM control indices
@@ -300 +300 @@
             ENDIF
          ENDIF
-         WRITE(numid_trc(j_id),*) 
-         WRITE(numid_trc(j_id),*) 'prt_tra_ctl :  Sum control indices'
-         WRITE(numid_trc(j_id),*) '~~~~~~~'
-         WRITE(numid_trc(j_id),*)
-         WRITE(numid_trc(j_id),9000)'                                nlej   = ', nlejtl(js), '              '
-         WRITE(numid_trc(j_id),9000)'                  ------------- njctle = ', njctle, ' -------------'
-         WRITE(numid_trc(j_id),9001)'                  |                                       |'
-         WRITE(numid_trc(j_id),9001)'                  |                                       |'
-         WRITE(numid_trc(j_id),9001)'                  |                                       |'
-         WRITE(numid_trc(j_id),9002)'           nictls = ', nictls,  '                           nictle = ', nictle
-         WRITE(numid_trc(j_id),9002)'           nldi   = ', nlditl(js),  '                           nlei   = ', nleitl(js)
-         WRITE(numid_trc(j_id),9001)'                  |                                       |'
-         WRITE(numid_trc(j_id),9001)'                  |                                       |'
-         WRITE(numid_trc(j_id),9001)'                  |                                       |'
-         WRITE(numid_trc(j_id),9004)'  njmpp  = ',njmpptl(js),'   ------------- njctls = ', njctls, ' -------------'
-         WRITE(numid_trc(j_id),9003)'           nimpp  = ', nimpptl(js), '        nldj   = ', nldjtl(js), '              '
-         WRITE(numid_trc(j_id),*)
-         WRITE(numid_trc(j_id),*)
+         WRITE(j_id,*) 
+         WRITE(j_id,*) 'prt_tra_ctl :  Sum control indices'
+         WRITE(j_id,*) '~~~~~~~'
+         WRITE(j_id,*)
+         WRITE(j_id,9000)'                                nlej   = ', nlejtl(js), '              '
+         WRITE(j_id,9000)'                  ------------- njctle = ', njctle, ' -------------'
+         WRITE(j_id,9001)'                  |                                       |'
+         WRITE(j_id,9001)'                  |                                       |'
+         WRITE(j_id,9001)'                  |                                       |'
+         WRITE(j_id,9002)'           nictls = ', nictls,  '                           nictle = ', nictle
+         WRITE(j_id,9002)'           nldi   = ', nlditl(js),  '                           nlei   = ', nleitl(js)
+         WRITE(j_id,9001)'                  |                                       |'
+         WRITE(j_id,9001)'                  |                                       |'
+         WRITE(j_id,9001)'                  |                                       |'
+         WRITE(j_id,9004)'  njmpp  = ',njmpptl(js),'   ------------- njctls = ', njctls, ' -------------'
+         WRITE(j_id,9003)'           nimpp  = ', nimpptl(js), '        nldj   = ', nldjtl(js), '              '
+         WRITE(j_id,*)
+         WRITE(j_id,*)
 
 9000     FORMAT(a41,i4.4,a14)