Changeset 11451

Timestamp:

2019-08-19T12:51:16+02:00 (5 years ago)

Author:

kingr

Message:

Modified sheck on min obs deth to use 3D depth array.

File:

• branches/UKMO/dev_r5518_obs_oper_update_NearSurfDepths/NEMOGCM/NEMO/OPA_SRC/OBS/obs_prep.F90 (modified) (3 diffs)

branches/UKMO/dev_r5518_obs_oper_update_NearSurfDepths/NEMOGCM/NEMO/OPA_SRC/OBS/obs_prep.F90

@@ -1159 +1159 @@
 #endif 
       REAL(KIND=wp), DIMENSION(2,2,kpk,kprofno) :: &
+         & zgdept, &
          & zgdepw         
       REAL(KIND=wp), DIMENSION(2,2,kprofno) :: &
@@ -1223 +1224 @@
       CALL obs_int_comm_3d( 2, 2, kprofno, kpi, kpj, kpk, igrdi, igrdj, fsdepw(:,:,:), &
          &                  zgdepw )
+      CALL obs_int_comm_3d( 2, 2, kprofno, kpi, kpj, kpk, igrdi, igrdj, fsdept(:,:,:), &
+         &                  zgdept )
 
       DO jobs = 1, kprofno
@@ -1328 +1331 @@
             
             ! Set observation depth equal to that of the first model depth
-            IF ( pobsdep(jobsp) <= pdep(1) ) THEN
-               pobsdep(jobsp) = pdep(1)  
+            IF ( pobsdep(jobsp) < MINVAL(zgdept(1:2,1:2,1,jobs) ) ) THEN
+               pobsdep(jobsp) = MINVAL(zgdept(1:2,1:2,1,jobs))
             ENDIF