Changeset 625 for trunk/NEMO/TOP_SRC/prtctl_trc.F90
- Timestamp:
- 2007-02-27T16:51:00+01:00 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/prtctl_trc.F90
r426 r625 15 15 16 16 !! * Module declaration 17 INTEGER, DIMENSION(:), ALLOCATABLE :: numid_trc ! logical unit 17 18 INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: & !: 18 19 nlditl , nldjtl , & !: first, last indoor index for each i-domain … … 81 82 82 83 !! * Local declarations 83 INTEGER :: overlap, numid, jn, js, sind, eind, kdir84 INTEGER :: overlap, jn, js, sind, eind, kdir, j_id 84 85 REAL(wp) :: zsum, zvctl 85 86 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zmask, ztab3d … … 119 120 DO js = sind, eind 120 121 121 numid = 90 + js122 122 ! Set logical unit 123 j_id = js - narea + 1 123 124 ! Set indices for the SUM control 124 125 IF( .NOT. lsp_area ) THEN … … 152 153 zsum = SUM( ztab3d(nictls:nictle,njctls:njctle,1:kdir) & 153 154 & *zmask(nictls:nictle,njctls:njctle,1:kdir) ) 154 WRITE(numid ,FMT="(3x,a,' : ',D23.16)") cl(jn), zsum-zvctl155 WRITE(numid_trc(j_id),FMT="(3x,a,' : ',D23.16)") cl(jn), zsum-zvctl 155 156 tra_ctl(jn,js) = zsum 156 157 ENDDO … … 160 161 zsum = SUM( ztab3d(nictls:nictle,njctls:njctle,1:kdir) & 161 162 & * zmask(nictls:nictle,njctls:njctle,1:kdir) ) 162 WRITE(numid ,FMT="(3x,a,' : ',D23.16)") cl(jn), zsum163 WRITE(numid_trc(j_id),FMT="(3x,a,' : ',D23.16)") cl(jn), zsum 163 164 END DO 164 165 ENDIF … … 185 186 186 187 !! * Local declarations 187 INTEGER :: numid, js, sind, eind188 INTEGER :: js, sind, eind, j_id 188 189 !!---------------------------------------------------------------------- 189 190 … … 200 201 ! Loop over each sub-domain, i.e. number of processors ijsplt 201 202 DO js = sind, eind 202 numid = 90 + js203 WRITE(numid ,*)clinfo203 j_id = js - narea + 1 204 WRITE(numid_trc(j_id),*)clinfo 204 205 ENDDO 205 206 … … 218 219 !!---------------------------------------------------------------------- 219 220 !! * Local declarations 220 INTEGER :: js, numid, sind, eind221 INTEGER :: js, sind, eind, j_id 221 222 CHARACTER (len=31) :: clfile_out 222 223 CHARACTER (len=27) :: clb_name … … 268 269 ENDIF 269 270 271 ALLOCATE(numid_trc(eind-sind+1)) 272 270 273 DO js = sind, eind 271 numid = 90 + js272 274 WRITE(clfile_out,FMT=clb_name) js-1 273 OPEN ( UNIT=numid, FILE=TRIM(clfile_out),FORM='FORMATTED' ) 274 WRITE(numid,*) 275 WRITE(numid,*) ' L O D Y C - I P S L' 276 WRITE(numid,*) ' O P A model' 277 WRITE(numid,*) ' Ocean General Circulation Model' 278 WRITE(numid,*) ' version OPA 9.0 (2005) ' 279 WRITE(numid,*) 280 WRITE(numid,*) ' PROC number: ', js 281 WRITE(numid,*) 282 WRITE(numid,FMT="(19x,a20)")cl_run 275 j_id = js -narea + 1 276 CALL ctlopn( numid_trc(j_id), clfile_out, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', & 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: ', js 285 WRITE(numid_trc(j_id),*) 286 WRITE(numid_trc(j_id),FMT="(19x,a20)")cl_run 283 287 284 288 ! Print the SUM control indices … … 296 300 ENDIF 297 301 ENDIF 298 WRITE(numid ,*)299 WRITE(numid ,*) 'prt_tra_ctl : Sum control indices'300 WRITE(numid ,*) '~~~~~~~'301 WRITE(numid ,*)302 WRITE(numid ,9000)' nlej = ', nlejtl(js), ' '303 WRITE(numid ,9000)' ------------- njctle = ', njctle, ' -------------'304 WRITE(numid ,9001)' | |'305 WRITE(numid ,9001)' | |'306 WRITE(numid ,9001)' | |'307 WRITE(numid ,9002)' nictls = ', nictls, ' nictle = ', nictle308 WRITE(numid ,9002)' nldi = ', nlditl(js), ' nlei = ', nleitl(js)309 WRITE(numid ,9001)' | |'310 WRITE(numid ,9001)' | |'311 WRITE(numid ,9001)' | |'312 WRITE(numid ,9004)' njmpp = ',njmpptl(js),' ------------- njctls = ', njctls, ' -------------'313 WRITE(numid ,9003)' nimpp = ', nimpptl(js), ' nldj = ', nldjtl(js), ' '314 WRITE(numid ,*)315 WRITE(numid ,*)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 = ', nictle 312 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),*) 316 320 317 321 9000 FORMAT(a41,i4.4,a14)
Note: See TracChangeset
for help on using the changeset viewer.