Changeset 13018


Ignore:
Timestamp:
2020-06-03T15:29:46+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/AMM15_v3_6_STABLE_package_collate_PS44/NEMOGCM/NEMO/OPA_SRC/OBS/obs_prep.F90

    r8058 r13018  
    17471747      REAL(KIND=wp), DIMENSION(2,2,kpk,kprofno) :: & 
    17481748         & zgmsk              ! Grid mask 
     1749      REAL(KIND=wp), DIMENSION(2,2,kpk,kprofno) :: & 
     1750         & zgdept, & 
     1751         & zgdepw          
    17491752      REAL(KIND=wp), DIMENSION(2,2,kprofno) :: & 
    17501753         & zglam, &           ! Model longitude at grid points 
     
    17921795      CALL obs_int_comm_2d( 2, 2, kprofno, igrdi, igrdj, plam, zglam ) 
    17931796      CALL obs_int_comm_2d( 2, 2, kprofno, igrdi, igrdj, pphi, zgphi ) 
     1797      ! Need to know the bathy depth for each observation for sco 
     1798      CALL obs_int_comm_3d( 2, 2, kprofno, kpi, kpj, kpk, igrdi, igrdj, fsdepw(:,:,:), & 
     1799         &                  zgdepw ) 
     1800      CALL obs_int_comm_3d( 2, 2, kprofno, kpi, kpj, kpk, igrdi, igrdj, fsdept(:,:,:), & 
     1801         &                  zgdept ) 
    17941802 
    17951803      DO jobs = 1, kprofno 
     
    18591867 
    18601868            ! Set observation depth equal to that of the first model depth 
    1861             IF ( pobsdep(jobsp) <= pdep(1) ) THEN 
    1862                pobsdep(jobsp) = pdep(1)   
     1869            IF ( pobsdep(jobsp) < MINVAL(zgdept(1:2,1:2,1,jobs) ) ) THEN 
     1870               pobsdep(jobsp) = MINVAL(zgdept(1:2,1:2,1,jobs)) 
    18631871            ENDIF 
    18641872             
Note: See TracChangeset for help on using the changeset viewer.