Changeset 11634


Ignore:
Timestamp:
2019-10-02T10:50:10+02:00 (13 months ago)
Author:
rrenshaw
Message:

Add Rob's change ticket 246 to use surface data at coorect depth

File:
1 edited

Legend:

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

    r8058 r11634  
    17471747      REAL(KIND=wp), DIMENSION(2,2,kpk,kprofno) :: & 
    17481748         & zgmsk              ! Grid mask 
     1749      ! RJR: Rob's change https://code.metoffice.gov.uk/trac/utils/ticket/246 
     1750      REAL(KIND=wp), DIMENSION(2,2,kpk,kprofno) :: & 
     1751         & zgdept, & 
     1752         & zgdepw 
    17491753      REAL(KIND=wp), DIMENSION(2,2,kprofno) :: & 
    17501754         & zglam, &           ! Model longitude at grid points 
     
    17921796      CALL obs_int_comm_2d( 2, 2, kprofno, igrdi, igrdj, plam, zglam ) 
    17931797      CALL obs_int_comm_2d( 2, 2, kprofno, igrdi, igrdj, pphi, zgphi ) 
     1798      ! RJR: Rob's change https://code.metoffice.gov.uk/trac/utils/ticket/246 
     1799      ! Need to know the bathy depth for each observation for sco 
     1800      CALL obs_int_comm_3d( 2, 2, kprofno, kpi, kpj, kpk, igrdi, igrdj, fsdepw(:,:,:), & 
     1801                           zgdepw ) 
     1802      CALL obs_int_comm_3d( 2, 2, kprofno, kpi, kpj, kpk, igrdi, igrdj, fsdept(:,:,:), & 
     1803                           zgdept ) 
    17941804 
    17951805      DO jobs = 1, kprofno 
     
    18581868            ENDIF 
    18591869 
     1870            ! RJR: Rob's change https://code.metoffice.gov.uk/trac/utils/ticket/246 
    18601871            ! Set observation depth equal to that of the first model depth 
    1861             IF ( pobsdep(jobsp) <= pdep(1) ) THEN 
    1862                pobsdep(jobsp) = pdep(1)   
     1872            IF ( pobsdep(jobsp) < MINVAL(zgdept(1:2,1:2,1,jobs) ) ) THEN 
     1873               pobsdep(jobsp) = MINVAL(zgdept(1:2,1:2,1,jobs)) 
    18631874            ENDIF 
    18641875             
Note: See TracChangeset for help on using the changeset viewer.