Changeset 12476

Timestamp:

2020-02-27T12:31:33+01:00 (4 years ago)

Author:

jcastill

Message:

Fix bug in boundary vertical interpolation

File:

• NEMO/branches/UKMO/r12083_India_uncoupled/src/OCE/SBC/fldread.F90 (modified) (1 diff)

NEMO/branches/UKMO/r12083_India_uncoupled/src/OCE/SBC/fldread.F90

@@ -926 +926 @@
             IF(     zdepth(jk) < pdta_read_z(jb,1,          1) ) THEN                ! above the first level of external data
                pdta(jb,1,jk) =  pdta_read(jb,1,1)
-            ELSEIF( zdepth(jk) > pdta_read_z(jb,1,ipkb) ) THEN                       ! below the last level of external data 
+            ELSEIF( zdepth(jk) > MAXVAL(pdta_read_z(jb,1,:),1) ) THEN                ! below the last level of external data 
                pdta(jb,1,jk) =  pdta_read(jb,1,MAXLOC(pdta_read_z(jb,1,:),1))
             ELSE                                                             ! inbetween: vertical interpolation between jkb & jkb+1