Changeset 10843 for NEMO/branches/UKMO/dev_r10037_GPU/src/TOP/trcdta.F90
- Timestamp:
- 2019-04-05T16:01:32+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/dev_r10037_GPU/src/TOP/trcdta.F90
r9950 r10843 40 40 41 41 !!---------------------------------------------------------------------- 42 !! NEMO/ OCE4.0 , NEMO Consortium (2018)42 !! NEMO/TOP 4.0 , NEMO Consortium (2018) 43 43 !! $Id$ 44 !! Software governed by the CeCILL licen ce (./LICENSE)44 !! Software governed by the CeCILL license (see ./LICENSE) 45 45 !!---------------------------------------------------------------------- 46 46 CONTAINS … … 217 217 ELSE !== z- or zps- coordinate ==! 218 218 ! zps-coordinate (partial steps) interpolation at the last ocean level 219 IF( ln_zps ) THEN220 DO jj = 1, jpj221 DO ji = 1, jpi222 ik = mbkt(ji,jj)223 IF( ik > 1 ) THEN224 zl = ( gdept_1d(ik) - gdept_0(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) )225 ptrcdta(ji,jj,ik) = (1.-zl) * ptrcdta(ji,jj,ik) + zl * ptrcdta(ji,jj,ik-1)226 ENDIF227 ik = mikt(ji,jj)228 IF( ik > 1 ) THEN229 zl = ( gdept_0(ji,jj,ik) - gdept_1d(ik) ) / ( gdept_1d(ik+1) - gdept_1d(ik) )230 ptrcdta(ji,jj,ik) = (1.-zl) * ptrcdta(ji,jj,ik) + zl * ptrcdta(ji,jj,ik+1)231 ENDIF232 END DO233 END DO234 ENDIF219 ! IF( ln_zps ) THEN 220 ! DO jj = 1, jpj 221 ! DO ji = 1, jpi 222 ! ik = mbkt(ji,jj) 223 ! IF( ik > 1 ) THEN 224 ! zl = ( gdept_1d(ik) - gdept_0(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 225 ! ptrcdta(ji,jj,ik) = (1.-zl) * ptrcdta(ji,jj,ik) + zl * ptrcdta(ji,jj,ik-1) 226 ! ENDIF 227 ! ik = mikt(ji,jj) 228 ! IF( ik > 1 ) THEN 229 ! zl = ( gdept_0(ji,jj,ik) - gdept_1d(ik) ) / ( gdept_1d(ik+1) - gdept_1d(ik) ) 230 ! ptrcdta(ji,jj,ik) = (1.-zl) * ptrcdta(ji,jj,ik) + zl * ptrcdta(ji,jj,ik+1) 231 ! ENDIF 232 ! END DO 233 ! END DO 234 ! ENDIF 235 235 ! 236 236 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.