*IF DEF,SEAICE ONT3F304.6
C ******************************COPYRIGHT****************************** GTS2F400.3025
C (c) CROWN COPYRIGHT 1995, METEOROLOGICAL OFFICE, All Rights Reserved. GTS2F400.3026
C GTS2F400.3027
C Use, duplication or disclosure of this code is subject to the GTS2F400.3028
C restrictions as set forth in the contract. GTS2F400.3029
C GTS2F400.3030
C Meteorological Office GTS2F400.3031
C London Road GTS2F400.3032
C BRACKNELL GTS2F400.3033
C Berkshire UK GTS2F400.3034
C RG12 2SZ GTS2F400.3035
C GTS2F400.3036
C If no contract has been raised with this copy of the code, the use, GTS2F400.3037
C duplication or disclosure of it is strictly prohibited. Permission GTS2F400.3038
C to do so must first be obtained in writing from the Head of Numerical GTS2F400.3039
C Modelling at the above address. GTS2F400.3040
C ******************************COPYRIGHT****************************** GTS2F400.3041
C GTS2F400.3042
!+ Blanks out heat fluxes at sea-ice points ONT3F304.7
! ONT3F304.8
! Subroutine Interface: ONT3F304.9
SUBROUTINE FLXBLANK(ICY,QFLUX,QFLUXALT,OCEANHEATFLUX 1OJC2F400.145
+,WME_ICE ONT3F304.12
+,SOL,SOLALT ONT3F304.15
+,WME,WMEALT ONT3F304.18
+,AICE ONT3F304.20
+,IMT) ONT3F304.21
ONT3F304.22
IMPLICIT NONE ONT3F304.23
! ONT3F304.24
! THIS ROUTINE BLANKS OUT THE HEAT FORCING FLUXES AT SEA-ICE ONT3F304.25
! POINTS. IT IS CALLED FROM TRACER. ONT3F304.26
! IT CAN BE COMPILED BY ANY FORTRAN COMPILER WHICH TOLERATES ONT3F304.27
! IN-LINE COMMENTS. ONT3F304.28
! AUTHOR: D L ROBERTS ONT3F304.29
! DATE: 17/10/90 ONT3F304.30
! REVIEWED BY: J F THOMPSON ONT3F304.31
! VERSION 1.2 ONT3F304.32
! REVISED BY: J.F.THOMSON (29/10/92) ONT3F304.33
! SOL (PENETRATING SOLAR RADIATION) UNALTERED UNDER ONT3F304.34
! THERMODYNAMIC SEA ICE ALLOWING BAND 1 SOLAR FLUX TO ONT3F304.35
! PENETRATE LEADS. ONT3F304.36
! Version 1.3 ONT3F304.37
! Revised By: N.K. Taylor (3/8/94) ONT3F304.38
! SOL WEIGHTED BY LEADS FRACTION IF COUPLED MODEL ONT3F304.39
! NOT SELECTED. ONT3F304.40
! *IF DEF STRUCTURE TIDIED UP ONT3F304.41
! Comments changed to UMDP3 v6 standard ONT3F304.42
! ONT3F304.43
! WME allowed to penetrate leads if coupled ONT3F304.44
! model forcing is selected. ONT3F304.45
! NB. This may be extended to non-coupled models ONT3F304.46
! at version 3.5 ONT3F304.47
! Modification History: ORH1F305.3764
! Version Date Details ORH1F305.3765
! ------- ------- ------------------------------------------ ORH1F305.3766
! 3.5 16.01.95 Remove *IF dependency. R.Hill ORH1F305.3767
! 4.0 Addition of OCEANHEATFLUX (initialised in OJC2F400.146
! INIT_OCEAN_ICE) to ensure energy conservation OJC2F400.147
! with sea ice dynamics. J.F.Crossley OJC2F400.148
! ONT3F304.48
! 4.0 26.06.95 WME always penetrates leads with dynamic ice model ORH0F400.31
ORH0F400.32
! Code Description: ONT3F304.49
! Language: FORTRAN 77 + common extensions ONT3F304.50
! This code is written to UMDP3 v6 programming standards. ONT3F304.51
! ONT3F304.52
! System component covered: P4 ONT3F304.53
! ONT3F304.54
*CALL CNTLOCN 
ORH1F305.3768
*CALL OARRYSIZ ORH1F305.3769
INTEGER IMT ! IN NUMBER OF COLUMNS. ONT3F304.55
LOGICAL ICY(IMT) ! IN TRUE WHERE SEA-ICE EXISTS. ONT3F304.56
! ONT3F304.57
REAL ONT3F304.58
+ QFLUX(IMT) ! IN NON-PENETRATIVE SURFACE HEAT FLUX. ONT3F304.59
& ! (IF L_OHANEY - NOT INCL. THE 'HANEY' TYPE CMPONENT) ORH1F305.3770
&,OCEANHEATFLUX(IMT_ICE) !IN QFLUX PORTION FOR OCEAN BUDGET OJC2F400.149
&,SOL (IMT_ICE_SOSO) ! IN PENETRATIVE SOLAR HEAT FLUX ORH1F305.3771
&,AICE (IMT_ICE) ! IN SEA ICE FRACTION ORH1F305.3772
&,WME(IMT_ICE_MIX) ! IN WIND MIXING ENERGY. ORH1F305.3773
&,WME_ICE(IMT_idr_MIX) ! IN WIND MXING ENERGY UNDER DYNAMIC SEA ICE ODC1F405.207
&,QFLUXALT(IMT_ICE) ! OUT QFLUX WITH ZEROES INSERTED AT ICE POINTS ORH1F305.3775
&,SOLALT(IMT_ICE_SOSO) ! OUT SOL WITH ZEROES INSERTED AT ICE POINTS ORH1F305.3776
&,WMEALT(IMT_ICE_MIX) ! OUT WME ORH1F305.3777
+ ! WITH ZEROES INSERTED AT ICE POINTS. ONT3F304.80
+ ! (ocean-only mode) ONT3F304.81
+ ! weighted by leads fraction ONT3F304.82
+ ! (coupled ocean-atmos mode) ONT3F304.83
ONT3F304.85
INTEGER I ! LOOP COUNTER FOR COLUMNS. ONT3F304.86
ONT3F304.87
DO 100 I = 1,IMT ONT3F304.88
IF ( ICY(I) ) THEN ONT3F304.89
QFLUXALT(I) = OCEANHEATFLUX(I) OJC2F400.150
IF (L_OSOLAR.OR.L_OSOLARAL) THEN ORH1F305.3778
IF (.NOT.(L_OICECOUP)) THEN ORH1F305.3779
SOLALT(I) = SOL(I)*(1.0-AICE(I)) ONT3F304.93
ELSE ORH1F305.3780
SOLALT(I) = SOL(I) ONT3F304.95
ENDIF ORH1F305.3781
ENDIF ORH1F305.3782
IF (L_OMIXLAY) THEN ORH1F305.3783
IF (L_OICECOUP) THEN ORH1F305.3784
WMEALT(I) = WME(I) ONT3F304.100
ELSE ORH1F305.3785
WMEALT(I) = 0.0 ONT3F304.102
ENDIF ORH1F305.3786
IF (L_ICEFREEDR) THEN ODC1F405.208
WMEALT(I) = WME_ICE(I) ONT3F304.105
ENDIF ORH1F305.3788
ENDIF ORH1F305.3789
ELSE ONT3F304.108
QFLUXALT(I) = OCEANHEATFLUX(I) OJC2F400.151
IF (L_OSOLAR.OR.L_OSOLARAL) THEN ORH1F305.3790
SOLALT(I) = SOL(I) ONT3F304.111
ENDIF ORH1F305.3791
IF (L_OMIXLAY) THEN ORH1F305.3792
WMEALT(I) = WME(I) ONT3F304.114
ENDIF ORH1F305.3793
ENDIF ONT3F304.116
100 CONTINUE ONT3F304.117
ONT3F304.118
RETURN ONT3F304.119
END ONT3F304.120
*ENDIF ONT3F304.121