New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 14035 for NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/vremap.F90 – NEMO

Ignore:
Timestamp:
2020-12-03T12:18:31+01:00 (4 years ago)
Author:
jchanut
Message:

#2222, remove linear vertical extrapolation over partial cells + correct initial state interp for top

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/vremap.F90

    r13939 r14035  
    380380         IF     ( pzout(jkout) <=  pzin(  1    ) ) THEN ! Surface extrapolation   
    381381            DO jn = 1, kn_var  
    382                ptout(jkout,jn) = ptin(1 ,jn) + & 
    383                                & (pzout(jkout) - pzin(1)) / (pzin(2)    - pzin(1)) & 
    384                                &                          * (ptin(2,jn) - ptin(1,jn)) 
     382! linear 
     383!               ptout(jkout,jn) = ptin(1 ,jn) + & 
     384!                               & (pzout(jkout) - pzin(1)) / (pzin(2)    - pzin(1)) & 
     385!                               &                          * (ptin(2,jn) - ptin(1,jn)) 
     386               ptout(jkout,jn) = ptin(1,jn) 
    385387            END DO 
    386388         ELSEIF ( pzout(jkout) >= pzin(kjpk_in) ) THEN ! Bottom extrapolation  
    387389            DO jn = 1, kn_var  
    388                ptout(jkout,jn) = ptin(kjpk_in ,jn) + & 
    389                                & (pzout(jkout) - pzin(kjpk_in)) / (pzin(kjpk_in)    - pzin(kjpk_in-1)) & 
    390                                &                                * (ptin(kjpk_in,jn) - ptin(kjpk_in-1,jn)) 
     390! linear 
     391!               ptout(jkout,jn) = ptin(kjpk_in ,jn) + & 
     392!                               & (pzout(jkout) - pzin(kjpk_in)) / (pzin(kjpk_in)    - pzin(kjpk_in-1)) & 
     393!                               &                                * (ptin(kjpk_in,jn) - ptin(kjpk_in-1,jn)) 
     394               ptout(jkout,jn) = ptin(kjpk_in ,jn) 
    391395            END DO 
    392396         ELSEIF ( ( pzout(jkout) > pzin(1) ).AND.( pzout(jkout) < pzin(kjpk_in) )) THEN 
Note: See TracChangeset for help on using the changeset viewer.