Ignore:
Timestamp:
2019-12-10T15:36:47+01:00 (7 months ago)
Author:
jchanut
Message:

#2222: fixes linear vertical interpolation of viscosities

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11233_AGRIF-05_jchanut_vert_coord_interp/src/NST/vremap.F90

    r11802 r12152  
    374374      DO jkout = 1, kjpk_out !  Loop over destination grid 
    375375         ! 
    376          IF     ( pzout(jkout) < pzin(  1    ) ) THEN    ; ptout(jkout,1:kn_var) = ptin(    1    ,1:kn_var)          
    377          ELSEIF ( pzout(jkout) > pzin(kjpk_in) ) THEN    ; ptout(jkout,1:kn_var) = ptin( kjpk_in ,1:kn_var) 
    378          ELSEIF ( ( pzout(jkout) >= pzin(1) ).AND.( pzout(jkout) <= pzin(kjpk_in) )) THEN 
     376         IF     ( pzout(jkout) <= pzin(  1    ) ) THEN  ; ptout(jkout,1:kn_var) = ptin(    1    ,1:kn_var)          
     377         ELSEIF ( pzout(jkout) >= pzin(kjpk_in) ) THEN   ; ptout(jkout,1:kn_var) = ptin( kjpk_in ,1:kn_var) 
     378         ELSEIF ( ( pzout(jkout) > pzin(1) ).AND.( pzout(jkout) < pzin(kjpk_in) )) THEN 
    379379            DO jkin = 1, kjpk_in - 1 !  Loop over source grid 
    380                IF ( pzout(jkout) >=  pzin(jkin) ) THEN 
     380               IF ( pzout(jkout) < pzin(jkin+1) ) THEN 
    381381                  DO jn = 1, kn_var 
    382382                     ptout(jkout,jn) =  ptin(jkin,jn) + & 
     
    390390         ! 
    391391      END DO 
    392              
     392 
    393393   END SUBROUTINE remap_linear 
    394394 
Note: See TracChangeset for help on using the changeset viewer.