Changeset 12476


Ignore:
Timestamp:
2020-02-27T12:31:33+01:00 (8 months ago)
Author:
jcastill
Message:

Fix bug in boundary vertical interpolation

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/r12083_India_uncoupled/src/OCE/SBC/fldread.F90

    r11715 r12476  
    926926            IF(     zdepth(jk) < pdta_read_z(jb,1,          1) ) THEN                ! above the first level of external data 
    927927               pdta(jb,1,jk) =  pdta_read(jb,1,1) 
    928             ELSEIF( zdepth(jk) > pdta_read_z(jb,1,ipkb) ) THEN                       ! below the last level of external data  
     928            ELSEIF( zdepth(jk) > MAXVAL(pdta_read_z(jb,1,:),1) ) THEN                ! below the last level of external data  
    929929               pdta(jb,1,jk) =  pdta_read(jb,1,MAXLOC(pdta_read_z(jb,1,:),1)) 
    930930            ELSE                                                             ! inbetween: vertical interpolation between jkb & jkb+1 
Note: See TracChangeset for help on using the changeset viewer.