Changeset 63
- Timestamp:
- 2004-04-22T13:57:46+02:00 (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/TRA/traldf_iso_zps.F90
r22 r63 14 14 USE oce ! ocean dynamics and active tracers 15 15 USE dom_oce ! ocean space and time domain 16 USE ldftra_oce ! ocean active tracers: lateral physics17 USE trdtra_oce ! ocean active tracers: trend16 USE ldftra_oce ! ocean active tracers: lateral physics 17 USE trdtra_oce ! ocean active tracers: trend 18 18 USE zdf_oce ! ocean vertical physics 19 19 USE in_out_manager ! I/O manager 20 20 USE ldfslp ! iso-neutral slopes 21 USE lbclnk 21 USE ptr ! poleward transport diagnostics 22 22 23 23 24 IMPLICIT NONE … … 90 91 REAL(wp) :: & 91 92 zabe1, zabe2, zcof1, zcof2, & ! temporary scalars 92 zmsku, zmskv, zbtr, zta, zsa, & 93 zcg1, zcg2, zuwk, zvwk, & 94 zuwk1, zvwk1, & 93 zmsku, zmskv, zbtr, zta, zsa, & ! " " 95 94 ztagu, ztagv, zsagu, zsagv 96 REAL(wp), DIMENSION(jpi,jpj) :: & 97 zdkt , zdk1t, zftug, zftvg, & ! temporary workspace 98 zdks , zdk1s, zfsug, zfsvg ! " " 95 #if defined key_traldf_eiv 96 REAL(wp) :: & 97 zcg1, zcg2, zuwk, zvwk, & ! temporary scalars 98 zuwk1, zvwk1 ! " " 99 #endif 100 REAL(wp), DIMENSION(jpi,jpj) :: & ! temporary workspace 101 zdkt , zdk1t, zdks , zdk1s ! " " 102 103 #if defined key_traldf_eiv 104 REAL(wp), DIMENSION(jpi,jpj) :: & ! temporary workspace 105 zftug, zftvg, zfsug, zfsvg ! " " 106 #endif 99 107 REAL(wp), DIMENSION(jpi,jpj,jpk) :: & 100 zftv, zgtbu, zgtbv, & ! temporary workspace101 zfsv, zgsbu, zgsbv ! " "108 zftv, zgtbu, zgtbv, & ! temporary workspace 109 zfsv, zgsbu, zgsbv ! " " 102 110 !!---------------------------------------------------------------------- 103 111 … … 192 200 END DO 193 201 194 ! ! ---------------------------------------! 195 IF( lk_traldf_eiv ) THEN ! Eddy induced vertical advective fluxes ! 196 ! ! ---------------------------------------! 202 #if defined key_traldf_eiv 203 ! ---------------------------------------! 204 ! Eddy induced vertical advective fluxes ! 205 ! ---------------------------------------! 197 206 DO jj = 1, jpjm1 198 207 DO ji = 1, fs_jpim1 ! vector opt. … … 220 229 END DO 221 230 END DO 222 ENDIF 231 #endif 223 232 224 233 ! II.4 Second derivative (divergence) and add to the general trend … … 238 247 END DO 239 248 END DO 240 #if defined key_trdtra || defined key_trdmld 241 IF( lk_traldf_eiv ) THEN 249 #if defined key_trdtra || defined key_trdmld && key_traldf_eiv 242 250 DO jj = 2 , jpjm1 243 251 DO ji = fs_2, fs_jpim1 ! vector opt. … … 255 263 END DO 256 264 END DO 257 ENDIF258 265 #endif 259 266 ! ! =============== … … 261 268 ! ! =============== 262 269 263 IF( l_ctl .AND. lwp) THEN ! print mean trends (used for debugging)270 IF( l_ctl ) THEN ! print mean trends (used for debugging) 264 271 zta = SUM( ta(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) 265 272 zsa = SUM( sa(2:jpim1,2:jpjm1,1:jpkm1) * tmask(2:jpim1,2:jpjm1,1:jpkm1) ) … … 273 280 IF( MOD( kt, nf_ptr ) == 0 ) THEN 274 281 ! "zonal" mean lateral diffusive heat and salt transports 275 pht_ldf(: ,:) = prt_vj( zftv(:,:,:) )276 pst_ldf(: ,:) = prt_vj( zfsv(:,:,:) )282 pht_ldf(:) = prt_vj( zftv(:,:,:) ) 283 pst_ldf(:) = prt_vj( zfsv(:,:,:) ) 277 284 ! "zonal" mean lateral eddy induced velocity heat and salt transports 278 pht_eiv(: ,:) = prt_vj( zftv(:,:,:) )279 pst_eiv(: ,:) = prt_vj( zfsv(:,:,:) )285 pht_eiv(:) = prt_vj( zftv(:,:,:) ) 286 pst_eiv(:) = prt_vj( zfsv(:,:,:) ) 280 287 ENDIF 281 288 #endif
Note: See TracChangeset
for help on using the changeset viewer.