Changes between Version 40 and Version 41 of user/ctlod/NewSurfaceModule

Timestamp:

2008-02-22T17:21:50+01:00 (16 years ago)

Author:

ctlod

Comment:

--

user/ctlod/NewSurfaceModule

@@ -13 +13 @@
 
 '''EXPERIENCES DESCRIPTION:'''[[BR]]
+== '''FULL_SDP7''': ==
+'''* FULL_SDP7''':[[BR]]
+Only for REF_CORE: following FULL_SDP6 exp., modifications in REF_CORE in the wind stress management over sea-ice and ocean.[[BR]]
+- '''in flx_core.h90''', splitting of the ocean/sea-ice wind stress:
+{{{
+          tauxt(:,:) = rhoa*dUnormt(:,:)*Cd(:,:)*(flxnow(:,:,2) - zut(:,:))
+          taux_it(:,:) = rhoa*dUnormt(:,:)*Cice*flxnow(:,:,2)
+
+!CT          tauxt(:,:) = rhoa*dUnormt(:,:)*( (1. - freeze(:,:))*Cd(:,:)*(flxnow(:,:,2) - zut(:,:))    &
+!CT               &                          + freeze(:,:)*Cice*flxnow(:,:,2) ) !lb correct pour glace
+         ! Tau_y at T-point
+         tauyt(:,:) = rhoa*dUnormt(:,:)*Cd(:,:)*(flxnow(:,:,3) - zvt(:,:))
+         tauy_it(:,:) = rhoa*dUnormt(:,:)*Cice*flxnow(:,:,3)
+}}}
+- the ocean wind stress is masked using tmask(:,:)
+{{{
+         !Tau_x at U-point
+         DO jj = 1, jpjm1
+            DO ji = 1, jpi
+               taux(ji,jj) = 0.5*(tauxt(ji,jj) + tauxt(ji+1,jj))*( 2. - tmask(ji,jj,1) )
+               tauy(ji,jj) = 0.5*(tauyt(ji,jj) + tauyt(ji,jj+1))*( 2. - tmask(ji,jj,1) )
+            END DO
+         END DO
+}}}
+- '''in icestp.F90''': compute wind stress gtau[xy] over sea-ice at I-point location using tau[xy]_it(:,:) fileds.[[BR]]
+- '''in limdyn.F90:''' compute the ice/ocean stress the same way it is done in NEW_SBC, i.e. only based on the wind speed [uv]_ice - [uv]_oce differences [[BR]]
+- '''in ocesbc.F90:''' the final wind stress over the ocean is the result of both ice/ocean & "pure" ocean stress mixing weighted with leads fraction.[[BR]]
+
+
 == '''FULL_SDP6''': ==
 '''* FULL_SDP6''':[[BR]]