Ignore:
Timestamp:
2018-01-30T17:58:41+01:00 (3 years ago)
Author:
dford
Message:

For log10(chlorophyll) profiles, mask out the rest of the profile below any excluded values, to avoid interpolation issues.

File:
1 edited

Legend:

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

    r9229 r9289  
    920920      INTEGER :: jtype             ! Data loop variable 
    921921      INTEGER :: jvar              ! Variable number 
    922       INTEGER :: ji, jj            ! Loop counters 
     922      INTEGER :: ji, jj, jk        ! Loop counters 
    923923      REAL(wp) :: tiny             ! small number 
    924924      REAL(wp), POINTER, DIMENSION(:,:,:,:) :: & 
     
    10041004                  zprofmask(:,:,:,:) = 0 
    10051005               END WHERE 
     1006               ! Mask out model below any excluded values, 
     1007               ! to avoid interpolation issues 
     1008               DO jvar = 1, profdataqc(jtype)%nvar 
     1009                 DO jj = 1, jpj 
     1010                    DO ji = 1, jpi 
     1011                       depth_loop: DO jk = 1, jpk 
     1012                          IF ( zprofmask(ji,jj,jk,jvar) == 0 ) THEN 
     1013                             zprofmask(ji,jj,jk:jpk,jvar) = 0 
     1014                             EXIT depth_loop 
     1015                          ENDIF 
     1016                       END DO depth_loop 
     1017                    END DO 
     1018                 END DO 
     1019              END DO 
    10061020 
    10071021            CASE('pchltot') 
Note: See TracChangeset for help on using the changeset viewer.