- Timestamp:
- 2010-11-25T20:03:49+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/DTA/dtatem.F90
r2392 r2436 21 21 USE in_out_manager ! I/O manager 22 22 USE phycst ! physical constants 23 #if defined key_orca_lev10 24 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 25 #endif 23 26 24 IMPLICIT NONE 27 25 PRIVATE … … 68 66 #endif 69 67 REAL(wp):: zl 70 #if defined key_orca_lev1071 INTEGER :: ikr, ikw, ikt, jjk72 REAL(wp):: zfac73 #endif74 68 REAL(wp), DIMENSION(jpk) :: ztemdta ! auxiliary array for interpolation 69 ! 75 70 CHARACTER(len=100) :: cn_dir ! Root directory for location of ssr files 76 71 TYPE(FLD_N) :: sn_tem … … 104 99 CALL ctl_stop( 'dta_tem: unable to allocate sf_tem structure' ) ; RETURN 105 100 ENDIF 106 107 #if defined key_orca_lev10108 ALLOCATE( sf_tem(1)%fnow(jpi,jpj,jpkdta) )109 IF( sn_tem%ln_tint ) ALLOCATE( sf_tem(1)%fdta(jpi,jpj,jpkdta,2) )110 #else111 101 ALLOCATE( sf_tem(1)%fnow(jpi,jpj,jpk) ) 112 102 IF( sn_tem%ln_tint ) ALLOCATE( sf_tem(1)%fdta(jpi,jpj,jpk,2) ) 113 #endif114 103 ! ! fill sf_tem with sn_tem and control print 115 104 CALL fld_fill( sf_tem, (/ sn_tem /), cn_dir, 'dta_tem', 'Temperature data', 'namdta_tem' ) … … 177 166 #endif 178 167 179 #if defined key_orca_lev10180 DO jjk = 1, 5181 t_dta(:,:,jjk) = sf_tem(1)%fnow(:,:,1)182 END DO183 DO jk = 1, jpk-20,10184 ik = jk+5185 ikr = INT(jk/10) + 1186 ikw = (ikr-1) *10 + 1187 ikt = ikw + 5188 DO jjk=ikt,ikt+9189 zfac = ( gdept_0(jjk ) - gdepw_0(ikt) ) / ( gdepw_0(ikt+10) - gdepw_0(ikt) )190 t_dta(:,:,jjk) = sf_tem(1)%fnow(:,:,ikr) + ( sf_tem(1)%fnow(:,:,ikr+1) - sf_tem(1)%fnow(:,:,ikr) ) * zfac191 END DO192 END DO193 DO jjk = jpk-5, jpk194 t_dta(:,:,jjk) = sf_tem(1)%fnow(:,:,jpkdta-1)195 END DO196 ! fill the overlap areas197 CALL lbc_lnk (t_dta(:,:,:),'Z',-999.,'no0')198 #else199 168 t_dta(:,:,:) = sf_tem(1)%fnow(:,:,:) 200 #endif201 169 202 170 IF( ln_sco ) THEN
Note: See TracChangeset
for help on using the changeset viewer.