Changeset 3691 for LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

May 30, 2020, 5:37:19 PM (4 years ago)

Author:

oboucher

Message:

Addition of two diagnostics for load factors for onshore and offshore wind mills

File:

• LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -40 +40 @@
          o_q2m, o_ustar, o_u10m, o_v10m, &
          o_wind10m, o_wind10max, o_wind100m, o_gusts, o_sicf, &
+         o_loadfactor_wind_onshore, o_loadfactor_wind_offshore, &
          o_psol, o_mass, o_qsurf, o_qsol, &
          o_precip, o_rain_fall, o_rain_con, o_ndayrain, o_plul, o_pluc, o_plun, &
@@ -411 +412 @@
     INTEGER :: itau_w
     INTEGER :: i, iinit, iinitend=1, iff, iq, iiq, nsrf, k, ll, naero
-    REAL, DIMENSION (klon) :: zx_tmp_fi2d, zpt_conv2d
+    REAL, DIMENSION (klon) :: zx_tmp_fi2d, zpt_conv2d, wind100m
     REAL, DIMENSION (klon,klev) :: zx_tmp_fi3d, zpt_conv
     REAL, DIMENSION (klon,klev+1) :: zx_tmp_fi3d1
@@ -425 +426 @@
     INTEGER, SAVE :: kmax_100m
 !$OMP THREADPRIVATE(kmax_100m)
+    REAL :: x
 #ifndef CPP_XIOS
     REAL :: missing_val
@@ -713 +715 @@
             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 )
+                wind100m(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)
+       CALL histwrite_phy(o_wind100m, wind100m)
+
+       IF (vars_defined) THEN
+         !--polynomial fit for 14,Vestas,1074,V136/3450 kW windmill - Olivier
+         DO i=1,klon
+           IF (pctsrf(i,is_ter).GT.0.05) THEN
+             x=wind100m(i)
+             IF (x.LE.3.0 .OR. x.GE.22.5) THEN
+               zx_tmp_fi2d(i)=0.0
+             ELSE IF (x.GE.10.0) THEN
+               zx_tmp_fi2d(i)=1.0
+             ELSE
+               zx_tmp_fi2d(i)= 1.059e-09*x**10. - 1.351e-07*x**9. + 7.478e-06*x**8. - 0.0002352*x**7. + 0.004627*x**6. & 
+                             - 0.05898*x**5. + 0.4893*x**4. - 2.59*x**3. + 8.339*x**2. - 14.69*x + 10.73
+               zx_tmp_fi2d(i)=MIN(MAX(zx_tmp_fi2d(i),0.0),1.0)
+             ENDIF
+           ELSE
+             zx_tmp_fi2d(i)=missing_val_nf90
+           ENDIF
+         ENDDO
+       ENDIF
+       CALL histwrite_phy(o_loadfactor_wind_onshore, zx_tmp_fi2d)
+
+       IF (vars_defined) THEN
+         !--polynomial fit for 14,Vestas,867,V164/8000 kW - Olivier
+         DO i=1,klon
+           IF (pctsrf(i,is_oce).GT.0.05) THEN
+             x=wind100m(i)
+             IF (x.LE.3.0 .OR. x.GE.25.5) THEN
+               zx_tmp_fi2d(i)=0.0
+             ELSE IF (x.GE.12.5) THEN
+               zx_tmp_fi2d(i)=1.0
+             ELSE
+               zx_tmp_fi2d(i)= 3.352e-10*x**10. - 4.959e-08*x**9. + 3.195e-06*x**8. - 0.0001175*x**7. + 0.002716*x**6. &
+                             - 0.04099*x**5. + 0.4065*x**4. - 2.601*x**3. + 10.25*x**2. - 22.39*x + 20.59
+               zx_tmp_fi2d(i)=MIN(MAX(zx_tmp_fi2d(i),0.0),1.0)
+             ENDIF
+           ELSE
+             zx_tmp_fi2d(i)=missing_val_nf90
+           ENDIF
+         ENDDO
+       ENDIF
+       CALL histwrite_phy(o_loadfactor_wind_offshore, zx_tmp_fi2d)
 
        IF (vars_defined) THEN