Changeset 3692

Timestamp:

Jun 2, 2020, 4:57:54 PM (4 years ago)

Author:

oboucher

Message:

correcting the missing_val for wind100m diagnostic and its impact on load factors for windmills

File:

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

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -711 +711 @@
             ENDIF
           ENDDO
-          zx_tmp_fi2d(:)=missing_val
+          wind100m(:)=missing_val
           DO k=1, kmax_100m-1                                      !--we could stop much lower
             DO i=1,klon
@@ -726 +726 @@
          !--polynomial fit for 14,Vestas,1074,V136/3450 kW windmill - Olivier
          DO i=1,klon
-           IF (pctsrf(i,is_ter).GT.0.05) THEN
+           IF (pctsrf(i,is_ter).GT.0.05 .AND. wind100m(i).NE.missing_val) THEN
              x=wind100m(i)
              IF (x.LE.3.0 .OR. x.GE.22.5) THEN
@@ -738 +738 @@
              ENDIF
            ELSE
-             zx_tmp_fi2d(i)=missing_val_nf90
+             zx_tmp_fi2d(i)=missing_val
            ENDIF
          ENDDO
@@ -747 +747 @@
          !--polynomial fit for 14,Vestas,867,V164/8000 kW - Olivier
          DO i=1,klon
-           IF (pctsrf(i,is_oce).GT.0.05) THEN
+           IF (pctsrf(i,is_oce).GT.0.05 .AND. wind100m(i).NE.missing_val) THEN
              x=wind100m(i)
              IF (x.LE.3.0 .OR. x.GE.25.5) THEN
@@ -759 +759 @@
              ENDIF
            ELSE
-             zx_tmp_fi2d(i)=missing_val_nf90
+             zx_tmp_fi2d(i)=missing_val
            ENDIF
          ENDDO