Changeset 911 for trunk/NEMO/OPA_SRC/DYN/wzvmod.F90
- Timestamp:
- 2008-04-28T11:31:32+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DYN/wzvmod.F90
r888 r911 7 7 !! 7.0 ! 96-01 (G. Madec) Statement function for e3 8 8 !! 8.5 ! 02-07 (G. Madec) Free form, F90 9 !! " ! 07-07 (D. Storkey) Zero zhdiv at open boundary (BDY) 9 10 !!---------------------------------------------------------------------- 10 11 !! wzv : Compute the vertical velocity … … 18 19 USE prtctl ! Print control 19 20 USE phycst 21 USE bdy_oce ! unstructured open boundaries 20 22 USE lbclnk ! ocean lateral boundary condition (or mpp link) 21 23 … … 54 56 55 57 !! * Local declarations 58 INTEGER :: jgrd, jb ! temporary scalars 56 59 INTEGER :: jk ! dummy loop indices 57 60 !! Variable volume … … 88 91 ! Horizontal divergence of barotropic transports 89 92 !-------------------------------------------------- 93 zhdiv(:,:) = 0.e0 90 94 DO jj = 2, jpjm1 91 95 DO ji = 2, jpim1 ! vector opt. … … 105 109 IF( lp_obc_north ) zhdiv(nin0 :nin1 ,njn0p1:njn1p1) = 0.e0 ! north 106 110 IF( lp_obc_south ) zhdiv(nis0 :nis1 ,njs0 :njs1) = 0.e0 ! south 111 #endif 112 113 #if defined key_bdy || defined key_bdy_tides 114 jgrd=1 !: tracer grid. 115 DO jb = 1, nblenrim(jgrd) 116 ji = nbi(jb,jgrd) 117 jj = nbj(jb,jgrd) 118 zhdiv(ji, jj) = 0.e0 119 zhdiv(ji, jj) = 0.e0 120 zhdiv(ji, jj) = 0.e0 121 END DO 107 122 #endif 108 123
Note: See TracChangeset
for help on using the changeset viewer.