Changeset 3616 for LMDZ6/branches/IPSLCM6.0.15

Timestamp:

Dec 16, 2019, 2:10:03 PM (5 years ago)

Author:

oboucher

Message:

Addition of diagnostics for the 100 meter wind

Location:

LMDZ6/branches/IPSLCM6.0.15

Files:

• DefLists/field_def_lmdz.xml (modified) (1 diff)
• libf/phylmd/phys_output_ctrlout_mod.F90 (modified) (1 diff)
• libf/phylmd/phys_output_write_mod.F90 (modified) (7 diffs)

LMDZ6/branches/IPSLCM6.0.15/DefLists/field_def_lmdz.xml

@@ -49 +49 @@
         <field id="wind10max"  long_name="10m wind speed max"   unit="m/s" />
         <field id="sfcWindmax" field_ref="wind10m" operation="maximum" />
+        <field id="wind100m"   long_name="100-m wind speed"      unit="m/s" />
         <field id="sicf"    long_name="Sea-ice fraction"        unit="-" />
         <field id="q2m"     long_name="Specific humidity 2m"    unit="kg/kg" />

LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -279 +279 @@
   TYPE(ctrl_out), SAVE :: o_wind10m = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11, 11/), &
     'wind10m', '10-m wind speed', 'm/s', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_wind100m = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'wind100m', '100-m wind speed', 'm/s', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_wind10max = ctrl_out((/ 10, 1, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'wind10max', '10m wind speed max', 'm/s', &

LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/phys_output_write_mod.F90

@@ -39 +39 @@
          o_t2m_min_mon, o_t2m_max_mon, &
          o_q2m, o_ustar, o_u10m, o_v10m, &
-         o_wind10m, o_wind10max, o_gusts, o_sicf, &
+         o_wind10m, o_wind10max, o_wind100m, o_gusts, o_sicf, &
          o_psol, o_mass, o_qsurf, o_qsol, &
          o_precip, o_rain_fall, o_rain_con, o_ndayrain, o_plul, o_pluc, o_plun, &
@@ -367 +367 @@
     USE print_control_mod, ONLY: prt_level,lunout
 
-
 #ifdef CPP_XIOS
     ! ug Pour les sorties XIOS
@@ -379 +378 @@
 #endif
     USE tracinca_mod, ONLY: config_inca
+
+    USE vertical_layers_mod, ONLY: presnivs
 
     IMPLICIT NONE
@@ -421 +422 @@
     REAL, DIMENSION(klev,2) :: Ahyb_mid_bounds, Bhyb_mid_bounds
     INTEGER :: ilev
+    INTEGER, SAVE :: kmax_100m
+!$OMP THREADPRIVATE(kmax_100m)
 #ifndef CPP_XIOS
     REAL :: missing_val
@@ -436 +439 @@
     LOGICAL, PARAMETER :: debug_strataer=.FALSE.
 #endif
+    REAL,DIMENSION(klon,klev) :: z, dz
+    REAL,DIMENSION(klon)      :: zrho, zt
 
     ! On calcul le nouveau tau:
@@ -446 +451 @@
     ELSE
        iinitend = 1
+    ENDIF
+
+    IF (.NOT.vars_defined) THEN 
+      kmax_100m=1
+      DO k=1, klev-1 
+        IF (presnivs(k).GT.0.97*101325.) kmax_100m = k !--finding out max level for 100 m with a good margin
+      ENDDO
     ENDIF
 
@@ -678 +690 @@
 
        CALL histwrite_phy(o_gusts, gustiness)
+
+       IF (vars_defined) THEN
+          DO k = 1, kmax_100m                                      !--we could stop much lower
+            zrho(:) = pplay(:,k)/t_seri(:,k)/RD                    ! air density in kg/m3
+            dz(:,k) = (paprs(:,k)-paprs(:,k+1))/zrho(:)/RG         ! layer thickness in m
+            IF (k==1) THEN
+              z(:,1) = (paprs(:,1)-pplay(:,1))/zrho(:)/RG          ! altitude middle of first layer in m
+              zt(:)  = dz(:,1)                                     ! altitude top of first layer in m
+            ELSE
+              z(:,k) = zt(:) + (paprs(:,k)-pplay(:,k))/zrho(:)/RG  ! altitude middle of layer k in m
+              zt(:)  = zt(:) + dz(:,k)                             ! altitude top of layer k in m
+            ENDIF
+          ENDDO
+          zx_tmp_fi2d(:)=missing_val
+          DO k=1, kmax_100m-1                                      !--we could stop much lower 
+            DO i=1,klon 
+              IF (z(i,k).LT.100..AND.z(i,k+1).GE.100.) THEN
+                zx_tmp_fi2d(i)=SQRT( (u_seri(i,k)+(100.-z(i,k))/(z(i,k+1)-z(i,k))*(u_seri(i,k+1)-u_seri(i,k)))**2.0 +  &
+                                     (v_seri(i,k)+(100.-z(i,k))/(z(i,k+1)-z(i,k))*(v_seri(i,k+1)-v_seri(i,k)))**2.0 )
+              ENDIF
+            ENDDO
+          ENDDO
+       ENDIF
+       CALL histwrite_phy(o_wind100m, zx_tmp_fi2d)
 
        IF (vars_defined) THEN