New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7841 for branches/UKMO/dev_rev5518_OBS_DoNotAssim/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_vel.F90 – NEMO

Ignore:
Timestamp:
2017-03-29T12:14:30+02:00 (7 years ago)
Author:
jwhile
Message:

Added "Do not Assimlate" funtionality to OBS code

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/UKMO/dev_rev5518_OBS_DoNotAssim/NEMOGCM/NEMO/OPA_SRC/OBS/obs_read_vel.F90

    r4990 r7841  
    326326                        IF ( inpfiles(jj)%pdep(ij,ji) >= 6000. ) & 
    327327                           & CYCLE 
    328                         IF ( ( inpfiles(jj)%ivlqc(ij,ji,1) <= 2 ) .AND. & 
    329                            & ( inpfiles(jj)%ivlqc(ij,ji,2) <= 2 ) .AND. & 
    330                            & ( inpfiles(jj)%idqc(ij,ji) <= 2 ) ) THEN 
     328                        IF ( .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,1),2) .AND. & 
     329                           & .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,2),2) .AND. & 
     330                           & .NOT. BTEST(inpfiles(jj)%idqc(ij,ji),2) ) THEN 
    331331                           iuv3dt = iuv3dt + 1 
    332332                           llvalprof = .TRUE. 
    … …  
    413413                  & CYCLE 
    414414                
    415                IF ( ( inpfiles(jj)%ivlqc(ij,ji,1) <= 2 ) .AND. & 
    416                   & ( inpfiles(jj)%ivlqc(ij,ji,2) <= 2 ) .AND. & 
    417                   & ( inpfiles(jj)%idqc(ij,ji) <= 2 ) ) THEN 
     415               IF ( .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,1),2) .AND. & 
     416                  & .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,2),2) .AND. & 
     417                  & .NOT. BTEST(inpfiles(jj)%idqc(ij,ji),2) ) THEN 
    418418 
    419419                  llvalprof = .TRUE. 
    … …  
    492492                     & CYCLE 
    493493 
    494                   IF ( ( inpfiles(jj)%ivlqc(ij,ji,1) <= 2 ) .AND. & 
    495                      & ( inpfiles(jj)%ivlqc(ij,ji,2) <= 2 ) .AND. & 
    496                      & ( inpfiles(jj)%idqc(ij,ji) <= 2 ) ) THEN 
     494                  IF ( .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,1),2) .AND. & 
     495                     & .NOT. BTEST(inpfiles(jj)%ivlqc(ij,ji,2),2) .AND. & 
     496                     & .NOT. BTEST(inpfiles(jj)%idqc(ij,ji),2) ) THEN 
    497497                     iuv3dt = iuv3dt + 1 
    498498                  ELSE 
Note: See TracChangeset for help on using the changeset viewer.