Changeset 32 for trunk/NEMO/OPA_SRC/TRA/traadv_cen2.F90
- Timestamp:
- 2004-02-17T10:20:15+01:00 (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/TRA/traadv_cen2.F90
r3 r32 63 63 !! Part I : horizontal advection 64 64 !! * centered flux: 65 !! * s-coordinate ( 'key_s_coord') or66 !! * z-coordinate with partial steps ( 'key_partial_steps'),65 !! * s-coordinate (lk_sco=T) or 66 !! * z-coordinate with partial steps (lk_zps=T), 67 67 !! the vertical scale factors e3. are inside the derivatives: 68 68 !! zcenu = e2u*e3u un mi(tn) … … 72 72 !! zcenv = e1v vn mj(tn) 73 73 !! * upstream flux: 74 !! * s-coordinate ( 'key_s_coord') or75 !! * z-coordinate with partial steps ( 'key_partial_steps')74 !! * s-coordinate (lk_sco=T) or 75 !! * z-coordinate with partial steps (lk_zps=T) 76 76 !! zupsu = e2u*e3u un (tb(i) or tb(i-1) ) [un>0 or <0] 77 77 !! zupsv = e1v*e3v vn (tb(j) or tb(j-1) ) [vn>0 or <0] … … 85 85 !! zwy = zcofj * zupsv + (1-zcofj) * zcenv 86 86 !! * horizontal advective trend (divergence of the fluxes) 87 !! * s-coordinate ( 'key_s_coord') or88 !! * z-coordinate with partial steps ( 'key_partial_steps')87 !! * s-coordinate (lk_sco=T) or 88 !! * z-coordinate with partial steps (lk_zps=T) 89 89 !! zta = 1/(e1t*e2t*e3t) { di-1[zwx] + dj-1[zwy] } 90 90 !! * z-coordinate (default key), e3t=e3u=e3v: … … 180 180 ! Advective bottom boundary layer 181 181 ! ------------------------------- 182 zun(:,:,:) = un 183 zvn(:,:,:) = vn 184 zwn(:,:,:) = wn 182 zun(:,:,:) = un(:,:,:) - u_bbl(:,:,:) 183 zvn(:,:,:) = vn(:,:,:) - v_bbl(:,:,:) 184 zwn(:,:,:) = wn(:,:,:) + w_bbl(:,:,:) 185 185 #endif 186 186
Note: See TracChangeset
for help on using the changeset viewer.