Changeset 9778 for NEMO/releases/release-3.6/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.F90

Timestamp:

2018-06-12T11:26:43+02:00 (6 years ago)

Author:

mathiot

Message:

v3.6 correction of ticket #2096

File:

• NEMO/releases/release-3.6/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.F90 (modified) (1 diff)

NEMO/releases/release-3.6/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.F90

@@ -407 +407 @@
             DO ji = fs_2, fs_jpim1
                !                                            ! ocean area at u- & v-points
-               zfrldu  = 0.5_wp * ( frld(ji,jj) + frld(ji+1,jj) )
-               zfrldv  = 0.5_wp * ( frld(ji,jj) + frld(ji,jj+1) )
+               zfrldu  = ( frld(ji,jj) * tmask(ji,jj) + frld(ji+1,jj) * tmask(ji+1,jj  ) ) / MAX(1.0_wp,tmask(ji,jj)+tmask(ji+1,jj  ))
+               zfrldv  = ( frld(ji,jj) * tmask(ji,jj) + frld(ji,jj+1) * tmask(ji  ,jj+1) ) / MAX(1.0_wp,tmask(ji,jj)+tmask(ji  ,jj+1))
                !                                            ! quadratic drag formulation without rotation
                !                                            ! using instantaneous surface ocean current