- Timestamp:
- 2016-11-28T17:04:10+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/OPA_SRC/FLO/flodom.F90
r5836 r7351 37 37 REAL(wp), ALLOCATABLE, DIMENSION(:) :: zgifl, zgjfl, zgkfl ! distances in indexes 38 38 39 !! * Substitutions40 # include "domzgr_substitute.h90"41 39 !!---------------------------------------------------------------------- 42 40 !! NEMO/OPA 3.3 , NEMO Consortium (2010) … … 176 174 ihtest(jfl) = ihtest(jfl)+1 177 175 DO jk = 1, jpk-1 178 IF( ( fsdepw(ji,jj,jk) <= flzz(jfl)) .AND. (fsdepw(ji,jj,jk+1) > flzz(jfl)) ) THEN176 IF( (gdepw_n(ji,jj,jk) <= flzz(jfl)) .AND. (gdepw_n(ji,jj,jk+1) > flzz(jfl)) ) THEN 179 177 ikmfl(jfl) = jk 180 178 ivtest(jfl) = ivtest(jfl) + 1 … … 238 236 zgifl(jfl)= (iimfl(jfl)-0.5) + zdxab/e1u(iimfl(jfl)-1,ijmfl(jfl)) + (mig(1)-jpizoom) 239 237 zgjfl(jfl)= (ijmfl(jfl)-0.5) + zdyad/e2v(iimfl(jfl),ijmfl(jfl)-1) + (mjg(1)-jpjzoom) 240 zgkfl(jfl) = (( fsdepw(iimfl(jfl),ijmfl(jfl),ikmfl(jfl)+1) - flzz(jfl) )* ikmfl(jfl)) &241 & / ( fsdepw(iimfl(jfl),ijmfl(jfl),ikmfl(jfl)+1) &242 & - fsdepw(iimfl(jfl),ijmfl(jfl),ikmfl(jfl) ) ) &243 & + (( flzz(jfl)- fsdepw(iimfl(jfl),ijmfl(jfl),ikmfl(jfl)) ) *(ikmfl(jfl)+1)) &244 & / ( fsdepw(iimfl(jfl),ijmfl(jfl),ikmfl(jfl)+1) &245 & - fsdepw(iimfl(jfl),ijmfl(jfl),ikmfl(jfl)) )238 zgkfl(jfl) = (( gdepw_n(iimfl(jfl),ijmfl(jfl),ikmfl(jfl)+1) - flzz(jfl) )* ikmfl(jfl)) & 239 & / ( gdepw_n(iimfl(jfl),ijmfl(jfl),ikmfl(jfl)+1) & 240 & - gdepw_n(iimfl(jfl),ijmfl(jfl),ikmfl(jfl) ) ) & 241 & + (( flzz(jfl)-gdepw_n(iimfl(jfl),ijmfl(jfl),ikmfl(jfl)) ) *(ikmfl(jfl)+1)) & 242 & / ( gdepw_n(iimfl(jfl),ijmfl(jfl),ikmfl(jfl)+1) & 243 & - gdepw_n(iimfl(jfl),ijmfl(jfl),ikmfl(jfl)) ) 246 244 ELSE 247 245 zgifl(jfl) = 0.e0
Note: See TracChangeset
for help on using the changeset viewer.