- Timestamp:
- 2018-11-21T11:25:53+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/trcdta.F90
r10068 r10345 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.