Changeset 13018
- Timestamp:
- 2020-06-03T15:29:46+02:00 (4 years ago)
- 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 1747 1747 REAL(KIND=wp), DIMENSION(2,2,kpk,kprofno) :: & 1748 1748 & zgmsk ! Grid mask 1749 REAL(KIND=wp), DIMENSION(2,2,kpk,kprofno) :: & 1750 & zgdept, & 1751 & zgdepw 1749 1752 REAL(KIND=wp), DIMENSION(2,2,kprofno) :: & 1750 1753 & zglam, & ! Model longitude at grid points … … 1792 1795 CALL obs_int_comm_2d( 2, 2, kprofno, igrdi, igrdj, plam, zglam ) 1793 1796 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 ) 1794 1802 1795 1803 DO jobs = 1, kprofno … … 1859 1867 1860 1868 ! Set observation depth equal to that of the first model depth 1861 IF ( pobsdep(jobsp) < = pdep(1) ) THEN1862 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)) 1863 1871 ENDIF 1864 1872
Note: See TracChangeset
for help on using the changeset viewer.