- Timestamp:
- 2017-12-14T11:29:10+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/TOOLS/DOMAINcfg/src/domzgr.f90
r8593 r9033 1957 1957 WRITE(numout,*) ' Bottom cell (Zb) (m) = H*rn_zb_a + rn_zb_b' 1958 1958 ENDIF 1959 1959 ! CEOD 1960 ! Add on some reference level to make nagative bathy positive to create 1961 ! sensible e3t_0's 1962 ! Can used commented code here to force bathy near open bdy not to have any WAD 1963 ! points 1964 1965 !CEOD MAKE points on the perimeter assuming a bondary is here on shallower than 1966 !10m so that tide can be specified here. 1967 ! DO jj = 1, jpj 1968 ! DO ji = 1, jpi 1969 ! !CEODIF( mig(ji) <= 2 .or. mjg(jj) <= 2 .or. mig(ji) >= jpiglo-2 .or. mjg(jj) >=jpjglo-2) then 1970 ! IF( mig(ji) <= 10 .or. mjg(jj) <=10 .or. mjg(ji) >=jpiglo-10+1 .or. mjg(jj) >=jpjglo-10+1 ) then 1971 ! IF( bathy(ji,jj) + rn_wd_ref_depth > 0.0) then 1972 ! bathy(ji,jj) = max(5.0,bathy(ji,jj)) ! 1973 ! ENDIF 1974 ! IF( bathy(ji,jj) + rn_wd_ref_depth <= 0.0) then 1975 ! bathy(ji,jj) = -rn_wd_ref_depth*2._wp ! Make sure it is land 1976 ! ENDIF 1977 ! ENDIF 1978 ! END DO 1979 ! END DO 1980 1981 IF( ln_wd) THEN 1982 bathy(:,:) = bathy(:,:) + rn_wd_ref_depth ! Arbitrary reference level for WAD, This needs 1983 ELSE 1984 bathy(:,:) = bathy(:,:) ! Arbitrary reference level for WAD, This needs 1985 ENDIF 1986 ! to be accounted for in istate in the actual model and bdys etc. 1960 1987 hift(:,:) = rn_sbot_min ! set the minimum depth for the s-coordinate 1961 1988 hifu(:,:) = rn_sbot_min
Note: See TracChangeset
for help on using the changeset viewer.