Changeset 13021


Ignore:
Timestamp:
2020-06-03T16:09:48+02:00 (4 months ago)
Author:
petesykes
Message:

Correct bug where obs above 3m are reset to 3m

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_obs_oper_update_PS44/NEMOGCM/NEMO/OPA_SRC/OBS/obs_prep.F90

    r9306 r13021  
    11591159#endif  
    11601160      REAL(KIND=wp), DIMENSION(2,2,kpk,kprofno) :: & 
     1161         & zgdept, & 
    11611162         & zgdepw          
    11621163      REAL(KIND=wp), DIMENSION(2,2,kprofno) :: & 
     
    12231224      CALL obs_int_comm_3d( 2, 2, kprofno, kpi, kpj, kpk, igrdi, igrdj, fsdepw(:,:,:), & 
    12241225         &                  zgdepw ) 
     1226      CALL obs_int_comm_3d( 2, 2, kprofno, kpi, kpj, kpk, igrdi, igrdj, fsdept(:,:,:), & 
     1227         &                  zgdept ) 
    12251228 
    12261229      DO jobs = 1, kprofno 
     
    13281331             
    13291332            ! Set observation depth equal to that of the first model depth 
    1330             IF ( pobsdep(jobsp) <= pdep(1) ) THEN 
    1331                pobsdep(jobsp) = pdep(1)   
     1333            IF ( pobsdep(jobsp) < MINVAL(zgdept(1:2,1:2,1,jobs) ) ) THEN 
     1334               pobsdep(jobsp) = MINVAL(zgdept(1:2,1:2,1,jobs)) 
    13321335            ENDIF 
    13331336             
Note: See TracChangeset for help on using the changeset viewer.