- Timestamp:
- 2013-11-20T17:28:04+01:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM/dtatsd.F90
r4245 r4292 221 221 DO ji = 1, jpi 222 222 DO jk = 1, jpk ! determines the intepolated T-S profiles at each (i,j) points 223 zl = fsdept_0(ji,jj,jk)224 IF( zl < gdept_ 0(1 ) ) THEN ! above the first level of data223 zl = gdept_0(ji,jj,jk) 224 IF( zl < gdept_1d(1 ) ) THEN ! above the first level of data 225 225 ztp(jk) = ptsd(ji,jj,1 ,jp_tem) 226 226 zsp(jk) = ptsd(ji,jj,1 ,jp_sal) 227 ELSEIF( zl > gdept_ 0(jpk) ) THEN ! below the last level of data227 ELSEIF( zl > gdept_1d(jpk) ) THEN ! below the last level of data 228 228 ztp(jk) = ptsd(ji,jj,jpkm1,jp_tem) 229 229 zsp(jk) = ptsd(ji,jj,jpkm1,jp_sal) 230 230 ELSE ! inbetween : vertical interpolation between jkk & jkk+1 231 231 DO jkk = 1, jpkm1 ! when gdept(jkk) < zl < gdept(jkk+1) 232 IF( (zl-gdept_ 0(jkk)) * (zl-gdept_0(jkk+1)) <= 0._wp ) THEN233 zi = ( zl - gdept_ 0(jkk) ) / (gdept_0(jkk+1)-gdept_0(jkk))232 IF( (zl-gdept_1d(jkk)) * (zl-gdept_1d(jkk+1)) <= 0._wp ) THEN 233 zi = ( zl - gdept_1d(jkk) ) / (gdept_1d(jkk+1)-gdept_1d(jkk)) 234 234 ztp(jk) = ptsd(ji,jj,jkk,jp_tem) + ( ptsd(ji,jj,jkk+1,jp_tem) - ptsd(ji,jj,jkk,jp_tem) ) * zi 235 235 zsp(jk) = ptsd(ji,jj,jkk,jp_sal) + ( ptsd(ji,jj,jkk+1,jp_sal) - ptsd(ji,jj,jkk,jp_sal) ) * zi … … 259 259 ik = mbkt(ji,jj) 260 260 IF( ik > 1 ) THEN 261 zl = ( gdept_ 0(ik) - fsdept_0(ji,jj,ik) ) / ( gdept_0(ik) - gdept_0(ik-1) )261 zl = ( gdept_1d(ik) - gdept_0(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 262 262 ptsd(ji,jj,ik,jp_tem) = (1.-zl) * ptsd(ji,jj,ik,jp_tem) + zl * ptsd(ji,jj,ik-1,jp_tem) 263 263 ptsd(ji,jj,ik,jp_sal) = (1.-zl) * ptsd(ji,jj,ik,jp_sal) + zl * ptsd(ji,jj,ik-1,jp_sal)
Note: See TracChangeset
for help on using the changeset viewer.