Changeset 627
- Timestamp:
- 2007-02-27T20:02:43+01:00 (17 years ago)
- Location:
- trunk/NEMO
- Files:
-
- 9 edited
-
OFF_SRC/IOM/in_out_manager.F90 (modified) (2 diffs)
-
OFF_SRC/ctlopn.f90 (modified) (2 diffs)
-
OFF_SRC/mppini.F90 (modified) (1 diff)
-
OFF_SRC/mppini_2.h90 (modified) (1 diff)
-
OFF_SRC/opa.F90 (modified) (2 diffs)
-
OFF_SRC/prtctl.F90 (modified) (10 diffs)
-
OPA_SRC/opa.F90 (modified) (2 diffs)
-
OPA_SRC/prtctl.F90 (modified) (6 diffs)
-
TOP_SRC/prtctl_trc.F90 (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OFF_SRC/IOM/in_out_manager.F90
r569 r627 56 56 !! logical units 57 57 !!---------------------------------------------------------------------- 58 INTEGER :: numstp = 1!: logical unit for time step59 INTEGER :: numout = 2!: logical unit for output print60 INTEGER :: numnam = 3!: logical unit for namelist61 INTEGER :: numnam_ice = 4!: logical unit for ice namelist62 INTEGER :: numevo_ice = 17!: logical unit for ice variables (temp. evolution)63 INTEGER :: numsol = 25!: logical unit for solver statistics64 INTEGER :: numwri = 40!: logical unit for output write65 INTEGER :: numisp = 41!: logical unit for island statistics66 INTEGER :: numgap = 45!: logical unit for differences diagnostic67 INTEGER :: numbol = 67!: logical unit for "bol" diagnostics68 INTEGER :: numptr = 68!: logical unit for Poleward TRansports69 INTEGER :: numflo = 69!: logical unit for drifting floats58 INTEGER :: numstp !: logical unit for time step 59 INTEGER :: numout !: logical unit for output print 60 INTEGER :: numnam !: logical unit for namelist 61 INTEGER :: numnam_ice !: logical unit for ice namelist 62 INTEGER :: numevo_ice !: logical unit for ice variables (temp. evolution) 63 INTEGER :: numsol !: logical unit for solver statistics 64 INTEGER :: numwri !: logical unit for output write 65 INTEGER :: numisp !: logical unit for island statistics 66 INTEGER :: numgap !: logical unit for differences diagnostic 67 INTEGER :: numbol !: logical unit for "bol" diagnostics 68 INTEGER :: numptr !: logical unit for Poleward TRansports 69 INTEGER :: numflo !: logical unit for drifting floats 70 70 71 71 !!---------------------------------------------------------------------- … … 157 157 LOGICAL :: llopn 158 158 !--------------------------------------------------------------------- 159 getunit = 50! choose a unit that is big enough then it is159 getunit = 15 ! choose a unit that is big enough then it is 160 160 ! not already used in OPA 161 161 llopn = .TRUE. -
trunk/NEMO/OFF_SRC/ctlopn.f90
r325 r627 12 12 !! 8.5 ! 02-06 (G. Madec) F90: Free form and module 13 13 !!---------------------------------------------------------------------- 14 USE in_out_manager 15 14 16 IMPLICIT NONE 15 17 18 INTEGER, INTENT( out ) :: & 19 knum ! logical unit to open 16 20 INTEGER, INTENT( in ) :: & 17 knum, & ! logical unit to open18 21 krequ, & ! =1 file required (stop if not exist) 19 22 ! ! =0 file not required (create the file if does not exist) … … 36 39 !!---------------------------------------------------------------------- 37 40 41 knum=getunit() 38 42 39 43 ! 1. Required file -
trunk/NEMO/OFF_SRC/mppini.F90
r482 r627 348 348 ! Save processor layout in layout.dat file 349 349 IF (lwp) THEN 350 inum = 11 351 352 OPEN(inum,FILE='layout.dat') 350 CALL ctlopn( inum, 'layout.dat', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', & 351 & 1, numout, .FALSE., 1 ) 353 352 WRITE(inum,'(6i8)') jpnij,jpi,jpj,jpk,jpiglo,jpjglo 354 353 WRITE(inum,'(a)') 'NAREA nlci nlcj nldi nldj nlei nlej nimpp njmpp' -
trunk/NEMO/OFF_SRC/mppini_2.h90
r495 r627 466 466 ! Save processor layout in ascii file 467 467 IF (lwp) THEN 468 inum = 11 ! how do we know that 11 is ok???469 OPEN(inum,FILE='layout.dat')468 CALL ctlopn( inum, 'layout.dat', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', & 469 & 1, numout, .FALSE., 1 ) 470 470 WRITE(inum,'(6i8)') jpnij,jpi,jpj,jpk,jpiglo,jpjglo 471 471 WRITE(inum,'(a)') 'NAREA nlci nlcj nldi nldj nlei nlej nimpp njmpp' -
trunk/NEMO/OFF_SRC/opa.F90
r496 r627 92 92 93 93 ! open listing and namelist units 94 IF ( numout /= 0 .AND. numout /= 6 ) THEN 95 CALL ctlopn( numout, file_out, 'UNKNOWN', 'FORMATTED', & 96 & 'SEQUENTIAL', 1, numout, .FALSE., 1 ) 97 ENDIF 94 CALL ctlopn( numout, file_out, 'UNKNOWN', 'FORMATTED', & 95 & 'SEQUENTIAL', 1, numout, .FALSE., 1 ) 98 96 99 97 namelistname = 'namelist' 100 98 CALL ctlopn( numnam, namelistname, 'OLD', 'FORMATTED', 'SEQUENTIAL', & 101 & 99 & 1, numout, .FALSE., 1 ) 102 100 103 101 WRITE(numout,*) … … 117 115 CLOSE( numout ) ! standard model output file 118 116 WRITE(file_out,FMT="('ocean.output_',I4.4)") narea-1 119 IF ( numout /= 0 .AND. numout /= 6 ) THEN 120 CALL ctlopn( numout, file_out, 'UNKNOWN', 'FORMATTED', & 121 & 'SEQUENTIAL', 1, numout, .FALSE., 1 ) 122 ENDIF 117 CALL ctlopn( numout, file_out, 'UNKNOWN', 'FORMATTED', & 118 & 'SEQUENTIAL', 1, numout, .FALSE., 1 ) 123 119 ! 124 120 WRITE(numout,*) -
trunk/NEMO/OFF_SRC/prtctl.F90
r568 r627 12 12 13 13 !! * Module declaration 14 INTEGER, DIMENSION(:), ALLOCATABLE :: numid 14 15 INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: & !: 15 16 nlditl , nldjtl , & !: first, last indoor index for each i-domain … … 91 92 92 93 !! * Local declarations 93 INTEGER :: overlap, numid, jn, sind, eind, kdir94 INTEGER :: overlap, jn, sind, eind, kdir,j_id 94 95 CHARACTER (len=15) :: cl2 95 96 REAL(wp) :: zsum1, zsum2, zvctl1, zvctl2 … … 136 137 ! Loop over each sub-domain, i.e. the total number of processors ijsplt 137 138 DO jn = sind, eind 138 139 numid = 90 + jn 140 139 ! Set logical unit 140 j_id = numid(jn - narea + 1) 141 141 ! Set indices for the SUM control 142 142 IF( .NOT. lsp_area ) THEN … … 187 187 ! Print the result 188 188 IF( PRESENT(clinfo3) ) THEN 189 WRITE( numid,FMT='(a,D23.16,3x,a,D23.16)')clinfo1, zsum1-zvctl1, cl2, zsum2-zvctl2189 WRITE(j_id,FMT='(a,D23.16,3x,a,D23.16)')clinfo1, zsum1-zvctl1, cl2, zsum2-zvctl2 190 190 SELECT CASE( clinfo3 ) 191 191 CASE ( 'tra-ta' ) … … 199 199 END SELECT 200 200 ELSEIF ( PRESENT(clinfo2) .OR. PRESENT(tab2d_2) .OR. PRESENT(tab3d_2) ) THEN 201 WRITE( numid,FMT='(a,D23.16,3x,a,D23.16)')clinfo1, zsum1, cl2, zsum2201 WRITE(j_id,FMT='(a,D23.16,3x,a,D23.16)')clinfo1, zsum1, cl2, zsum2 202 202 ELSE 203 WRITE( numid,FMT='(a,D23.16)')clinfo1, zsum1203 WRITE(j_id,FMT='(a,D23.16)')clinfo1, zsum1 204 204 ENDIF 205 205 … … 232 232 233 233 !! * Local declarations 234 INTEGER :: numid, jn, sind, eind, iltime234 INTEGER :: jn, sind, eind, iltime, j_id 235 235 !!---------------------------------------------------------------------- 236 236 … … 257 257 ! Loop over each sub-domain, i.e. number of processors ijsplt 258 258 DO jn = sind, eind 259 260 numid = 90 + jn 259 260 ! Set logical unit 261 j_id = numid(jn - narea + 1) 261 262 262 263 IF( PRESENT(ivar1) .AND. PRESENT(clinfo2) .AND. PRESENT(ivar2) ) THEN 263 WRITE( numid,*)clinfo1, ivar1, clinfo2, ivar2264 WRITE(j_id,*)clinfo1, ivar1, clinfo2, ivar2 264 265 ELSEIF ( PRESENT(ivar1) .AND. PRESENT(clinfo2) .AND. .NOT. PRESENT(ivar2) ) THEN 265 WRITE( numid,*)clinfo1, ivar1, clinfo2266 WRITE(j_id,*)clinfo1, ivar1, clinfo2 266 267 ELSEIF ( PRESENT(ivar1) .AND. .NOT. PRESENT(clinfo2) .AND. PRESENT(ivar2) ) THEN 267 WRITE( numid,*)clinfo1, ivar1, ivar2268 WRITE(j_id,*)clinfo1, ivar1, ivar2 268 269 ELSEIF ( PRESENT(ivar1) .AND. .NOT. PRESENT(clinfo2) .AND. .NOT. PRESENT(ivar2) ) THEN 269 WRITE( numid,*)clinfo1, ivar1270 WRITE(j_id,*)clinfo1, ivar1 270 271 ELSE 271 WRITE( numid,*)clinfo1272 WRITE(j_id,*)clinfo1 272 273 ENDIF 273 274 … … 287 288 !!---------------------------------------------------------------------- 288 289 !! * Local declarations 289 INTEGER :: jn, numid, sind, eind290 INTEGER :: jn, sind, eind, j_id 290 291 CHARACTER (len=28) :: clfile_out 291 292 CHARACTER (len=23) :: clb_name … … 344 345 ENDIF 345 346 347 ALLOCATE(numid(eind-sind+1)) 348 346 349 DO jn = sind, eind 347 numid = 90 + jn348 350 WRITE(clfile_out,FMT=clb_name) jn-1 349 OPEN ( UNIT=numid, FILE=TRIM(clfile_out),FORM='FORMATTED' ) 350 WRITE(numid,*) 351 WRITE(numid,*) ' L O D Y C - I P S L' 352 WRITE(numid,*) ' O P A model' 353 WRITE(numid,*) ' Ocean General Circulation Model' 354 WRITE(numid,*) ' version OPA 9.0 (2005) ' 355 WRITE(numid,*) 356 WRITE(numid,*) ' PROC number: ', jn 357 WRITE(numid,*) 358 WRITE(numid,FMT="(19x,a20)")cl_run 351 j_id = numid(jn -narea + 1) 352 CALL ctlopn( j_id, clfile_out, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', & 353 & 1, numout, .FALSE., 1 ) 354 WRITE(j_id,*) 355 WRITE(j_id,*) ' L O D Y C - I P S L' 356 WRITE(j_id,*) ' O P A model' 357 WRITE(j_id,*) ' Ocean General Circulation Model' 358 WRITE(j_id,*) ' version OPA 9.0 (2005) ' 359 WRITE(j_id,*) 360 WRITE(j_id,*) ' PROC number: ', jn 361 WRITE(j_id,*) 362 WRITE(j_id,FMT="(19x,a20)")cl_run 359 363 360 364 ! Print the SUM control indices … … 365 369 njctle = njmpptl(jn) + nlejtl(jn) - 1 366 370 ENDIF 367 WRITE( numid,*)368 WRITE( numid,*) 'prt_ctl : Sum control indices'369 WRITE( numid,*) '~~~~~~~'370 WRITE( numid,*)371 WRITE( numid,9000)' nlej = ', nlejtl(jn), ' '372 WRITE( numid,9000)' ------------- njctle = ', njctle, ' -------------'373 WRITE( numid,9001)' | |'374 WRITE( numid,9001)' | |'375 WRITE( numid,9001)' | |'376 WRITE( numid,9002)' nictls = ', nictls, ' nictle = ', nictle377 WRITE( numid,9002)' nldi = ', nlditl(jn), ' nlei = ', nleitl(jn)378 WRITE( numid,9001)' | |'379 WRITE( numid,9001)' | |'380 WRITE( numid,9001)' | |'381 WRITE( numid,9004)' njmpp = ',njmpptl(jn),' ------------- njctls = ', njctls, ' -------------'382 WRITE( numid,9003)' nimpp = ', nimpptl(jn), ' nldj = ', nldjtl(jn), ' '383 WRITE( numid,*)384 WRITE( numid,*)371 WRITE(j_id,*) 372 WRITE(j_id,*) 'prt_ctl : Sum control indices' 373 WRITE(j_id,*) '~~~~~~~' 374 WRITE(j_id,*) 375 WRITE(j_id,9000)' nlej = ', nlejtl(jn), ' ' 376 WRITE(j_id,9000)' ------------- njctle = ', njctle, ' -------------' 377 WRITE(j_id,9001)' | |' 378 WRITE(j_id,9001)' | |' 379 WRITE(j_id,9001)' | |' 380 WRITE(j_id,9002)' nictls = ', nictls, ' nictle = ', nictle 381 WRITE(j_id,9002)' nldi = ', nlditl(jn), ' nlei = ', nleitl(jn) 382 WRITE(j_id,9001)' | |' 383 WRITE(j_id,9001)' | |' 384 WRITE(j_id,9001)' | |' 385 WRITE(j_id,9004)' njmpp = ',njmpptl(jn),' ------------- njctls = ', njctls, ' -------------' 386 WRITE(j_id,9003)' nimpp = ', nimpptl(jn), ' nldj = ', nldjtl(jn), ' ' 387 WRITE(j_id,*) 388 WRITE(j_id,*) 385 389 386 390 9000 FORMAT(a41,i4.4,a14) -
trunk/NEMO/OPA_SRC/opa.F90
r624 r627 203 203 namelistname = 'namelist' 204 204 CALL ctlopn( numnam, namelistname, 'OLD', 'FORMATTED', 'SEQUENTIAL', & 205 & 205 & 1, numout, .FALSE., 1 ) 206 206 207 207 ! Namelist namctl : Control prints & Benchmark … … 224 224 IF( narea-1 > 0 ) THEN 225 225 WRITE(file_out,FMT="('ocean.output_',I4.4)") narea-1 226 IF( numout /= 0 .AND. numout /= 6 ) THEN 227 CALL ctlopn( numout, file_out, 'UNKNOWN', 'FORMATTED', & 228 & 'SEQUENTIAL', 1, numout, .FALSE., 1 ) 229 ENDIF 226 CALL ctlopn( numout, file_out, 'UNKNOWN', 'FORMATTED', & 227 & 'SEQUENTIAL', 1, numout, .FALSE., 1 ) 230 228 lwp = .TRUE. 231 229 ! -
trunk/NEMO/OPA_SRC/prtctl.F90
r624 r627 138 138 DO jn = sind, eind 139 139 ! Set logical unit 140 j_id = jn - narea + 1140 j_id = numid(jn - narea + 1) 141 141 ! Set indices for the SUM control 142 142 IF( .NOT. lsp_area ) THEN … … 187 187 ! Print the result 188 188 IF( PRESENT(clinfo3) ) THEN 189 WRITE( numid(j_id),FMT='(a,D23.16,3x,a,D23.16)')clinfo1, zsum1-zvctl1, cl2, zsum2-zvctl2189 WRITE(j_id,FMT='(a,D23.16,3x,a,D23.16)')clinfo1, zsum1-zvctl1, cl2, zsum2-zvctl2 190 190 SELECT CASE( clinfo3 ) 191 191 CASE ( 'tra-ta' ) … … 199 199 END SELECT 200 200 ELSEIF ( PRESENT(clinfo2) .OR. PRESENT(tab2d_2) .OR. PRESENT(tab3d_2) ) THEN 201 WRITE( numid(j_id),FMT='(a,D23.16,3x,a,D23.16)')clinfo1, zsum1, cl2, zsum2201 WRITE(j_id,FMT='(a,D23.16,3x,a,D23.16)')clinfo1, zsum1, cl2, zsum2 202 202 ELSE 203 WRITE( numid(j_id),FMT='(a,D23.16)')clinfo1, zsum1203 WRITE(j_id,FMT='(a,D23.16)')clinfo1, zsum1 204 204 ENDIF 205 205 … … 259 259 260 260 ! Set logical unit 261 j_id = jn - narea +1261 j_id = numid(jn - narea + 1) 262 262 263 263 IF( PRESENT(ivar1) .AND. PRESENT(clinfo2) .AND. PRESENT(ivar2) ) THEN 264 WRITE( numid(j_id),*)clinfo1, ivar1, clinfo2, ivar2264 WRITE(j_id,*)clinfo1, ivar1, clinfo2, ivar2 265 265 ELSEIF ( PRESENT(ivar1) .AND. PRESENT(clinfo2) .AND. .NOT. PRESENT(ivar2) ) THEN 266 WRITE( numid(j_id),*)clinfo1, ivar1, clinfo2266 WRITE(j_id,*)clinfo1, ivar1, clinfo2 267 267 ELSEIF ( PRESENT(ivar1) .AND. .NOT. PRESENT(clinfo2) .AND. PRESENT(ivar2) ) THEN 268 WRITE( numid(jn-narea+1),*)clinfo1, ivar1, ivar2268 WRITE(j_id,*)clinfo1, ivar1, ivar2 269 269 ELSEIF ( PRESENT(ivar1) .AND. .NOT. PRESENT(clinfo2) .AND. .NOT. PRESENT(ivar2) ) THEN 270 WRITE( numid(j_id),*)clinfo1, ivar1270 WRITE(j_id,*)clinfo1, ivar1 271 271 ELSE 272 WRITE( numid(j_id),*)clinfo1272 WRITE(j_id,*)clinfo1 273 273 ENDIF 274 274 … … 349 349 DO jn = sind, eind 350 350 WRITE(clfile_out,FMT=clb_name) jn-1 351 j_id = jn -narea + 1352 CALL ctlopn( numid(j_id), clfile_out, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', &351 j_id = numid(jn -narea + 1) 352 CALL ctlopn( j_id, clfile_out, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', & 353 353 & 1, numout, .FALSE., 1 ) 354 WRITE( numid(j_id),*)355 WRITE( numid(j_id),*) ' L O D Y C - I P S L'356 WRITE( numid(j_id),*) ' O P A model'357 WRITE( numid(j_id),*) ' Ocean General Circulation Model'358 WRITE( numid(j_id),*) ' version OPA 9.0 (2005) '359 WRITE( numid(j_id),*)360 WRITE( numid(j_id),*) ' PROC number: ', jn361 WRITE( numid(j_id),*)362 WRITE( numid(j_id),FMT="(19x,a20)")cl_run354 WRITE(j_id,*) 355 WRITE(j_id,*) ' L O D Y C - I P S L' 356 WRITE(j_id,*) ' O P A model' 357 WRITE(j_id,*) ' Ocean General Circulation Model' 358 WRITE(j_id,*) ' version OPA 9.0 (2005) ' 359 WRITE(j_id,*) 360 WRITE(j_id,*) ' PROC number: ', jn 361 WRITE(j_id,*) 362 WRITE(j_id,FMT="(19x,a20)")cl_run 363 363 364 364 ! Print the SUM control indices … … 369 369 njctle = njmpptl(jn) + nlejtl(jn) - 1 370 370 ENDIF 371 WRITE( numid(j_id),*)372 WRITE( numid(j_id),*) 'prt_ctl : Sum control indices'373 WRITE( numid(j_id),*) '~~~~~~~'374 WRITE( numid(j_id),*)375 WRITE( numid(j_id),9000)' nlej = ', nlejtl(jn), ' '376 WRITE( numid(j_id),9000)' ------------- njctle = ', njctle, ' -------------'377 WRITE( numid(j_id),9001)' | |'378 WRITE( numid(j_id),9001)' | |'379 WRITE( numid(j_id),9001)' | |'380 WRITE( numid(j_id),9002)' nictls = ', nictls, ' nictle = ', nictle381 WRITE( numid(j_id),9002)' nldi = ', nlditl(jn), ' nlei = ', nleitl(jn)382 WRITE( numid(j_id),9001)' | |'383 WRITE( numid(j_id),9001)' | |'384 WRITE( numid(j_id),9001)' | |'385 WRITE( numid(j_id),9004)' njmpp = ',njmpptl(jn),' ------------- njctls = ', njctls, ' -------------'386 WRITE( numid(j_id),9003)' nimpp = ', nimpptl(jn), ' nldj = ', nldjtl(jn), ' '387 WRITE( numid(j_id),*)388 WRITE( numid(j_id),*)371 WRITE(j_id,*) 372 WRITE(j_id,*) 'prt_ctl : Sum control indices' 373 WRITE(j_id,*) '~~~~~~~' 374 WRITE(j_id,*) 375 WRITE(j_id,9000)' nlej = ', nlejtl(jn), ' ' 376 WRITE(j_id,9000)' ------------- njctle = ', njctle, ' -------------' 377 WRITE(j_id,9001)' | |' 378 WRITE(j_id,9001)' | |' 379 WRITE(j_id,9001)' | |' 380 WRITE(j_id,9002)' nictls = ', nictls, ' nictle = ', nictle 381 WRITE(j_id,9002)' nldi = ', nlditl(jn), ' nlei = ', nleitl(jn) 382 WRITE(j_id,9001)' | |' 383 WRITE(j_id,9001)' | |' 384 WRITE(j_id,9001)' | |' 385 WRITE(j_id,9004)' njmpp = ',njmpptl(jn),' ------------- njctls = ', njctls, ' -------------' 386 WRITE(j_id,9003)' nimpp = ', nimpptl(jn), ' nldj = ', nldjtl(jn), ' ' 387 WRITE(j_id,*) 388 WRITE(j_id,*) 389 389 390 390 9000 FORMAT(a41,i4.4,a14) -
trunk/NEMO/TOP_SRC/prtctl_trc.F90
r625 r627 121 121 122 122 ! Set logical unit 123 j_id = js - narea + 1123 j_id = numid_trc(js - narea + 1) 124 124 ! Set indices for the SUM control 125 125 IF( .NOT. lsp_area ) THEN … … 153 153 zsum = SUM( ztab3d(nictls:nictle,njctls:njctle,1:kdir) & 154 154 & *zmask(nictls:nictle,njctls:njctle,1:kdir) ) 155 WRITE( numid_trc(j_id),FMT="(3x,a,' : ',D23.16)") cl(jn), zsum-zvctl155 WRITE(j_id,FMT="(3x,a,' : ',D23.16)") cl(jn), zsum-zvctl 156 156 tra_ctl(jn,js) = zsum 157 157 ENDDO … … 161 161 zsum = SUM( ztab3d(nictls:nictle,njctls:njctle,1:kdir) & 162 162 & * zmask(nictls:nictle,njctls:njctle,1:kdir) ) 163 WRITE( numid_trc(j_id),FMT="(3x,a,' : ',D23.16)") cl(jn), zsum163 WRITE(j_id,FMT="(3x,a,' : ',D23.16)") cl(jn), zsum 164 164 END DO 165 165 ENDIF … … 201 201 ! Loop over each sub-domain, i.e. number of processors ijsplt 202 202 DO js = sind, eind 203 j_id = js - narea + 1204 WRITE( numid_trc(j_id),*)clinfo203 j_id = numid_trc(js - narea + 1) 204 WRITE(j_id,*)clinfo 205 205 ENDDO 206 206 … … 273 273 DO js = sind, eind 274 274 WRITE(clfile_out,FMT=clb_name) js-1 275 j_id = js -narea + 1276 CALL ctlopn( numid_trc(j_id), clfile_out, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', &275 j_id = numid_trc(js -narea + 1) 276 CALL ctlopn( j_id, clfile_out, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', & 277 277 & 1, numout, .FALSE., 1 ) 278 WRITE( numid_trc(j_id),*)279 WRITE( numid_trc(j_id),*) ' L O D Y C - I P S L'280 WRITE( numid_trc(j_id),*) ' O P A model'281 WRITE( numid_trc(j_id),*) ' Ocean General Circulation Model'282 WRITE( numid_trc(j_id),*) ' version OPA 9.0 (2005) '283 WRITE( numid_trc(j_id),*)284 WRITE( numid_trc(j_id),*) ' PROC number: ', js285 WRITE( numid_trc(j_id),*)286 WRITE( numid_trc(j_id),FMT="(19x,a20)")cl_run278 WRITE(j_id,*) 279 WRITE(j_id,*) ' L O D Y C - I P S L' 280 WRITE(j_id,*) ' O P A model' 281 WRITE(j_id,*) ' Ocean General Circulation Model' 282 WRITE(j_id,*) ' version OPA 9.0 (2005) ' 283 WRITE(j_id,*) 284 WRITE(j_id,*) ' PROC number: ', js 285 WRITE(j_id,*) 286 WRITE(j_id,FMT="(19x,a20)")cl_run 287 287 288 288 ! Print the SUM control indices … … 300 300 ENDIF 301 301 ENDIF 302 WRITE( numid_trc(j_id),*)303 WRITE( numid_trc(j_id),*) 'prt_tra_ctl : Sum control indices'304 WRITE( numid_trc(j_id),*) '~~~~~~~'305 WRITE( numid_trc(j_id),*)306 WRITE( numid_trc(j_id),9000)' nlej = ', nlejtl(js), ' '307 WRITE( numid_trc(j_id),9000)' ------------- njctle = ', njctle, ' -------------'308 WRITE( numid_trc(j_id),9001)' | |'309 WRITE( numid_trc(j_id),9001)' | |'310 WRITE( numid_trc(j_id),9001)' | |'311 WRITE( numid_trc(j_id),9002)' nictls = ', nictls, ' nictle = ', nictle312 WRITE( numid_trc(j_id),9002)' nldi = ', nlditl(js), ' nlei = ', nleitl(js)313 WRITE( numid_trc(j_id),9001)' | |'314 WRITE( numid_trc(j_id),9001)' | |'315 WRITE( numid_trc(j_id),9001)' | |'316 WRITE( numid_trc(j_id),9004)' njmpp = ',njmpptl(js),' ------------- njctls = ', njctls, ' -------------'317 WRITE( numid_trc(j_id),9003)' nimpp = ', nimpptl(js), ' nldj = ', nldjtl(js), ' '318 WRITE( numid_trc(j_id),*)319 WRITE( numid_trc(j_id),*)302 WRITE(j_id,*) 303 WRITE(j_id,*) 'prt_tra_ctl : Sum control indices' 304 WRITE(j_id,*) '~~~~~~~' 305 WRITE(j_id,*) 306 WRITE(j_id,9000)' nlej = ', nlejtl(js), ' ' 307 WRITE(j_id,9000)' ------------- njctle = ', njctle, ' -------------' 308 WRITE(j_id,9001)' | |' 309 WRITE(j_id,9001)' | |' 310 WRITE(j_id,9001)' | |' 311 WRITE(j_id,9002)' nictls = ', nictls, ' nictle = ', nictle 312 WRITE(j_id,9002)' nldi = ', nlditl(js), ' nlei = ', nleitl(js) 313 WRITE(j_id,9001)' | |' 314 WRITE(j_id,9001)' | |' 315 WRITE(j_id,9001)' | |' 316 WRITE(j_id,9004)' njmpp = ',njmpptl(js),' ------------- njctls = ', njctls, ' -------------' 317 WRITE(j_id,9003)' nimpp = ', nimpptl(js), ' nldj = ', nldjtl(js), ' ' 318 WRITE(j_id,*) 319 WRITE(j_id,*) 320 320 321 321 9000 FORMAT(a41,i4.4,a14)
Note: See TracChangeset
for help on using the changeset viewer.
