- Timestamp:
- 2015-12-16T10:25:22+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_merge_2015/NEMOGCM/NEMO/TOP_SRC/trcdta.F90
r5385 r6060 36 36 TYPE(FLD), SAVE, PUBLIC, ALLOCATABLE, DIMENSION(:) :: sf_trcdta ! structure of input SST (file informations, fields read) 37 37 !$AGRIF_END_DO_NOT_TREAT 38 !! * Substitutions 39 # include "domzgr_substitute.h90" 38 40 39 !!---------------------------------------------------------------------- 41 40 !! NEMO/OPA 3.3 , NEMO Consortium (2010) … … 189 188 DO ji = 1, jpi 190 189 DO jk = 1, jpk ! determines the intepolated T-S profiles at each (i,j) points 191 zl = fsdept_n(ji,jj,jk)190 zl = gdept_n(ji,jj,jk) 192 191 IF( zl < gdept_1d(1 ) ) THEN ! above the first level of data 193 192 ztp(jk) = sf_dta(1)%fnow(ji,jj,1) … … 220 219 ik = mbkt(ji,jj) 221 220 IF( ik > 1 ) THEN 222 zl = ( gdept_1d(ik) - fsdept_n(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) )221 zl = ( gdept_1d(ik) - gdept_n(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 223 222 sf_dta(1)%fnow(ji,jj,ik) = (1.-zl) * sf_dta(1)%fnow(ji,jj,ik) + zl * sf_dta(1)%fnow(ji,jj,ik-1) 224 223 ENDIF … … 236 235 sf_dta(1)%fnow(:,:,:) = sf_dta(1)%fnow(:,:,:) * zrf_trfac ! multiplicative factor 237 236 ! 238 IF( lwp .AND. kt == nit000 ) THEN239 clndta = TRIM( sf_dta(1)%clvar )240 WRITE(numout,*) ''//clndta//' data '241 WRITE(numout,*)242 WRITE(numout,*)' level = 1'243 CALL prihre( sf_dta(1)%fnow(:,:,1), jpi, jpj, 1, jpi, 20, 1, jpj, 20, 1., numout )244 WRITE(numout,*)' level = ', jpk/2245 CALL prihre( sf_dta(1)%fnow(:,:,jpk/2), jpi, jpj, 1, jpi, 20, 1, jpj, 20, 1., numout )246 WRITE(numout,*)' level = ', jpkm1247 CALL prihre( sf_dta(1)%fnow(:,:,jpkm1), jpi, jpj, 1, jpi, 20, 1, jpj, 20, 1., numout )248 WRITE(numout,*)249 ENDIF250 237 ENDIF 251 238 ! … … 253 240 ! 254 241 END SUBROUTINE trc_dta 242 255 243 #else 256 244 !!---------------------------------------------------------------------- … … 262 250 END SUBROUTINE trc_dta 263 251 #endif 252 264 253 !!====================================================================== 265 254 END MODULE trcdta
Note: See TracChangeset
for help on using the changeset viewer.