Changeset 12789

Timestamp:

2020-04-21T14:01:44+02:00 (4 years ago)

Author:

rrenshaw

Message:

adopt one of Rob's profile fixes

File:

• branches/UKMO/AMM15_v3_6_STABLE_package_reanalysis4/NEMOGCM/NEMO/OPA_SRC/OBS/obs_level_search.h90 (modified) (3 diffs)

branches/UKMO/AMM15_v3_6_STABLE_package_reanalysis4/NEMOGCM/NEMO/OPA_SRC/OBS/obs_level_search.h90

@@ -13 +13 @@
       !! ** Method  : Straightforward search
       !!
-      !! ** Action  : 
+      !! ** Action  : Will return level associated with T-point below the obs
+      !!              depth, except when observation is in the top box will
+      !!              return level 2. Also, if obs depth greater than depth
+      !!              of last wet T-point (kpk-1) will return level kpk.
       !!
       !! History :
@@ -21 +24 @@
       !!        !  2006-10  (A. Weaver) Cleanup
       !!        !  2008-10  (K. Mogensen) Remove assumptions on grid.
-      !!        !  2019-07  (R. Renshaw) Don't exceed model bottom.
+      !!        !  2020-04  (R. King/RJR) Add Rob's fix
+      !!        !            https://code.metoffice.gov.uk/trac/utils/ticket/335
       !!----------------------------------------------------------------------
 
@@ -44 +48 @@
       DO ji = 1, kobs 
          kobsk(ji) = 1
-         depk: DO jk = 2, kgrd
+         depk: DO jk = 2, kgrd-1
             IF ( pgrddep(jk) >= pobsdep(ji) ) EXIT depk
          END DO depk
-         kobsk(ji) = MIN( jk, kgrd )
+         kobsk(ji) = jk
       END DO