- Timestamp:
- 2014-04-06T17:28:25+02:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_CNRS0_NOC1_LDF/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_ubs.F90
r4499 r4616 120 120 DO jj = 1, jpjm1 ! First derivative (gradient) 121 121 DO ji = 1, fs_jpim1 ! vector opt. 122 zeeu = e2 u(ji,jj) * fse3u(ji,jj,jk) / e1u(ji,jj) * umask(ji,jj,jk)123 zeev = e1 v(ji,jj) * fse3v(ji,jj,jk) / e2v(ji,jj) * vmask(ji,jj,jk)122 zeeu = e2_e1u(ji,jj) * fse3u(ji,jj,jk) * umask(ji,jj,jk) 123 zeev = e1_e2v(ji,jj) * fse3v(ji,jj,jk) * vmask(ji,jj,jk) 124 124 ztu(ji,jj,jk) = zeeu * ( ptb(ji+1,jj ,jk,jn) - ptb(ji,jj,jk,jn) ) 125 125 ztv(ji,jj,jk) = zeev * ( ptb(ji ,jj+1,jk,jn) - ptb(ji,jj,jk,jn) ) … … 164 164 DO jj = 2, jpjm1 165 165 DO ji = fs_2, fs_jpim1 ! vector opt. 166 zbtr = 1. / ( e1 t(ji,jj) *e2t(ji,jj) * fse3t(ji,jj,jk) )166 zbtr = 1. / ( e1e2t(ji,jj) * fse3t(ji,jj,jk) ) 167 167 ! horizontal advective 168 168 ztra = - zbtr * ( zwx(ji,jj,jk) - zwx(ji-1,jj ,jk) & … … 219 219 DO jj = 2, jpjm1 220 220 DO ji = fs_2, fs_jpim1 ! vector opt. 221 zbtr = 1. / ( e1 t(ji,jj) *e2t(ji,jj) * fse3t(ji,jj,jk) )221 zbtr = 1. / ( e1e2t(ji,jj) * fse3t(ji,jj,jk) ) 222 222 ztak = - ( ztw(ji,jj,jk) - ztw(ji,jj,jk+1) ) * zbtr 223 223 pta(ji,jj,jk,jn) = pta(ji,jj,jk,jn) + ztak … … 245 245 DO jj = 2, jpjm1 246 246 DO ji = fs_2, fs_jpim1 ! vector opt. 247 zbtr = 1. / ( e1 t(ji,jj) *e2t(ji,jj) * fse3t(ji,jj,jk) )247 zbtr = 1. / ( e1e2t(ji,jj) * fse3t(ji,jj,jk) ) 248 248 ! k- vertical advective trends 249 249 ztra = - zbtr * ( ztw(ji,jj,jk) - ztw(ji,jj,jk+1) ) … … 259 259 DO jj = 2, jpjm1 260 260 DO ji = fs_2, fs_jpim1 ! vector opt. 261 zbtr = 1.e0 / ( e1 t(ji,jj) *e2t(ji,jj) * fse3t(ji,jj,jk) )261 zbtr = 1.e0 / ( e1e2t(ji,jj) * fse3t(ji,jj,jk) ) 262 262 z_hdivn = ( pwn(ji,jj,jk) - pwn(ji,jj,jk+1) ) * zbtr 263 263 zltv(ji,jj,jk) = pta(ji,jj,jk,jn) - zltv(ji,jj,jk) + ptn(ji,jj,jk,jn) * z_hdivn … … 358 358 zneg = MAX( 0., pcc(ji ,jj ,jk ) ) - MIN( 0., pcc(ji ,jj ,jk+1) ) 359 359 ! up & down beta terms 360 zbt = e1 t(ji,jj) *e2t(ji,jj) * fse3t(ji,jj,jk) / z2dtt360 zbt = e1e2t(ji,jj) * fse3t(ji,jj,jk) / z2dtt 361 361 zbetup(ji,jj,jk) = ( zbetup(ji,jj,jk) - paft(ji,jj,jk) ) / (zpos+zrtrn) * zbt 362 362 zbetdo(ji,jj,jk) = ( paft(ji,jj,jk) - zbetdo(ji,jj,jk) ) / (zneg+zrtrn) * zbt
Note: See TracChangeset
for help on using the changeset viewer.