Changeset 6027
- Timestamp:
- 2015-12-09T12:40:24+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r4650_general_vert_coord_obsoper/NEMOGCM/NEMO/OPA_SRC/OBS/obs_prep.F90
r5838 r6027 1757 1757 REAL(KIND=wp), DIMENSION(2,2,kpk,kprofno) :: & 1758 1758 & zgmsk ! Grid mask 1759 REAL(KIND=wp), DIMENSION(2,2,kpk,kprofno) :: & 1760 & zgdepw 1759 1761 REAL(KIND=wp), DIMENSION(2,2,kprofno) :: & 1760 1762 & zglam, & ! Model longitude at grid points … … 1803 1805 CALL obs_int_comm_2d( 2, 2, kprofno, igrdi, igrdj, plam, zglam ) 1804 1806 CALL obs_int_comm_2d( 2, 2, kprofno, igrdi, igrdj, pphi, zgphi ) 1807 ! Need to know the bathy depth for each observation for sco 1808 CALL obs_int_comm_3d( 2, 2, kprofno, kpk, igrdi, igrdj, fsdepw(:,:,:), & 1809 & zgdepw ) 1805 1810 1806 1811 DO jobs = 1, kprofno … … 1878 1883 ! Flag if the observation is deeper than the bathymetry 1879 1884 ! Or if it is within the mask 1880 IF ( ALL( fsdepw(iig-1:iig+1,ijg-1:ijg+1,kpk) < pobsdep(jobsp) ) &1885 IF ( ALL( zgdepw(1:2,1:2,kpk,jobs) < pobsdep(jobsp) ) & 1881 1886 & .OR. & 1882 1887 & ( SUM( zgmsk(1:2,1:2,kobsk(jobsp)-1:kobsk(jobsp),jobs) ) & … … 1908 1913 ENDIF 1909 1914 1915 ! Flag if the observation falls is close to land 1916 IF ( MINVAL( zgmsk(1:2,1:2,kobsk(jobsp)-1:kobsk(jobsp),jobs) ) == & 1917 & 0.0_wp) THEN 1918 IF (ld_nea) kobsqc(jobsp) = kobsqc(jobsp) + 14 1919 knlaobs = knlaobs + 1 1920 ENDIF 1910 1921 1911 1922 ! Set observation depth equal to that of the first model depth
Note: See TracChangeset
for help on using the changeset viewer.