Changeset 8932

Timestamp:

2017-12-07T11:57:04+01:00 (7 years ago)

Author:

jcastill

Message:

Fix relative wind calculation: the winds are received on the U,V grid, not on the T grid

File:

• branches/UKMO/CO6_KD490/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx.F90 (modified) (1 diff)

branches/UKMO/CO6_KD490/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx.F90

@@ -157 +157 @@
 
             IF( ln_rel_wind ) THEN
-               DO jj = 2, jpjm1
-                  DO ji = fs_2, fs_jpim1   ! vect. opt.
-                     zwnd_i(ji,jj) = ( sf(jp_utau)%fnow(ji,jj,1) - rn_wfac * 0.5 * ( ssu_m(ji-1,jj  ) + ssu_m(ji,jj) ))
-                     zwnd_j(ji,jj) = ( sf(jp_vtau)%fnow(ji,jj,1) - rn_wfac * 0.5 * ( ssv_m(ji  ,jj-1) + ssv_m(ji,jj) ))
+               DO jj = 1, jpj
+                  DO ji = 1, jpi
+                     zwnd_i(ji,jj) = sf(jp_utau)%fnow(ji,jj,1) - rn_wfac * ssu_m(ji,jj)
+                     zwnd_j(ji,jj) = sf(jp_vtau)%fnow(ji,jj,1) - rn_wfac * ssv_m(ji,jj)
                   END DO
                END DO
-               CALL lbc_lnk( zwnd_i(:,:) , 'T', -1. )
-               CALL lbc_lnk( zwnd_j(:,:) , 'T', -1. )
             ELSE
                zwnd_i(:,:) = sf(jp_utau)%fnow(:,:,1)