Changeset 14275 for NEMO/trunk/src/OCE/OBS/obs_grid.F90
- Timestamp:
- 2021-01-07T13:13:16+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/OBS/obs_grid.F90
r13286 r14275 130 130 CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 131 131 & 1, jpi, 1, jpj, & 132 & n proc, jpnij,&132 & narea-1, jpnij, & 133 133 & glamt, gphit, tmask, & 134 134 & kobsin, plam, pphi, & … … 137 137 CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 138 138 & 1, jpi, 1, jpj, & 139 & n proc, jpnij,&139 & narea-1, jpnij, & 140 140 & glamu, gphiu, umask, & 141 141 & kobsin, plam, pphi, & … … 144 144 CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 145 145 & 1, jpi, 1, jpj, & 146 & n proc, jpnij,&146 & narea-1, jpnij, & 147 147 & glamv, gphiv, vmask, & 148 148 & kobsin, plam, pphi, & … … 151 151 CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 152 152 & 1, jpi, 1, jpj, & 153 & n proc, jpnij,&153 & narea-1, jpnij, & 154 154 & glamf, gphif, fmask, & 155 155 & kobsin, plam, pphi, & … … 176 176 !! 177 177 !! ** Action : Return kproc holding the observation and kiobsi,kobsj 178 !! valid on kproc=n procprocessor only.178 !! valid on kproc=narea-1 processor only. 179 179 !! 180 180 !! History : … … 248 248 jlon = jpiglo 249 249 jlat = jpjglo 250 joffset = n proc250 joffset = narea-1 251 251 jostride = jpnij 252 252 ELSE … … 513 513 IF ( ABS( zlam - zplam(jo) ) < 1e-6 ) THEN 514 514 IF ( llinvalidcell(ji,jj) ) THEN 515 kproc(jo) = n proc+ 1000000515 kproc(jo) = narea-1 + 1000000 516 516 kobsi(jo) = ji + 1 517 517 kobsj(jo) = jj + 1 518 518 CYCLE 519 519 ELSE 520 kproc(jo) = n proc520 kproc(jo) = narea-1 521 521 kobsi(jo) = ji + 1 522 522 kobsj(jo) = jj + 1 … … 552 552 & zlamtm(:,ji,jj), zphitm(:,ji,jj) ) ) THEN 553 553 IF ( llinvalidcell(ji,jj) ) THEN 554 kproc(jo) = n proc+ 1000000554 kproc(jo) = narea-1 + 1000000 555 555 kobsi(jo) = ji + 1 556 556 kobsj(jo) = jj + 1 557 557 CYCLE 558 558 ELSE 559 kproc(jo) = n proc559 kproc(jo) = narea-1 560 560 kobsi(jo) = ji + 1 561 561 kobsj(jo) = jj + 1 … … 584 584 & zlamtm(:,ji,jj), zphitm(:,ji,jj) ) ) THEN 585 585 IF ( llinvalidcell(ji,jj) ) THEN 586 kproc(jo) = n proc+ 1000000586 kproc(jo) = narea-1 + 1000000 587 587 kobsi(jo) = ji + 1 588 588 kobsj(jo) = jj + 1 589 589 CYCLE 590 590 ELSE 591 kproc(jo) = n proc591 kproc(jo) = narea-1 592 592 kobsi(jo) = ji + 1 593 593 kobsj(jo) = jj + 1 … … 716 716 ! define the following format: "(a,a,ix.x,a,ix.x,a,ix.x,a)" 717 717 WRITE(clfmt, "('(a,a,i', i1, '.', i1',a,i', i1, '.', i1',a,i', i1, '.', i1',a)')") idg, idg, idg, idg, idg, idg 718 WRITE(cfname, clfmt ) TRIM(cn_gridsearchfile),'_', n proc,'of', jpni,'by', jpnj,'.nc'718 WRITE(cfname, clfmt ) TRIM(cn_gridsearchfile),'_', narea-1,'of', jpni,'by', jpnj,'.nc' 719 719 ENDIF 720 720 … … 820 820 CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 821 821 & 1, jpi, 1, jpj, & 822 & n proc, jpnij,&822 & narea-1, jpnij, & 823 823 & glamt, gphit, tmask, & 824 824 & nlons*nlats, lonsi, latsi, & … … 1070 1070 1071 1071 IF ( ( .NOT. ln_grid_global ) .OR. & 1072 & ( ( ln_grid_global ) .AND. ( n proc==0 ) ) ) THEN1072 & ( ( ln_grid_global ) .AND. ( narea-1==0 ) ) ) THEN 1073 1073 1074 1074 CALL chkerr( nf90_create (TRIM(cfname), nf90_clobber, idfile), &
Note: See TracChangeset
for help on using the changeset viewer.