Changeset 1639 for trunk/NEMO
- Timestamp:
- 2009-10-05T19:06:17+02:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DOM/domzgr.F90
r1632 r1639 1139 1139 ! ! Define the envelop bathymetry (hbatt) 1140 1140 ! ! ============================= 1141 ! Smooth the bathymetry (if required)1142 scosrf(:,:) = 0.e0 ! ocean surface depth (here zero: no under ice-shelf sea)1143 scobot(:,:) = bathy(:,:) ! ocean bottom depth1144 !1145 1141 ! use r-value to create hybrid coordinates 1146 1142 DO jj = 1, jpj … … 1149 1145 END DO 1150 1146 END DO 1147 ! 1148 ! Smooth the bathymetry (if required) 1149 scosrf(:,:) = 0.e0 ! ocean surface depth (here zero: no under ice-shelf sea) 1150 scobot(:,:) = bathy(:,:) ! ocean bottom depth 1151 ! 1151 1152 jl = 0 1152 1153 zrmax = 1.e0 … … 1337 1338 esigw3(ji,jj,jk+1) = gsigt3(ji,jj,jk+1) - gsigt3(ji,jj,jk) 1338 1339 END DO 1339 esigw3(ji,jj,1 ) = esigw3(ji,jj,2)1340 esigt3(ji,jj,jpk) = esigt3(ji,jj,jpkm1)1340 esigw3(ji,jj,1 ) = 2.0_wp * (gsigt3(ji,jj,1 ) - gsigw3(ji,jj,1 )) 1341 esigt3(ji,jj,jpk) = 2.0_wp * (gsigt3(ji,jj,jpk) - gsigw3(ji,jj,jpk)) 1341 1342 1342 1343 ! Coefficients for vertical depth as the sum of e3w scale factors … … 1351 1352 gdept (ji,jj,jk) = (scosrf(ji,jj)+(hbatt(ji,jj)-rn_hc)*gsigt3(ji,jj,jk)+rn_hc*zcoeft) 1352 1353 gdepw (ji,jj,jk) = (scosrf(ji,jj)+(hbatt(ji,jj)-rn_hc)*gsigw3(ji,jj,jk)+rn_hc*zcoefw) 1353 gdep3w(ji,jj,jk) = (scosrf(ji,jj)+(hbatt(ji,jj)-rn_hc)*gsi3w3(ji,jj,jk)+rn_hc*zcoef w)1354 gdep3w(ji,jj,jk) = (scosrf(ji,jj)+(hbatt(ji,jj)-rn_hc)*gsi3w3(ji,jj,jk)+rn_hc*zcoeft) 1354 1355 END DO 1355 1356 … … 1402 1403 esigw(jk+1) = gsigt(jk+1) - gsigt(jk) 1403 1404 END DO 1404 esigw( 1 ) = esigw( 2 )1405 esigt(jpk) = esigt(jpkm1)1405 esigw( 1 ) = 2.0_wp * (gsigt(1) - gsigw(1)) 1406 esigt(jpk) = 2.0_wp * (gsigt(jpk) - gsigw(jpk)) 1406 1407 1407 1408 !!gm original form … … 1423 1424 gdept (:,:,jk) = (scosrf(:,:)+(hbatt(:,:)-hift(:,:))*gsigt(jk)+hift(:,:)*zcoeft) 1424 1425 gdepw (:,:,jk) = (scosrf(:,:)+(hbatt(:,:)-hift(:,:))*gsigw(jk)+hift(:,:)*zcoefw) 1425 gdep3w(:,:,jk) = (scosrf(:,:)+(hbatt(:,:)-hift(:,:))*gsi3w(jk)+hift(:,:)*zcoef w)1426 gdep3w(:,:,jk) = (scosrf(:,:)+(hbatt(:,:)-hift(:,:))*gsi3w(jk)+hift(:,:)*zcoeft) 1426 1427 END DO 1427 1428 !!gm: e3uw, e3vw can be suppressed (modif in dynzdf, dynzdf_iso, zdfbfr) (save 2 3D arrays)
Note: See TracChangeset
for help on using the changeset viewer.