- Timestamp:
- 2019-10-02T10:50:10+02:00 (5 years ago)
- 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 1747 1747 REAL(KIND=wp), DIMENSION(2,2,kpk,kprofno) :: & 1748 1748 & 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 1749 1753 REAL(KIND=wp), DIMENSION(2,2,kprofno) :: & 1750 1754 & zglam, & ! Model longitude at grid points … … 1792 1796 CALL obs_int_comm_2d( 2, 2, kprofno, igrdi, igrdj, plam, zglam ) 1793 1797 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 ) 1794 1804 1795 1805 DO jobs = 1, kprofno … … 1858 1868 ENDIF 1859 1869 1870 ! RJR: Rob's change https://code.metoffice.gov.uk/trac/utils/ticket/246 1860 1871 ! Set observation depth equal to that of the first model depth 1861 IF ( pobsdep(jobsp) < = pdep(1) ) THEN1862 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)) 1863 1874 ENDIF 1864 1875
Note: See TracChangeset
for help on using the changeset viewer.