Changes between Version 34 and Version 35 of user/ctlod/NewSurfaceModule

Timestamp:

2008-02-14T16:19:45+01:00 (16 years ago)

Author:

ctlod

Comment:

--

user/ctlod/NewSurfaceModule

@@ -12 +12 @@
 
 '''EXPERIENCES DESCRIPTION:'''[[BR]]
-
-'''* T2X''':[[BR]]
-In doing the following modifications in both experiments, we are waiting for almost exactly the same results.[[BR]]
-# The wind stress ([uv]tau/gtau[xy]) and wind module (zwind_speed_t(:,:)/dUnormt(:,:) at the ocean surface (even [[BR]]
-above sea-ice) is forced to be the same in both experiments. (sbcblk_core.F90/flx_core.h90)[[BR]]
-# The friction velocity ust2s(:,:) between sea-ice and ocean is also forced to be the same (limdyn.F90)[[BR]]
-
-'''* FULL''': done for REF_CORE and SBC_NEW: [[BR]]
-We fully relax constrains imposed in T2X experiences, i.e. on stress, frictional velocity and so on ...
-
-'''* FULL_SDP''':[[BR]]
-similar to FULL but with a strong restoring SSS damping of ~ 787 mm/day correponding to 12 days of time restoring [[BR]]
-It is the same time damping as the SST which corresponds to -40 W/m2/K used for SST over a 10 meters ocean water column.[[BR]]
-Modifications have been necessary in sbcssr.F90 module such as using the array sdta(:,:) computed in dtasal.F90 to ensure same [[BR]]
-SSS damping field in both experiments SBC_NEW and REF_CORE.[[BR]]
-
-'''* FULL_SDP2''':[[BR]]
-Still SSS damping term without the ice-dynamics, i.e. ln_limdyn=.FALSE. in the sea-ice namelist.
-
-'''* FULL_SDP3''':[[BR]]
-Only for the SBC_NEW code, this time this run is identical to FULL_SDP but in sbcblk_core.F90, the sea-ice
-velocity ([uv]i_ice()) has been removed in the computation of the wind stress and wind module over sea-ice.[[BR]]
-Since the computation is clearly done in a better way now, it could implies important modifications which should be identify. [[BR]]
-'''Small impact O(±2/3 cm) in the Arctic area but shoul be underlined.'''
-
-
+'''* FULL_SDP6''':[[BR]]
+Correction of an other bug in limdyn.F90:[[BR]]
+The friction velocity ust2s() is computed at T-point. To do so, we must compute a mean
+of speed velocities using velocities around T-point at i,j location and not i-1,j-1 location as it is done in the reference.[[BR]]
+[[Color(white,red,'''NOTICE''')]]: this run includes the modification done for SDP5. See above.[[BR]]
+{{{
+        ! computation of friction velocity
+         DO jj = 2, jpjm1
+            DO ji = 2, jpim1
+
+!CT start SDP6
+               zu_ice   = u_ice(ji+1,jj+1) - u_oce(ji+1,jj+1)
+               zv_ice   = v_ice(ji+1,jj+1) - v_oce(ji+1,jj+1)
+!CT SDP6              zu_ice   = u_ice(ji-1,jj-1) - u_oce(ji-1,jj-1)
+!CT SDP6              zv_ice   = v_ice(ji-1,jj-1) - v_oce(ji-1,jj-1)
+               zt11  = rhoco * ( zu_ice * zu_ice + zv_ice * zv_ice )
+
+               zu_ice   = u_ice(ji,jj+1) - u_oce(ji,jj+1)
+               zv_ice   = v_ice(ji,jj+1) - v_oce(ji,jj+1)
+!CT SDP6              zu_ice   = u_ice(ji-1,jj) - u_oce(ji-1,jj)
+!CT SDP6              zv_ice   = v_ice(ji-1,jj) - v_oce(ji-1,jj)
+               zt12  = rhoco * ( zu_ice * zu_ice + zv_ice * zv_ice )
+
+               zu_ice   = u_ice(ji+1,jj) - u_oce(ji+1,jj)
+               zv_ice   = v_ice(ji+1,jj) - v_oce(ji+1,jj)
+!CT SDP6              zu_ice   = u_ice(ji,jj-1) - u_oce(ji,jj-1)
+!CT SDP6              zv_ice   = v_ice(ji,jj-1) - v_oce(ji,jj-1)
+               zt21  = rhoco * ( zu_ice * zu_ice + zv_ice * zv_ice )
+
+               zu_ice   = u_ice(ji,jj) - u_oce(ji,jj)
+               zv_ice   = v_ice(ji,jj) - v_oce(ji,jj)
+               zt22  = rhoco * ( zu_ice * zu_ice + zv_ice * zv_ice )
+!CT end SDP6
+
+               ztair2 = gtaux(ji,jj) * gtaux(ji,jj) + gtauy(ji,jj) * gtauy(ji,jj)
+
+!CT start SDP5
+               zustm =  0.25 * ( zt11 + zt12 + zt21 + zt22 )
+!CT start               zustm =  ( 1 - frld(ji,jj) ) * 0.25 * ( zt11 + zt12 + zt21 + zt22 )        &
+!CT start                  &  +        frld(ji,jj)   * SQRT( ztair2 )
+!CT start SDP5
+
+               ust2s(ji,jj) = ( zustm / rau0 ) * ( rone + sdvt(ji,jj) ) * tms(ji,jj)
+            END DO
+         END DO
+}}}
+
+----
+== '''FULL_SDP5 Experiences''': ==
+'''* FULL_SDP5''':[[BR]]
+Only for the REF_CORE code, we remove the fraction leads frld in the computation of 
+the friction velocity between sea-ice & ocean in limdyn.F90. It is not necessary to mix since 
+the friction velocity is only used to compute flux at Ice/ocean interface? So:[[BR]]
+{{{
+               ztair2 = gtaux(ji,jj) * gtaux(ji,jj) + gtauy(ji,jj) * gtauy(ji,jj)
+!CT start SDP5
+               zustm =  0.25 * ( zt11 + zt12 + zt21 + zt22 )
+!CT start               zustm =  ( 1 - frld(ji,jj) ) * 0.25 * ( zt11 + zt12 + zt21 + zt22 )        &
+!CT start                  &  +        frld(ji,jj)   * SQRT( ztair2 )
+!CT start SDP5
+               ust2s(ji,jj) = ( zustm / rau0 ) * ( rone + sdvt(ji,jj) ) * tms(ji,jj)
+}}}
+Above, gtau[xy] are wind stress components over the ocean at I-point.[[BR]]
+In doing that we are closer to what is done in limdyn.F90 of the NEW SBC.[[BR]]
+It results in a sea-ice surface temperature colder O(-1/2°C) than in SDP Ref mainly in the limit sea-ice extent.[[BR]]
+And a systematic decrease of the icet thickness by O(-2/4 cm) over all the Arctic reaching in some area O(-10 cm) Hudson bay or near 
+the Bering strait.[[BR]]
+So its impact is significant and must be take into account in the rewritten of this part in the NEW _SBC.
+[[Color(white,red,'''BUT''')]] it remains a bug in the comptutation of zt11, zt12 ..etc see SDP6 experience below[[BR]]
+
+'''Ice temperature diffs. between REF_CORE SDP5 - REF_CORE SDP / Year 10 / Winter March''':[[BR]]
+[[Image(SDP5_Ref_SDP_Ref_icetemp_March.png, 90%)]]
+----
+'''Ice thickness diffs. between REF_CORE SDP5 - REF_CORE SDP / Year 10 / Winter March''':[[BR]]
+[[Image(SDP5_Ref_SDP_Ref_icethic_March.png, 90%)]]
+----
+
+----
+== '''FULL_SDP4 Experiences''': ==
 '''* FULL_SDP4''':[[BR]]
 Only for the REF_CORE code, correction of a bug in the computation of the final wind stress in ocesbc.F90 module.[[BR]]
@@ -59 +112 @@
 The impact of this modification is weak as seen on plots below change in ice thickness are O(±1/2 cm)) 
 in Arctic during winter (March). Same behaviour in summer (Sept).
-
-'''* FULL_SDP5''':[[BR]]
-Only for the REF_CORE code, we remove the fraction leads frld in the computation of 
-the friction velocity between sea-ice & ocean in limdyn.F90. It is not necessary to mix since 
-the friction velocity is only used to compute flux at Ice/ocean interface? So:[[BR]]
-{{{
-               ztair2 = gtaux(ji,jj) * gtaux(ji,jj) + gtauy(ji,jj) * gtauy(ji,jj)
-!CT start SDP5
-               zustm =  0.25 * ( zt11 + zt12 + zt21 + zt22 )
-!CT start               zustm =  ( 1 - frld(ji,jj) ) * 0.25 * ( zt11 + zt12 + zt21 + zt22 )        &
-!CT start                  &  +        frld(ji,jj)   * SQRT( ztair2 )
-!CT start SDP5
-               ust2s(ji,jj) = ( zustm / rau0 ) * ( rone + sdvt(ji,jj) ) * tms(ji,jj)
-}}}
-Above, gtau[xy] are wind stress components over the ocean at I-point.[[BR]]
-In doing that we are closer to what is done in limdyn.F90 of the NEW SBC.[[BR]]
-It results in a sea-ice surface temperature colder O(-1/2°C) than in SDP Ref mainly in the limit sea-ice extent.[[BR]]
-And a systematic decrease of the icet thickness by O(-2/4 cm) over all the Arctic reaching in some area O(-10 cm) Hudson bay or near 
-the Bering strait.[[BR]]
-So its impact is significant and must be take into account in the rewritten of this part in the NEW _SBC.
-[[Color(white,red,'''BUT''')]] it remains a bug in the comptutation of zt11, zt12 ..etc see SDP6 experience below[[BR]]
-
-
-'''* FULL_SDP6''':[[BR]]
-Correction of an other bug in limdyn.F90:[[BR]]
-The friction velocity ust2s() is computed at T-point. To do so, we must compute a mean
-of speed velocities using velocities around T-point at i,j location and not i-1,j-1 location as it is done in the reference.[[BR]]
-[[Color(white,red,'''NOTICE''')]]: this run includes the modification done for SDP5. See above.[[BR]]
-{{{
-        ! computation of friction velocity
-         DO jj = 2, jpjm1
-            DO ji = 2, jpim1
-
-!CT start SDP6
-               zu_ice   = u_ice(ji+1,jj+1) - u_oce(ji+1,jj+1)
-               zv_ice   = v_ice(ji+1,jj+1) - v_oce(ji+1,jj+1)
-!CT SDP6              zu_ice   = u_ice(ji-1,jj-1) - u_oce(ji-1,jj-1)
-!CT SDP6              zv_ice   = v_ice(ji-1,jj-1) - v_oce(ji-1,jj-1)
-               zt11  = rhoco * ( zu_ice * zu_ice + zv_ice * zv_ice )
-
-               zu_ice   = u_ice(ji,jj+1) - u_oce(ji,jj+1)
-               zv_ice   = v_ice(ji,jj+1) - v_oce(ji,jj+1)
-!CT SDP6              zu_ice   = u_ice(ji-1,jj) - u_oce(ji-1,jj)
-!CT SDP6              zv_ice   = v_ice(ji-1,jj) - v_oce(ji-1,jj)
-               zt12  = rhoco * ( zu_ice * zu_ice + zv_ice * zv_ice )
-
-               zu_ice   = u_ice(ji+1,jj) - u_oce(ji+1,jj)
-               zv_ice   = v_ice(ji+1,jj) - v_oce(ji+1,jj)
-!CT SDP6              zu_ice   = u_ice(ji,jj-1) - u_oce(ji,jj-1)
-!CT SDP6              zv_ice   = v_ice(ji,jj-1) - v_oce(ji,jj-1)
-               zt21  = rhoco * ( zu_ice * zu_ice + zv_ice * zv_ice )
-
-               zu_ice   = u_ice(ji,jj) - u_oce(ji,jj)
-               zv_ice   = v_ice(ji,jj) - v_oce(ji,jj)
-               zt22  = rhoco * ( zu_ice * zu_ice + zv_ice * zv_ice )
-!CT end SDP6
-
-               ztair2 = gtaux(ji,jj) * gtaux(ji,jj) + gtauy(ji,jj) * gtauy(ji,jj)
-
-!CT start SDP5
-               zustm =  0.25 * ( zt11 + zt12 + zt21 + zt22 )
-!CT start               zustm =  ( 1 - frld(ji,jj) ) * 0.25 * ( zt11 + zt12 + zt21 + zt22 )        &
-!CT start                  &  +        frld(ji,jj)   * SQRT( ztair2 )
-!CT start SDP5
-
-               ust2s(ji,jj) = ( zustm / rau0 ) * ( rone + sdvt(ji,jj) ) * tms(ji,jj)
-            END DO
-         END DO
-}}}
-
-----
-== '''FULL_SDP5 Experiences''': ==
-'''Ice temperature diffs. between REF_CORE SDP5 - REF_CORE SDP / Year 10 / Winter March''':[[BR]]
-[[Image(SDP5_Ref_SDP_Ref_icetemp_March.png, 90%)]]
-----
-'''Ice thickness diffs. between REF_CORE SDP5 - REF_CORE SDP / Year 10 / Winter March''':[[BR]]
-[[Image(SDP5_Ref_SDP_Ref_icethic_March.png, 90%)]]
-----
-
-----
-== '''FULL_SDP4 Experiences''': ==
 '''Ice thickness diffs. between REF_CORE SDP4 - REF_CORE SDP / Year 10 / Winter March''':[[BR]]
 [[Image(SDP4_Ref_SDP_Ref_icethic_March.png, 90%)]]
@@ -149 +121 @@
 ----
 == '''FULL_SDP3 Experiences''': ==
+'''* FULL_SDP3''':[[BR]]
+Only for the SBC_NEW code, this time this run is identical to FULL_SDP but in sbcblk_core.F90, the sea-ice
+velocity ([uv]i_ice()) has been removed in the computation of the wind stress and wind module over sea-ice.[[BR]]
+Since the computation is clearly done in a better way now, it could implies important modifications which should be identify. [[BR]]
+'''Small impact O(±2/3 cm) in the Arctic area but shoul be underlined.'''
+
 '''Ice thickness in NEW_SBC experience / Year 10 / Winter March''':[[BR]]
 [[Image(SDP3_New_icethic_March.png, 90%)]]
@@ -174 +152 @@
 
 == '''FULL_SDP2 Experiences''': ==
+'''* FULL_SDP2''':[[BR]]
+Still SSS damping term without the ice-dynamics, i.e. ln_limdyn=.FALSE. in the sea-ice namelist.
+
 '''Ice thickness in NEW_SBC experience / Year 10''': NO SEA-ICE DYNAMICS = ln_limdyn=FALSE[[BR]]
 [[Image(FULL_SDP2_SBC_NEW_icethic_March_Arctic.png, 90%)]]
@@ -182 +163 @@
 
 == '''FULL_SDP Experiences''': ==
+'''* FULL_SDP''':[[BR]]
+similar to FULL but with a strong restoring SSS damping of ~ 787 mm/day correponding to 12 days of time restoring [[BR]]
+It is the same time damping as the SST which corresponds to -40 W/m2/K used for SST over a 10 meters ocean water column.[[BR]]
+Modifications have been necessary in sbcssr.F90 module such as using the array sdta(:,:) computed in dtasal.F90 to ensure same [[BR]]
+SSS damping field in both experiments SBC_NEW and REF_CORE.[[BR]]
+
 '''Ice thickness in NEW_SBC experience / Year 10'''[[BR]]
 The sea-ice extension in winter seems to be larger in this run than in the REF_CORE one.[[BR]]
@@ -190 +177 @@
 ----
 
+== '''FULL Experiences''': ==
+
+'''* FULL''': done for REF_CORE and SBC_NEW: [[BR]]
+We fully relax constrains imposed in T2X experiences, i.e. on stress, frictional velocity and so on ...
+
 
    
 == '''T24 Experience''': ==
+'''* T2X''':[[BR]]
+In doing the following modifications in both experiments, we are waiting for almost exactly the same results.[[BR]]
+# The wind stress ([uv]tau/gtau[xy]) and wind module (zwind_speed_t(:,:)/dUnormt(:,:) at the ocean surface (even [[BR]]
+above sea-ice) is forced to be the same in both experiments. (sbcblk_core.F90/flx_core.h90)[[BR]]
+# The friction velocity ust2s(:,:) between sea-ice and ocean is also forced to be the same (limdyn.F90)[[BR]]
 '''Ice thickness in NEW_SBC experience / Date 11/2000'''[[BR]]
 Image(T24_Icethic_New_Polaire_112000.png, 90%)