- Timestamp:
- 2015-10-31T08:40:45+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/SOL/solmat.F90
r5836 r5845 89 89 DO ji = 2, jpim1 90 90 zcoef = z2dt * z2dt * grav * bmask(ji,jj) 91 zcoefs = -zcoef * hv (ji ,jj-1) * e1_e2v(ji ,jj-1) ! south coefficient92 zcoefw = -zcoef * hu (ji-1,jj ) * e2_e1u(ji-1,jj ) ! west coefficient93 zcoefe = -zcoef * hu (ji ,jj ) * e2_e1u(ji ,jj ) ! east coefficient94 zcoefn = -zcoef * hv (ji ,jj ) * e1_e2v(ji ,jj ) ! north coefficient91 zcoefs = -zcoef * hv_n(ji ,jj-1) * e1_e2v(ji ,jj-1) ! south coefficient 92 zcoefw = -zcoef * hu_n(ji-1,jj ) * e2_e1u(ji-1,jj ) ! west coefficient 93 zcoefe = -zcoef * hu_n(ji ,jj ) * e2_e1u(ji ,jj ) ! east coefficient 94 zcoefn = -zcoef * hv_n(ji ,jj ) * e1_e2v(ji ,jj ) ! north coefficient 95 95 gcp(ji,jj,1) = zcoefs 96 96 gcp(ji,jj,2) = zcoefw … … 110 110 111 111 ! south coefficient 112 zcoefs = -zcoef * hv (ji,jj-1) * e1_e2v(ji,jj-1)112 zcoefs = -zcoef * hv_n(ji,jj-1) * e1_e2v(ji,jj-1) 113 113 zcoefs = zcoefs * bdyvmask(ji,jj-1) 114 114 gcp(ji,jj,1) = zcoefs 115 115 116 116 ! west coefficient 117 zcoefw = -zcoef * hu (ji-1,jj) * e2_e1u(ji-1,jj)117 zcoefw = -zcoef * hu_n(ji-1,jj) * e2_e1u(ji-1,jj) 118 118 zcoefw = zcoefw * bdyumask(ji-1,jj) 119 119 gcp(ji,jj,2) = zcoefw 120 120 121 121 ! east coefficient 122 zcoefe = -zcoef * hu (ji,jj) * e2_e1u(ji,jj)122 zcoefe = -zcoef * hu_n(ji,jj) * e2_e1u(ji,jj) 123 123 zcoefe = zcoefe * bdyumask(ji,jj) 124 124 gcp(ji,jj,3) = zcoefe 125 125 126 126 ! north coefficient 127 zcoefn = -zcoef * hv (ji,jj) * e1_e2v(ji,jj)127 zcoefn = -zcoef * hv_n(ji,jj) * e1_e2v(ji,jj) 128 128 zcoefn = zcoefn * bdyvmask(ji,jj) 129 129 gcp(ji,jj,4) = zcoefn … … 148 148 ! south coefficient 149 149 IF( ( nbondj == -1 .OR. nbondj == 2 ) .AND. ( jj == 3 ) ) THEN 150 zcoefs = -zcoef * hv (ji,jj-1) * e1_e2v(ji,jj-1)*(1.-vmask(ji,jj-1,1))150 zcoefs = -zcoef * hv_n(ji,jj-1) * e1_e2v(ji,jj-1)*(1.-vmask(ji,jj-1,1)) 151 151 ELSE 152 zcoefs = -zcoef * hv (ji,jj-1) * e1_e2v(ji,jj-1)152 zcoefs = -zcoef * hv_n(ji,jj-1) * e1_e2v(ji,jj-1) 153 153 END IF 154 154 gcp(ji,jj,1) = zcoefs … … 156 156 ! west coefficient 157 157 IF( ( nbondi == -1 .OR. nbondi == 2 ) .AND. ( ji == 3 ) ) THEN 158 zcoefw = -zcoef * hu (ji-1,jj) * e2_e1u(ji-1,jj)*(1.-umask(ji-1,jj,1))158 zcoefw = -zcoef * hu_n(ji-1,jj) * e2_e1u(ji-1,jj)*(1.-umask(ji-1,jj,1)) 159 159 ELSE 160 zcoefw = -zcoef * hu (ji-1,jj) * e2_e1u(ji-1,jj)160 zcoefw = -zcoef * hu_n(ji-1,jj) * e2_e1u(ji-1,jj) 161 161 END IF 162 162 gcp(ji,jj,2) = zcoefw … … 164 164 ! east coefficient 165 165 IF( ( nbondi == 1 .OR. nbondi == 2 ) .AND. ( ji == nlci-2 ) ) THEN 166 zcoefe = -zcoef * hu (ji,jj) * e2_e1u(ji,jj)*(1.-umask(ji,jj,1))166 zcoefe = -zcoef * hu_n(ji,jj) * e2_e1u(ji,jj)*(1.-umask(ji,jj,1)) 167 167 ELSE 168 zcoefe = -zcoef * hu (ji,jj) * e2_e1u(ji,jj)168 zcoefe = -zcoef * hu_n(ji,jj) * e2_e1u(ji,jj) 169 169 END IF 170 170 gcp(ji,jj,3) = zcoefe … … 172 172 ! north coefficient 173 173 IF( ( nbondj == 1 .OR. nbondj == 2 ) .AND. ( jj == nlcj-2 ) ) THEN 174 zcoefn = -zcoef * hv (ji,jj) * e1_e2v(ji,jj)*(1.-vmask(ji,jj,1))174 zcoefn = -zcoef * hv_n(ji,jj) * e1_e2v(ji,jj)*(1.-vmask(ji,jj,1)) 175 175 ELSE 176 zcoefn = -zcoef * hv (ji,jj) * e1_e2v(ji,jj)176 zcoefn = -zcoef * hv_n(ji,jj) * e1_e2v(ji,jj) 177 177 END IF 178 178 gcp(ji,jj,4) = zcoefn
Note: See TracChangeset
for help on using the changeset viewer.