Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

    r3691 r3692  
    711711            ENDIF
    712712          ENDDO
    713           zx_tmp_fi2d(:)=missing_val
     713          wind100m(:)=missing_val
    714714          DO k=1, kmax_100m-1                                      !--we could stop much lower
    715715            DO i=1,klon
     
    726726         !--polynomial fit for 14,Vestas,1074,V136/3450 kW windmill - Olivier
    727727         DO i=1,klon
    728            IF (pctsrf(i,is_ter).GT.0.05) THEN
     728           IF (pctsrf(i,is_ter).GT.0.05 .AND. wind100m(i).NE.missing_val) THEN
    729729             x=wind100m(i)
    730730             IF (x.LE.3.0 .OR. x.GE.22.5) THEN
     
    738738             ENDIF
    739739           ELSE
    740              zx_tmp_fi2d(i)=missing_val_nf90
     740             zx_tmp_fi2d(i)=missing_val
    741741           ENDIF
    742742         ENDDO
     
    747747         !--polynomial fit for 14,Vestas,867,V164/8000 kW - Olivier
    748748         DO i=1,klon
    749            IF (pctsrf(i,is_oce).GT.0.05) THEN
     749           IF (pctsrf(i,is_oce).GT.0.05 .AND. wind100m(i).NE.missing_val) THEN
    750750             x=wind100m(i)
    751751             IF (x.LE.3.0 .OR. x.GE.25.5) THEN
     
    759759             ENDIF
    760760           ELSE
    761              zx_tmp_fi2d(i)=missing_val_nf90
     761             zx_tmp_fi2d(i)=missing_val
    762762           ENDIF
    763763         ENDDO
Note: See TracChangeset for help on using the changeset viewer.