Changes between Version 53 and Version 54 of user/ctlod/NewSurfaceModule

Timestamp:

2008-02-29T11:45:06+01:00 (16 years ago)

Author:

ctlod

Comment:

--

user/ctlod/NewSurfaceModule

@@ -13 +13 @@
 
 '''EXPERIENCES DESCRIPTION:'''[[BR]]
+== '''FULL_SDP10''': ==
+'''* FULL_SDP10''':[[BR]]
+Both for NEW_SBC & REF_CORE: follow SDP9 BUT with ln_limdyn=.FALSE.[[BR]]
+Modifications in the REF_CORE limdyn.F90 and ocesbc.F90:[[BR]]
+limdyn.F90:[[BR]]
+{{{
+!CT start SDP10
+!CT gtau[xy] fields are computed at I-points, whereas ust2s() must be at T-point
+!CT so we take here directly the wind stress over the ocean tau[xy] calculated in flx_core.h90
+!CT the computation is exactly the same than in NEW_SBC
+                ztair2       = taux(ji,jj) * taux(ji,jj) + taux(ji-1,jj) * taux(ji-1,jj)   &
+                     &       + tauy(ji,jj) * tauy(ji,jj) + tauy(ji,jj-1) * tauy(ji,jj-1)
+                zustm        = SQRT( 0.5 * ztair2  )
+
+                ust2s(ji,jj) = ( zustm / rau0 ) * tms(ji,jj)
+
+!CT SDP10               ztair2       = gtaux(ji,jj) * gtaux(ji,jj) + gtauy(ji,jj) * gtauy(ji,jj)
+!CT SDP10               zustm        = SQRT( ztair2  )
+!CT
+!CT SDP10               ust2s(ji,jj) = ( zustm / rau0 ) * ( rone + sdvt(ji,jj) ) * tms(ji,jj)
+!CT end SDP10
+}}}
+ocesbc.F90: add the logical ln_limdyn to avoid modifications of the ocean stress tau[xy].[[BR]]
+
+
+
 == '''FULL_SDP9''': ==
 '''* FULL_SDP9''':[[BR]]