*IF DEF,S40_1A TRNS2A1A.2
C ******************************COPYRIGHT****************************** GTS2F400.10639
C (c) CROWN COPYRIGHT 1995, METEOROLOGICAL OFFICE, All Rights Reserved. GTS2F400.10640
C GTS2F400.10641
C Use, duplication or disclosure of this code is subject to the GTS2F400.10642
C restrictions as set forth in the contract. GTS2F400.10643
C GTS2F400.10644
C Meteorological Office GTS2F400.10645
C London Road GTS2F400.10646
C BRACKNELL GTS2F400.10647
C Berkshire UK GTS2F400.10648
C RG12 2SZ GTS2F400.10649
C GTS2F400.10650
C If no contract has been raised with this copy of the code, the use, GTS2F400.10651
C duplication or disclosure of it is strictly prohibited. Permission GTS2F400.10652
C to do so must first be obtained in writing from the Head of Numerical GTS2F400.10653
C Modelling at the above address. GTS2F400.10654
C ******************************COPYRIGHT****************************** GTS2F400.10655
C GTS2F400.10656
C*LL TRNS2A1A.3
CLL SUBROUTINE TRANSS2A TRNS2A1A.4
CLL ------------------- TRNS2A1A.5
CLL TRNS2A1A.6
CLL THIS ROUTINE TRANSFERS DATA NEEDED FOR COUPLING TRNS2A1A.7
CLL FROM THE SLAB OCEAN TO THE ATMOSPHERE,PERFORMING VARIOUS TRNS2A1A.8
CLL MANIPULATIONS ON THE WAY. SJT1F304.749
CLL THE SNOW DEPTH FIELD OVER LAND POINTS CONTAINS LAND SNOW SJT1F304.750
CLL DEPTHS USED BY THE ATMOSPHERE MODEL AND SO IS NOT RESET SJT1F304.751
CLL IT CAN BE COMPILED BY CFT77, BUT DOES SJT1F304.752
CLL NOT CONFORM TO THE ANSI FORTRAN77 STANDARDS TRNS2A1A.10
CLL BECAUSE OF IN_LINE COMMENTS AND THE USE OF ENDDO SJT1F304.753
CLL CALLED BY: SLABCNTL TRNS2A1A.11
CLL VERSION NUMBER 1.1 TRNS2A1A.12
CLL WRITTEN BY D L ROBERTS (14/1/91) TRNS2A1A.13
CLL MODIFIED BY A.B.KEEN (02/02/93) TRNS2A1A.14
CLL MODIFIED BY C.A.SENIOR (25/02/94) SJT1F304.754
CLL REVIEWED BY W.INGRAM (01/03/93) TRNS2A1A.15
CLL FOLLOWS DOCUMENTATION PAPER 3, VERSION 5 FOR STANDARDS. TRNS2A1A.16
CLL TRNS2A1A.17
CLLEND TRNS2A1A.18
C*L TRNS2A1A.19
C----------------------------------------------------------------- TRNS2A1A.20
SUBROUTINE TRANSS2A(L1,L2, 1SJT1F304.755
+ LAND, SJT1F304.756
+ SLABTEMP, SJT1F304.757
+ TSTARATM, SJT1F304.758
+ AICESLB, SJT1F304.759
+ AICEATM, SJT1F304.760
+ HICESLB, SJT1F304.761
+ HICEATM, SJT1F304.762
+ HICEMIN, SJT1F304.763
+ HSNOWSLB, SJT1F304.764
+ HSNOWATM, SJT1F304.765
+ AICEMIN) SJT1F304.766
C TRNS2A1A.23
C THE FLOW OF CONTROL IS STRAIGHTFORWARD. TRNS2A1A.24
C TRNS2A1A.27
INTEGER TRNS2A1A.28
+ L1, ! IN SIZE OF DATA VECTORS TRNS2A1A.29
+ L2 ! IN AMOUNT OF DATA TO BE PROCESSED TRNS2A1A.30
C TRNS2A1A.31
LOGICAL LAND(L1) ! IN ATMOSPHERIC MODEL LAND-SEA SJT1F304.767
+ ! MASK (FALSE AT OCEAN POINTS). TRNS2A1A.33
C TRNS2A1A.34
REAL TRNS2A1A.35
+ SLABTEMP(L1) ! IN TEMPERATURE OF OCEAN SURFACE LAYER SJT1F304.768
+,TSTARATM(L1) ! INOUT SURFACE TEMPERATURE OF ATMOSPHERIC MODEL SJT1F304.769
+,AICESLB(L1) ! IN ICE CONCENTRATION FROM SLAB SJT1F304.770
+,AICEATM(L1) ! INOUT ICE CONCENTRATION IN ATMOS MODEL SJT1F304.771
+,HICESLB(L1) ! IN ICE DEPTH FROM SLAB SJT1F304.772
+,HICEATM(L1) ! OUT ICE DEPTH IN ATMOSPHERIC MODEL SJT1F304.773
+,HSNOWSLB(L1) ! IN SNOW DEPTH FROM SLAB SJT1F304.774
+,HSNOWATM(L1) ! INOUT SNOW DEPTH IN ATMOSPHERIC MODEL SJT1F304.775
C* TRNS2A1A.44
REAL TRNS2A1A.45
+ AICEMIN ! IN MIN ICE CONCENTRATION IF ICE PRESENT TRNS2A1A.46
+,HICEMIN ! IN MINIMUM DEPTH OF ICE IF ICE PRESENT SJT1F304.776
+ ! PREVENTS SMALL ICE DEPTHS CAUSING FAILURE SJT1F304.777
C TRNS2A1A.47
C Include COMDECKS TRNS2A1A.48
C TRNS2A1A.49
*CALL C_SLAB 
TRNS2A1A.50
*CALL C_0_DG_C TRNS2A1A.51
*CALL C_MDI TRNS2A1A.52
C TRNS2A1A.53
C SJT1F304.778
C LOCAL VARIABLES TRNS2A1A.54
C TRNS2A1A.55
INTEGER SJT1F304.779
+ J ! LOOP COUNTER SJT1F304.780
REAL SJT1F304.781
+ ONEEM8 ! SMALL +VE VALUE TO ELIMINATE ROUNDING SJT1F304.782
C TRNS2A1A.57
PARAMETER(ONEEM8 = 1.0E-08 ) SJT1F304.783
C TRNS2A1A.60
C TRNS2A1A.61
C TRNS2A1A.62
C ---------------------------------------------------------- TRNS2A1A.63
C TRNS2A1A.64
DO J = 1,L2 SJT1F304.784
C TRNS2A1A.66
IF ( .NOT. LAND(J) ) THEN SJT1F304.785
C TRNS2A1A.70
CL 1. ICE DEPTH. SJT1F304.786
C TRNS2A1A.83
C CONVERT FROM THE GRID BOX MEAN ACTUAL ICE DEPTH TO THE TRNS2A1A.84
C EQUIVALENT ICE DEPTH AVERAGED OVER ICE AREA. SJT1F304.787
C THIS PROCESS USES THE ICE CONCENTRATION AND SNOW DEPTH FIELDS. TRNS2A1A.86
C NOTE THAT AN EXTRA PIECE OF ICE OF DEPTH HICEMIN IS ADDED SJT1F304.788
C TO PREVENT VERY SMALL ICE DEPTHS OCCURING WHICH CAN CAUSE SJT1F304.789
C FAILURE IN THE ATMOSPHERE MODEL SJT1F304.790
C TRNS2A1A.87
IF ( AICESLB(J) .LT. ( AICEMIN - ONEEM8 ) ) THEN SJT1F304.791
HICEATM(J) = 0.0 SJT1F304.792
ELSE TRNS2A1A.92
HICEATM(J) = ( HICESLB(J) + HICEMIN )/ AICESLB(J) SJT1F304.793
+ + CONRATIO * HSNOWSLB(J) SJT1F304.794
ENDIF TRNS2A1A.95
C TRNS2A1A.100
CL 2. SNOW DEPTH. SJT1F304.795
C TRNS2A1A.103
C NOTE THAT THIS HAS TO BE CONVERTED FROM M TO KG/M**2. TRNS2A1A.104
C TRNS2A1A.105
C THE SNOW DEPTH OVER LAND POINTS IS NOT SET TO 'RMDI' TRNS2A1A.106
C BECAUSE THIS FIELD CONTAINS LAND SNOW DEPTHS USED BY THE SJT1F304.796
C ATMOSPHERE MODEL. THIS IS NOT THE CASE FOR THE OTHER TRNS2A1A.108
C VARIABLES PASSED TO THE ATMOSPHERE, WHICH ARE INITIALISED FOR TRNS2A1A.109
C SAFETY! TRNS2A1A.110
C TRNS2A1A.111
HSNOWATM(J) = HSNOWSLB(J) * RHOSNOW SJT1F304.797
C TRNS2A1A.117
CL 3. SEA SURFACE TEMPERATURE. SJT1F304.798
C TRNS2A1A.120
C NOTE THAT THIS HAS TO BE CONVERTED FROM CELSIUS TO KELVIN. TRNS2A1A.121
C TRNS2A1A.122
C AT SEA-ICE POINTS, THE GRID BOX MEAN SURFACE TEMPERATURE IS TRNS2A1A.123
C ALTERED IN SUCH A WAY THAT THE SURFACE TEMPERATURE OF THE ICY TRNS2A1A.124
C PORTION OF THE BOX IS THE SAME AS IT WAS AT THE END OF THE LAST TRNS2A1A.125
C ATMOSPHERIC PHASE. (NOTE - THIS IS NON-CONSERVATIVE) SJT1F304.799
C HOWEVER, IF ICE APPEARED DURING THE SJT1F304.800
C MOST RECENT OCEAN PHASE, ITS TEMPERATURE IS INITIALISED AT THE TRNS2A1A.127
C FREEZING POINT OF SEAWATER. SJT1F304.801
C THIS CODE USES THE OLD VALUES OF ICE CONCENTRATION. SJT1F304.802
C TRNS2A1A.131
IF ( AICESLB(J) .EQ. 0.0 ) THEN SJT1F304.803
TSTARATM(J) = SLABTEMP(J) + ZERODEGC SJT1F304.804
ELSEIF ( AICEATM(J) .GE. SJT1F304.805
+ ( AICEMIN - ONEEM8 ) ) THEN SJT1F304.806
TSTARATM(J) = TFS + ( AICESLB(J) / AICEATM(J) ) SJT1F304.807
+ *( TSTARATM(J) - TFS ) SJT1F304.808
ELSE TRNS2A1A.139
TSTARATM(J) = TFS SJT1F304.809
ENDIF TRNS2A1A.141
C SJT1F304.810
CL SECTION 4: ICE CONCENTRATION. SJT1F304.811
C SJT1F304.812
C FINALLY, UPDATE THE ICE CONCENTRATION THAT IS PASSED TO SJT1F304.813
C THE ATMOSPHERE MODEL. SJT1F304.814
C SJT1F304.815
C SJT1F304.816
AICEATM(J) = AICESLB(J) SJT1F304.817
C SJT1F304.818
ELSE SJT1F304.819
C SJT1F304.820
HICEATM(J) = RMDI SJT1F304.821
AICEATM(J) = RMDI SJT1F304.822
ENDIF TRNS2A1A.142
END DO TRNS2A1A.143
C TRNS2A1A.144
RETURN TRNS2A1A.146
END TRNS2A1A.147
*ENDIF TRNS2A1A.148