Changeset 1695 for trunk/NEMO/LIM_SRC_2/limsbc_2.F90

Timestamp:

2009-10-30T17:42:38+01:00 (14 years ago)

Author:

smasson

Message:

wind stress module directly at T-point, see ticket:577

File:

• trunk/NEMO/LIM_SRC_2/limsbc_2.F90 (modified) (3 diffs)

trunk/NEMO/LIM_SRC_2/limsbc_2.F90

@@ -90 +90 @@
       REAL(wp), DIMENSION(jpi,jpj,1) ::   zalb     ! albedo of ice under overcast sky
       REAL(wp), DIMENSION(jpi,jpj,1) ::   zalbp    ! albedo of ice under clear sky
-      REAL(wp) ::   zsang, zmod, zfm
-      REAL(wp), DIMENSION(jpi,jpj) ::   ztio_u, ztio_v   ! ocean stress below sea-ice
+      REAL(wp) ::   zsang, zmod, zztmp, zfm
+      REAL(wp), DIMENSION(jpi,jpj) ::   ztio_u, ztio_v   ! component of ocean stress below sea-ice at I-point
+      REAL(wp), DIMENSION(jpi,jpj) ::   ztiomi           ! module    of ocean stress below sea-ice at I-point
       REAL(wp), DIMENSION(jpi,jpj) ::   zqnsoce          ! save qns before its modification by ice model
 
@@ -249 +250 @@
                zu_io  = u_ice(ji,jj) - u_oce(ji,jj)
                zv_io  = v_ice(ji,jj) - v_oce(ji,jj)
-               zmod   = rhoco * SQRT( zu_io * zu_io + zv_io * zv_io )
-               ! ... ice stress over ocean with a ice-ocean rotation angle (at I-point)
-               ztio_u(ji,jj) = zmod * ( cangvg * zu_io - zsang * zv_io )
-               ztio_v(ji,jj) = zmod * ( cangvg * zv_io + zsang * zu_io )
+               zmod   = SQRT( zu_io * zu_io + zv_io * zv_io )
+               zztmp  = rhoco * zmod
+               ! ... components of ice stress over ocean with a ice-ocean rotation angle (at I-point)
+               ztio_u(ji,jj) = zztmp * ( cangvg * zu_io - zsang * zv_io )
+               ztio_v(ji,jj) = zztmp * ( cangvg * zv_io + zsang * zu_io )
+               ! ... module of ice stress over ocean (at I-point)
+               ztiomi(ji,jj) = zztmp * zmod
                ! 
             END DO
@@ -259 +263 @@
          DO jj = 2, jpjm1
             DO ji = 2, jpim1   ! NO vector opt.
-               ! ... ice-cover wheighted ice-ocean stress at U and V-points  (from I-point values)
+               ! ... components of ice-ocean stress at U and V-points  (from I-point values)
                zutau  = 0.5 * ( ztio_u(ji+1,jj) + ztio_u(ji+1,jj+1) )
                zvtau  = 0.5 * ( ztio_v(ji,jj+1) + ztio_v(ji+1,jj+1) )
                ! ... open-ocean (lead) fraction at U- & V-points (from T-point values)
-               zfrldu = 0.5 * ( frld (ji,jj) + frld (ji+1,jj  ) )
-               zfrldv = 0.5 * ( frld (ji,jj) + frld (ji  ,jj+1) )
-               ! update surface ocean stress
+               zfrldu = 0.5 * ( frld(ji,jj) + frld(ji+1,jj  ) )
+               zfrldv = 0.5 * ( frld(ji,jj) + frld(ji  ,jj+1) )
+               ! ... update components of surface ocean stress (ice-cover wheighted)
                utau(ji,jj) = zfrldu * utau(ji,jj) + ( 1. - zfrldu ) * zutau
                vtau(ji,jj) = zfrldv * vtau(ji,jj) + ( 1. - zfrldv ) * zvtau
+               ! ... module of ice-ocean stress at T-points (from I-point values)
+               zztmp = 0.25 * ( ztiomi(ji,jj) + ztiomi(ji+1,jj) + ztiomi(ji,jj+1) + ztiomi(ji+1,jj+1) )
+               ! ... update module of surface ocean stress (ice-cover wheighted)
+               taum(ji,jj) = frld(ji,jj) * taum(ji,jj) + ( 1. - frld(ji,jj) ) * zztmp
                !
             END DO
          END DO
 
-         ! boundary condition on the stress (utau,vtau)
+         ! boundary condition on the stress (utau,vtau,taum)
          CALL lbc_lnk( utau, 'U', -1. )
          CALL lbc_lnk( vtau, 'V', -1. )
+         CALL lbc_lnk( taum, 'T',  1. )
 
       ENDIF