Changeset 10970 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/FLO/flodom.F90
- Timestamp:
- 2019-05-13T14:02:19+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/FLO/flodom.F90
r10425 r10970 44 44 CONTAINS 45 45 46 SUBROUTINE flo_dom 46 SUBROUTINE flo_dom( Kmm ) 47 47 !! --------------------------------------------------------------------- 48 48 !! *** ROUTINE flo_dom *** … … 53 53 !! the longitude (degree) and the depth (m). 54 54 !!---------------------------------------------------------------------- 55 INTEGER, INTENT(in) :: Kmm ! ocean time level index 56 ! 55 57 INTEGER :: jfl ! dummy loop 56 58 INTEGER :: inum ! logical unit for file read … … 94 96 CALL flo_add_new_ariane_floats(jpnrstflo+1,jpnfl) 95 97 ELSE !Add new floats with long/lat convention 96 CALL flo_add_new_floats( jpnrstflo+1,jpnfl)98 CALL flo_add_new_floats(Kmm,jpnrstflo+1,jpnfl) 97 99 ENDIF 98 100 ENDIF … … 106 108 CALL flo_add_new_ariane_floats(1,jpnfl) 107 109 ELSE !Add new floats with long/lat convention 108 CALL flo_add_new_floats( 1,jpnfl)110 CALL flo_add_new_floats(Kmm,1,jpnfl) 109 111 ENDIF 110 112 … … 113 115 END SUBROUTINE flo_dom 114 116 115 SUBROUTINE flo_add_new_floats( kfl_start, kfl_end)117 SUBROUTINE flo_add_new_floats(Kmm, kfl_start, kfl_end) 116 118 !! ------------------------------------------------------------- 117 119 !! *** SUBROUTINE add_new_arianefloats *** … … 128 130 !! ** Method : 129 131 !!---------------------------------------------------------------------- 132 INTEGER, INTENT(in) :: Kmm 130 133 INTEGER, INTENT(in) :: kfl_start, kfl_end 131 134 !! … … 174 177 ihtest(jfl) = ihtest(jfl)+1 175 178 DO jk = 1, jpk-1 176 IF( (gdepw _n(ji,jj,jk) <= flzz(jfl)) .AND. (gdepw_n(ji,jj,jk+1) > flzz(jfl)) ) THEN179 IF( (gdepw(ji,jj,jk,Kmm) <= flzz(jfl)) .AND. (gdepw(ji,jj,jk+1,Kmm) > flzz(jfl)) ) THEN 177 180 ikmfl(jfl) = jk 178 181 ivtest(jfl) = ivtest(jfl) + 1 … … 236 239 zgifl(jfl)= (iimfl(jfl)-0.5) + zdxab/e1u(iimfl(jfl)-1,ijmfl(jfl)) + (mig(1)-1) 237 240 zgjfl(jfl)= (ijmfl(jfl)-0.5) + zdyad/e2v(iimfl(jfl),ijmfl(jfl)-1) + (mjg(1)-1) 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)) )241 zgkfl(jfl) = (( gdepw(iimfl(jfl),ijmfl(jfl),ikmfl(jfl)+1,Kmm) - flzz(jfl) )* ikmfl(jfl)) & 242 & / ( gdepw(iimfl(jfl),ijmfl(jfl),ikmfl(jfl)+1,Kmm) & 243 & - gdepw(iimfl(jfl),ijmfl(jfl),ikmfl(jfl) ,Kmm) ) & 244 & + (( flzz(jfl)-gdepw(iimfl(jfl),ijmfl(jfl),ikmfl(jfl),Kmm) ) *(ikmfl(jfl)+1)) & 245 & / ( gdepw(iimfl(jfl),ijmfl(jfl),ikmfl(jfl)+1,Kmm) & 246 & - gdepw(iimfl(jfl),ijmfl(jfl),ikmfl(jfl),Kmm) ) 244 247 ELSE 245 248 zgifl(jfl) = 0.e0
Note: See TracChangeset
for help on using the changeset viewer.