Changeset 7646 for trunk/NEMOGCM/NEMO/OPA_SRC/LBC/mppini.F90
- Timestamp:
- 2017-02-06T10:25:03+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/LBC/mppini.F90
r6412 r7646 65 65 WRITE(numout,*) 66 66 WRITE(numout,*) 'mpp_init(2) : NO massively parallel processing' 67 WRITE(numout,*) '~~~~~~~~~~~ :'67 WRITE(numout,*) '~~~~~~~~~~~ ' 68 68 WRITE(numout,*) ' nperio = ', nperio 69 69 WRITE(numout,*) ' npolj = ', npolj … … 76 76 & 'the domain is lay out for distributed memory computing! ' ) 77 77 78 IF( jperio == 7 ) CALL ctl_stop( ' jperio = 7 needs distributed memory computing ', & 79 & ' with 1 process. Add key_mpp_mpi in the list of active cpp keys ' ) 78 80 END SUBROUTINE mpp_init 79 81 … … 265 267 266 268 IF(lwp) WRITE(numout,*) 267 IF(lwp) WRITE(numout,*) ' mpp_init: defines mpp subdomains' 268 IF(lwp) WRITE(numout,*) ' ~~~~~~ ----------------------' 269 IF(lwp) WRITE(numout,*) 270 IF(lwp) WRITE(numout,*) 'iresti=',iresti,' irestj=',irestj 271 IF(lwp) WRITE(numout,*) 272 IF(lwp) WRITE(numout,*) 'jpni=',jpni,' jpnj=',jpnj 269 IF(lwp) WRITE(numout,*) ' defines mpp subdomains' 270 IF(lwp) WRITE(numout,*) ' jpni=', jpni, ' iresti=', iresti 271 IF(lwp) WRITE(numout,*) ' jpnj=', jpnj, ' irestj=', irestj 273 272 zidom = nreci 274 273 DO ji = 1, jpni … … 276 275 END DO 277 276 IF(lwp) WRITE(numout,*) 278 IF(lwp) WRITE(numout,*)' sum ilcit(i,1)=', zidom, ' jpiglo=', jpiglo277 IF(lwp) WRITE(numout,*)' sum ilcit(i,1)=', zidom, ' jpiglo=', jpiglo 279 278 280 279 zjdom = nrecj … … 282 281 zjdom = zjdom + ilcjt(1,jj) - nrecj 283 282 END DO 284 IF(lwp) WRITE(numout,*)' sum ilcit(1,j)=', zjdom, ' jpjglo=', jpjglo 285 IF(lwp) WRITE(numout,*) 283 IF(lwp) WRITE(numout,*)' sum ilcit(1,j)=', zjdom, ' jpjglo=', jpjglo 286 284 287 285 IF(lwp) THEN … … 360 358 njmpp = njmppt(narea) 361 359 362 ! Save processor layout in layout.dat file363 IF(lwp) THEN364 CALL ctl_opn( inum, 'layout.dat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE., narea )365 WRITE(inum,'(a)') ' jpnij jpi jpj jpk jpiglo jpjglo'366 WRITE(inum,'(6i8)') jpnij,jpi,jpj,jpk,jpiglo,jpjglo367 WRITE(inum,'(a)') 'NAREA nlci nlcj nldi nldj nlei nlej nimpp njmpp'368 369 DOjn = 1, jpnij370 WRITE(inum,'(9i5)') jn, nlcit(jn), nlcjt(jn), &371 372 373 374 END DO375 CLOSE(inum)360 ! Save processor layout in layout.dat file 361 IF(lwp) THEN 362 CALL ctl_opn( inum, 'layout.dat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE., narea ) 363 WRITE(inum,'(a)') ' jpnij jpi jpj jpk jpiglo jpjglo' 364 WRITE(inum,'(6i8)') jpnij,jpi,jpj,jpk,jpiglo,jpjglo 365 WRITE(inum,'(a)') 'NAREA nlci nlcj nldi nldj nlei nlej nimpp njmpp' 366 ! 367 DO jn = 1, jpnij 368 WRITE(inum,'(9i5)') jn, nlcit(jn), nlcjt(jn), & 369 & nldit(jn), nldjt(jn), & 370 & nleit(jn), nlejt(jn), & 371 & nimppt(jn), njmppt(jn) 372 END DO 373 CLOSE(inum) 376 374 END IF 377 375 378 379 376 ! w a r n i n g narea (zone) /= nproc (processors)! 380 377 381 IF( jperio == 1 .OR. jperio == 4 .OR. jperio == 6 ) THEN378 IF( jperio == 1 .OR. jperio == 4 .OR. jperio == 6 .OR. jperio == 7 ) THEN 382 379 IF( jpni == 1 )THEN 383 380 nbondi = 2 … … 428 425 429 426 IF(lwp) THEN 430 WRITE(numout,*) ' nproc = ', nproc 431 WRITE(numout,*) ' nowe = ', nowe , ' noea = ', noea 432 WRITE(numout,*) ' nono = ', nono , ' noso = ', noso 433 WRITE(numout,*) ' nbondi = ', nbondi 434 WRITE(numout,*) ' nbondj = ', nbondj 435 WRITE(numout,*) ' npolj = ', npolj 436 WRITE(numout,*) ' nperio = ', nperio 437 WRITE(numout,*) ' nlci = ', nlci 438 WRITE(numout,*) ' nlcj = ', nlcj 439 WRITE(numout,*) ' nimpp = ', nimpp 440 WRITE(numout,*) ' njmpp = ', njmpp 441 WRITE(numout,*) ' nreci = ', nreci , ' npse = ', npse 442 WRITE(numout,*) ' nrecj = ', nrecj , ' npsw = ', npsw 443 WRITE(numout,*) ' jpreci = ', jpreci , ' npne = ', npne 444 WRITE(numout,*) ' jprecj = ', jprecj , ' npnw = ', npnw 427 WRITE(numout,*) ' nproc = ', nproc 428 WRITE(numout,*) ' nowe = ', nowe , ' noea = ', noea 429 WRITE(numout,*) ' nono = ', nono , ' noso = ', noso 430 WRITE(numout,*) ' nbondi = ', nbondi, ' nbondj = ', nbondj 431 WRITE(numout,*) ' npolj = ', npolj 432 WRITE(numout,*) ' nperio = ', nperio 433 WRITE(numout,*) ' nlci = ', nlci , ' nlcj = ', nlcj 434 WRITE(numout,*) ' nimpp = ', nimpp , ' njmpp = ', njmpp 435 WRITE(numout,*) ' nreci = ', nreci , ' npse = ', npse 436 WRITE(numout,*) ' nrecj = ', nrecj , ' npsw = ', npsw 437 WRITE(numout,*) ' jpreci = ', jpreci, ' npne = ', npne 438 WRITE(numout,*) ' jprecj = ', jprecj, ' npnw = ', npnw 445 439 WRITE(numout,*) 446 440 ENDIF 447 441 442 IF( jperio == 7 .AND. ( jpni /= 1 .OR. jpnj /= 1 ) ) & 443 & CALL ctl_stop( ' mpp_init: error jperio = 7 works only with jpni = jpnj = 1' ) 448 444 IF( nperio == 1 .AND. jpni /= 1 ) CALL ctl_stop( ' mpp_init: error on cyclicity' ) 449 445
Note: See TracChangeset
for help on using the changeset viewer.