Ignore:
Timestamp:
2019-02-21T16:44:54+01:00 (3 years ago)
Author:
emmafiedler
Message:

Correct flagging of ice data at zeroth timestep, QC for sea ice thickness

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90

    r10276 r10712  
    191191         & cn_altbiasfile        ! Altimeter bias input filename 
    192192 
     193      LOGICAL :: ln_seaicetypes = .FALSE.  ! Logical switch indicating data type is sea ice 
    193194 
    194195      LOGICAL :: ln_t3d          ! Logical switch for temperature profiles 
     
    844845               &               rn_dobsini, rn_dobsend, ln_ignmis, .FALSE., llnightav(jtype) ) 
    845846 
    846             CALL obs_pre_surf( surfdata(jtype), surfdataqc(jtype), ln_nea, ln_bound_reject ) 
     847 
     848            IF ( TRIM(cobstypessurf(jtype)) == 'sit' ) THEN 
     849              ln_seaicetypes = .TRUE. 
     850            ELSE 
     851              ln_seaicetypes = .FALSE. 
     852            ENDIF            
     853            CALL obs_pre_surf( surfdata(jtype), surfdataqc(jtype), ln_nea, ln_bound_reject, ln_seaicetypes ) 
    847854 
    848855            IF ( TRIM(cobstypessurf(jtype)) == 'sla' ) THEN 
     
    12331240                        &           'time-step but some obs are valid then.' ) 
    12341241                     WRITE(numout,*)surfdataqc(jtype)%nsstpmpp(1), & 
    1235                         &           ' sea-ice concentration obs will be missed' 
    1236                   ENDIF 
    1237                   surfdataqc(jtype)%nsurfup = surfdataqc(jtype)%nsurfup + & 
    1238                      &                        surfdataqc(jtype)%nsstp(1) 
    1239                   CYCLE 
     1242                        &           ' sea-ice concentration obs will be missed'                   
     1243                  ENDIF                                 
    12401244               ELSE 
    12411245#if defined key_cice 
     
    12491253               ENDIF 
    12501254            CASE('sit') 
    1251                IF ( kstp == 0 ) THEN ! **Copied from SIC, check applies to SIT!** 
     1255               IF ( kstp == 0 ) THEN 
    12521256                  IF ( lwp .AND. surfdataqc(jtype)%nsstpmpp(1) > 0 ) THEN 
    12531257                     CALL ctl_warn( 'Sea-ice not initialised on zeroth '// & 
    12541258                        &           'time-step but some obs are valid then.' ) 
    12551259                     WRITE(numout,*)surfdataqc(jtype)%nsstpmpp(1), & 
    1256                         &           ' sea-ice thickness obs will be missed' 
    1257                   ENDIF 
    1258                   surfdataqc(jtype)%nsurfup = surfdataqc(jtype)%nsurfup + & 
    1259                      &                        surfdataqc(jtype)%nsstp(1) 
    1260                   CYCLE 
    1261                ELSE             
     1260                        &           ' sea-ice thickness obs will be missed and QC flag set to 4' 
     1261                  ENDIF                                   
     1262               ELSE        
    12621263#if defined key_cice 
    12631264                  zsurfvar(:,:) = thick_i(:,:) 
Note: See TracChangeset for help on using the changeset viewer.