SUBROUTINE ZONMCTL ( 1,4@DYALLOC.3882
*CALL ARGSIZE 
@DYALLOC.3883
*CALL ARGD1 @DYALLOC.3884
*CALL ARGDUMA @DYALLOC.3885
*CALL ARGPTRA @DYALLOC.3886
*CALL ARGCONA @DYALLOC.3887
* BANNER, NF171193.91
* P_FIELDDA, NF171193.92
* P_LEVELSDA, NF171193.93
& ST_LEVELSDA, AJS1F401.1601
& SM_LEVELSDA, AJS1F401.1602
* U_FIELDDA, NF171193.95
* Q_LEVELSDA, NF171193.96
* ICODE, CMESSAGE) NF171193.97
C ZONMCTL1.37
IMPLICIT NONE ZONMCTL1.38
@DYALLOC.3889
*CALL CMAXSIZE @DYALLOC.3890
*CALL TYPSIZE @DYALLOC.3891
*CALL TYPD1 @DYALLOC.3892
*CALL TYPDUMA @DYALLOC.3893
*CALL TYPPTRA @DYALLOC.3894
*CALL TYPCONA @DYALLOC.3895
@DYALLOC.3896
CHARACTER*80 BANNER ! IN - Description of data fields @DYALLOC.3897
INTEGER NF171193.98
& P_FIELDDA, ! IN : Copy of P_FIELD NF171193.99
& P_LEVELSDA, ! IN : Copy of P_LEVELS NF171193.100
& ST_LEVELSDA, ! IN : Copy of ST_LEVELS AJS1F401.1603
& SM_LEVELSDA, ! IN : Copy of SM_LEVELS AJS1F401.1604
& U_FIELDDA, ! IN : Copy of U_FIELD NF171193.102
& Q_LEVELSDA ! IN : Copy of Q_LEVELS NF171193.103
INTEGER ICODE ! OUT - Error exit code ZONMCTL1.42
CHARACTER*(80) CMESSAGE ! OUT - Error message TS150793.227
C ZONMCTL1.44
C*---------------------------------------------------------------------- ZONMCTL1.45
C Common blocks ZONMCTL1.46
C ZONMCTL1.47
*CALL CZONMPTR ZONMCTL1.51
*CALL CZONPRIT ZONMCTL1.52
*CALL C_LHEAT ZONMCTL1.53
*CALL C_R_CP ZONMCTL1.54
C ZONMCTL1.55
C Subroutines called ZONMCTL1.56
C ZONMCTL1.57
EXTERNAL FROM_LAND_POINTS,ZONM_ATM ZONMCTL1.58
C ZONMCTL1.59
C Local variables ZONMCTL1.60
C ZONMCTL1.61
INTEGER ZONMCTL1.62
* I, ! Field index ZONMCTL1.63
* LEVEL, ! Vertical level ZONMCTL1.64
* NLANDPT ! No of land pts ZONMCTL1.65
REAL ZONMCTL1.66
& T(P_FIELDDA,P_LEVELSDA) ! Temperature NF171193.104
& ,SOILT(P_FIELDDA,ST_LEVELSDA) ! Soil temperature AJS1F401.1605
& ,SOILM(P_FIELDDA) ! Soil moisture NF171193.106
& ,CANOPYW(P_FIELDDA) ! Canopy water NF171193.107
& ,PMSL(P_FIELDDA) ! PMSL NF171193.108
& ,SH(P_FIELDDA) ! Sensible Heat NF171193.109
& ,EVAP(P_FIELDDA) ! Evaporation NF171193.110
& ,PPTN(P_FIELDDA) ! Total Precipitation NF171193.111
& ,LSRN(P_FIELDDA) ! Dynamic Rain NF171193.112
& ,CVRN(P_FIELDDA) ! Convective Rain NF171193.113
& ,SNOW(P_FIELDDA) ! Snowfall NF171193.114
& ,TAUX(U_FIELDDA) ! Windstress (X) NF171193.115
& ,TAUY(U_FIELDDA) ! Windstress (Y) NF171193.116
& ,SUBL(P_FIELDDA) ! Sublimation NF171193.117
& ,SOEV(P_FIELDDA) ! Soil evaporation NF171193.118
& ,SOHF(P_FIELDDA) ! Soil Heat flux NF171193.119
& ,THRF(P_FIELDDA) ! Throughfall NF171193.120
& ,SNML(P_FIELDDA) ! Snowmelt NF171193.121
& ,SFRU(P_FIELDDA) ! Surface runoff NF171193.122
& ,SBRU(P_FIELDDA) ! Subsurface runoff NF171193.123
& ,SDTR(P_FIELDDA) ! Surface downward rad total NF171193.124
& ,SDSR(P_FIELDDA) ! Surface downward rad solar NF171193.125
& ,TDTR(P_FIELDDA) ! TOA net downward rad NF171193.126
& ,TOLR(P_FIELDDA) ! TOA net outgoing rad(LW) NF171193.127
& ,TOSW(P_FIELDDA) ! total outgoing solar rad NF171193.128
& ,TISW(P_FIELDDA) ! total incoming solar rad NF171193.129
& ,WFCA(P_FIELDDA) ! net water flux into canopy NF171193.130
& ,CLLIQ(P_FIELDDA,Q_LEVELSDA) ! cloud liquid water NF171193.131
& ,CLICE(P_FIELDDA,Q_LEVELSDA) ! cloud ice water NF171193.132
& ,SIHF(P_FIELDDA) ! Heat flux through sea-ice * ice fraction NF171193.133
& ,SIMH(P_FIELDDA) ! Sea-ice top melting heat flux NF171193.134
& ,SISH(P_FIELDDA) ! sea-ice sensible heat flux NF171193.135
& ,SIST(P_FIELDDA) ! Sea-ice surface total downward radiation NF171193.136
& ,SISS(P_FIELDDA) ! Sea-ice surface solar downward radiation NF171193.137
& ,SIEF(P_FIELDDA) ! Net energy flux into Sea-ice NF171193.138
& ,P_EXNER_FULL ! Exner Pressure at full model level @DYALLOC.3933
& ,P_LEV1 ! Pressure at half level (LEVEL+1) @DYALLOC.3934
& ,P_LEV ! Pressure at half level (LEVEL) @DYALLOC.3935
ZONMCTL1.105
CL---------------------------------------------------------------------- ZONMCTL1.106
*CALL P_EXNERC ZONMCTL1.107
ZONMCTL1.108
CL 1. Gather together the primary fields in the right format, ZONMCTL1.109
CL ie. convert theta to temp, get the area weights on the right ZONMCTL1.110
CL grid, etc. ZONMCTL1.111
CL ZONMCTL1.112
CL ZONMCTL1.113
CL 1.1 Convert THETA to temperature at model levels ZONMCTL1.114
CL ZONMCTL1.115
DO LEVEL=1,P_LEVELS ZONMCTL1.116
DO I=1,P_FIELD ZONMCTL1.117
P_LEV1 = AKH(LEVEL+1) + BKH(LEVEL+1)*D1(JPSTAR+I-1) ZONMCTL1.118
P_LEV = AKH(LEVEL) + BKH(LEVEL) *D1(JPSTAR+I-1) ZONMCTL1.119
P_EXNER_FULL = P_EXNER_C ZONMCTL1.120
* (D1(JP_EXNER(LEVEL+1)+I-1),D1(JP_EXNER(LEVEL)+I-1), ZONMCTL1.121
* P_LEV1,P_LEV,KAPPA) ZONMCTL1.122
T(I,LEVEL) = D1(JTHETA(LEVEL)+I-1) * P_EXNER_FULL ZONMCTL1.123
ENDDO ZONMCTL1.124
ENDDO ZONMCTL1.125
ZONMCTL1.126
CL Expand soil temperature fields to full field ZONMCTL1.127
DO LEVEL=1,ST_LEVELS AJS1F401.1606
CALL FROM_LAND_POINTS(SOILT(1,LEVEL), ZONMCTL1.129
& D1(J_DEEP_SOIL_TEMP(LEVEL)),D1(JLAND),P_FIELD,NLANDPT) @DYALLOC.3936
ENDDO ZONMCTL1.131
ZONMCTL1.132
CL Expand soil moisture and canopy water fields to full field ZONMCTL1.133
CALL FROM_LAND_POINTS(SOILM,D1(JSMC), ZONMCTL1.134
& D1(JLAND),P_FIELD,NLANDPT) @DYALLOC.3937
CALL FROM_LAND_POINTS(CANOPYW,D1(JCANOPY_WATER), ZONMCTL1.136
& D1(JLAND),P_FIELD,NLANDPT) @DYALLOC.3938
CL ZONMCTL1.138
CL 1.2 Copy STASH fields from D1 if present in dump; otherwise set to ZONMCTL1.139
CL large value (printed as *****) ZONMCTL1.140
CL ZONMCTL1.141
IF (JPMSL.EQ.0) THEN ZONMCTL1.142
DO I=1,P_FIELD ZONMCTL1.143
PMSL(I) = 1.E6 ZONMCTL1.144
ENDDO ZONMCTL1.145
ELSE ZONMCTL1.146
DO I=1,P_FIELD ZONMCTL1.147
PMSL(I) = D1(JPMSL+I-1) ZONMCTL1.148
ENDDO ZONMCTL1.149
ENDIF ZONMCTL1.150
ZONMCTL1.151
IF (JEVAP.EQ.0) THEN ZONMCTL1.152
DO I=1,P_FIELD ZONMCTL1.153
EVAP(I) = 1.E6 ZONMCTL1.154
ENDDO ZONMCTL1.155
ELSE ZONMCTL1.156
DO I=1,P_FIELD ZONMCTL1.157
EVAP(I) = D1(JEVAP+I-1) ZONMCTL1.158
ENDDO ZONMCTL1.159
ENDIF ZONMCTL1.160
ZONMCTL1.161
IF (JLSRN.EQ.0) THEN ZONMCTL1.162
DO I=1,P_FIELD ZONMCTL1.163
LSRN(I) = 1.E6 ZONMCTL1.164
ENDDO ZONMCTL1.165
ELSE ZONMCTL1.166
IF (IPPTN.EQ.2) THEN ! convert to accumulation per day ARS2F304.96
DO I=1,P_FIELD ARS2F304.97
LSRN(I) = D1(JLSRN+I-1)*86400.0 ARS2F304.98
ENDDO ARS2F304.99
ELSE ARS2F304.100
DO I=1,P_FIELD ZONMCTL1.167
LSRN(I) = D1(JLSRN+I-1) ZONMCTL1.168
ENDDO ZONMCTL1.169
ENDIF ARS2F304.101
ENDIF ZONMCTL1.170
ZONMCTL1.171
IF (JCVRN.EQ.0) THEN ZONMCTL1.172
DO I=1,P_FIELD ZONMCTL1.173
CVRN(I) = 1.E6 ZONMCTL1.174
ENDDO ZONMCTL1.175
ELSE ZONMCTL1.176
IF (IPPTN.EQ.2) THEN ! convert to accumulation per day ARS2F304.102
DO I=1,P_FIELD ARS2F304.103
CVRN(I) = D1(JCVRN+I-1)*86400.0 ARS2F304.104
ENDDO ARS2F304.105
ELSE ARS2F304.106
DO I=1,P_FIELD ZONMCTL1.177
CVRN(I) = D1(JCVRN+I-1) ZONMCTL1.178
ENDDO ZONMCTL1.179
ENDIF ARS2F304.107
ENDIF ZONMCTL1.180
ZONMCTL1.181
IF (ISNOW.EQ.0) THEN ZONMCTL1.182
DO I=1,P_FIELD ZONMCTL1.183
SNOW(I) = 1.E6 ZONMCTL1.184
ENDDO ZONMCTL1.185
ELSE ZONMCTL1.186
IF (IPPTN.EQ.2) THEN ! convert to accumulation per day ARS2F304.108
DO I=1,P_FIELD ARS2F304.109
SNOW(I) = (D1(JLSSN+I-1)+D1(JCVSN+I-1))*86400.0 ARS2F304.110
ENDDO ARS2F304.111
ELSE ! already accumulation ARS2F304.112
DO I=1,P_FIELD ZONMCTL1.187
SNOW(I) = D1(JLSSN+I-1)+D1(JCVSN+I-1) ZONMCTL1.188
ENDDO ZONMCTL1.189
ENDIF ARS2F304.113
ENDIF ZONMCTL1.190
ZONMCTL1.191
IF (IPPTN.EQ.0) THEN ZONMCTL1.192
DO I=1,P_FIELD ZONMCTL1.193
PPTN(I) = 1.E6 ZONMCTL1.194
ENDDO ZONMCTL1.195
ELSE ! note components already converted to rates ARS2F304.114
DO I=1,P_FIELD ZONMCTL1.197
PPTN(I) = CVRN(I)+LSRN(I)+SNOW(I) ARS2F304.115
ENDDO ZONMCTL1.199
ENDIF ZONMCTL1.200
ZONMCTL1.201
IF (JSH .EQ.0) THEN ZONMCTL1.202
DO I=1,P_FIELD ZONMCTL1.203
SH(I) = 1.E6 ZONMCTL1.204
ENDDO ZONMCTL1.205
ELSE ZONMCTL1.206
DO I=1,P_FIELD ZONMCTL1.207
SH(I) = D1(JSH+I-1) ZONMCTL1.208
ENDDO ZONMCTL1.209
ENDIF ZONMCTL1.210
ZONMCTL1.211
IF (JTAUX.EQ.0) THEN ZONMCTL1.212
DO I=1,U_FIELD ZONMCTL1.213
TAUX(I) = 1.E6 ZONMCTL1.214
ENDDO ZONMCTL1.215
ELSE ZONMCTL1.216
DO I=1,U_FIELD ZONMCTL1.217
TAUX(I) = D1(JTAUX+I-1) ZONMCTL1.218
ENDDO ZONMCTL1.219
ENDIF ZONMCTL1.220
ZONMCTL1.221
IF (JTAUY.EQ.0) THEN ZONMCTL1.222
DO I=1,U_FIELD ZONMCTL1.223
TAUY(I) = 1.E6 ZONMCTL1.224
ENDDO ZONMCTL1.225
ELSE ZONMCTL1.226
DO I=1,U_FIELD ZONMCTL1.227
TAUY(I) = D1(JTAUY+I-1) ZONMCTL1.228
ENDDO ZONMCTL1.229
ENDIF ZONMCTL1.230
ZONMCTL1.231
IF (JSUBL.EQ.0) THEN ZONMCTL1.232
DO I=1,P_FIELD ZONMCTL1.233
SUBL(I) = 1.E6 ZONMCTL1.234
ENDDO ZONMCTL1.235
ELSE ZONMCTL1.236
DO I=1,P_FIELD ZONMCTL1.237
SUBL(I) = D1(JSUBL+I-1) ZONMCTL1.238
ENDDO ZONMCTL1.239
ENDIF ZONMCTL1.240
ZONMCTL1.241
IF (JSOEV.EQ.0) THEN ZONMCTL1.242
DO I=1,P_FIELD ZONMCTL1.243
SOEV(I) = 1.E6 ZONMCTL1.244
ENDDO ZONMCTL1.245
ELSE ZONMCTL1.246
DO I=1,P_FIELD ZONMCTL1.247
SOEV(I) = D1(JSOEV+I-1) ZONMCTL1.248
ENDDO ZONMCTL1.249
ENDIF ZONMCTL1.250
CL ZONMCTL1.251
CL Extra fields only initialised if requested ZONMCTL1.252
IF (IPRTEXTRA.EQ.1) THEN ZONMCTL1.253
CL Extra fields from Hydrology ONLY initialised if all fields requested ZONMCTL1.254
CL ZONMCTL1.255
IF (JSFRU.NE.0.AND.JSBRU.NE.0) THEN ARS2F304.116
DO I=1,P_FIELD ZONMCTL1.257
SFRU(I) = D1(JSFRU+I-1) ZONMCTL1.258
SBRU(I) = D1(JSBRU+I-1) ZONMCTL1.259
ENDDO ARS2F304.117
ENDIF ARS2F304.118
IF (IHYDRO.GT.1) THEN ARS2F304.119
DO I=1,P_FIELD ARS2F304.120
THRF(I) = D1(JTHRF+I-1) ZONMCTL1.260
SNML(I) = D1(JSNML+I-1) ZONMCTL1.261
ENDDO ZONMCTL1.262
ENDIF ZONMCTL1.263
CL ZONMCTL1.264
CL Radiation ONLY initialised if all fields available ZONMCTL1.265
CL ZONMCTL1.266
IF (IRAD.GT.1) THEN ARS2F304.121
DO I=1,P_FIELD ZONMCTL1.268
TDTR(I) = D1(JTISW+I-1) - D1(JTOSW+I-1) ZONMCTL1.269
* - D1(JTOLR+I-1) ZONMCTL1.270
SDTR(I) = D1(JSDSR+I-1)+D1(JSDLR+I-1) ZONMCTL1.271
SDSR(I) = D1(JSDSR+I-1) ZONMCTL1.272
TOLR(I) = D1(JTOLR+I-1) ZONMCTL1.273
TOSW(I) = D1(JTOSW+I-1) ZONMCTL1.274
TISW(I) = D1(JTISW+I-1) ZONMCTL1.275
ENDDO ZONMCTL1.276
ENDIF ZONMCTL1.277
CL ZONMCTL1.278
CL Cloud water ONLY initialised if available ZONMCTL1.279
CL ZONMCTL1.280
IF (ICLOUD.GT.1) THEN ARS2F304.122
DO LEVEL=1,Q_LEVELS ZONMCTL1.282
DO I=1,P_FIELD ZONMCTL1.283
CLLIQ(I,LEVEL) = D1(JCLLIQ+P_FIELD*(LEVEL-1)+I-1) ZONMCTL1.284
CLICE(I,LEVEL) = D1(JCLICE+P_FIELD*(LEVEL-1)+I-1) ZONMCTL1.285
ENDDO ZONMCTL1.286
ENDDO ZONMCTL1.287
ENDIF ZONMCTL1.288
C ZONMCTL1.289
CL Fluxes over land + net energy flux into atmosphere ZONMCTL1.290
CL only initialised if all can be calculated ZONMCTL1.291
C ZONMCTL1.292
IF (IFLUXL.GT.1) THEN ARS2F304.123
DO I=1,P_FIELD ZONMCTL1.294
SOHF(I) = D1(JSOHF+I-1) ! Soil heat flux between top 2 layers ZONMCTL1.295
WFCA(I)= PPTN(I)-D1(JCAEV+I-1)-THRF(I) ! water flux into canopy ZONMCTL1.296
ENDDO ZONMCTL1.297
ENDIF ZONMCTL1.298
CL ZONMCTL1.299
CL Sea-Ice fluxes Only initialised if ZONMCTL1.300
CL only initialised if all can be calculated ZONMCTL1.301
C ZONMCTL1.302
IF (ISICE.GT.1) THEN ARS2F304.124
DO I=1,P_FIELD ZONMCTL1.304
SIHF(I) = D1(JSIHF+I-1) ZONMCTL1.305
SIMH(I) = D1(JSIMH+I-1) ZONMCTL1.306
SISH(I) = D1(JSH+I-1)*D1(JICE_FRACTION+I-1) ZONMCTL1.307
SIST(I) = SDTR(I)*D1(JICE_FRACTION+I-1) ZONMCTL1.308
SISS(I) = SDSR(I)*D1(JICE_FRACTION+I-1) ZONMCTL1.309
SIEF(I) = (SDTR(I)-D1(JSH+I-1)-(LC+LF)*SUBL(I)/84600.)* ZONMCTL1.310
& D1(JICE_FRACTION+I-1) -D1(JSIHF+I-1) - D1(JSIMH+I-1) ZONMCTL1.311
ENDDO ZONMCTL1.312
ENDIF ZONMCTL1.313
ENDIF ZONMCTL1.314
CL ZONMCTL1.315
CL 1.3 Call zonal mean routine (initially, only allowed for up to ZONMCTL1.316
CL 20 vertical levels - general version available later) ZONMCTL1.317
CL ZONMCTL1.318
CALL ZONM_ATM( PMSL, ZONMCTL1.319
1 D1(JPSTAR), T, D1(JQ(1)), D1(JU(1)), D1(JV(1)), ZONMCTL1.320
2 D1(JTSTAR), SOILT, ZONMCTL1.321
3 SOILM, D1(JSNODEP), CANOPYW, ZONMCTL1.322
4 SH ,EVAP , PPTN, ZONMCTL1.323
5 LSRN ,CVRN , SNOW, ZONMCTL1.324
& D1(JICE_FRACTION), D1(JICE_THICKNESS), ZONMCTL1.325
6 TAUX ,TAUY , ZONMCTL1.326
* SUBL, SOEV, SOHF, THRF,SNML, SFRU,SBRU, ZONMCTL1.327
* SDTR, SDSR, TDTR, TOLR, TOSW, TISW, ZONMCTL1.328
* WFCA, ZONMCTL1.329
* CLLIQ,CLICE, ZONMCTL1.330
* SIHF,SIMH,SISH,SIST,SISS,SIEF, ZONMCTL1.331
7 D1(JLAND), @DYALLOC.3939
8 A_LEVDEPC(JDELTA_AK),A_LEVDEPC(JDELTA_BK), ZONMCTL1.333
9 COS_P_LATITUDE, COS_U_LATITUDE, ZONMCTL1.334
A ROW_LENGTH, P_ROWS, U_ROWS, P_LEVELS, Q_LEVELS, ZONMCTL1.335
& ST_LEVELS, SM_LEVELS, P_FIELD, U_FIELD, BANNER, AJS1F401.1607
& IHYDRO,IRAD,ICLOUD,IFLUXL,ISICE, ZONMCTL1.337
& IPRTWIND,IPRTKE,IPRTQ,IPRTTEMP,IPRTVAR,IPRTEXTRA, ZONMCTL1.338
C ICODE, CMESSAGE) ZONMCTL1.339
C ZONMCTL1.340
RETURN ZONMCTL1.341
CL---------------------------------------------------------------------- ZONMCTL1.342
END ZONMCTL1.343
*ENDIF ZONMCTL1.344