Changeset 2255 for LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90

Timestamp:

Apr 3, 2015, 3:51:31 PM (9 years ago)

Author:

jyg

Message:

Changes to pbl_surface and other routines concerning split/no-split.
+ pbl_surface_mod.F90: call cdrag for (w) region.
+ phyredem.F90: write wake_delta_pbl_TKE.
+ phys_output_write_mod.F90: control output of wake_delta_pbl_TKE by
IF(iflag_pbl_split>=1).
+ lmdz1d.F90: initialize wake_delta_pbl_TKE=0.
+ phys_output_ctrlout_mod.F90: suppression of accents in some variable
attributes.
+ cva_driver.F90: suppression of a print introduced in version 2253.

File:

• LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -1351 +1351 @@
         IF (prt_level >=10) print *,'clcdrag -> ycdragh_x ', ycdragh_x
 !
-        CALL clcdrag( knon, nsrf, ypaprs, ypplay, &
-            yu_w(:,1), yv_w(:,1), yt_w(:,1), yq_w(:,1), &
-            yts_w, yqsurf, yz0m, &
-            ycdragm_w, ycdragh_w )
+! Faire disparaitre les lignes commentees fin 2015 (le temps des tests)
+!        CALL clcdrag( knon, nsrf, ypaprs, ypplay, &
+!            yu_w(:,1), yv_w(:,1), yt_w(:,1), yq_w(:,1), &
+!            yts_w, yqsurf, yz0m, &
+!            ycdragm_w, ycdragh_w )
+! Fuxing WANG, 04/03/2015, replace the clcdrag by the merged version: cdrag
+        DO i = 1, knon
+           zgeo1_w(i) = RD * yt_w(i,1) / (0.5*(ypaprs(i,1)+ypplay(i,1))) &
+                * (ypaprs(i,1)-ypplay(i,1))
+           speed_w(i) = SQRT(yu_w(i,1)**2+yv_w(i,1)**2)
+        END DO
+        CALL cdrag(knon, nsrf, &
+            speed_w, yt_w(:,1), yq_w(:,1), zgeo1_w, ypaprs(:,1),&
+            yts_w, yqsurf, yz0m, yz0h, &
+            ycdragm_w, ycdragh_w, zri1_w, pref_w )
 
 ! --- special Dice. JYG+MPL 25112013