Changeset 5082 for LMDZ6/branches/Amaury_dev/libf/phylmd/iophy.F90
- Timestamp:
- Jul 19, 2024, 5:41:58 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/iophy.F90
r4848 r5082 328 328 plon_bounds(i,1)=rlon_glo(tabij(i)-1) 329 329 plon_bounds(i,2)=rlon_glo(tabij(i)+1) 330 IF (plon_bounds(i,2) .LE.0..AND.plon_bounds(i,1).GE.0.) THEN331 IF (rlon_glo(tabij(i)) .GE.0.) THEN330 IF (plon_bounds(i,2)<=0..AND.plon_bounds(i,1)>=0.) THEN 331 IF (rlon_glo(tabij(i))>=0.) THEN 332 332 plon_bounds(i,2)=-1*plon_bounds(i,2) 333 333 ENDIF 334 334 ENDIF 335 IF (plon_bounds(i,2) .GE.0..AND.plon_bounds(i,1).LE.0.) THEN336 IF (rlon_glo(tabij(i)) .LE.0.) THEN335 IF (plon_bounds(i,2)>=0..AND.plon_bounds(i,1)<=0.) THEN 336 IF (rlon_glo(tabij(i))<=0.) THEN 337 337 plon_bounds(i,2)=-1*plon_bounds(i,2) 338 338 ENDIF 339 339 ENDIF 340 340 ! 341 IF ( tabij(i) .LE.nbp_lon) THEN341 IF ( tabij(i)<=nbp_lon) THEN 342 342 plat_bounds(i,1)=rlat_glo(tabij(i)) 343 343 ELSE … … 361 361 362 362 CALL grid1dTo2d_glo(rlon_glo,zx_lon) 363 IF ((nbp_lon*nbp_lat) .GT.1) THEN363 IF ((nbp_lon*nbp_lat)>1) THEN 364 364 DO i = 1, nbp_lon 365 365 zx_lon(i,1) = rlon_glo(i+1) … … 375 375 plon_bounds(i,2)=zx_lon(ipt(i)+1,jpt(i)) 376 376 377 IF (ipt(i) .EQ.1) THEN377 IF (ipt(i)==1) THEN 378 378 plon_bounds(i,1)=zx_lon(nbp_lon,jpt(i)) 379 379 plon_bounds(i,2)=360.+zx_lon(ipt(i)+1,jpt(i)) 380 380 ENDIF 381 381 382 IF (ipt(i) .EQ.nbp_lon) THEN382 IF (ipt(i)==nbp_lon) THEN 383 383 plon_bounds(i,2)=360.+zx_lon(1,jpt(i)) 384 384 ENDIF … … 387 387 plat_bounds(i,2)=zx_lat(ipt(i),jpt(i)+1) 388 388 389 IF (jpt(i) .EQ.1) THEN389 IF (jpt(i)==1) THEN 390 390 plat_bounds(i,1)=zx_lat(ipt(i),1)+0.001 391 391 plat_bounds(i,2)=zx_lat(ipt(i),1)-0.001 392 392 ENDIF 393 393 394 IF (jpt(i) .EQ.nbp_lat) THEN394 IF (jpt(i)==nbp_lat) THEN 395 395 plat_bounds(i,1)=zx_lat(ipt(i),nbp_lat)+0.001 396 396 plat_bounds(i,2)=zx_lat(ipt(i),nbp_lat)-0.001 … … 412 412 tabprocbeg=klon_mpi_begin 413 413 tabprocend=klon_mpi_end 414 IF(tabij(ip) .GE.tabprocbeg.AND.tabij(ip).LE.tabprocend) THEN414 IF(tabij(ip)>=tabprocbeg.AND.tabij(ip)<=tabprocend) THEN 415 415 npproc=npproc+1 416 416 npstn=npproc … … 425 425 npproc=0 426 426 DO ip=1, pim 427 IF(tabij(ip) .GE.tabprocbeg.AND.tabij(ip).LE.tabprocend) THEN427 IF(tabij(ip)>=tabprocbeg.AND.tabij(ip)<=tabprocend) THEN 428 428 npproc=npproc+1 429 429 nptabij(npproc)=tabij(ip) … … 852 852 DO ip=1, npstn 853 853 ! print*,'histwrite2d is_sequential npstn ip name nptabij',npstn,ip,name,nptabij(ip) 854 IF(nptabij(ip) .GE.klon_mpi_begin.AND. &855 nptabij(ip) .LE.klon_mpi_end) THEN854 IF(nptabij(ip)>=klon_mpi_begin.AND. & 855 nptabij(ip)<=klon_mpi_end) THEN 856 856 fieldok(ip)=buffer_omp(nptabij(ip)-klon_mpi_begin+1) 857 857 ENDIF … … 923 923 DO n=1, nlev 924 924 DO ip=1, npstn 925 IF(nptabij(ip) .GE.klon_mpi_begin.AND. &926 nptabij(ip) .LE.klon_mpi_end) THEN925 IF(nptabij(ip)>=klon_mpi_begin.AND. & 926 nptabij(ip)<=klon_mpi_end) THEN 927 927 fieldok(ip,n)=buffer_omp(nptabij(ip)-klon_mpi_begin+1,n) 928 928 ENDIF … … 1127 1127 DO ip=1, npstn 1128 1128 write(lunout,*)'histwrite2d_phy is_sequential npstn ip namenptabij',npstn,ip,var%name,nptabij(ip) 1129 IF(nptabij(ip) .GE.klon_mpi_begin.AND. &1130 nptabij(ip) .LE.klon_mpi_end) THEN1129 IF(nptabij(ip)>=klon_mpi_begin.AND. & 1130 nptabij(ip)<=klon_mpi_end) THEN 1131 1131 fieldok(ip)=buffer_omp(nptabij(ip)-klon_mpi_begin+1) 1132 1132 ENDIF … … 1331 1331 DO n=1, nlev 1332 1332 DO ip=1, npstn 1333 IF(nptabij(ip) .GE.klon_mpi_begin.AND. &1334 nptabij(ip) .LE.klon_mpi_end) THEN1333 IF(nptabij(ip)>=klon_mpi_begin.AND. & 1334 nptabij(ip)<=klon_mpi_end) THEN 1335 1335 fieldok(ip,n)=buffer_omp(nptabij(ip)-klon_mpi_begin+1,n) 1336 1336 ENDIF … … 1420 1420 DO ip=1, npstn 1421 1421 PRINT*,'histwrite2d_xios is_sequential npstn ip namenptabij',npstn,ip,field_name,nptabij(ip) 1422 IF(nptabij(ip) .GE.klon_mpi_begin.AND. &1423 nptabij(ip) .LE.klon_mpi_end) THEN1422 IF(nptabij(ip)>=klon_mpi_begin.AND. & 1423 nptabij(ip)<=klon_mpi_end) THEN 1424 1424 fieldok(ip)=buffer_omp(nptabij(ip)-klon_mpi_begin+1) 1425 1425 ENDIF … … 1510 1510 DO n=1, nlev 1511 1511 DO ip=1, npstn 1512 IF(nptabij(ip) .GE.klon_mpi_begin.AND. &1513 nptabij(ip) .LE.klon_mpi_end) THEN1512 IF(nptabij(ip)>=klon_mpi_begin.AND. & 1513 nptabij(ip)<=klon_mpi_end) THEN 1514 1514 fieldok(ip,n)=buffer_omp(nptabij(ip)-klon_mpi_begin+1,n) 1515 1515 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.