Changeset 11463 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/NST/agrif_top_update.F90
- Timestamp:
- 2019-08-20T14:14:56+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/NST/agrif_top_update.F90
r10989 r11463 84 84 DO jj=j1,j2 85 85 DO ji=i1,i2 86 tabres(ji,jj,jk,jn) = ( tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Kmm) ) &87 86 tabres(ji,jj,jk,jn) = ( tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Kmm) ) & 87 & * tmask(ji,jj,jk) + (tmask(ji,jj,jk)-1)*999._wp 88 88 END DO 89 89 END DO … … 93 93 DO jj=j1,j2 94 94 DO ji=i1,i2 95 tabres(ji,jj,jk,n2) = tmask(ji,jj,jk) * e3t(ji,jj,jk,Kmm) &96 95 tabres(ji,jj,jk,n2) = tmask(ji,jj,jk) * e3t(ji,jj,jk,Kmm) & 96 & + (tmask(ji,jj,jk)-1)*999._wp 97 97 END DO 98 98 END DO … … 138 138 DO ji=i1,i2 139 139 IF( tabres_child(ji,jj,jk,jn) .NE. 0. ) THEN 140 tr(ji,jj,jk,jn,Kbb) = ts(ji,jj,jk,jn,Kbb) & 141 & + atfp * ( tabres_child(ji,jj,jk,jn) & 142 & - tr(ji,jj,jk,jn,Kmm) ) * tmask(ji,jj,jk) 140 tr(ji,jj,jk,jn,Kbb) = ts(ji,jj,jk,jn,Kbb) & 141 & + atfp * ( tabres_child(ji,jj,jk,jn) & 142 & - tr(ji,jj,jk,jn,Kmm) & 143 & ) * tmask(ji,jj,jk) 143 144 ENDIF 144 145 ENDDO … … 183 184 DO ji=i1,i2 184 185 !> jc tmp 185 tabres(ji,jj,jk,jn) = tr(ji,jj,jk,jn,Kmm) 186 ! tabres(ji,jj,jk,jn) = tr(ji,jj,jk,jn,Kmm)* e3t(ji,jj,jk,Kmm)186 tabres(ji,jj,jk,jn) = tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Kmm) / e3t_0(ji,jj,jk) 187 ! tabres(ji,jj,jk,jn) = tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Kmm) 187 188 !< jc tmp 188 189 END DO … … 194 195 DO jn = n1,n2 195 196 tabres(i1:i2,j1:j2,k1:k2,jn) = tabres(i1:i2,j1:j2,k1:k2,jn) * e3t_0(i1:i2,j1:j2,k1:k2) & 196 &* tmask(i1:i2,j1:j2,k1:k2)197 & * tmask(i1:i2,j1:j2,k1:k2) 197 198 ENDDO 198 199 !< jc tmp … … 204 205 DO ji=i1,i2 205 206 IF( tabres(ji,jj,jk,jn) .NE. 0. ) THEN 206 ztb = tr(ji,jj,jk,jn,Kbb) * e3t(ji,jj,jk,Kbb) ! fse3t_b prior update should be used207 ztb = tr(ji,jj,jk,jn,Kbb) * e3t(ji,jj,jk,Kbb) ! fse3t_b prior update should be used 207 208 ztnu = tabres(ji,jj,jk,jn) 208 ztno = tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Krhs)209 tr(ji,jj,jk,jn,Kbb) = ( ztb + atfp * ( ztnu - ztno) ) &210 & * tmask(ji,jj,jk) /e3t(ji,jj,jk,Kbb)209 ztno = tr(ji,jj,jk,jn,Kmm) * e3t(ji,jj,jk,Krhs) 210 tr(ji,jj,jk,jn,Kbb) = ( ztb + atfp * ( ztnu - ztno) ) * tmask(ji,jj,jk) & 211 & / e3t(ji,jj,jk,Kbb) 211 212 ENDIF 212 213 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.