Changeset 7841 for branches/UKMO/dev_rev5518_OBS_DoNotAssim/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_prof.F90
- Timestamp:
- 2017-03-29T12:14:30+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_rev5518_OBS_DoNotAssim/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_prof.F90
r4990 r7841 308 308 inowin = 0 309 309 DO ji = 1, inpfiles(jj)%nobs 310 IF ( inpfiles(jj)%ioqc(ji) > 2) CYCLE311 IF ( ( inpfiles(jj)%ivqc(ji,1) > 2) .AND. &312 & ( inpfiles(jj)%ivqc(ji,2) > 2 )) CYCLE310 IF ( BTEST(inpfiles(jj)%ioqc(ji),2 ) ) CYCLE 311 IF ( BTEST(inpfiles(jj)%ivqc(ji,1),2) .AND. & 312 & BTEST(inpfiles(jj)%ivqc(ji,2),2) ) CYCLE 313 313 IF ( ( inpfiles(jj)%ptim(ji) > djulini(jj) ) .AND. & 314 314 & ( inpfiles(jj)%ptim(ji) <= djulend(jj) ) ) THEN … … 323 323 inowin = 0 324 324 DO ji = 1, inpfiles(jj)%nobs 325 IF ( inpfiles(jj)%ioqc(ji) > 2) CYCLE326 IF ( ( inpfiles(jj)%ivqc(ji,1) > 2) .AND. &327 & ( inpfiles(jj)%ivqc(ji,2) > 2 )) CYCLE325 IF ( BTEST(inpfiles(jj)%ioqc(ji),2 ) ) CYCLE 326 IF ( BTEST(inpfiles(jj)%ivqc(ji,1),2) .AND. & 327 & BTEST(inpfiles(jj)%ivqc(ji,2),2) ) CYCLE 328 328 IF ( ( inpfiles(jj)%ptim(ji) > djulini(jj) ) .AND. & 329 329 & ( inpfiles(jj)%ptim(ji) <= djulend(jj) ) ) THEN … … 338 338 inowin = 0 339 339 DO ji = 1, inpfiles(jj)%nobs 340 IF ( inpfiles(jj)%ioqc(ji) > 2) CYCLE341 IF ( ( inpfiles(jj)%ivqc(ji,1) > 2) .AND. &342 & ( inpfiles(jj)%ivqc(ji,2) > 2 )) CYCLE340 IF ( BTEST(inpfiles(jj)%ioqc(ji),2 ) ) CYCLE 341 IF ( BTEST(inpfiles(jj)%ivqc(ji,1),2) .AND. & 342 & BTEST(inpfiles(jj)%ivqc(ji,2),2) ) CYCLE 343 343 IF ( ( inpfiles(jj)%ptim(ji) > djulini(jj) ) .AND. & 344 344 & ( inpfiles(jj)%ptim(ji) <= djulend(jj) ) ) THEN … … 352 352 353 353 DO ji = 1, inpfiles(jj)%nobs 354 IF ( inpfiles(jj)%ioqc(ji) > 2) CYCLE355 IF ( ( inpfiles(jj)%ivqc(ji,1) > 2) .AND. &356 & ( inpfiles(jj)%ivqc(ji,2) > 2 )) CYCLE354 IF ( BTEST(inpfiles(jj)%ioqc(ji),2 ) ) CYCLE 355 IF ( BTEST(inpfiles(jj)%ivqc(ji,1),2) .AND. & 356 & BTEST(inpfiles(jj)%ivqc(ji,2),2) ) CYCLE 357 357 IF ( ( inpfiles(jj)%ptim(ji) > djulini(jj) ) .AND. & 358 358 & ( inpfiles(jj)%ptim(ji) <= djulend(jj) ) ) THEN … … 367 367 IF ( inpfiles(jj)%pdep(ij,ji) >= 6000. ) & 368 368 & CYCLE 369 IF ( ( inpfiles(jj)%ivlqc(ij,ji,1) <= 2) .AND. &370 & ( inpfiles(jj)%idqc(ij,ji) <= 2) ) THEN369 IF ( .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,1),2) .AND. & 370 & .NOT. BTEST(inpfiles(jj)%idqc(ij,ji),2) ) THEN 371 371 it3dt0 = it3dt0 + 1 372 372 ENDIF … … 377 377 IF ( inpfiles(jj)%pdep(ij,ji) >= 6000. ) & 378 378 & CYCLE 379 IF ( ( inpfiles(jj)%ivlqc(ij,ji,2) <= 2) .AND. &380 & ( inpfiles(jj)%idqc(ij,ji) <= 2) ) THEN379 IF ( .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,1),2) .AND. & 380 & .NOT. BTEST(inpfiles(jj)%idqc(ij,ji),2) ) THEN 381 381 is3dt0 = is3dt0 + 1 382 382 ENDIF … … 386 386 IF ( inpfiles(jj)%pdep(ij,ji) >= 6000. ) & 387 387 & CYCLE 388 IF ( ( ( ( inpfiles(jj)%ivlqc(ij,ji,1) <= 2) .AND. &389 & ( inpfiles(jj)%idqc(ij,ji) <= 2 )) .AND. &390 & 391 & ( ( ( inpfiles(jj)%ivlqc(ij,ji,2) <= 2) .AND. &392 & ( inpfiles(jj)%idqc(ij,ji) <= 2 )) .AND. &393 & 388 IF ( ( .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,1),2) .AND. & 389 & .NOT. BTEST(inpfiles(jj)%idqc(ij,ji),2) .AND. & 390 & ldt3d ) .OR. & 391 & ( .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,2),2) .AND. & 392 & .NOT. BTEST(inpfiles(jj)%idqc(ij,ji),2) .AND. & 393 & lds3d ) ) THEN 394 394 ip3dt = ip3dt + 1 395 395 llvalprof = .TRUE. … … 416 416 DO jj = 1, inobf 417 417 DO ji = 1, inpfiles(jj)%nobs 418 IF ( inpfiles(jj)%ioqc(ji) > 2) CYCLE419 IF ( ( inpfiles(jj)%ivqc(ji,1) > 2) .AND. &420 & ( inpfiles(jj)%ivqc(ji,2) > 2 )) CYCLE418 IF ( BTEST(inpfiles(jj)%ioqc(ji),2 ) ) CYCLE 419 IF ( BTEST(inpfiles(jj)%ivqc(ji,1),2) .AND. & 420 & BTEST(inpfiles(jj)%ivqc(ji,2),2) ) CYCLE 421 421 IF ( ( inpfiles(jj)%ptim(ji) > djulini(jj) ) .AND. & 422 422 & ( inpfiles(jj)%ptim(ji) <= djulend(jj) ) ) THEN … … 431 431 DO jj = 1, inobf 432 432 DO ji = 1, inpfiles(jj)%nobs 433 IF ( inpfiles(jj)%ioqc(ji) > 2) CYCLE434 IF ( ( inpfiles(jj)%ivqc(ji,1) > 2) .AND. &435 & ( inpfiles(jj)%ivqc(ji,2) > 2 )) CYCLE433 IF ( BTEST(inpfiles(jj)%ioqc(ji),2 ) ) CYCLE 434 IF ( BTEST(inpfiles(jj)%ivqc(ji,1),2) .AND. & 435 & BTEST(inpfiles(jj)%ivqc(ji,2),2) ) CYCLE 436 436 IF ( ( inpfiles(jj)%ptim(ji) > djulini(jj) ) .AND. & 437 437 & ( inpfiles(jj)%ptim(ji) <= djulend(jj) ) ) THEN … … 480 480 ji = iprofidx(iindx(jk)) 481 481 482 IF ( inpfiles(jj)%ioqc(ji) > 2) CYCLE483 IF ( ( inpfiles(jj)%ivqc(ji,1) > 2) .AND. &484 & ( inpfiles(jj)%ivqc(ji,2) > 2 )) CYCLE482 IF ( BTEST(inpfiles(jj)%ioqc(ji),2 ) ) CYCLE 483 IF ( BTEST(inpfiles(jj)%ivqc(ji,1),2) .AND. & 484 & BTEST(inpfiles(jj)%ivqc(ji,2),2) ) CYCLE 485 485 486 486 IF ( ( inpfiles(jj)%ptim(ji) > djulini(jj) ) .AND. & … … 495 495 llvalprof = .FALSE. 496 496 497 IF ( inpfiles(jj)%ioqc(ji) > 2 ) CYCLE 498 499 IF ( ( inpfiles(jj)%ivqc(ji,1) > 2 ) .AND. & 500 & ( inpfiles(jj)%ivqc(ji,2) > 2 ) ) CYCLE 497 IF ( BTEST(inpfiles(jj)%ioqc(ji),2 ) ) CYCLE 498 IF ( BTEST(inpfiles(jj)%ivqc(ji,1),2) .AND. & 499 & BTEST(inpfiles(jj)%ivqc(ji,2),2) ) CYCLE 501 500 502 501 loop_prof : DO ij = 1, inpfiles(jj)%nlev … … 505 504 & CYCLE 506 505 507 IF ( ( inpfiles(jj)%ivlqc(ij,ji,1) <= 2) .AND. &508 & ( inpfiles(jj)%idqc(ij,ji) <= 2) ) THEN506 IF ( .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,1),2) .AND. & 507 & .NOT. BTEST(inpfiles(jj)%idqc(ij,ji),2) ) THEN 509 508 510 509 llvalprof = .TRUE. … … 513 512 ENDIF 514 513 515 IF ( ( inpfiles(jj)%ivlqc(ij,ji,2) <= 2) .AND. &516 & ( inpfiles(jj)%idqc(ij,ji) <= 2) ) THEN514 IF ( .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,2),2) .AND. & 515 & .NOT. BTEST(inpfiles(jj)%idqc(ij,ji),2) ) THEN 517 516 518 517 llvalprof = .TRUE. … … 592 591 IF (ldsatt) THEN 593 592 594 IF ( ( ( ( inpfiles(jj)%ivlqc(ij,ji,1) <= 2) .AND. &595 & ( inpfiles(jj)%idqc(ij,ji) <= 2 )) .AND. &596 & 597 & ( ( ( inpfiles(jj)%ivlqc(ij,ji,2) <= 2) .AND. &598 & ( inpfiles(jj)%idqc(ij,ji) <= 2 )) .AND. &599 & 593 IF ( ( .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,1),2) .AND. & 594 & .NOT. BTEST(inpfiles(jj)%idqc(ij,ji),2) .AND. & 595 & ldt3d ) .OR. & 596 & ( .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,2),2) .AND. & 597 & .NOT. BTEST(inpfiles(jj)%idqc(ij,ji),2) .AND. & 598 & lds3d ) ) THEN 600 599 ip3dt = ip3dt + 1 601 600 ELSE … … 605 604 ENDIF 606 605 607 IF ( ( ( ( inpfiles(jj)%ivlqc(ij,ji,1) <= 2) .AND. &608 & ( inpfiles(jj)%idqc(ij,ji) <= 2 )) .AND. &609 &ldt3d ) .OR. ldsatt ) THEN606 IF ( ( .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,1),2) .AND. & 607 & .NOT. BTEST(inpfiles(jj)%idqc(ij,ji),2) .AND. & 608 & ldt3d ) .OR. ldsatt ) THEN 610 609 611 610 IF (ldsatt) THEN … … 667 666 ENDIF 668 667 669 IF ( ( ( ( inpfiles(jj)%ivlqc(ij,ji,2) <= 2) .AND. &670 & ( inpfiles(jj)%idqc(ij,ji) <= 2 )) .AND. &671 & lds3d ) .OR. ldsatt ) THEN668 IF ( ( .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,1),2) .AND. & 669 & .NOT. BTEST(inpfiles(jj)%idqc(ij,ji),2) .AND. & 670 & ldt3d ) .OR. ldsatt ) THEN 672 671 673 672 IF (ldsatt) THEN
Note: See TracChangeset
for help on using the changeset viewer.
