Changeset 14035 for NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/vremap.F90
- Timestamp:
- 2020-12-03T12:18:31+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/vremap.F90
r13939 r14035 380 380 IF ( pzout(jkout) <= pzin( 1 ) ) THEN ! Surface extrapolation 381 381 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) 385 387 END DO 386 388 ELSEIF ( pzout(jkout) >= pzin(kjpk_in) ) THEN ! Bottom extrapolation 387 389 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) 391 395 END DO 392 396 ELSEIF ( ( pzout(jkout) > pzin(1) ).AND.( pzout(jkout) < pzin(kjpk_in) )) THEN
Note: See TracChangeset
for help on using the changeset viewer.