Changeset 6025
- Timestamp:
- 2015-12-09T11:10:32+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_MetOffice_merge_2015/NEMOGCM/NEMO/OPA_SRC/OBS/obs_prep.F90
r5998 r6025 1100 1100 REAL(KIND=wp), DIMENSION(2,2,kpk,kprofno) :: & 1101 1101 & zgmsk ! Grid mask 1102 REAL(KIND=wp), DIMENSION(2,2,kpk,kprofno) :: & 1103 & zgdepw 1102 1104 REAL(KIND=wp), DIMENSION(2,2,kprofno) :: & 1103 1105 & zglam, & ! Model longitude at grid points … … 1146 1148 CALL obs_int_comm_2d( 2, 2, kprofno, kpi, kpj, igrdi, igrdj, plam, zglam ) 1147 1149 CALL obs_int_comm_2d( 2, 2, kprofno, kpi, kpj, igrdi, igrdj, pphi, zgphi ) 1150 IF ( .NOT.( ln_zps .OR. ln_zco ) ) THEN 1151 ! Need to know the bathy depth for each observation for sco 1152 CALL obs_int_comm_3d( 2, 2, kprofno, kpi, kpj, kpk, igrdi, igrdj, fsdepw(:,:,:), & 1153 & zgdepw ) 1154 ENDIF 1148 1155 1149 1156 DO jobs = 1, kprofno … … 1221 1228 ! Flag if the observation is deeper than the bathymetry 1222 1229 ! Or if it is within the mask 1223 IF ( ALL( fsdepw(iig-1:iig+1,ijg-1:ijg+1,kpk) < pobsdep(jobsp) ) &1230 IF ( ALL( zgdepw(1:2,1:2,kpk,jobs) < pobsdep(jobsp) ) & 1224 1231 & .OR. & 1225 1232 & ( SUM( zgmsk(1:2,1:2,kobsk(jobsp)-1:kobsk(jobsp),jobs) ) &
Note: See TracChangeset
for help on using the changeset viewer.