Changeset 2159 for LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

Timestamp:

Nov 27, 2014, 4:48:31 PM (10 years ago)

Author:

jyg

Message:

1/ Splitting of the boundary layer : the climbing down and up of Pbl_surface is
split between the off-wake and wake regions ; the thermal scheme is applied
only to the off-wake region.
2/ Elimination of wake_scal and calwake_scal.

File:

• LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -477 +477 @@
   TYPE(ctrl_out), SAVE :: o_alp_wk = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
     'alp_wk', 'ALP WK', 'm2/s2', (/ ('', i=1, 9) /))
+!!!
+!nrlmd+jyg<
+  type(ctrl_out),save :: o_dtvdf_x        = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
+    'dtvdf_x', ' dtvdf off_wake','K/s', (/ ('', i=1, 9) /))
+  type(ctrl_out),save :: o_dtvdf_w        = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
+    'dtvdf_w', ' dtvdf within_wake','K/s', (/ ('', i=1, 9) /))
+  type(ctrl_out),save :: o_dqvdf_x        = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
+    'dqvdf_x', ' dqvdf off_wake','kg/kg/s', (/ ('', i=1, 9) /))
+  type(ctrl_out),save :: o_dqvdf_w        = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
+    'dqvdf_w', ' dqvdf within_wake','kg/kg/s', (/ ('', i=1, 9) /))
+!!
+  type(ctrl_out),save :: o_sens_x        = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+'sens_x', 'ALP WK', 'm2/s2', (/ ('', i=1, 9) /))
+  type(ctrl_out),save :: o_sens_w        = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+'sens_w', 'ALP WK', 'm2/s2', (/ ('', i=1, 9) /))                                                                                   
+  type(ctrl_out),save :: o_flat_x        = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+'flat_x', 'ALP WK', 'm2/s2', (/ ('', i=1, 9) /))                                                                                   
+  type(ctrl_out),save :: o_flat_w        = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+'flat_w', 'ALP WK', 'm2/s2', (/ ('', i=1, 9) /))
+!!
+  type(ctrl_out),save :: o_delta_tsurf    = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+'delta_tsurf', 'Temperature difference (w-x)', 'K', (/ ('', i=1, 9) /))                                                                               
+  type(ctrl_out),save :: o_cdragh_x       = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+'cdragh_x', 'cdragh off-wake', '', (/ ('', i=1, 9) /))
+  type(ctrl_out),save :: o_cdragh_w       = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+'cdragh_w', 'cdragh within-wake', '', (/ ('', i=1, 9) /))                                                                                  
+  type(ctrl_out),save :: o_cdragm_x       = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+'cdragm_x', 'cdragm off-wake', '', (/ ('', i=1, 9) /))
+  type(ctrl_out),save :: o_cdragm_w       = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+'cdragm_w', 'cdrgam within-wake', '', (/ ('', i=1, 9) /))                                                                                  
+  type(ctrl_out),save :: o_kh             = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+'kh', 'Kh', 'kg/s/m2', (/ ('', i=1, 9) /))                                                                                        
+  type(ctrl_out),save :: o_kh_x           = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+'kh_x', 'Kh off-wake', 'kg/s/m2', (/ ('', i=1, 9) /))                                                                                      
+  type(ctrl_out),save :: o_kh_w           = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+'kh_w', 'Kh within-wake', 'kg/s/m2', (/ ('', i=1, 9) /))
+!>nrlmd+jyg
+!!!
   TYPE(ctrl_out), SAVE :: o_ale = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
     'ale', 'ALE', 'm2/s2', (/ ('', i=1, 9) /))
@@ -693 +731 @@
       (/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
          "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /)) /)
+
+  TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_dltpbltke_srf      = (/             &
+      ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'dltpbltke_ter',       &
+      "TKE difference (w - x) "//clnsurf(1),"-", (/ ('', i=1, 9) /)), &
+      ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'dltpbltke_lic',       &
+      "TKE difference (w - x) "//clnsurf(2),"-", (/ ('', i=1, 9) /)), &
+      ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'dltpbltke_oce',       &
+      "TKE difference (w - x) "//clnsurf(3),"-", (/ ('', i=1, 9) /)), &
+      ctrl_out((/ 10, 4, 10, 10, 10, 10, 11, 11, 11 /),'dltpbltke_sic',       &
+      "TKE difference (w - x) "//clnsurf(4),"-", (/ ('', i=1, 9) /)) /)
 
   TYPE(ctrl_out), SAVE :: o_kz = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &