Changeset 15675

Timestamp:

2022-01-27T14:03:44+01:00 (2 years ago)

Author:

hadjt

Message:

Bugfix SBC/sbcflx.F90:

Enda emailed (27 January 2022 09:28) to say a bug had been spotted in SBC/sbcflx.F90, where the wind stress was halved.

I have removed this 0.5 factor in (my) line 250.

File:

• NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate/src/OCE/SBC/sbcflx.F90 (modified) (1 diff)

NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate/src/OCE/SBC/sbcflx.F90

@@ -248 +248 @@
             ztx = ( utau(ji-1,jj  ) + utau(ji,jj) ) * 0.5_wp * ( 2._wp - MIN( umask(ji-1,jj  ,1), umask(ji,jj,1) ) )
             zty = ( vtau(ji  ,jj-1) + vtau(ji,jj) ) * 0.5_wp * ( 2._wp - MIN( vmask(ji  ,jj-1,1), vmask(ji,jj,1) ) ) 
-            zmod = 0.5_wp * SQRT( ztx * ztx + zty * zty ) * tmask(ji,jj,1)
+            zmod = SQRT( ztx * ztx + zty * zty ) * tmask(ji,jj,1)
             taum(ji,jj) = zmod
             wndm(ji,jj) = SQRT( zmod * zcoef )  !!clem: not used?