Changeset 2974 for branches/2011/dev_r2739_LOCEAN8_ZTC
- Timestamp:
- 2011-10-21T16:42:39+02:00 (13 years ago)
- Location:
- branches/2011/dev_r2739_LOCEAN8_ZTC/NEMOGCM/NEMO/OPA_SRC
- Files:
-
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_r2739_LOCEAN8_ZTC/NEMOGCM/NEMO/OPA_SRC/DOM/dom_oce.F90
r2905 r2974 157 157 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ht_0 !: refernce depth at t- points (meters) 158 158 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: hu_0 , hv_0 !: refernce depth at u- and v-points (meters) 159 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: e 1ur , e2vr!: scale factor coeffs at U--V points159 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: e2_1u , e1_2v !: scale factor coeffs at U--V points 160 160 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: e12t , e12t_1 !: horizontal cell surface and inverse at T points 161 161 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: e12u , e12u_1 !: horizontal cell surface and inverse at U points 162 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: e12v_1 , e12f_1 !: inverse horizontal cell surface at V--F points 162 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: e12v , e12v_1 !: horizontal cell surface and inverse at V points 163 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: e12f , e12f_1 !: inverse horizontal cell surface at V--F points 163 164 164 165 INTEGER, PUBLIC :: nla10 !: deepest W level Above ~10m (nlb10 - 1) … … 293 294 #endif 294 295 ! 295 ALLOCATE( hu (jpi,jpj) , hur (jpi,jpj) , hu_0 (jpi,jpj) , ht_0 (jpi,jpj) , & 296 & hv (jpi,jpj) , hvr (jpi,jpj) , hv_0 (jpi,jpj) , & 297 & e1ur(jpi,jpj) , e2vr (jpi,jpj) , e12t (jpi,jpj) , e12t_1(jpi,jpj) , & 298 & e12u(jpi,jpj) , e12u_1(jpi,jpj) , e12v_1(jpi,jpj) , e12f_1(jpi,jpj) , STAT=ierr(6) ) 296 ALLOCATE( hu (jpi,jpj) , hur (jpi,jpj) , hu_0(jpi,jpj) , ht_0 (jpi,jpj) , & 297 & hv (jpi,jpj) , hvr (jpi,jpj) , hv_0(jpi,jpj) , & 298 & e2_1u(jpi,jpj) , e1_2v (jpi,jpj) , e12t(jpi,jpj) , e12t_1(jpi,jpj) , & 299 & e12u (jpi,jpj) , e12u_1(jpi,jpj) , e12v(jpi,jpj) , e12v_1(jpi,jpj) , & 300 & e12f (jpi,jpj) , e12f_1(jpi,jpj) , STAT=ierr(6) ) 299 301 ! 300 302 ALLOCATE( gdept_0(jpk) , gdepw_0(jpk) , & -
branches/2011/dev_r2739_LOCEAN8_ZTC/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90
r2917 r2974 90 90 END IF 91 91 ! 92 ! - ML - Used in domvvl and traldf_(lab/bilap/iso) but could be usefull in many other modules 92 ! - ML - Used in dom_vvl_sf_nxt and lateral diffusion routines 93 ! but could be usefull in many other routines 93 94 e12t (:,:) = e1t(:,:) * e2t(:,:) 94 95 e12u (:,:) = e1u(:,:) * e2u(:,:) 95 e12u_1(:,:) = 0.5 / e12u(:,:) 96 e12v_1(:,:) = 0.5 / ( e1v(:,:) * e2v(:,:) ) 97 e12f_1(:,:) = 0.5 / ( e1f(:,:) * e2f(:,:) ) 98 e1ur (:,:) = e2u(:,:) / e1u(:,:) 99 e2vr (:,:) = e1v(:,:) / e2v(:,:) 96 e12v (:,:) = e1v(:,:) * e2v(:,:) 97 e12f (:,:) = e1f(:,:) * e2f(:,:) 98 e12t_1(:,:) = 1. / e12t(:,:) 99 e12u_1(:,:) = 1. / e12u(:,:) 100 e12v_1(:,:) = 1. / e12v(:,:) 101 e12f_1(:,:) = 1. / e12f(:,:) 102 e2_1u (:,:) = e2u(:,:) / e1u(:,:) 103 e1_2v (:,:) = e1v(:,:) / e2v(:,:) 100 104 ! 101 105 hu(:,:) = 0.e0 ! Ocean depth at U- and V-points -
branches/2011/dev_r2739_LOCEAN8_ZTC/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90
r2970 r2974 308 308 DO jj = 1, jpjm1 309 309 DO ji = 1, fs_jpim1 ! vector opt. 310 un_td(ji,jj,jk) = ahe3 * umask(ji,jj,jk) * e 1ur(ji,jj) * ( e3t_t_b(ji,jj,jk) - e3t_t_b(ji+1,jj ,jk) )311 vn_td(ji,jj,jk) = ahe3 * vmask(ji,jj,jk) * e 2vr(ji,jj) * ( e3t_t_b(ji,jj,jk) - e3t_t_b(ji ,jj+1,jk) )310 un_td(ji,jj,jk) = ahe3 * umask(ji,jj,jk) * e2_1u(ji,jj) * ( e3t_t_b(ji,jj,jk) - e3t_t_b(ji+1,jj ,jk) ) 311 vn_td(ji,jj,jk) = ahe3 * vmask(ji,jj,jk) * e1_2v(ji,jj) * ( e3t_t_b(ji,jj,jk) - e3t_t_b(ji ,jj+1,jk) ) 312 312 zwu(ji,jj) = zwu(ji,jj) + un_td(ji,jj,jk) 313 313 zwv(ji,jj) = zwv(ji,jj) + vn_td(ji,jj,jk) … … 332 332 END DO 333 333 END DO 334 ! d - thickness diffusion equivalenttransport: boundary conditions334 ! d - thickness diffusion transport: boundary conditions 335 335 ! (stored for tracer advction and continuity equation) 336 336 CALL lbc_lnk( un_td , 'U' , -1.) … … 557 557 DO jj = 2, jpjm1 558 558 DO ji = 1, fs_jpim1 ! vector opt. 559 pe3_out(ji,jj,jk) = umask(ji,jj,jk) * e12u_1(ji,jj)&560 & * ( e12t(ji ,jj) * ( pe3_in(ji ,jj,jk) - fse3t_0(ji ,jj,jk) ) &561 & + e12t(ji+1,jj) * ( pe3_in(ji+1,jj,jk) - fse3t_0(ji+1,jj,jk) ) )559 pe3_out(ji,jj,jk) = 0.5 * umask(ji,jj,jk) * e12u_1(ji,jj) & 560 & * ( e12t(ji ,jj) * ( pe3_in(ji ,jj,jk) - fse3t_0(ji ,jj,jk) ) & 561 & + e12t(ji+1,jj) * ( pe3_in(ji+1,jj,jk) - fse3t_0(ji+1,jj,jk) ) ) 562 562 END DO 563 563 END DO … … 573 573 DO jj = 1, jpjm1 574 574 DO ji = fs_2, fs_jpim1 ! vector opt. 575 pe3_out(ji,jj,jk) = umask(ji,jj,jk) * e12v_1(ji,jj)&576 & * ( e12t(ji,jj ) * ( pe3_in(ji,jj ,jk) - fse3t_0(ji,jj ,jk) ) &577 & + e12t(ji,jj+1) * ( pe3_in(ji,jj+1,jk) - fse3t_0(ji,jj+1,jk) ) )575 pe3_out(ji,jj,jk) = 0.5 * vmask(ji,jj,jk) * e12v_1(ji,jj) & 576 & * ( e12t(ji,jj ) * ( pe3_in(ji,jj ,jk) - fse3t_0(ji,jj ,jk) ) & 577 & + e12t(ji,jj+1) * ( pe3_in(ji,jj+1,jk) - fse3t_0(ji,jj+1,jk) ) ) 578 578 END DO 579 579 END DO … … 589 589 DO jj = 1, jpjm1 590 590 DO ji = 1, fs_jpim1 ! vector opt. 591 pe3_out(ji,jj,jk) = umask(ji,jj,jk) * umask(ji,jj+1,jk) * e12f_1(ji,jj)&592 & * ( e12u(ji,jj ) * ( pe3_in(ji,jj ,jk) - fse3u_0(ji,jj ,jk) ) &593 & + e12u(ji,jj+1) * ( pe3_in(ji,jj+1,jk) - fse3u_0(ji,jj+1,jk) ) )591 pe3_out(ji,jj,jk) = 0.5 * umask(ji,jj,jk) * umask(ji,jj+1,jk) * e12f_1(ji,jj) & 592 & * ( e12u(ji,jj ) * ( pe3_in(ji,jj ,jk) - fse3u_0(ji,jj ,jk) ) & 593 & + e12u(ji,jj+1) * ( pe3_in(ji,jj+1,jk) - fse3u_0(ji,jj+1,jk) ) ) 594 594 END DO 595 595 END DO -
branches/2011/dev_r2739_LOCEAN8_ZTC/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_bilap.F90
r2905 r2974 107 107 DO jj = 1, jpjm1 108 108 DO ji = 1, fs_jpim1 ! vector opt. 109 zeeu(ji,jj) = e 1ur(ji,jj) * fse3u_n(ji,jj,jk) * umask(ji,jj,jk)110 zeev(ji,jj) = e 2vr(ji,jj) * fse3v_n(ji,jj,jk) * vmask(ji,jj,jk)109 zeeu(ji,jj) = e2_1u(ji,jj) * fse3u_n(ji,jj,jk) * umask(ji,jj,jk) 110 zeev(ji,jj) = e1_2v(ji,jj) * fse3v_n(ji,jj,jk) * vmask(ji,jj,jk) 111 111 END DO 112 112 END DO … … 130 130 DO jj = 2, jpjm1 ! Second derivative (divergence) time the eddy diffusivity coefficient 131 131 DO ji = fs_2, fs_jpim1 ! vector opt. 132 zbtr = 1.0 / ( e1 t(ji,jj) * e2t(ji,jj) * fse3t_n(ji,jj,jk) )132 zbtr = 1.0 / ( e12t(ji,jj) * fse3t_n(ji,jj,jk) ) 133 133 zlt(ji,jj) = fsahtt(ji,jj,jk) * zbtr * ( ztu(ji,jj,jk) - ztu(ji-1,jj,jk) & 134 134 & + ztv(ji,jj,jk) - ztv(ji,jj-1,jk) ) … … 148 148 DO ji = fs_2, fs_jpim1 ! vector opt. 149 149 ! horizontal diffusive trends 150 zbtr = 1.0 / ( e1 t(ji,jj) * e2t(ji,jj) * fse3t_n(ji,jj,jk) )150 zbtr = 1.0 / ( e12t(ji,jj) * fse3t_n(ji,jj,jk) ) 151 151 ztra = zbtr * ( ztu(ji,jj,jk) - ztu(ji-1,jj,jk) + ztv(ji,jj,jk) - ztv(ji,jj-1,jk) ) 152 152 ! add it to the general tracer trends -
branches/2011/dev_r2739_LOCEAN8_ZTC/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_bilapg.F90
r2715 r2974 207 207 DO jj = 1, jpjm1 208 208 DO ji = 1, jpim1 209 zabe1 = e2 u(ji,jj) * fse3u(ji,jj,jk) / e1u(ji,jj)210 zabe2 = e1 v(ji,jj) * fse3v(ji,jj,jk) / e2v(ji,jj)209 zabe1 = e2_1u(ji,jj) * fse3u_n(ji,jj,jk) 210 zabe2 = e1_2v(ji,jj) * fse3v_n(ji,jj,jk) 211 211 212 212 zmku = 1./MAX( tmask(ji+1,jj,jk )+tmask(ji,jj,jk+1) & … … 275 275 DO jk = 2, jpkm1 276 276 DO ji = 2, jpim1 277 zcof0 = e1 t(ji,jj) * e2t(ji,jj) / fse3w(ji,jj,jk) &277 zcof0 = e12t(ji,jj) / fse3w_n(ji,jj,jk) & 278 278 & * ( wslpi(ji,jj,jk) * wslpi(ji,jj,jk) & 279 279 & + wslpj(ji,jj,jk) * wslpj(ji,jj,jk) ) … … 306 306 DO ji = 2, jpim1 307 307 ! eddy coef. divided by the volume element 308 zbtr = 1.0 / ( e1 t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) )308 zbtr = 1.0 / ( e12t(ji,jj) * fse3t_n(ji,jj,jk) ) 309 309 ! vertical divergence 310 310 ztav = fsahtt(ji,jj,jk) * ( zftw(ji,jk) - zftw(ji,jk+1) ) … … 318 318 DO ji = 2, jpim1 319 319 ! inverse of the volume element 320 zbtr = 1.0 / ( e1 t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) )320 zbtr = 1.0 / ( e12t(ji,jj) * fse3t_n(ji,jj,jk) ) 321 321 ! vertical divergence 322 322 ztav = zftw(ji,jk) - zftw(ji,jk+1) -
branches/2011/dev_r2739_LOCEAN8_ZTC/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_iso.F90
r2905 r2974 172 172 DO jj = 1 , jpjm1 173 173 DO ji = 1, fs_jpim1 ! vector opt. 174 zabe1 = ( fsahtu(ji,jj,jk) + pahtb0 ) * e 1ur(ji,jj) * fse3u_n(ji,jj,jk)175 zabe2 = ( fsahtv(ji,jj,jk) + pahtb0 ) * e 2vr(ji,jj) * fse3v_n(ji,jj,jk)174 zabe1 = ( fsahtu(ji,jj,jk) + pahtb0 ) * e2_1u(ji,jj) * fse3u_n(ji,jj,jk) 175 zabe2 = ( fsahtv(ji,jj,jk) + pahtb0 ) * e1_2v(ji,jj) * fse3v_n(ji,jj,jk) 176 176 ! 177 177 zmsku = 1. / MAX( tmask(ji+1,jj,jk ) + tmask(ji,jj,jk+1) & … … 197 197 DO jj = 2 , jpjm1 198 198 DO ji = fs_2, fs_jpim1 ! vector opt. 199 zbtr = 1.0 / ( e1 t(ji,jj) * e2t(ji,jj) * fse3t_n(ji,jj,jk) )199 zbtr = 1.0 / ( e12t(ji,jj) * fse3t_n(ji,jj,jk) ) 200 200 ztra = zbtr * ( zftu(ji,jj,jk) - zftu(ji-1,jj,jk) + zftv(ji,jj,jk) - zftv(ji,jj-1,jk) ) 201 201 pta(ji,jj,jk,jn) = pta(ji,jj,jk,jn) + ztra … … 282 282 DO jj = 2, jpjm1 283 283 DO ji = fs_2, fs_jpim1 ! vector opt. 284 zbtr = 1.0 / ( e1 t(ji,jj) * e2t(ji,jj) * fse3t_n(ji,jj,jk) )284 zbtr = 1.0 / ( e12t(ji,jj) * fse3t_n(ji,jj,jk) ) 285 285 ztra = ( ztfw(ji,jj,jk) - ztfw(ji,jj,jk+1) ) * zbtr 286 286 pta(ji,jj,jk,jn) = pta(ji,jj,jk,jn) + ztra -
branches/2011/dev_r2739_LOCEAN8_ZTC/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_iso_grif.F90
r2715 r2974 159 159 DO jj = 1, jpjm1 160 160 DO ji = 1, fs_jpim1 161 ze3wr = 1._wp / fse3w (ji+ip,jj,jk+kp)162 zbu = 0.25_wp * e1 u(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,jk)161 ze3wr = 1._wp / fse3w_n(ji+ip,jj,jk+kp) 162 zbu = 0.25_wp * e12u(ji,jj) * fse3u_n(ji,jj,jk) 163 163 zah = fsahtu(ji,jj,jk) ! fsaht(ji+ip,jj,jk) 164 164 zslope_skew = triadi_g(ji+ip,jj,jk,1-ip,kp) 165 zslope2 = zslope_skew - ( fsdept (ji+1,jj,jk) - fsdept(ji ,jj ,jk) ) * ze1ur * umask(ji,jj,jk+kp)165 zslope2 = zslope_skew - ( fsdept_n(ji+1,jj,jk) - fsdept_n(ji ,jj ,jk) ) * ze1ur * umask(ji,jj,jk+kp) 166 166 zslope2 = zslope2 *zslope2 167 167 ah_wslp2(ji+ip,jj,jk+kp) = ah_wslp2(ji+ip,jj,jk+kp) & 168 & + zah * ( zbu * ze3wr / ( e1 t(ji+ip,jj) * e2t(ji+ip,jj) ) ) * zslope2168 & + zah * ( zbu * ze3wr / ( e12t(ji+ip,jj) ) ) * zslope2 169 169 IF( ln_traldf_gdia ) THEN 170 170 zaei_slp = fsaeiw(ji+ip,jj,jk) * zslope_skew !fsaeit(ji+ip,jj,jk)*zslope_skew … … 182 182 DO jj = 1, jpjm1 183 183 DO ji=1,fs_jpim1 184 ze3wr = 1.0_wp / fse3w (ji,jj+jp,jk+kp)185 zbv = 0.25_wp * e1 v(ji,jj) * e2v(ji,jj) * fse3v(ji,jj,jk)184 ze3wr = 1.0_wp / fse3w_n(ji,jj+jp,jk+kp) 185 zbv = 0.25_wp * e12v(ji,jj) * fse3v_n(ji,jj,jk) 186 186 zah = fsahtu(ji,jj,jk) !fsaht(ji,jj+jp,jk) 187 187 zslope_skew = triadj_g(ji,jj+jp,jk,1-jp,kp) 188 zslope2 = zslope_skew - ( fsdept (ji,jj+1,jk) - fsdept(ji,jj,jk) ) * ze2vr * vmask(ji,jj,jk+kp)188 zslope2 = zslope_skew - ( fsdept_n(ji,jj+1,jk) - fsdept_n(ji,jj,jk) ) * ze2vr * vmask(ji,jj,jk+kp) 189 189 zslope2 = zslope2 * zslope2 190 190 ah_wslp2(ji,jj+jp,jk+kp) = ah_wslp2(ji,jj+jp,jk+kp) & 191 & + zah * ( zbv * ze3wr / ( e1 t(ji,jj+jp) * e2t(ji,jj+jp) ) ) * zslope2191 & + zah * ( zbv * ze3wr / ( e12t(ji,jj+jp) ) ) * zslope2 192 192 IF( ln_traldf_gdia ) THEN 193 193 zaei_slp = fsaeiw(ji,jj+jp,jk) * zslope_skew !fsaeit(ji,jj+jp,jk)*zslope_skew … … 255 255 ze1ur = 1._wp / e1u(ji,jj) 256 256 zdxt = zdit(ji,jj,jk) * ze1ur 257 ze3wr = 1._wp / fse3w (ji+ip,jj,jk+kp)257 ze3wr = 1._wp / fse3w_n(ji+ip,jj,jk+kp) 258 258 zdzt = zdkt(ji+ip,jj,kp) * ze3wr 259 259 zslope_skew = triadi_g(ji+ip,jj,jk,1-ip,kp) 260 260 zslope_iso = triadi(ji+ip,jj,jk,1-ip,kp) 261 261 262 zbu = 0.25_wp * e1 u(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,jk)262 zbu = 0.25_wp * e12u(ji,jj) * fse3u_n(ji,jj,jk) 263 263 zah = fsahtu(ji,jj,jk) !*umask(ji,jj,jk+kp) !fsaht(ji+ip,jj,jk) ===>> ???? 264 264 zah_slp = zah * zslope_iso … … 277 277 ze2vr = 1._wp / e2v(ji,jj) 278 278 zdyt = zdjt(ji,jj,jk) * ze2vr 279 ze3wr = 1._wp / fse3w (ji,jj+jp,jk+kp)279 ze3wr = 1._wp / fse3w_n(ji,jj+jp,jk+kp) 280 280 zdzt = zdkt(ji,jj+jp,kp) * ze3wr 281 281 zslope_skew = triadj_g(ji,jj+jp,jk,1-jp,kp) 282 282 zslope_iso = triadj(ji,jj+jp,jk,1-jp,kp) 283 zbv = 0.25_wp * e1 v(ji,jj) * e2v(ji,jj) * fse3v(ji,jj,jk)283 zbv = 0.25_wp * e12v(ji,jj) * fse3v_n(ji,jj,jk) 284 284 zah = fsahtv(ji,jj,jk) !*vmask(ji,jj,jk+kp) !fsaht(ji,jj+jp,jk) 285 285 zah_slp = zah * zslope_iso … … 295 295 DO jj = 2 , jpjm1 296 296 DO ji = fs_2, fs_jpim1 ! vector opt. 297 zbtr = 1._wp / ( e1 t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) )297 zbtr = 1._wp / ( e12t(ji,jj) * fse3t_n(ji,jj,jk) ) 298 298 pta(ji,jj,jk,jn) = pta(ji,jj,jk,jn) + zbtr * ( zftu(ji-1,jj,jk) - zftu(ji,jj,jk) & 299 299 & + zftv(ji,jj-1,jk) - zftv(ji,jj,jk) ) … … 307 307 DO ji = fs_2, fs_jpim1 ! vector opt. 308 308 pta(ji,jj,jk,jn) = pta(ji,jj,jk,jn) + ( ztfw(ji,jj,jk+1) - ztfw(ji,jj,jk) ) & 309 & / ( e1 t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) )309 & / ( e12t(ji,jj) * fse3t_n(ji,jj,jk) ) 310 310 END DO 311 311 END DO -
branches/2011/dev_r2739_LOCEAN8_ZTC/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_lap.F90
r2905 r2974 90 90 DO jj = 1, jpjm1 91 91 DO ji = 1, fs_jpim1 ! vector opt. 92 zabe1 = fsahtu(ji,jj,jk) * umask(ji,jj,jk) * e 1ur(ji,jj) * fse3u_n(ji,jj,jk)93 zabe2 = fsahtv(ji,jj,jk) * vmask(ji,jj,jk) * e 2vr(ji,jj) * fse3v_n(ji,jj,jk)92 zabe1 = fsahtu(ji,jj,jk) * umask(ji,jj,jk) * e2_1u(ji,jj) * fse3u_n(ji,jj,jk) 93 zabe2 = fsahtv(ji,jj,jk) * vmask(ji,jj,jk) * e1_2v(ji,jj) * fse3v_n(ji,jj,jk) 94 94 ztu(ji,jj,jk) = zabe1 * ( ptb(ji+1,jj ,jk,jn) - ptb(ji,jj,jk,jn) ) 95 95 ztv(ji,jj,jk) = zabe2 * ( ptb(ji ,jj+1,jk,jn) - ptb(ji,jj,jk,jn) ) … … 103 103 ikv = mbkv(ji,jj) 104 104 IF( iku == jk ) THEN 105 zabe1 = fsahtu(ji,jj,iku) * umask(ji,jj,iku) * e 1ur(ji,jj) * fse3u_n(ji,jj,iku)105 zabe1 = fsahtu(ji,jj,iku) * umask(ji,jj,iku) * e2_1u(ji,jj) * fse3u_n(ji,jj,iku) 106 106 ztu(ji,jj,jk) = zabe1 * pgu(ji,jj,jn) 107 107 ENDIF 108 108 IF( ikv == jk ) THEN 109 zabe2 = fsahtv(ji,jj,ikv) * vmask(ji,jj,ikv) * e 2vr(ji,jj) * fse3v_n(ji,jj,ikv)109 zabe2 = fsahtv(ji,jj,ikv) * vmask(ji,jj,ikv) * e1_2v(ji,jj) * fse3v_n(ji,jj,ikv) 110 110 ztv(ji,jj,jk) = zabe2 * pgv(ji,jj,jn) 111 111 ENDIF … … 119 119 DO jj = 2, jpjm1 120 120 DO ji = fs_2, fs_jpim1 ! vector opt. 121 zbtr = 1._wp / ( e1 t(ji,jj) *e2t(ji,jj) * fse3t_n(ji,jj,jk) )121 zbtr = 1._wp / ( e12t(ji,jj) * fse3t_n(ji,jj,jk) ) 122 122 ! horizontal diffusive trends added to the general tracer trends 123 123 pta(ji,jj,jk,jn) = pta(ji,jj,jk,jn) + zbtr * ( ztu(ji,jj,jk) - ztu(ji-1,jj,jk) &
Note: See TracChangeset
for help on using the changeset viewer.