*IF DEF,A87_1A ZONMAT1A.2
C ******************************COPYRIGHT****************************** GTS2F400.12223
C (c) CROWN COPYRIGHT 1995, METEOROLOGICAL OFFICE, All Rights Reserved. GTS2F400.12224
C GTS2F400.12225
C Use, duplication or disclosure of this code is subject to the GTS2F400.12226
C restrictions as set forth in the contract. GTS2F400.12227
C GTS2F400.12228
C Meteorological Office GTS2F400.12229
C London Road GTS2F400.12230
C BRACKNELL GTS2F400.12231
C Berkshire UK GTS2F400.12232
C RG12 2SZ GTS2F400.12233
C GTS2F400.12234
C If no contract has been raised with this copy of the code, the use, GTS2F400.12235
C duplication or disclosure of it is strictly prohibited. Permission GTS2F400.12236
C to do so must first be obtained in writing from the Head of Numerical GTS2F400.12237
C Modelling at the above address. GTS2F400.12238
C ******************************COPYRIGHT****************************** GTS2F400.12239
C GTS2F400.12240
CLL Routine: ZONM_ATM ------------------------------------------------- ZONMAT1A.3
CLL ZONMAT1A.4
CLL Purpose: Calculates zonal mean, quarter global mean and global ZONMAT1A.5
CLL mean values from atmospheric prognostic fields, and ZONMAT1A.6
CLL prints formatted output summary on UNIT 6. ZONMAT1A.7
CLL NB: zonal means are performed on p-grid or u-grid as ZONMAT1A.8
CLL appropriate. Quarter global means on p-grid include ZONMAT1A.9
CLL overlap row(s). Formatted output assumes up to 20 ZONMAT1A.10
CLL levels initially (general no. of levels to be allowed ZONMAT1A.11
CLL in future). ZONMAT1A.12
CLL ZONMAT1A.13
CLL Tested under compiler: cft77 ZONMAT1A.14
CLL Tested under OS version: UNICOS 5.0 ZONMAT1A.15
CLL ZONMAT1A.16
CLL Author N.Farnon ZONMAT1A.17
CLL ZONMAT1A.18
CLL Model Modification history from model version 3.0: ZONMAT1A.19
CLL version Date ZONMAT1A.20
CLL 3.2 13/07/93 Changed CHARACTER*(*) to CHARACTER*(80) for TS150793.222
CLL portability. Author Tracey Smith. TS150793.223
CLL 3.4 14/09/94 Reduce workspace required by altered definitions ARS2F304.125
CLL of IHYDRO, IRAD etc. ARS2F304.126
CLL 4.0 24/10/95 Simplify do loop structure to remove Cray ARR9F400.1
CLL compiler optimisation mis-translation and improve ARR9F400.2
CLL vectorisation. R.Rawlins. ARR9F400.3
!LL 4.3 5/3/97 Correct error in zonal mean print. R A Stratton. AAS1F403.1
CLL ZONMAT1A.21
CLL Programming standard: UM Doc Paper 3, version 1 (15/1/90) ZONMAT1A.22
CLL ZONMAT1A.23
CLL Logical components covered: D61,D64 ZONMAT1A.24
CLL ZONMAT1A.25
CLL Project task: C61 ZONMAT1A.26
CLL ZONMAT1A.27
CLL External documentation: UM document C61 - Zonal mean calculations. ZONMAT1A.28
CLL ZONMAT1A.29
CLLEND------------------------------------------------------------------ ZONMAT1A.30
C ZONMAT1A.31
C*L Interface and arguments: ------------------------------------------ ZONMAT1A.32
C ZONMAT1A.33
SUBROUTINE ZONM_ATM ( PMSL, 1,223ZONMAT1A.34
1 PSTAR, T, Q, ZONMAT1A.35
2 U, V, TSTAR, SOILT, SOILM, SNOWD, CANOPYW, ZONMAT1A.36
3 SH, EVAP, PPTN, LSRN, CVRN, SNOW, AICE, HICE, ZONMAT1A.37
& TAUX, TAUY, ZONMAT1A.38
* SUBL, SOEV, SOHF, THRF,SNML, SFRU,SBRU, ZONMAT1A.39
* SDTR, SDSR, TDTR, TOLR, TOSW, TISW, ZONMAT1A.40
* WFCA, ZONMAT1A.41
* CLLIQ,CLICE, ZONMAT1A.42
* SIHF,SIMH,SISH,SIST,SISS,SIEF, ZONMAT1A.43
4 LAND, ZONMAT1A.44
5 DELTA_AK, DELTA_BK, ZONMAT1A.45
6 COS_P_LATITUDE, COS_U_LATITUDE, ZONMAT1A.46
7 ROW_LENGTH, P_ROWS, U_ROWS, P_LEVELS, Q_LEVELS, ZONMAT1A.47
8 ST_LEVELS, SM_LEVELS, P_FIELD, U_FIELD, BANNER, AJS1F401.1533
& IHYDRO,IRAD,ICLOUD,IFLUXL,ISICE, ZONMAT1A.49
* IPRTWIND,IPRTKE,IPRTQ,IPRTTEMP,IPRTVAR,IPRTEXTRA, ZONMAT1A.50
9 ICODE, CMESSAGE) ZONMAT1A.51
C ZONMAT1A.52
IMPLICIT NONE ZONMAT1A.53
C ZONMAT1A.54
*CALL C_LHEAT 
ZONMAT1A.55
C ZONMAT1A.56
C Parameters set by user interface. ZONMAT1A.57
INTEGER ZONMAT1A.58
1 ROW_LENGTH, ! IN - Number of points per row ZONMAT1A.59
2 P_ROWS, ! IN - Number of rows (p-grid) ZONMAT1A.60
3 P_LEVELS, ! IN - Number of pressure levels ZONMAT1A.61
4 Q_LEVELS, ! IN - Number of wet levels ZONMAT1A.62
! levels in multilevel hydrology scheme: AJS1F401.1534
5 ST_LEVELS, ! IN - Number of soil layers AJS1F401.1535
6 SM_LEVELS ! IN - Number of soil moisture AJS1F401.1536
C Parameters derived from those above. ZONMAT1A.64
INTEGER ZONMAT1A.65
1 P_FIELD, ! IN - Size of p-grid 2D field ZONMAT1A.66
2 U_ROWS, ! IN - Number of rows (u-grid) ZONMAT1A.67
3 U_FIELD ! IN - Size of u-grid 2D field ZONMAT1A.68
INTEGER ARS2F304.127
& IHYDRO, ! IN - Hydrology printout ARS2F304.128
& IRAD, ! IN - Radiation printout ARS2F304.129
& ICLOUD, ! IN - Cloud water printout ARS2F304.130
& IFLUXL, ! IN - Fluxes over land ARS2F304.131
& ISICE ! IN - sea-ice ARS2F304.132
LOGICAL ZONMAT1A.69
1 LAND(P_FIELD) ! IN - Land-sea mask (p-grid) ZONMAT1A.70
REAL ZONMAT1A.71
1 PMSL(P_FIELD), ! IN - Pressure at mean sea-level ZONMAT1A.72
2 PSTAR(P_FIELD), ! IN - Surface pressure ZONMAT1A.73
3 T(P_FIELD,P_LEVELS), ! IN - Atmospheric temperature ZONMAT1A.74
4 Q(P_FIELD,Q_LEVELS), ! IN - Atmospheric humidity ZONMAT1A.75
5 U(U_FIELD,P_LEVELS), ! IN - Atmospheric u-velocity ZONMAT1A.76
6 V(U_FIELD,P_LEVELS) ! IN - Atmospheric v-velocity ZONMAT1A.77
REAL ZONMAT1A.78
1 TSTAR(P_FIELD), ! IN - Surface temperature ZONMAT1A.79
2 SOILT(P_FIELD,ST_LEVELS),! IN - Soil temperature AJS1F401.1537
3 SOILM(P_FIELD), ! IN - Soil moisture ZONMAT1A.81
4 SNOWD(P_FIELD), ! IN - Snow depth ZONMAT1A.82
5 CANOPYW(P_FIELD), ! IN - Canopy water ZONMAT1A.83
6 SH(P_FIELD), ! IN - Surface sensible heat ZONMAT1A.84
7 EVAP(P_FIELD), ! IN - Surface evaporation ZONMAT1A.85
8 PPTN(P_FIELD), ! IN - Surface total precipitation ZONMAT1A.86
9 LSRN(P_FIELD), ! IN - Surface dynamic rain ZONMAT1A.87
A CVRN(P_FIELD), ! IN - Surface convective rain ZONMAT1A.88
B SNOW(P_FIELD), ! IN - Surface snowfall ZONMAT1A.89
& AICE(P_FIELD), ! IN - Seaice fraction (sea points) ZONMAT1A.90
& HICE(P_FIELD), ! IN - Seaice thickness (sea points) ZONMAT1A.91
C TAUX(U_FIELD), ! IN - Surface windstress (x) ZONMAT1A.92
D TAUY(U_FIELD), ! IN - Surface windstress (y) ZONMAT1A.93
& SUBL(P_FIELD), ! IN - Sublimation ZONMAT1A.94
& SOEV(P_FIELD), ! IN - Soil evaporation ZONMAT1A.95
& THRF(IHYDRO), ! IN - Throughfall ARS2F304.133
& SNML(IHYDRO), ! IN - Snow melt ARS2F304.134
& SFRU(P_FIELD), ! IN - surface runoff ZONMAT1A.98
& SBRU(P_FIELD), ! IN - Sub surface runoff ZONMAT1A.99
& SOHF(IFLUXL), ! IN - Soil heat flux ARS2F304.135
& SDTR(IRAD), ! IN - net surface downward rad ARS2F304.136
& SDSR(IRAD), ! IN - surface solar downward rad ARS2F304.137
& TDTR(IRAD), ! IN - TOA total downward rad ARS2F304.138
& TOLR(IRAD), ! IN - TOA outgoing lw rad ARS2F304.139
& TOSW(IRAD), ! IN - outgoing solar radiation ARS2F304.140
& TISW(IRAD), ! IN - incoming solar radiation ARS2F304.141
& WFCA(IFLUXL), ! IN - water flux into canopy ARS2F304.142
& CLLIQ(ICLOUD,Q_LEVELS), ! IN - cloud liquid water ARS2F304.143
& CLICE(ICLOUD,Q_LEVELS), ! IN - cloud Ice water ARS2F304.144
& SIHF(ISICE), ! IN - sea-ice heat flux ARS2F304.145
& SIMH(ISICE), ! IN - sea-ice melting heat ARS2F304.146
& SISH(ISICE), ! IN - sea-ice sensible heat ARS2F304.147
& SIST(ISICE), ! IN - sea-ice surface total rad ARS2F304.148
& SISS(ISICE), ! IN - sea-ice surface solar ARS2F304.149
& SIEF(ISICE) ! IN - sea-ice energy flux ARS2F304.150
REAL ZONMAT1A.116
1 DELTA_AK(P_LEVELS), ! IN - Hybrid layer thickness A ZONMAT1A.117
2 DELTA_BK(P_LEVELS), ! IN - Hybrid layer thickness B ZONMAT1A.118
3 COS_P_LATITUDE(P_FIELD), ! IN - Cosine latitude (p-grid) ZONMAT1A.119
4 COS_U_LATITUDE(U_FIELD) ! IN - Cosine latitude (u-grid) ZONMAT1A.120
INTEGER ZONMAT1A.127
& IPRTWIND, ! IN - wind printout ZONMAT1A.128
& IPRTKE, ! IN - Kinetic energy printout ZONMAT1A.129
& IPRTQ, ! IN - Q printout ZONMAT1A.130
& IPRTTEMP, ! IN - temp printout ZONMAT1A.131
& IPRTVAR, ! IN - variance printout ZONMAT1A.132
& IPRTEXTRA ! IN - extra printouts ZONMAT1A.133
INTEGER ZONMAT1A.134
1 ICODE ! OUT - Error exit code ZONMAT1A.135
CHARACTER*80 ZONMAT1A.136
1 BANNER ! IN - Description of data fields ZONMAT1A.137
CHARACTER*(80) TS150793.224
1 CMESSAGE ! OUT - Error message ZONMAT1A.139
C ZONMAT1A.140
C*---------------------------------------------------------------------- ZONMAT1A.141
C ZONMAT1A.142
C Subroutines called ZONMAT1A.143
C ZONMAT1A.144
EXTERNAL ZONM,COLM,GLBM,P_TO_UV ZONMAT1A.145
C ZONMAT1A.146
ZONMAT1A.147
C Local variables ZONMAT1A.148
C ZONMAT1A.149
REAL ZONMAT1A.150
& AMASK(P_FIELD), ! LOC - Mask of all ones ZONMAT1A.151
& LMASK(P_FIELD), ! LOC - Land mask (p-grid) ZONMAT1A.152
& SMASK(P_FIELD) ! LOC - Sea mask (p-grid) ZONMAT1A.153
LOGICAL ZONMAT1A.154
1 LAPTS(P_ROWS,2), ! LOC - Marks all rows ZONMAT1A.155
2 LLPTS(P_ROWS,2), ! LOC - Marks rows with land pts ZONMAT1A.156
3 LSPTS(P_ROWS,2), ! LOC - Marks rows with sea pts ZONMAT1A.157
4 LQGAPTS(4,2), ! LOC - Marks all pts for each 1/4 ZONMAT1A.158
5 LQGLPTS(4,2), ! LOC - Marks land pts for each 1/4 ZONMAT1A.159
6 LQGSPTS(4,2), ! LOC - Marks sea pts for each 1/4 ZONMAT1A.160
7 LGAPTS(2), ! LOC - Marks all pts of globe ZONMAT1A.161
8 LGLPTS(2), ! LOC - Marks land pts of globe ZONMAT1A.162
9 LGSPTS(2) ! LOC - Marks sea pts of globe ZONMAT1A.163
INTEGER ZONMAT1A.164
1 I, ! LOC - Field index ZONMAT1A.165
2 ROW, ! LOC - Row number ZONMAT1A.166
3 J, ! LOC - Row length ZONMAT1A.167
4 LEVEL, ! LOC - Vertical level ZONMAT1A.168
8 LPTS(P_ROWS,2), ! LOC - No of land points/row ZONMAT1A.169
9 SPTS(P_ROWS,2), ! LOC - No of sea points/row ZONMAT1A.170
A QGLPTS(4,2), ! LOC - No of land pts/ 1/4 globe ZONMAT1A.171
B QGSPTS(4,2), ! LOC - No of sea pts/ 1/4 globe ZONMAT1A.172
C GLPTS(2), ! LOC - No of land pts/globe ZONMAT1A.173
D GSPTS(2), ! LOC - No of sea pts/globe ZONMAT1A.174
E QUART, ! LOC - Index for each 1/4 globe ZONMAT1A.175
F IGRID ! LOC - Index for p or u-grid ZONMAT1A.176
REAL ZONMAT1A.177
1 P_MASS(P_FIELD,P_LEVELS),! IN - Mass weighting (p-grid) ZONMAT1A.178
2 U_MASS(U_FIELD,P_LEVELS),! LOC - Mass weighting (u-grid) ZONMAT1A.179
3 S_PMASS(P_FIELD), ! LOC - Dummy wgt for surf(p-grid) ZONMAT1A.180
4 S_UMASS(U_FIELD), ! LOC - Dummy wgt for surf(u-grid) ZONMAT1A.181
5 PSTAR_U(U_FIELD) ! LOC - wgt for atmos (u-grid) ZONMAT1A.182
ZONMAT1A.183
C Additional local variables used in the calculation of Energy ZONMAT1A.184
C diagnostics ZONMAT1A.185
ZONMAT1A.186
REAL ZONMAT1A.187
& WORK1(P_FIELD,P_LEVELS), ! 3D work field ZONMAT1A.188
& WORK2(P_FIELD) ! 2D work field ZONMAT1A.189
& ,WORK3(U_FIELD,P_LEVELS) ! 3D work field ZONMAT1A.190
ZONMAT1A.191
C Values used in the calculation of Zonal means ZONMAT1A.192
ZONMAT1A.193
REAL ZONMAT1A.194
& Z_TKE(U_ROWS), ! Total Kinetic Energy ZONMAT1A.195
& Z_ZKE(U_ROWS), ! Zonal Kinetic Energy ZONMAT1A.196
& Z_EKE(U_ROWS), ! Eddy Kinetic Energy ZONMAT1A.197
& Z_VAR_T(P_ROWS), ! Variance of T ZONMAT1A.198
& Z_VAR_Q(P_ROWS) ! Variance of Q ZONMAT1A.199
REAL ZONMAT1A.200
& Z_VAR_T_LEV(P_ROWS,P_LEVELS), ! Variance of T at each level ZONMAT1A.201
& Z_VAR_Q_LEV(P_ROWS,Q_LEVELS), ! Variance of Q at each level ZONMAT1A.202
& Z_TKE_LEV(U_ROWS,P_LEVELS), ! Total Kinetic Energy per level ZONMAT1A.203
& Z_ZKE_LEV(U_ROWS,P_LEVELS), ! Zonal Kinetic Energy per level ZONMAT1A.204
& Z_EKE_LEV(U_ROWS,P_LEVELS) ! Eddy Kinetic Energy per level ZONMAT1A.205
ZONMAT1A.206
C Values used in the calculation of Global means ZONMAT1A.207
ZONMAT1A.208
REAL ZONMAT1A.209
& G_TKE, ! Total Kinetic Energy ZONMAT1A.210
& G_ZKE, ! Zonal Kinetic Energy ZONMAT1A.211
& G_EKE, ! Eddy Kinetic Energy ZONMAT1A.212
& G_VAR_T, ! Variance of T ZONMAT1A.213
& G_VAR_Q ! Variance of Q ZONMAT1A.214
REAL ZONMAT1A.215
& G_VAR_T_LEV(P_LEVELS), ! Variance of T at each level ZONMAT1A.216
& G_VAR_Q_LEV(Q_LEVELS), ! Variance of Q at each level ZONMAT1A.217
& G_TKE_LEV(P_LEVELS), ! Total Kinetic Energy per level ZONMAT1A.218
& G_ZKE_LEV(P_LEVELS), ! Zonal Kinetic Energy per level ZONMAT1A.219
& G_EKE_LEV(P_LEVELS) ! Eddy Kinetic Energy per level ZONMAT1A.220
REAL ZONMAT1A.221
1 TATMOS(P_FIELD), ! LOC - Col mean atmos temperature ZONMAT1A.222
2 QATMOS(P_FIELD), ! LOC - Col mean atmos humidity ZONMAT1A.223
3 UATMOS(U_FIELD), ! LOC - Col mean atmos temperature ZONMAT1A.224
4 VATMOS(U_FIELD), ! LOC - Col mean atmos humidity ZONMAT1A.225
5 CLLIQ_F(ICLOUD), ! LOC - Cloud liquid water ARS2F304.151
6 CLICE_F(ICLOUD) ! LOC - Cloud ice water ARS2F304.152
INTEGER ! NOTE -Start & end pts for each 1/4 ZONMAT1A.228
1 START_ROW(2,4), ! LOC - Start of 1/4 global section ZONMAT1A.229
& END_ROW(2,4), ! LOC - End of 1/4 global section ARR9F400.4
& LAST_ROW ! LOC - p_rows or u_rows ARR9F400.5
REAL ZONMAT1A.231
1 Z_TATMOS(P_ROWS), ! LOC - Zonm of col mn atmos temperat ZONMAT1A.232
2 Z_QATMOS(P_ROWS), ! LOC - Zonm of col mean atmos humidi ZONMAT1A.233
3 Z_UATMOS(U_ROWS), ! LOC - Zonal mn of col mean atmos u- ZONMAT1A.234
4 Z_VATMOS(U_ROWS) ! LOC - Zonal mn of col mean atmos v- ZONMAT1A.235
REAL ZONMAT1A.236
1 Z_T(P_ROWS,P_LEVELS), ! LOC - Atmospheric temperature ZONMAT1A.237
2 Z_Q(P_ROWS,Q_LEVELS), ! LOC - Atmospheric humidity ZONMAT1A.238
3 Z_U(U_ROWS,P_LEVELS), ! LOC - Atmospheric u-velocity ZONMAT1A.239
4 Z_V(U_ROWS,P_LEVELS) ! LOC - Atmospheric v-velocity ZONMAT1A.240
REAL ZONMAT1A.241
1 Z_PMSL(P_ROWS), ! LOC - Pressure at mean sea-level ZONMAT1A.242
2 Z_PSTAR(P_ROWS), ! LOC - Surface pressure ZONMAT1A.243
3 Z_TSTAR(P_ROWS), ! LOC - Surface temperature ZONMAT1A.244
4 Z_SOILT(P_ROWS,ST_LEVELS), ! LOC - Soil temperature AJS1F401.1538
5 Z_SOILM(P_ROWS), ! LOC - Soil moisture ZONMAT1A.246
6 Z_SNOWD(P_ROWS), ! LOC - Snow depth ZONMAT1A.247
7 Z_CANOPYW(P_ROWS), ! LOC - Canopy water ZONMAT1A.248
8 Z_SH(P_ROWS), ! LOC - Surface sensible heat ZONMAT1A.249
9 Z_EVAP(P_ROWS), ! LOC - Surface evaporation ZONMAT1A.250
A Z_PPTN(P_ROWS), ! LOC - Surface total precipitation ZONMAT1A.251
B Z_LSRN(P_ROWS), ! LOC - Surface dynamic rain ZONMAT1A.252
C Z_CVRN(P_ROWS), ! LOC - Surface convective rain ZONMAT1A.253
D Z_SNOW(P_ROWS), ! LOC - Surface snowfall ZONMAT1A.254
E Z_TAUX(U_ROWS), ! LOC - Surface windstress (x) ZONMAT1A.255
F Z_TAUY(U_ROWS), ! LOC - Surface windstress (y) ZONMAT1A.256
G Z_CLLIQ(P_ROWS), ! LOC - Cloud liquid water ZONMAT1A.257
H Z_CLICE(P_ROWS), ! LOC - Cloud ice water ZONMAT1A.258
I Z_TDTR(P_ROWS), ! LOC - TOA net total down rad ZONMAT1A.259
J Z_TOLR(P_ROWS), ! LOC - TOA outgoing LW rad ZONMAT1A.260
K Z_TOSW(P_ROWS), ! LOC - outgoing solar rad ZONMAT1A.261
L Z_TISW(P_ROWS), ! LOC - incoming solar rad ZONMAT1A.262
M Z_ALBEDO(P_ROWS) ! LOC - planetary albedo ZONMAT1A.263
C land only ZONMAT1A.264
REAL ZONMAT1A.265
1 Z_L_TSTAR(P_ROWS), ! LOC - Surface temperature ZONMAT1A.266
2 Z_L_SH(P_ROWS), ! LOC - Surface sensible heat ZONMAT1A.267
3 Z_L_EVAP(P_ROWS), ! LOC - Surface evaporation ZONMAT1A.268
4 Z_L_PPTN(P_ROWS), ! LOC - Surface total precipitation ZONMAT1A.269
5 Z_L_SNOW(P_ROWS), ! LOC - Surface snowfall ZONMAT1A.270
6 Z_L_TAUX(U_ROWS), ! LOC - Surface windstress (x) ZONMAT1A.271
7 Z_L_TAUY(U_ROWS), ! LOC - Surface windstress (y) ZONMAT1A.272
8 Z_L_SUBL(P_ROWS), ! LOC - sublimation ZONMAT1A.273
9 Z_L_SFRU(P_ROWS), ! LOC - surface runoff ZONMAT1A.274
A Z_L_SBRU(P_ROWS), ! LOC - Subsurface runoff ZONMAT1A.275
B Z_L_SDTR(P_ROWS), ! LOC - surface downward total rad ZONMAT1A.276
C Z_L_SDSR(P_ROWS) ! LOC - surface solar downward rad ZONMAT1A.277
C Sea only ZONMAT1A.278
REAL ZONMAT1A.279
1 Z_S_TSTAR(P_ROWS), ! LOC - Surface temperature ZONMAT1A.280
2 Z_S_SH(P_ROWS), ! LOC - Surface sensible heat ZONMAT1A.281
3 Z_S_EVAP(P_ROWS), ! LOC - Surface evaporation ZONMAT1A.282
4 Z_S_PPTN(P_ROWS), ! LOC - Surface total precipitation ZONMAT1A.283
5 Z_S_SNOW(P_ROWS), ! LOC - Surface snowfall ZONMAT1A.284
& Z_S_AICE(P_ROWS), ! LOC - Seaice fraction ZONMAT1A.285
& Z_S_HICE(P_ROWS), ! LOC - Seaice thickness ZONMAT1A.286
6 Z_S_TAUX(U_ROWS), ! LOC - Surface windstress (x) ZONMAT1A.287
7 Z_S_TAUY(U_ROWS), ! LOC - Surface windstress (y) ZONMAT1A.288
8 Z_S_SUBL(P_ROWS), ! LOC - Sublimation ZONMAT1A.289
9 Z_S_SDTR(P_ROWS), ! LOC - surface downward total rad ZONMAT1A.290
A Z_S_SDSR(P_ROWS) ! LOC - surface solar downward rad ZONMAT1A.291
REAL ZONMAT1A.292
1 G_Q_PMSL(4), ! LOC - 1/4 Global mean PMSL ZONMAT1A.293
2 G_Q_PSTAR(4), ! LOC - 1/4 Global mean surface press ZONMAT1A.294
3 G_Q_TSTAR(4), ! LOC - Surface temperature ZONMAT1A.295
4 G_Q_L_TSTAR(4), ! LOC - Surface temperature ZONMAT1A.296
5 G_Q_S_TSTAR(4), ! LOC - Surface temperature ZONMAT1A.297
6 G_Q_SOILT(4,ST_LEVELS), ! LOC - Soil temperature AJS1F401.1539
7 G_Q_SOILM(4), ! LOC - Soil moisture ZONMAT1A.299
8 G_Q_SNOWD(4), ! LOC - Snow depth ZONMAT1A.300
9 G_Q_CANOPYW(4), ! LOC - Canopy water ZONMAT1A.301
A G_Q_SH(4), ! LOC - Surface sensible heat ZONMAT1A.302
B G_Q_EVAP(4), ! LOC - Surface evaporation ZONMAT1A.303
C G_Q_PPTN(4), ! LOC - Surface total precipitation ZONMAT1A.304
D G_Q_LSRN(4), ! LOC - Surface dynamic rain ZONMAT1A.305
E G_Q_CVRN(4), ! LOC - Surface convective rain ZONMAT1A.306
F G_Q_SNOW(4), ! LOC - Surface snowfall ZONMAT1A.307
& G_Q_SUBL(4), ! LOC -Sublimation ZONMAT1A.308
& G_Q_SDTR(4), ! LOC -Net surface downward total rad ZONMAT1A.309
& G_Q_SDSR(4), ! LOC -Net surface downward total rad ZONMAT1A.310
& G_Q_TDTR(4), ! LOC -TOA downward total rad ZONMAT1A.311
& G_Q_TOLR(4), ! LOC -TOA outgoing lw rad ZONMAT1A.312
& G_Q_TOSW(4), ! LOC -outgoing radiation ZONMAT1A.313
& G_Q_TISW(4), ! LOC -incoming radiation ZONMAT1A.314
& G_Q_ALBEDO(4), ! LOC -Planetary albedo ZONMAT1A.315
& G_Q_CLLIQL(4,Q_LEVELS), ! LOC - Cloud liquid water ZONMAT1A.316
& G_Q_CLICEL(4,Q_LEVELS), ! LOC - Cloud Ice water ZONMAT1A.317
& G_Q_CLLIQ(4), ! LOC - Cloud liquid water ZONMAT1A.318
& G_Q_CLICE(4), ! LOC - Cloud Ice water ZONMAT1A.319
& G_Q_SIHF(4), ! LOC - Sea-ice heat flux ZONMAT1A.320
& G_Q_SIMH(4), ! LOC - Sea-ice melting heat ZONMAT1A.321
& G_Q_SISH(4), ! LOC - Sea-ice sensible heat ZONMAT1A.322
& G_Q_SIST(4), ! LOC - Sea-ice surface total rad ZONMAT1A.323
& G_Q_SISS(4), ! LOC - Sea-ice surface solar rad ZONMAT1A.324
& G_Q_SIEF(4) ! LOC - Sea-ice energy flux ZONMAT1A.325
C Land only quarter globe means ZONMAT1A.326
REAL ZONMAT1A.327
& G_Q_L_SH(4), ! LOC - Surface sensible heat ZONMAT1A.328
& G_Q_L_EVAP(4), ! LOC - Surface evaporation ZONMAT1A.329
& G_Q_L_PPTN(4), ! LOC - Surface total precipitation ZONMAT1A.330
& G_Q_L_SNOW(4), ! LOC - Surface snowfall ZONMAT1A.331
& G_Q_L_SUBL(4), ! LOC - Sublimation ZONMAT1A.332
& G_Q_L_SOEV(4), ! LOC - Soil evaporation ZONMAT1A.333
& G_Q_L_THRF(4), ! LOC - throughfall ZONMAT1A.334
& G_Q_L_SNML(4), ! LOC - Snowmelt ZONMAT1A.335
& G_Q_L_SFRU(4), ! LOC - Surface runoff ZONMAT1A.336
& G_Q_L_SBRU(4), ! LOC - Subsurface runoff ZONMAT1A.337
& G_Q_L_SOHF(4), ! LOC -Soil heat flux (top 2 layers) ZONMAT1A.338
& G_Q_L_SDTR(4), ! LOC -Net surface downward total rad ZONMAT1A.339
& G_Q_L_SDSR(4) ! LOC -Net surface downward total rad ZONMAT1A.340
C Land only quarter globe means derived from above ZONMAT1A.341
REAL ZONMAT1A.342
& G_Q_L_WFSS(4), ! LOC - Net water flux surface snow ZONMAT1A.343
& G_Q_L_WFCA(4), ! LOC - Net water flux canopy ZONMAT1A.344
& G_Q_L_WFSO(4), ! LOC - Net water flux into soil ZONMAT1A.345
& G_Q_WAFL(4), ! LOC - Net water flux atmosphere ZONMAT1A.346
& G_Q_L_ENFS(4), ! LOC - Net energy flux into soil ZONMAT1A.347
& G_Q_ENFL(4) ! LOC - Net energy flux atmosphere ZONMAT1A.348
C Sea only quarter globe means ZONMAT1A.349
REAL ZONMAT1A.350
& G_Q_S_SH(4), ! LOC - Surface sensible heat ZONMAT1A.351
& G_Q_S_EVAP(4), ! LOC - Surface evaporation ZONMAT1A.352
& G_Q_S_PPTN(4), ! LOC - Surface total precipitation ZONMAT1A.353
& G_Q_S_SNOW(4), ! LOC - Surface snowfall ZONMAT1A.354
& G_Q_S_SUBL(4), ! LOC - Sublimation ZONMAT1A.355
& G_Q_S_SDTR(4), ! LOC -Net surface downward total rad ZONMAT1A.356
& G_Q_S_SDSR(4) ! LOC -Net surface downward total rad ZONMAT1A.357
C Sea only quarter globe means derived from above ZONMAT1A.358
C REAL ZONMAT1A.359
C & G_Q_S_SH(4) ! LOC - Surface sensible heat ZONMAT1A.360
REAL ZONMAT1A.361
1 G_Q_T(4,P_LEVELS), ! LOC - 1/4Global mean atmos temperat ZONMAT1A.362
2 G_Q_Q(4,Q_LEVELS), ! LOC - 1/4 Global mean atmos humidit ZONMAT1A.363
3 G_Q_U(4,P_LEVELS), ! LOC - Atmospheric u-velocity ZONMAT1A.364
4 G_Q_V(4,P_LEVELS), ! LOC - Atmospheric v-velocity ZONMAT1A.365
5 G_Q_TAUX(4), ! LOC - Surface windstress (x) ZONMAT1A.366
6 G_Q_TAUY(4), ! LOC - Surface windstress (y) ZONMAT1A.367
7 G_Q_L_TAUX(4), ! LOC - Surface windstress (x) ZONMAT1A.368
8 G_Q_L_TAUY(4), ! LOC - Surface windstress (y) ZONMAT1A.369
9 G_Q_S_TAUX(4), ! LOC - Surface windstress (x) ZONMAT1A.370
A G_Q_S_TAUY(4), ! LOC - Surface windstress (y) ZONMAT1A.371
B G_Q_TKE(4), ! LOC - Total Kinetic Energy ZONMAT1A.372
C G_Q_ZKE(4), ! LOC - Zonal Kinetic Energy ZONMAT1A.373
D G_Q_EKE(4) ! LOC - Eddy Kinetic Energy ZONMAT1A.374
REAL ZONMAT1A.375
1 G_Q_TATMOS(4), ! LOC - 1/4 global mn atmos temperatu ZONMAT1A.376
2 G_Q_QATMOS(4), ! LOC - 1/4 global mean atmos humidit ZONMAT1A.377
3 G_Q_UATMOS(4), ! LOC - 1/4 global mean atmos u-winds ZONMAT1A.378
4 G_Q_VATMOS(4) ! LOC - 1/4 global mean atmos v-winds ZONMAT1A.379
REAL ZONMAT1A.380
1 G_PMSL, ! LOC - 1/4 Global mean PMSL ZONMAT1A.381
2 G_PSTAR, ! LOC - 1/4 Global mean surface press ZONMAT1A.382
3 G_TSTAR, ! LOC - Surface temperature ZONMAT1A.383
4 G_L_TSTAR, ! LOC - Surface temperature ZONMAT1A.384
5 G_S_TSTAR, ! LOC - Surface temperature ZONMAT1A.385
6 G_SOILT(ST_LEVELS), ! LOC - Soil temperature AJS1F401.1540
7 G_SOILM, ! LOC - Soil moisture ZONMAT1A.387
8 G_SNOWD, ! LOC - Snow depth ZONMAT1A.388
9 G_CANOPYW, ! LOC - Canopy water ZONMAT1A.389
A G_SH, ! LOC - Surface sensible heat ZONMAT1A.390
B G_EVAP, ! LOC - Surface evaporation ZONMAT1A.391
C G_PPTN, ! LOC - Surface total precipitation ZONMAT1A.392
D G_LSRN, ! LOC - Surface dynamic rain ZONMAT1A.393
E G_CVRN, ! LOC - Surface convective rain ZONMAT1A.394
F G_SNOW ! LOC - Surface snowfall ZONMAT1A.395
REAL ZONMAT1A.396
1 G_T(P_LEVELS), ! LOC - 1/4Global mean atmos temperat ZONMAT1A.397
2 G_Q(Q_LEVELS), ! LOC - 1/4 Global mean atmos humidit ZONMAT1A.398
3 G_U(P_LEVELS), ! LOC - Atmospheric u-velocity ZONMAT1A.399
4 G_V(P_LEVELS), ! LOC - Atmospheric v-velocity ZONMAT1A.400
5 G_TAUX, ! LOC - Surface windstress (x) ZONMAT1A.401
6 G_TAUY, ! LOC - Surface windstress (y) ZONMAT1A.402
7 G_L_TAUX, ! LOC - Surface windstress (x) ZONMAT1A.403
8 G_L_TAUY, ! LOC - Surface windstress (y) ZONMAT1A.404
9 G_S_TAUX, ! LOC - Surface windstress (x) ZONMAT1A.405
A G_S_TAUY, ! LOC - Surface windstress (y) ZONMAT1A.406
& G_SUBL, ! LOC - Sublimation ZONMAT1A.407
& G_SDTR, ! LOC - Surface net total down rad ZONMAT1A.408
& G_SDSR, ! LOC - Surface net solar down rad ZONMAT1A.409
& G_TDTR, ! LOC - TOA net total down rad ZONMAT1A.410
& G_TOLR, ! LOC - TOA out going LW rad ZONMAT1A.411
& G_TOSW, ! LOC - outgoing solar radiation ZONMAT1A.412
& G_TISW, ! LOC - incoming solar radiation ZONMAT1A.413
& G_ALBEDO, ! LOC - Planetary albedo ZONMAT1A.414
& G_CLLIQL(Q_LEVELS), ! LOC - Cloud liquid water ZONMAT1A.415
& G_CLICEL(Q_LEVELS), ! LOC - Cloud Ice water ZONMAT1A.416
& G_CLLIQ, ! LOC - Cloud liquid water ZONMAT1A.417
& G_CLICE, ! LOC - Cloud Ice water ZONMAT1A.418
& G_SIHF, ! LOC - Sea-ice heat flux ZONMAT1A.419
& G_SIMH, ! LOC - Sea-ice melting heat ZONMAT1A.420
& G_SISH, ! LOC - Sea-ice sensible heat ZONMAT1A.421
& G_SIST, ! LOC - Sea-ice surface total rad ZONMAT1A.422
& G_SISS, ! LOC - Sea-ice surface solar rad ZONMAT1A.423
& G_SIEF ! LOC - Sea-ice energy flux ZONMAT1A.424
REAL ZONMAT1A.425
1 G_TATMOS, ! LOC - Global mn atmos temperature ZONMAT1A.426
2 G_QATMOS, ! LOC - Global mean atmos humidity ZONMAT1A.427
3 G_UATMOS, ! LOC - Global mean atmos u-winds ZONMAT1A.428
4 G_VATMOS ! LOC - Global mean atmos v-winds ZONMAT1A.429
C Land only global means ZONMAT1A.430
REAL ZONMAT1A.431
& G_L_SH, ! LOC - Surface sensible heat ZONMAT1A.432
& G_L_EVAP, ! LOC - Surface evaporation ZONMAT1A.433
& G_L_PPTN, ! LOC - Surface total precipitation ZONMAT1A.434
& G_L_SNOW, ! LOC - Surface snowfall ZONMAT1A.435
& G_L_SUBL, ! LOC - Sublimation ZONMAT1A.436
& G_L_SOEV, ! LOC - Soil evaporation ZONMAT1A.437
& G_L_THRF, ! LOC - throughfall ZONMAT1A.438
& G_L_SNML, ! LOC - Snowmelt ZONMAT1A.439
& G_L_SFRU, ! LOC - Surface runoff ZONMAT1A.440
& G_L_SBRU, ! LOC - Subsurface runoff ZONMAT1A.441
& G_L_SOHF, ! LOC -Soil heat flux (top 2 layers) ZONMAT1A.442
& G_L_SDTR, ! LOC -Net surface downward total rad ZONMAT1A.443
& G_L_SDSR ! LOC -Net surface downward total rad ZONMAT1A.444
C Land only global means derived from above ZONMAT1A.445
REAL ZONMAT1A.446
& G_L_WFSS, ! LOC - Net water flux surface snow ZONMAT1A.447
& G_L_WFCA, ! LOC - Net water flux canopy ZONMAT1A.448
& G_L_WFSO, ! LOC - Net water flux into soil ZONMAT1A.449
& G_WAFL, ! LOC - Net water flux atmosphere ZONMAT1A.450
& G_ENFL, ! LOC - Net energy flux atmosphere ZONMAT1A.451
& G_L_ENFS ! LOC - Net energy flux into soil ZONMAT1A.452
C Sea only global means ZONMAT1A.453
REAL ZONMAT1A.454
& G_S_SH, ! LOC - Surface sensible heat ZONMAT1A.455
& G_S_EVAP, ! LOC - Surface evaporation ZONMAT1A.456
& G_S_PPTN, ! LOC - Surface total precipitation ZONMAT1A.457
& G_S_SNOW, ! LOC - Surface snowfall ZONMAT1A.458
& G_S_SUBL, ! LOC - Sublimation ZONMAT1A.459
& G_S_SDTR, ! LOC -Net surface downward total rad ZONMAT1A.460
& G_S_SDSR ! LOC -Net surface downward total rad ZONMAT1A.461
C ZONMAT1A.462
CHARACTER*130 ZONMAT1A.463
& CTITLE ! title for print ZONMAT1A.464
& ,CHEAD1 ! header for table ZONMAT1A.465
& ,CHEAD2 ! second header for table ZONMAT1A.466
CL---------------------------------------------------------------------- ZONMAT1A.467
ZONMAT1A.468
CL 1. Initializing the start and end pt for each 1/4 of globe ZONMAT1A.469
CL ZONMAT1A.470
CL 1.1 For p-grid ZONMAT1A.471
CL ZONMAT1A.472
START_ROW(1,1)=1 ! 1st quarter start ZONMAT1A.473
END_ROW(1,1)=(P_ROWS-1)/3 + 1 ! 1st quarter end ZONMAT1A.474
START_ROW(1,2)=(P_ROWS-1)/3 + 1 ! 2nd quarter start ZONMAT1A.475
END_ROW(1,2)=(P_ROWS-1)/2 + 1 ! 2nd quarter end ZONMAT1A.476
START_ROW(1,3)=(P_ROWS-1)/2 + 1 ! 3rd quarter start ZONMAT1A.477
END_ROW(1,3)=2*(P_ROWS-1)/3 + 1 ! 3rd quarter end ZONMAT1A.478
START_ROW(1,4)=2*(P_ROWS-1)/3 + 1 ! 4th quarter start ZONMAT1A.479
END_ROW(1,4)=P_ROWS ! 4th quarter end ZONMAT1A.480
CL ZONMAT1A.481
CL 1.2 For u-grid ZONMAT1A.482
CL ZONMAT1A.483
START_ROW(2,1)=1 ! 1st quarter start ZONMAT1A.484
END_ROW(2,1)=(U_ROWS)/3 ! 1st quarter end ZONMAT1A.485
START_ROW(2,2)=(U_ROWS)/3 + 1 ! 2nd quarter start ZONMAT1A.486
END_ROW(2,2)=(U_ROWS)/2 ! 2nd quarter end ZONMAT1A.487
START_ROW(2,3)=(U_ROWS)/2 + 1 ! 3rd quarter start ZONMAT1A.488
END_ROW(2,3)=2*(U_ROWS)/3 ! 3rd quarter end ZONMAT1A.489
START_ROW(2,4)=2*(U_ROWS)/3 + 1 ! 4th quarter start ZONMAT1A.490
END_ROW(2,4)=U_ROWS ! 4th quarter end ZONMAT1A.491
CL---------------------------------------------------------------------- ZONMAT1A.492
CL 2. Calculate mass weights and set up masks for weighted sums ZONMAT1A.493
CL ZONMAT1A.494
DO I=1,P_FIELD ZONMAT1A.495
AMASK(I) = 1.0 ! Set mask of all land & sea pts =1 ZONMAT1A.496
IF (LAND(I)) THEN ZONMAT1A.497
LMASK(I) = 1.0 ! Set mask of all land pts =1 ZONMAT1A.498
ELSE ZONMAT1A.499
LMASK(I) = 0.0 ! Set mask of all other pts =0 ZONMAT1A.500
ENDIF ZONMAT1A.501
SMASK(I) = 1.0-LMASK(I) ! Set mask of all sea pts =1 or 0 ZONMAT1A.502
END DO ZONMAT1A.503
CL---------------------------------------------------------------------- ZONMAT1A.504
CL 3. Calculate no of land/sea points on row-by-row basis ZONMAT1A.505
CL and set logical arrays to denote active land/sea rows ZONMAT1A.506
ZONMAT1A.507
DO IGRID=1,2 ! p-grid, then u-grid ARR9F400.6
ARR9F400.7
IF(IGRID.EQ.1) THEN ARR9F400.8
LAST_ROW=P_ROWS ARR9F400.9
ELSE ARR9F400.10
LAST_ROW=U_ROWS ARR9F400.11
ENDIF ARR9F400.12
ARR9F400.13
DO ROW=1,LAST_ROW ! iterate over all rows in grid ARR9F400.14
ARR9F400.15
LPTS(ROW,IGRID) = 0 ! Initialize no. of land pts =0 ARR9F400.16
DO J=1,ROW_LENGTH ! Loop through each pt/row ARR9F400.17
LPTS(ROW,IGRID) = LPTS(ROW,IGRID)+ ARR9F400.18
& LMASK((ROW-1)*ROW_LENGTH+J) ! Sum land pts for grid ARR9F400.19
ENDDO ! J 1,ROW_LENGTH ARR9F400.20
SPTS(ROW,IGRID) = ROW_LENGTH-LPTS(ROW,IGRID) ! No. of sea pts ARR9F400.21
LAPTS(ROW,IGRID) = .TRUE. ! land or sea pts/row?= always true ARR9F400.22
LLPTS(ROW,IGRID) = LPTS(ROW,IGRID).GT.0 ! Any land pts/row=t ARR9F400.23
LSPTS(ROW,IGRID) = SPTS(ROW,IGRID).GT.0 ! Any sea pts/row=t ARR9F400.24
ARR9F400.25
ENDDO ! ROW 1, last row ARR9F400.26
ENDDO ! IGRID 1,2 ARR9F400.27
ARR9F400.28
DO 21 IGRID=1,2 ZONMAT1A.508
GLPTS(IGRID) = 0 ! Initialize no. of global pts =0 ZONMAT1A.509
DO 20 I=1,4 ! Loop through each 1/4 of globe ZONMAT1A.510
QGLPTS(I,IGRID) = 0 ! Init. no. of 1/4 global pts =0 ZONMAT1A.511
DO 30 ROW=START_ROW(IGRID,I),END_ROW(IGRID,I) ZONMAT1A.512
& ! Loop from start to end pt for each ZONMAT1A.513
QGLPTS(I,IGRID) = QGLPTS(I,IGRID)+LPTS(ROW,IGRID)! Land pts/qu ZONMAT1A.523
30 CONTINUE ZONMAT1A.524
QGSPTS(I,IGRID) = ROW_LENGTH * (END_ROW(IGRID,I) ZONMAT1A.525
& - START_ROW(IGRID,I)) - QGLPTS(I,IGRID) ! No sea pts ZONMAT1A.526
LQGAPTS(I,IGRID) = .TRUE. ! All land&seapts/quar ZONMAT1A.527
LQGLPTS(I,IGRID) = QGLPTS(I,IGRID).GT.0 ! All land pts/quart=t ZONMAT1A.528
LQGSPTS(I,IGRID) = QGSPTS(I,IGRID).GT.0 ! All sea pts/quart=tr ZONMAT1A.529
GLPTS(IGRID) = GLPTS(IGRID) + QGLPTS(I,IGRID) ! No. of global la ZONMAT1A.530
20 CONTINUE ZONMAT1A.531
ZONMAT1A.532
GSPTS(IGRID) = ROW_LENGTH * P_ROWS - GLPTS(IGRID) ! No. of global ZONMAT1A.533
LGAPTS(IGRID) = .TRUE. ! All land&sea pts=true fo ZONMAT1A.534
LGLPTS(IGRID) = GLPTS(IGRID) .GT. 0 ! Logical land pts=true fo ZONMAT1A.535
LGSPTS(IGRID) = GSPTS(IGRID) .GT. 0 ! Logical sea pts=true for ZONMAT1A.536
21 CONTINUE ZONMAT1A.537
ZONMAT1A.538
CL---------------------------------------------------------------------- ZONMAT1A.539
CL 4. Compute mass weights from P* and hybrid coordinates ZONMAT1A.540
CL Allow for negative sign of delta_ak and delta_bk ZONMAT1A.541
CL ZONMAT1A.542
DO LEVEL=1,P_LEVELS ZONMAT1A.543
DO I=1,P_FIELD ZONMAT1A.544
P_MASS(I,LEVEL) = -DELTA_AK(LEVEL)-DELTA_BK(LEVEL)*PSTAR(I) ZONMAT1A.545
END DO ZONMAT1A.546
CL ZONMAT1A.547
CL 4.1 Interpolate onto U-Grid ZONMAT1A.548
CL ZONMAT1A.549
CALL P_TO_UV(P_MASS(1,LEVEL),U_MASS(1,LEVEL),P_FIELD,U_FIELD, ZONMAT1A.550
& ROW_LENGTH,P_ROWS) ZONMAT1A.551
END DO ZONMAT1A.552
ZONMAT1A.553
CL ZONMAT1A.554
CL 4.2 Interpolate PSTAR onto U-Grid ZONMAT1A.555
CL ZONMAT1A.556
CALL P_TO_UV(PSTAR,PSTAR_U,P_FIELD,U_FIELD,ROW_LENGTH,P_ROWS) ZONMAT1A.557
CL ZONMAT1A.558
CL 4.3 Set dummy weighting for surface variables to one ZONMAT1A.559
CL Following two DO loops labelled due to fpp translation problem ZONMAT1A.560
CL ZONMAT1A.561
DO 430 I=1,P_FIELD ZONMAT1A.562
S_PMASS(I)=1.0 ! Dummy weight for surface vars(p- ZONMAT1A.563
430 CONTINUE ZONMAT1A.564
ZONMAT1A.565
DO 431 I=1,U_FIELD ZONMAT1A.566
S_UMASS(I)=1.0 ! Dummy weight for surface vars(u- ZONMAT1A.567
431 CONTINUE ZONMAT1A.568
ZONMAT1A.569
CL---------------------------------------------------------------------- ZONMAT1A.570
CL 5. Compute column means of primary atmospheric variables ZONMAT1A.571
CL ZONMAT1A.572
CALL COLM(T,TATMOS,P_MASS,P_ROWS,ROW_LENGTH,P_LEVELS) ZONMAT1A.573
CALL COLM(Q,QATMOS,P_MASS,P_ROWS,ROW_LENGTH,Q_LEVELS) ZONMAT1A.574
CALL COLM(U,UATMOS,U_MASS,U_ROWS,ROW_LENGTH,P_LEVELS) ZONMAT1A.575
CALL COLM(V,VATMOS,U_MASS,U_ROWS,ROW_LENGTH,P_LEVELS) ZONMAT1A.576
CL---------------------------------------------------------------------- ZONMAT1A.577
CL 6. Compute zonal means ZONMAT1A.578
CL ZONMAT1A.579
CL 6.1 Mass weighted on p-grid for surface variables (all points) ZONMAT1A.580
CL ZONMAT1A.581
CALL ZONM(PMSL,Z_PMSL,AMASK,S_PMASS,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.582
& P_ROWS) ZONMAT1A.583
CALL ZONM(PSTAR,Z_PSTAR,AMASK,S_PMASS,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.584
& P_ROWS) ZONMAT1A.585
CALL ZONM(TSTAR,Z_TSTAR,AMASK,S_PMASS,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.586
& P_ROWS) ZONMAT1A.587
CALL ZONM(SNOWD,Z_SNOWD,AMASK,S_PMASS,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.588
& P_ROWS) ZONMAT1A.589
CALL ZONM(SH,Z_SH,AMASK,S_PMASS,LAPTS(1,1),ROW_LENGTH,P_ROWS) ZONMAT1A.590
ZONMAT1A.591
CALL ZONM(EVAP,Z_EVAP,AMASK,S_PMASS,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.592
& P_ROWS) ZONMAT1A.593
CALL ZONM(PPTN,Z_PPTN,AMASK,S_PMASS,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.594
& P_ROWS) ZONMAT1A.595
CALL ZONM(LSRN,Z_LSRN,AMASK,S_PMASS,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.596
& P_ROWS) ZONMAT1A.597
CALL ZONM(CVRN,Z_CVRN,AMASK,S_PMASS,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.598
& P_ROWS) ZONMAT1A.599
CALL ZONM(SNOW,Z_SNOW,AMASK,S_PMASS,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.600
& P_ROWS) ZONMAT1A.601
CL ZONMAT1A.602
CL 6.2 Mass weighted on p-grid for column meaned atmos variables (all po ZONMAT1A.603
CL ZONMAT1A.604
CALL ZONM(TATMOS,Z_TATMOS,AMASK,PSTAR,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.605
& P_ROWS) ZONMAT1A.606
CALL ZONM(QATMOS,Z_QATMOS,AMASK,PSTAR,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.607
& P_ROWS) ZONMAT1A.608
CL ZONMAT1A.609
CL 6.3 Mass weighted on u-grid for column meaned atmos variables (all po ZONMAT1A.610
CL ZONMAT1A.611
CALL ZONM(UATMOS,Z_UATMOS,AMASK,PSTAR_U,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.612
& U_ROWS) ZONMAT1A.613
CALL ZONM(VATMOS,Z_VATMOS,AMASK,PSTAR_U,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.614
& U_ROWS) ZONMAT1A.615
CL ZONMAT1A.616
CL 6.4 Mass weighted on u-grid for surface variables (all points) ZONMAT1A.617
CL ZONMAT1A.618
CALL ZONM(TAUX,Z_TAUX,AMASK,S_UMASS,LAPTS(1,2),ROW_LENGTH, ZONMAT1A.619
& U_ROWS) ZONMAT1A.620
CALL ZONM(TAUY,Z_TAUY,AMASK,S_UMASS,LAPTS(1,2),ROW_LENGTH, ZONMAT1A.621
& U_ROWS) ZONMAT1A.622
CL ZONMAT1A.623
CL 6.5 Mass weighted on p-levels for atmospheric variables (all points) ZONMAT1A.624
CL ZONMAT1A.625
DO LEVEL=1,P_LEVELS ZONMAT1A.626
CALL ZONM(U(1,LEVEL),Z_U(1,LEVEL),AMASK,U_MASS(1,LEVEL), ZONMAT1A.627
& LAPTS(1,2),ROW_LENGTH,U_ROWS) ZONMAT1A.628
CALL ZONM(V(1,LEVEL),Z_V(1,LEVEL),AMASK,U_MASS(1,LEVEL), ZONMAT1A.629
& LAPTS(1,2),ROW_LENGTH,U_ROWS) ZONMAT1A.630
CALL ZONM(T(1,LEVEL),Z_T(1,LEVEL),AMASK,P_MASS(1,LEVEL), ZONMAT1A.631
& LAPTS(1,1),ROW_LENGTH,P_ROWS) ZONMAT1A.632
ZONMAT1A.633
ZONMAT1A.634
CL ZONMAT1A.635
CL 6.6 Mass weighted on p_levels for kinetic energy ZONMAT1A.636
CL ZONMAT1A.637
DO I=1,U_FIELD ZONMAT1A.638
WORK3(I,LEVEL)= ZONMAT1A.639
& U(I,LEVEL)*U(I,LEVEL)+V(I,LEVEL)*V(I,LEVEL) ZONMAT1A.640
END DO ZONMAT1A.641
ZONMAT1A.642
CALL ZONM(WORK3(1,LEVEL),Z_TKE_LEV(1,LEVEL),AMASK, ZONMAT1A.643
& U_MASS(1,LEVEL),LAPTS(1,2),ROW_LENGTH,U_ROWS) ZONMAT1A.644
ZONMAT1A.645
DO I=1,U_ROWS ZONMAT1A.646
Z_ZKE_LEV(I,LEVEL)=Z_U(I,LEVEL)*Z_U(I,LEVEL)+ ZONMAT1A.647
& Z_V(I,LEVEL)*Z_V(I,LEVEL) ZONMAT1A.648
Z_EKE_LEV(I,LEVEL)=Z_TKE_LEV(I,LEVEL)-Z_ZKE_LEV(I,LEVEL) ZONMAT1A.649
END DO ZONMAT1A.650
ZONMAT1A.651
END DO ZONMAT1A.652
ZONMAT1A.653
C Total KE in vertical mean sense, mass weighted ZONMAT1A.654
ZONMAT1A.655
CALL COLM(WORK3,WORK2,U_MASS,U_ROWS,ROW_LENGTH,P_LEVELS) ZONMAT1A.656
CALL ZONM(WORK2,Z_TKE, ZONMAT1A.657
& AMASK,PSTAR_U,LAPTS(1,2),ROW_LENGTH,U_ROWS) ZONMAT1A.658
ZONMAT1A.659
DO LEVEL=1,P_LEVELS ZONMAT1A.660
DO I=1,U_FIELD ZONMAT1A.661
WORK3(I,LEVEL) = Z_ZKE_LEV((I-1)/ROW_LENGTH+1,LEVEL) ZONMAT1A.662
ENDDO ZONMAT1A.663
ENDDO ZONMAT1A.664
ZONMAT1A.665
CALL COLM(WORK3,WORK2,U_MASS,U_ROWS,ROW_LENGTH,P_LEVELS) ZONMAT1A.666
CALL ZONM(WORK2,Z_ZKE,AMASK,PSTAR_U,LAPTS(1,2),ROW_LENGTH, ZONMAT1A.667
& U_ROWS) ZONMAT1A.668
DO I=1,U_ROWS ZONMAT1A.669
Z_EKE(I)=Z_TKE(I)-Z_ZKE(I) ZONMAT1A.670
END DO ZONMAT1A.671
CL ZONMAT1A.672
CL 6.7 Mass weighted on p-levels for temperature variance ZONMAT1A.673
CL ZONMAT1A.674
DO LEVEL=1,P_LEVELS ZONMAT1A.675
ZONMAT1A.676
DO I=1,P_FIELD ZONMAT1A.677
WORK1(I,LEVEL)=T(I,LEVEL)*T(I,LEVEL) ZONMAT1A.678
END DO ZONMAT1A.679
ZONMAT1A.680
CALL ZONM(WORK1(1,LEVEL),WORK2,AMASK,P_MASS(1,LEVEL), ZONMAT1A.681
& LAPTS(1,1),ROW_LENGTH,P_ROWS) ZONMAT1A.682
ZONMAT1A.683
DO I=1,P_ROWS ZONMAT1A.684
Z_VAR_T_LEV(I,LEVEL)=WORK2(I)-Z_T(I,LEVEL)*Z_T(I,LEVEL) ZONMAT1A.685
END DO ZONMAT1A.686
ZONMAT1A.687
END DO ZONMAT1A.688
ZONMAT1A.689
C Vertical mean temperature variance, mass weighted ZONMAT1A.690
ZONMAT1A.691
CALL COLM(WORK1,WORK2,P_MASS,P_ROWS,ROW_LENGTH,P_LEVELS) ZONMAT1A.692
CALL ZONM(WORK2,Z_VAR_T,AMASK,PSTAR,LAPTS(1,1),ROW_LENGTH,P_ROWS) ZONMAT1A.693
ZONMAT1A.694
DO I=1,P_ROWS ZONMAT1A.695
Z_VAR_T(I)=Z_VAR_T(I)-Z_TATMOS(I)*Z_TATMOS(I) ZONMAT1A.696
END DO ZONMAT1A.697
CL ZONMAT1A.698
CL 6.8 Mass weighted on q-levels for moisture and moisture variance ZONMAT1A.699
CL ZONMAT1A.700
DO LEVEL=1,Q_LEVELS ZONMAT1A.701
CALL ZONM(Q(1,LEVEL),Z_Q(1,LEVEL),AMASK,P_MASS(1,LEVEL), ZONMAT1A.702
& LAPTS(1,1),ROW_LENGTH,P_ROWS) ZONMAT1A.703
ZONMAT1A.704
DO I=1,P_FIELD ZONMAT1A.705
WORK1(I,LEVEL)=Q(I,LEVEL)*Q(I,LEVEL) ZONMAT1A.706
END DO ZONMAT1A.707
ZONMAT1A.708
CALL ZONM(WORK1(1,LEVEL),WORK2, ZONMAT1A.709
& AMASK,P_MASS(1,LEVEL),LAPTS(1,1),ROW_LENGTH,P_ROWS) ZONMAT1A.710
ZONMAT1A.711
DO I=1,P_ROWS ZONMAT1A.712
Z_VAR_Q_LEV(I,LEVEL)=WORK2(I)-Z_Q(I,LEVEL)*Z_Q(I,LEVEL) ZONMAT1A.713
END DO ZONMAT1A.714
ZONMAT1A.715
END DO ZONMAT1A.716
ZONMAT1A.717
C Vertical mean moisture variance, mass weighted ZONMAT1A.718
ZONMAT1A.719
CALL COLM(WORK1,WORK2,P_MASS,P_ROWS,ROW_LENGTH,Q_LEVELS) ZONMAT1A.720
CALL ZONM(WORK2,Z_VAR_Q,AMASK,PSTAR,LAPTS(1,1),ROW_LENGTH,P_ROWS) ZONMAT1A.721
ZONMAT1A.722
DO I=1,P_ROWS ZONMAT1A.723
Z_VAR_Q(I)=Z_VAR_Q(I)-Z_QATMOS(I)*Z_QATMOS(I) ZONMAT1A.724
END DO ZONMAT1A.725
CL ZONMAT1A.726
CL 6.9 Mass weighted on p/u-grid for surface variables (land points) ZONMAT1A.727
CL ZONMAT1A.728
DO LEVEL=1,ST_LEVELS AJS1F401.1541
CALL ZONM(SOILT(1,LEVEL),Z_SOILT(1,LEVEL),LMASK,S_PMASS, ZONMAT1A.730
& LLPTS(1,1),ROW_LENGTH,P_ROWS) ZONMAT1A.731
END DO ZONMAT1A.732
ZONMAT1A.733
CALL ZONM(SOILM,Z_SOILM,LMASK,S_PMASS,LLPTS(1,1),ROW_LENGTH, ZONMAT1A.734
& P_ROWS) ZONMAT1A.735
CALL ZONM(TSTAR,Z_L_TSTAR,LMASK,S_PMASS,LLPTS(1,1),ROW_LENGTH, ZONMAT1A.736
& P_ROWS) ZONMAT1A.737
CALL ZONM(CANOPYW,Z_CANOPYW,LMASK,S_PMASS,LLPTS(1,1),ROW_LENGTH, ZONMAT1A.738
& P_ROWS) ZONMAT1A.739
CALL ZONM(SH,Z_L_SH,LMASK,S_PMASS,LLPTS(1,1),ROW_LENGTH,P_ROWS) ZONMAT1A.740
CALL ZONM(EVAP,Z_L_EVAP,LMASK,S_PMASS,LLPTS(1,1),ROW_LENGTH, ZONMAT1A.741
& P_ROWS) ZONMAT1A.742
CALL ZONM(PPTN,Z_L_PPTN,LMASK,S_PMASS,LLPTS(1,1),ROW_LENGTH, ZONMAT1A.743
& P_ROWS) ZONMAT1A.744
CALL ZONM(SNOW,Z_L_SNOW,LMASK,S_PMASS,LLPTS(1,1),ROW_LENGTH, ZONMAT1A.745
& P_ROWS) ZONMAT1A.746
CALL ZONM(TAUX,Z_L_TAUX,LMASK,S_UMASS,LLPTS(1,2),ROW_LENGTH, ZONMAT1A.747
& U_ROWS) ZONMAT1A.748
CALL ZONM(TAUY,Z_L_TAUY,LMASK,S_UMASS,LLPTS(1,2),ROW_LENGTH, ZONMAT1A.749
& U_ROWS) ZONMAT1A.750
CL ZONMAT1A.751
CL 6.10 Mass weighted on p/u-grid for surface variables (sea points) ZONMAT1A.752
CL ZONMAT1A.753
CALL ZONM(TSTAR,Z_S_TSTAR,SMASK,S_PMASS,LSPTS(1,1),ROW_LENGTH, ZONMAT1A.754
& P_ROWS) ZONMAT1A.755
CALL ZONM(SH,Z_S_SH,SMASK,S_PMASS,LSPTS(1,1),ROW_LENGTH,P_ROWS) ZONMAT1A.756
CALL ZONM(EVAP,Z_S_EVAP,SMASK,S_PMASS,LSPTS(1,1),ROW_LENGTH, ZONMAT1A.757
& P_ROWS) ZONMAT1A.758
CALL ZONM(PPTN,Z_S_PPTN,SMASK,S_PMASS,LSPTS(1,1),ROW_LENGTH, ZONMAT1A.759
& P_ROWS) ZONMAT1A.760
CALL ZONM(SNOW,Z_S_SNOW,SMASK,S_PMASS,LSPTS(1,1),ROW_LENGTH, ZONMAT1A.761
& P_ROWS) ZONMAT1A.762
CALL ZONM(AICE,Z_S_AICE,SMASK,S_PMASS,LSPTS(1,1),ROW_LENGTH, ZONMAT1A.763
& P_ROWS) ZONMAT1A.764
CALL ZONM(HICE,Z_S_HICE,SMASK,S_PMASS,LSPTS(1,1),ROW_LENGTH, ZONMAT1A.765
& P_ROWS) ZONMAT1A.766
CALL ZONM(TAUX,Z_S_TAUX,SMASK,S_UMASS,LSPTS(1,2),ROW_LENGTH, ZONMAT1A.767
& U_ROWS) ZONMAT1A.768
CALL ZONM(TAUY,Z_S_TAUY,SMASK,S_UMASS,LSPTS(1,2),ROW_LENGTH, ZONMAT1A.769
& U_ROWS) ZONMAT1A.770
CL---------------------------------------------------------------------- ZONMAT1A.771
CL 7. Compute quarter global and global means ZONMAT1A.772
CL ZONMAT1A.773
DO 90,I=1,4 ZONMAT1A.774
CALL GLBM(PMSL,G_Q_PMSL(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.775
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.776
CALL GLBM(PSTAR,G_Q_PSTAR(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.777
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.778
CALL GLBM(TSTAR,G_Q_TSTAR(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.779
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.780
CALL GLBM(TSTAR,G_Q_L_TSTAR(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.781
& COS_P_LATITUDE,LMASK,LQGLPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.782
CALL GLBM(TSTAR,G_Q_S_TSTAR(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.783
& COS_P_LATITUDE,SMASK,LQGSPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.784
CALL GLBM(SOILM,G_Q_SOILM(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.785
& COS_P_LATITUDE,LMASK,LQGLPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.786
CALL GLBM(SNOWD,G_Q_SNOWD(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.787
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.788
CALL GLBM(CANOPYW,G_Q_CANOPYW(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.789
& COS_P_LATITUDE,LMASK,LQGLPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.790
CALL GLBM(SH,G_Q_SH(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.791
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.792
CALL GLBM(SH,G_Q_L_SH(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.793
& COS_P_LATITUDE,LMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.794
CALL GLBM(SH,G_Q_S_SH(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.795
& COS_P_LATITUDE,SMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.796
CALL GLBM(EVAP,G_Q_EVAP(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.797
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.798
CALL GLBM(EVAP,G_Q_L_EVAP(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.799
& COS_P_LATITUDE,LMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.800
CALL GLBM(EVAP,G_Q_S_EVAP(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.801
& COS_P_LATITUDE,SMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.802
CALL GLBM(PPTN,G_Q_PPTN(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.803
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.804
CALL GLBM(PPTN,G_Q_L_PPTN(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.805
& COS_P_LATITUDE,LMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.806
CALL GLBM(PPTN,G_Q_S_PPTN(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.807
& COS_P_LATITUDE,SMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.808
CALL GLBM(LSRN,G_Q_LSRN(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.809
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.810
CALL GLBM(CVRN,G_Q_CVRN(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.811
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.812
CALL GLBM(SNOW,G_Q_SNOW(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.813
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.814
CALL GLBM(SNOW,G_Q_L_SNOW(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.815
& COS_P_LATITUDE,LMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.816
CALL GLBM(SNOW,G_Q_S_SNOW(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.817
& COS_P_LATITUDE,SMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.818
CALL GLBM(TATMOS,G_Q_TATMOS(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.819
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.820
CALL GLBM(QATMOS,G_Q_QATMOS(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.821
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.822
CALL GLBM(UATMOS,G_Q_UATMOS(I),START_ROW(2,I),END_ROW(2,I), ZONMAT1A.823
& COS_U_LATITUDE,AMASK,LQGAPTS(I,1),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.824
CALL GLBM(VATMOS,G_Q_VATMOS(I),START_ROW(2,I),END_ROW(2,I), ZONMAT1A.825
& COS_U_LATITUDE,AMASK,LQGAPTS(I,1),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.826
CALL GLBM(TAUX,G_Q_TAUX(I),START_ROW(2,I),END_ROW(2,I), ZONMAT1A.827
& COS_U_LATITUDE,AMASK,LQGAPTS(I,2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.828
CALL GLBM(TAUY,G_Q_TAUY(I),START_ROW(2,I),END_ROW(2,I), ZONMAT1A.829
& COS_U_LATITUDE,AMASK,LQGAPTS(I,2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.830
CALL GLBM(TAUX,G_Q_L_TAUX(I),START_ROW(2,I),END_ROW(2,I), ZONMAT1A.831
& COS_U_LATITUDE,LMASK,LQGLPTS(I,2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.832
CALL GLBM(TAUY,G_Q_L_TAUY(I),START_ROW(2,I),END_ROW(2,I), ZONMAT1A.833
& COS_U_LATITUDE,LMASK,LQGLPTS(I,2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.834
CALL GLBM(TAUX,G_Q_S_TAUX(I),START_ROW(2,I),END_ROW(2,I), ZONMAT1A.835
& COS_U_LATITUDE,SMASK,LQGSPTS(I,2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.836
CALL GLBM(TAUY,G_Q_S_TAUY(I),START_ROW(2,I),END_ROW(2,I), ZONMAT1A.837
& COS_U_LATITUDE,SMASK,LQGSPTS(I,2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.838
90 CONTINUE ZONMAT1A.839
ZONMAT1A.840
DO LEVEL=1,ST_LEVELS AJS1F401.1542
DO I=1,4 ZONMAT1A.842
CALL GLBM(SOILT(1,LEVEL),G_Q_SOILT(I,LEVEL), START_ROW(1,I), ZONMAT1A.843
& END_ROW(1,I),COS_P_LATITUDE,LMASK, LQGLPTS(I,1),P_ROWS, ZONMAT1A.844
& ROW_LENGTH,.TRUE.) ZONMAT1A.845
END DO ZONMAT1A.846
END DO ZONMAT1A.847
ZONMAT1A.848
CL 7.1 Winds and Kinetic energy at each level, mass weighted ZONMAT1A.849
ZONMAT1A.850
DO 120 LEVEL=1,P_LEVELS ZONMAT1A.851
ZONMAT1A.852
DO I=1,U_FIELD ZONMAT1A.853
WORK3(I,LEVEL)= ZONMAT1A.854
& U(I,LEVEL)*U(I,LEVEL)+V(I,LEVEL)*V(I,LEVEL) ZONMAT1A.855
WORK2(I)=Z_ZKE_LEV((I-1)/ROW_LENGTH+1,LEVEL) ZONMAT1A.856
END DO ZONMAT1A.857
ZONMAT1A.858
CALL GLBM(WORK3(1,LEVEL),G_TKE_LEV(LEVEL), ZONMAT1A.859
& 1,U_ROWS,COS_U_LATITUDE, ZONMAT1A.860
& U_MASS(1,LEVEL),LGAPTS(2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.861
CALL GLBM(WORK2,G_ZKE_LEV(LEVEL),1,U_ROWS,COS_U_LATITUDE, ZONMAT1A.862
& U_MASS(1,LEVEL),LGAPTS(2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.863
G_EKE_LEV(LEVEL)=G_TKE_LEV(LEVEL)-G_ZKE_LEV(LEVEL) ZONMAT1A.864
ZONMAT1A.865
DO 130,I=1,4 ZONMAT1A.866
CALL GLBM(U(1,LEVEL),G_Q_U(I,LEVEL),START_ROW(2,I), ZONMAT1A.867
& END_ROW(2,I),COS_U_LATITUDE,U_MASS(1,LEVEL),LQGAPTS(I,2), ZONMAT1A.868
& U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.869
CALL GLBM(V(1,LEVEL),G_Q_V(I,LEVEL),START_ROW(2,I), ZONMAT1A.870
& END_ROW(2,I),COS_U_LATITUDE,U_MASS(1,LEVEL),LQGAPTS(I,2), ZONMAT1A.871
& U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.872
130 CONTINUE ZONMAT1A.873
120 CONTINUE ZONMAT1A.874
ZONMAT1A.875
CL 7.2 Column mean energies, mass weighted ZONMAT1A.876
ZONMAT1A.877
CALL COLM(WORK3,WORK2,U_MASS,U_ROWS,ROW_LENGTH,P_LEVELS) ZONMAT1A.878
DO I=1,4 ZONMAT1A.879
CALL GLBM(WORK2,G_Q_TKE(I),START_ROW(2,I),END_ROW(2,I), ZONMAT1A.880
& COS_U_LATITUDE,PSTAR_U,LQGAPTS(I,2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.881
END DO ZONMAT1A.882
CALL GLBM(WORK2,G_TKE,1,U_ROWS, ZONMAT1A.883
& COS_U_LATITUDE,PSTAR_U,LGAPTS(2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.884
DO I=1,U_FIELD ZONMAT1A.885
WORK2(I)=Z_ZKE((I-1)/ROW_LENGTH+1) ZONMAT1A.886
END DO ZONMAT1A.887
DO I=1,4 ZONMAT1A.888
CALL GLBM(WORK2,G_Q_ZKE(I),START_ROW(2,I),END_ROW(2,I), ZONMAT1A.889
& COS_U_LATITUDE,PSTAR_U,LQGAPTS(I,2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.890
G_Q_EKE(I)=G_Q_TKE(I)-G_Q_ZKE(I) ZONMAT1A.891
END DO ZONMAT1A.892
CALL GLBM(WORK2,G_ZKE,1,U_ROWS, ZONMAT1A.893
& COS_U_LATITUDE,PSTAR_U,LGAPTS(2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.894
G_EKE=G_TKE-G_ZKE ZONMAT1A.895
ZONMAT1A.896
CL 7.3 Temperature and temperature variance at each level, mass weighted ZONMAT1A.897
ZONMAT1A.898
DO 140 LEVEL=1,P_LEVELS ZONMAT1A.899
ZONMAT1A.900
DO I=1,P_FIELD ZONMAT1A.901
WORK1(I,LEVEL)=T(I,LEVEL)*T(I,LEVEL) ZONMAT1A.902
END DO ZONMAT1A.903
ZONMAT1A.904
DO I=1,4 ZONMAT1A.905
CALL GLBM(T(1,LEVEL),G_Q_T(I,LEVEL), ZONMAT1A.906
& START_ROW(1,I),END_ROW(1,I),COS_P_LATITUDE, ZONMAT1A.907
& P_MASS(1,LEVEL),LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.908
ENDDO ZONMAT1A.909
ZONMAT1A.910
CALL GLBM(T(1,LEVEL),G_T(LEVEL),1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.911
& P_MASS(1,LEVEL),LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.912
CALL GLBM(WORK1(1,LEVEL),G_VAR_T_LEV(LEVEL),1,P_ROWS, ZONMAT1A.913
& COS_P_LATITUDE,P_MASS(1,LEVEL),LGAPTS(1),P_ROWS,ROW_LENGTH, ZONMAT1A.914
& .TRUE.) ZONMAT1A.915
G_VAR_T_LEV(LEVEL)=G_VAR_T_LEV(LEVEL)-G_T(LEVEL)*G_T(LEVEL) ZONMAT1A.916
ZONMAT1A.917
140 CONTINUE ZONMAT1A.918
ZONMAT1A.919
CL 7.4 Column mean temperature and temperature variance, mass weighted ZONMAT1A.920
ZONMAT1A.921
CALL COLM(WORK1,WORK2,P_MASS,P_ROWS,ROW_LENGTH,P_LEVELS) ZONMAT1A.922
CALL GLBM(TATMOS,G_TATMOS,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.923
& PSTAR,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.924
CALL GLBM(WORK2,G_VAR_T,1,P_ROWS, ZONMAT1A.925
& COS_P_LATITUDE,PSTAR,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.926
G_VAR_T=G_VAR_T-G_TATMOS*G_TATMOS ZONMAT1A.927
ZONMAT1A.928
CL 7.5 Moisture and moisture variance at each level, mass weighted ZONMAT1A.929
ZONMAT1A.930
DO LEVEL=1,Q_LEVELS ZONMAT1A.931
ZONMAT1A.932
DO I=1,P_FIELD ZONMAT1A.933
WORK1(I,LEVEL)=Q(I,LEVEL)*Q(I,LEVEL) ZONMAT1A.934
END DO ZONMAT1A.935
ZONMAT1A.936
DO I=1,4 ZONMAT1A.937
CALL GLBM(Q(1,LEVEL),G_Q_Q(I,LEVEL), ZONMAT1A.938
& START_ROW(1,I),END_ROW(1,I),COS_P_LATITUDE, ZONMAT1A.939
& P_MASS(1,LEVEL),LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.940
ENDDO ZONMAT1A.941
ZONMAT1A.942
CALL GLBM(Q(1,LEVEL),G_Q(LEVEL),1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.943
& P_MASS(1,LEVEL),LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.944
CALL GLBM(WORK1(1,LEVEL),G_VAR_Q_LEV(LEVEL),1,P_ROWS, ZONMAT1A.945
& COS_P_LATITUDE,P_MASS(1,LEVEL),LGAPTS(1),P_ROWS,ROW_LENGTH, ZONMAT1A.946
& .TRUE.) ZONMAT1A.947
G_VAR_Q_LEV(LEVEL)=G_VAR_Q_LEV(LEVEL)-G_Q(LEVEL)*G_Q(LEVEL) ZONMAT1A.948
ZONMAT1A.949
END DO ZONMAT1A.950
ZONMAT1A.951
CL 7.6 Column mean moisture and moisture variance, mass weighted ZONMAT1A.952
ZONMAT1A.953
CALL COLM(WORK1,WORK2,P_MASS,P_ROWS,ROW_LENGTH,Q_LEVELS) ZONMAT1A.954
CALL GLBM(QATMOS,G_QATMOS,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.955
& PSTAR,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.956
CALL GLBM(WORK2,G_VAR_Q,1,P_ROWS, ZONMAT1A.957
& COS_P_LATITUDE,PSTAR,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.958
G_VAR_Q=G_VAR_Q-G_QATMOS*G_QATMOS ZONMAT1A.959
CL ZONMAT1A.960
CL 7.7 Cloud water liquid and ice over whole atmosphere, put in rows ZONMAT1A.961
CL ZONMAT1A.962
IF ((IPRTEXTRA.EQ.1) .AND. (ICLOUD.GT.1)) THEN ARS2F304.153
CL ZONMAT1A.964
CALL COLM(CLLIQ,CLLIQ_F,P_MASS,P_ROWS,ROW_LENGTH,Q_LEVELS) ZONMAT1A.965
CALL ZONM(CLLIQ_F,Z_CLLIQ,AMASK,PSTAR,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.966
& P_ROWS) ZONMAT1A.967
CALL COLM(CLICE,CLICE_F,P_MASS,P_ROWS,ROW_LENGTH,Q_LEVELS) ZONMAT1A.968
CALL ZONM(CLICE_F,Z_CLICE,AMASK,PSTAR,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.969
& P_ROWS) ZONMAT1A.970
ZONMAT1A.971
END IF ZONMAT1A.972
CL ZONMAT1A.973
CL---------------------------------------------------------------------- ZONMAT1A.974
CL 8. More global means ZONMAT1A.975
CL ZONMAT1A.976
CALL GLBM(PMSL,G_PMSL,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.977
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.978
CALL GLBM(PSTAR,G_PSTAR,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.979
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.980
CALL GLBM(TSTAR,G_TSTAR,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.981
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.982
CALL GLBM(TSTAR,G_L_TSTAR,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.983
& LMASK,LGLPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.984
CALL GLBM(TSTAR,G_S_TSTAR,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.985
& SMASK,LGSPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.986
CALL GLBM(SOILM,G_SOILM,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.987
& LMASK,LGLPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.988
CALL GLBM(SNOWD,G_SNOWD,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.989
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.990
CALL GLBM(CANOPYW,G_CANOPYW,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.991
& LMASK,LGLPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.992
CALL GLBM(SH,G_SH,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.993
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.994
CALL GLBM(SH,G_L_SH,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.995
& LMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.996
CALL GLBM(SH,G_S_SH,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.997
& SMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.998
CALL GLBM(EVAP,G_EVAP,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.999
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1000
CALL GLBM(EVAP,G_L_EVAP,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1001
& LMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1002
CALL GLBM(EVAP,G_S_EVAP,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1003
& SMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1004
CALL GLBM(PPTN,G_PPTN,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1005
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1006
CALL GLBM(PPTN,G_L_PPTN,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1007
& LMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1008
CALL GLBM(PPTN,G_S_PPTN,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1009
& SMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1010
CALL GLBM(LSRN,G_LSRN,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1011
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1012
CALL GLBM(CVRN,G_CVRN,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1013
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1014
CALL GLBM(SNOW,G_SNOW,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1015
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1016
CALL GLBM(SNOW,G_L_SNOW,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1017
& LMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1018
CALL GLBM(SNOW,G_S_SNOW,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1019
& SMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1020
CALL GLBM(UATMOS,G_UATMOS,1,U_ROWS,COS_U_LATITUDE, ZONMAT1A.1021
& PSTAR_U,LGAPTS(2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.1022
CALL GLBM(VATMOS,G_VATMOS,1,U_ROWS,COS_U_LATITUDE, ZONMAT1A.1023
& PSTAR_U,LGAPTS(2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.1024
CALL GLBM(TAUX,G_TAUX,1,U_ROWS,COS_U_LATITUDE, ZONMAT1A.1025
& AMASK,LGAPTS(2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.1026
CALL GLBM(TAUY,G_TAUY,1,U_ROWS,COS_U_LATITUDE, ZONMAT1A.1027
& AMASK,LGAPTS(2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.1028
CALL GLBM(TAUX,G_L_TAUX,1,U_ROWS,COS_U_LATITUDE, ZONMAT1A.1029
& LMASK,LGLPTS(2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.1030
CALL GLBM(TAUY,G_L_TAUY,1,U_ROWS,COS_U_LATITUDE, ZONMAT1A.1031
& LMASK,LGLPTS(2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.1032
CALL GLBM(TAUX,G_S_TAUX,1,U_ROWS,COS_U_LATITUDE, ZONMAT1A.1033
& SMASK,LGSPTS(2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.1034
CALL GLBM(TAUY,G_S_TAUY,1,U_ROWS,COS_U_LATITUDE, ZONMAT1A.1035
& SMASK,LGSPTS(2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.1036
ZONMAT1A.1037
DO LEVEL=1,P_LEVELS ZONMAT1A.1038
CALL GLBM(U(1,LEVEL),G_U(LEVEL),1,U_ROWS,COS_U_LATITUDE, ZONMAT1A.1039
& U_MASS(1,LEVEL),LGAPTS(2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.1040
CALL GLBM(V(1,LEVEL),G_V(LEVEL),1,U_ROWS,COS_U_LATITUDE, ZONMAT1A.1041
& U_MASS(1,LEVEL),LGAPTS(2),U_ROWS,ROW_LENGTH,.FALSE.) ZONMAT1A.1042
END DO ZONMAT1A.1043
C ZONMAT1A.1044
DO LEVEL=1,ST_LEVELS AJS1F401.1543
CALL GLBM(SOILT(1,LEVEL),G_SOILT(LEVEL),1,P_ROWS, ZONMAT1A.1046
& COS_P_LATITUDE,LMASK,LGLPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1047
END DO ZONMAT1A.1048
C ZONMAT1A.1049
CALL ZONM(SUBL,Z_L_SUBL,LMASK,S_PMASS,LLPTS(1,1),ROW_LENGTH, ZONMAT1A.1050
& P_ROWS) ZONMAT1A.1051
CALL ZONM(SUBL,Z_S_SUBL,SMASK,S_PMASS,LSPTS(1,1),ROW_LENGTH, ZONMAT1A.1052
& P_ROWS) ZONMAT1A.1053
DO I=1,4 ZONMAT1A.1054
CALL GLBM(SUBL,G_Q_L_SUBL(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1055
& COS_P_LATITUDE,LMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1056
CALL GLBM(SUBL,G_Q_SUBL(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1057
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1058
CALL GLBM(SUBL,G_Q_S_SUBL(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1059
& COS_P_LATITUDE,SMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1060
ENDDO ZONMAT1A.1061
CALL GLBM(SUBL,G_SUBL,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1062
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1063
CALL GLBM(SUBL,G_L_SUBL,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1064
& LMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1065
CALL GLBM(SUBL,G_S_SUBL,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1066
& SMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1067
DO I=1,4 ZONMAT1A.1068
CALL GLBM(SOEV,G_Q_L_SOEV(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1069
& COS_P_LATITUDE,LMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1070
ENDDO ZONMAT1A.1071
CALL GLBM(SOEV,G_L_SOEV,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1072
& LMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1073
C ZONMAT1A.1074
CL---------------------------------------------------------------------- ZONMAT1A.1075
CL ZONMAT1A.1076
CL 8.a Additional zonal means ONLY calculated if field available ZONMAT1A.1077
CL ZONMAT1A.1078
IF (IPRTEXTRA.EQ.1) THEN ZONMAT1A.1079
IF (IHYDRO.GT.1) THEN ARS2F304.154
CALL ZONM(SFRU,Z_L_SFRU,LMASK,S_PMASS,LLPTS(1,1),ROW_LENGTH, ZONMAT1A.1081
& P_ROWS) ZONMAT1A.1082
CALL ZONM(SBRU,Z_L_SBRU,LMASK,S_PMASS,LLPTS(1,1),ROW_LENGTH, ZONMAT1A.1083
& P_ROWS) ZONMAT1A.1084
DO I=1,4 ZONMAT1A.1085
CALL GLBM(THRF,G_Q_L_THRF(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1086
& COS_P_LATITUDE,LMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1087
CALL GLBM(SNML,G_Q_L_SNML(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1088
& COS_P_LATITUDE,LMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1089
CALL GLBM(SFRU,G_Q_L_SFRU(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1090
& COS_P_LATITUDE,LMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1091
CALL GLBM(SBRU,G_Q_L_SBRU(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1092
& COS_P_LATITUDE,LMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1093
ENDDO ZONMAT1A.1094
CALL GLBM(SFRU,G_L_SFRU,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1095
& LMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1096
CALL GLBM(SBRU,G_L_SBRU,1,P_ROWS,COS_P_LATITUDE, ARS2F304.155
& LMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1098
CALL GLBM(THRF,G_L_THRF,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1099
& LMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1100
CALL GLBM(SNML,G_L_SNML,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1101
& LMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1102
ENDIF ZONMAT1A.1103
IF (IRAD.GT.1) THEN ! RADIATION ARS2F304.156
CALL ZONM(SDTR,Z_L_SDTR,LMASK,S_PMASS,LLPTS(1,1),ROW_LENGTH, ZONMAT1A.1105
& P_ROWS) ZONMAT1A.1106
CALL ZONM(SDTR,Z_S_SDTR,SMASK,S_PMASS,LSPTS(1,1),ROW_LENGTH, ZONMAT1A.1107
& P_ROWS) ZONMAT1A.1108
CALL ZONM(SDSR,Z_L_SDSR,LMASK,S_PMASS,LLPTS(1,1),ROW_LENGTH, ZONMAT1A.1109
& P_ROWS) ZONMAT1A.1110
CALL ZONM(SDSR,Z_S_SDSR,SMASK,S_PMASS,LSPTS(1,1),ROW_LENGTH, ZONMAT1A.1111
& P_ROWS) ZONMAT1A.1112
CALL ZONM(TDTR,Z_TDTR,AMASK,S_PMASS,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.1113
& P_ROWS) ZONMAT1A.1114
CALL ZONM(TOLR,Z_TOLR,AMASK,S_PMASS,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.1115
& P_ROWS) ZONMAT1A.1116
CALL ZONM(TOSW,Z_TOSW,AMASK,S_PMASS,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.1117
& P_ROWS) ZONMAT1A.1118
CALL ZONM(TISW,Z_TISW,AMASK,S_PMASS,LAPTS(1,1),ROW_LENGTH, ZONMAT1A.1119
& P_ROWS) ZONMAT1A.1120
DO I=1,P_ROWS ZONMAT1A.1121
IF (Z_TISW(I) .GT. 10.0E-10) THEN ZONMAT1A.1122
Z_ALBEDO(I)=Z_TOSW(I)/Z_TISW(I) ZONMAT1A.1123
ELSE ZONMAT1A.1124
Z_ALBEDO(I) = 0.0 ZONMAT1A.1125
END IF ZONMAT1A.1126
END DO ZONMAT1A.1127
DO I=1,4 ZONMAT1A.1128
CALL GLBM(SDTR,G_Q_L_SDTR(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1129
& COS_P_LATITUDE,LMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1130
CALL GLBM(SDTR,G_Q_SDTR(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1131
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1132
CALL GLBM(SDTR,G_Q_L_SDTR(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1133
& COS_P_LATITUDE,LMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1134
CALL GLBM(SDTR,G_Q_S_SDTR(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1135
& COS_P_LATITUDE,SMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1136
CALL GLBM(SDSR,G_Q_SDSR(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1137
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1138
CALL GLBM(SDSR,G_Q_L_SDSR(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1139
& COS_P_LATITUDE,LMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1140
CALL GLBM(SDSR,G_Q_S_SDSR(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1141
& COS_P_LATITUDE,SMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1142
CALL GLBM(TDTR,G_Q_TDTR(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1143
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1144
CALL GLBM(TOLR,G_Q_TOLR(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1145
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1146
CALL GLBM(TOSW,G_Q_TOSW(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1147
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1148
CALL GLBM(TISW,G_Q_TISW(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1149
& COS_P_LATITUDE,AMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1150
G_Q_ALBEDO(I)=G_Q_TOSW(I)/G_Q_TISW(I) ZONMAT1A.1151
ENDDO ZONMAT1A.1152
CALL GLBM(SDTR,G_L_SDTR,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1153
& LMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1154
CALL GLBM(SDTR,G_SDTR,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1155
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1156
CALL GLBM(SDTR,G_L_SDTR,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1157
& LMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1158
CALL GLBM(SDTR,G_S_SDTR,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1159
& SMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1160
CALL GLBM(SDSR,G_SDSR,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1161
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1162
CALL GLBM(SDSR,G_L_SDSR,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1163
& LMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1164
CALL GLBM(SDSR,G_S_SDSR,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1165
& SMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1166
CALL GLBM(TDTR,G_TDTR,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1167
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1168
CALL GLBM(TOLR,G_TOLR,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1169
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1170
CALL GLBM(TOSW,G_TOSW,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1171
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1172
CALL GLBM(TISW,G_TISW,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1173
& AMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1174
G_ALBEDO=G_TOSW/G_TISW ZONMAT1A.1175
ENDIF ZONMAT1A.1176
CL ZONMAT1A.1177
CL Fluxes over land - except for energy flux ENFL ZONMAT1A.1178
CL ZONMAT1A.1179
IF (IFLUXL.GT.1) THEN ARS2F304.157
DO I=1,4 ZONMAT1A.1181
CALL GLBM(SOHF,G_Q_L_SOHF(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1182
& COS_P_LATITUDE,LMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1183
G_Q_L_WFSS(I)=G_Q_L_SNOW(I)-G_Q_L_SUBL(I)-G_Q_L_SNML(I) ZONMAT1A.1184
CALL GLBM(WFCA,G_Q_L_WFCA(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1185
& COS_P_LATITUDE,LMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1186
G_Q_L_WFSO(I)=G_Q_L_THRF(I)+G_Q_L_SNML(I)-G_Q_L_SOEV(I) ZONMAT1A.1187
& -G_Q_L_SFRU(I)-G_Q_L_SBRU(I) ZONMAT1A.1188
G_Q_WAFL(I)=G_Q_EVAP(I)*86400.-G_Q_PPTN(I) ZONMAT1A.1189
G_Q_L_ENFS(I)=G_Q_L_SDTR(I)-G_Q_L_SH(I) - LC*G_Q_L_EVAP(I) ZONMAT1A.1190
& -LF*(G_Q_L_SUBL(I)+G_Q_L_SNML(I))/86400. ZONMAT1A.1191
G_Q_ENFL(I)=G_Q_TDTR(I)-G_Q_SDTR(I)+G_Q_SH(I) ZONMAT1A.1192
& +(LC*G_Q_PPTN(I)+LF*G_Q_SNOW(I))/86400. ZONMAT1A.1193
ENDDO ZONMAT1A.1194
CALL GLBM(SOHF,G_L_SOHF,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1195
& LMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1196
G_L_WFSS=G_L_SNOW-G_L_SUBL-G_L_SNML ZONMAT1A.1197
CALL GLBM(WFCA,G_L_WFCA,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1198
& LMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1199
G_L_WFSO=G_L_THRF+G_L_SNML-G_L_SOEV-G_L_SFRU-G_L_SBRU ZONMAT1A.1200
G_WAFL=G_EVAP*86400.-G_PPTN ZONMAT1A.1201
G_L_ENFS=G_L_SDTR-G_L_SH - LC*G_L_EVAP ZONMAT1A.1202
& -LF*(G_L_SUBL+G_L_SNML)/86400. ZONMAT1A.1203
G_ENFL=G_TDTR-G_SDTR+G_SH+(LC*G_PPTN+LF*G_SNOW)/86400. ZONMAT1A.1204
ENDIF ZONMAT1A.1205
CL ZONMAT1A.1206
CL Cloud water quarter globe means ZONMAT1A.1207
CL ZONMAT1A.1208
IF (ICLOUD.GT.1) THEN ARS2F304.158
DO LEVEL=1,Q_LEVELS ZONMAT1A.1210
DO I=1,4 ZONMAT1A.1211
CALL GLBM(CLLIQ(1,LEVEL),G_Q_CLLIQL(I,LEVEL), ZONMAT1A.1212
& START_ROW(1,I),END_ROW(1,I),COS_P_LATITUDE, ZONMAT1A.1213
& P_MASS(1,LEVEL),LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1214
CALL GLBM(CLICE(1,LEVEL),G_Q_CLICEL(I,LEVEL), ZONMAT1A.1215
& START_ROW(1,I),END_ROW(1,I),COS_P_LATITUDE, ZONMAT1A.1216
& P_MASS(1,LEVEL),LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1217
ENDDO ZONMAT1A.1218
ZONMAT1A.1219
CALL GLBM(CLLIQ(1,LEVEL),G_CLLIQL(LEVEL),1,P_ROWS, ZONMAT1A.1220
& COS_P_LATITUDE,P_MASS(1,LEVEL),LGAPTS(1),P_ROWS,ROW_LENGTH, ZONMAT1A.1221
& .TRUE.) ZONMAT1A.1222
CALL GLBM(CLICE(1,LEVEL),G_CLICEL(LEVEL),1,P_ROWS, ZONMAT1A.1223
& COS_P_LATITUDE,P_MASS(1,LEVEL),LGAPTS(1),P_ROWS,ROW_LENGTH, ZONMAT1A.1224
& .TRUE.) ZONMAT1A.1225
ZONMAT1A.1226
END DO ZONMAT1A.1227
CALL COLM(CLLIQ,WORK2,P_MASS,P_ROWS,ROW_LENGTH,Q_LEVELS) ZONMAT1A.1228
DO I=1,4 ZONMAT1A.1229
CALL GLBM(WORK2,G_Q_CLLIQ(I), ZONMAT1A.1230
& START_ROW(1,I),END_ROW(1,I),COS_P_LATITUDE, ZONMAT1A.1231
& PSTAR,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1232
ENDDO ZONMAT1A.1233
CALL GLBM(WORK2,G_CLLIQ,1,P_ROWS, ZONMAT1A.1234
& COS_P_LATITUDE,PSTAR,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1235
CALL COLM(CLICE,WORK2,P_MASS,P_ROWS,ROW_LENGTH,Q_LEVELS) ZONMAT1A.1236
DO I=1,4 ZONMAT1A.1237
CALL GLBM(WORK2,G_Q_CLICE(I), ZONMAT1A.1238
& START_ROW(1,I),END_ROW(1,I),COS_P_LATITUDE, ZONMAT1A.1239
& PSTAR,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1240
ENDDO ZONMAT1A.1241
CALL GLBM(WORK2,G_CLICE,1,P_ROWS, ZONMAT1A.1242
& COS_P_LATITUDE,PSTAR,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1243
ENDIF ZONMAT1A.1244
CL ZONMAT1A.1245
CL Fluxes over Sea-ice ZONMAT1A.1246
CL ZONMAT1A.1247
IF (ISICE.GT.1) THEN ARS2F304.159
DO I=1,4 ZONMAT1A.1249
CALL GLBM(SIHF,G_Q_SIHF(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1250
& COS_P_LATITUDE,SMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1251
CALL GLBM(SIMH,G_Q_SIMH(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1252
& COS_P_LATITUDE,SMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1253
CALL GLBM(SISS,G_Q_SISS(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1254
& COS_P_LATITUDE,SMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1255
CALL GLBM(SIST,G_Q_SIST(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1256
& COS_P_LATITUDE,SMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1257
CALL GLBM(SIEF,G_Q_SIEF(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1258
& COS_P_LATITUDE,SMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1259
CALL GLBM(SISH,G_Q_SISH(I),START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1260
& COS_P_LATITUDE,SMASK,LQGAPTS(I,1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1261
ENDDO ZONMAT1A.1262
CALL GLBM(SIHF,G_SIHF,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1263
& SMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1264
CALL GLBM(SIMH,G_SIMH,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1265
& SMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1266
CALL GLBM(SISS,G_SISS,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1267
& SMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1268
CALL GLBM(SIST,G_SIST,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1269
& SMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1270
CALL GLBM(SIEF,G_SIEF,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1271
& SMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1272
CALL GLBM(SISH,G_SISH,1,P_ROWS,COS_P_LATITUDE, ZONMAT1A.1273
& SMASK,LGAPTS(1),P_ROWS,ROW_LENGTH,.TRUE.) ZONMAT1A.1274
ENDIF ZONMAT1A.1275
ENDIF ZONMAT1A.1276
ZONMAT1A.1277
CL---------------------------------------------------------------------- ZONMAT1A.1278
CL 9. Print results neatly formatted ZONMAT1A.1279
CL ZONMAT1A.1280
WRITE(6,'(1X)') ZONMAT1A.1281
WRITE(6,'(A80/)') BANNER ZONMAT1A.1282
WRITE(6,115) ZONMAT1A.1283
115 FORMAT(T2,'ROW SEA LAND PMSL PSTAR SNOWD CANOPYW LSRN', ZONMAT1A.1284
&T54,'CVRN TATMOS QATMOS SOILM T*SEA T*LAND SOILT1 SOILT2 ', AJS1F401.1544
&T108,'SOILT3 SOILT4',T122,'SRUN SUBRUN') ZONMAT1A.1286
WRITE(6,125) ZONMAT1A.1287
125 FORMAT(T10,'ONLY',T16,'MB',T24,'MB',T34,'MM',T42,'MM',T50,'MM/DAY' ZONMAT1A.1288
&,T59,'DEG.K',T66,'GM/KG',T74,'MM',T80,'DEG.K',T87,'DEG.K',T94, ZONMAT1A.1289
&'DEG.K',T101,'DEG.K',T108,'DEG.K',T115,'DEG.K',T123,'MM/DAY'/) ZONMAT1A.1290
ZONMAT1A.1291
DO ROW=1,P_ROWS ZONMAT1A.1292
IF (IHYDRO.GT.1) THEN ARS2F304.160
WRITE(6,135)ROW,SPTS(ROW,1),LPTS(ROW,1),(Z_PMSL(ROW)*0.01), ZONMAT1A.1294
&(Z_PSTAR(ROW)*0.01),Z_SNOWD(ROW),Z_CANOPYW(ROW),Z_LSRN(ROW), ZONMAT1A.1295
&Z_CVRN(ROW),Z_TATMOS(ROW),(Z_QATMOS(ROW)*1000.),Z_SOILM(ROW), ZONMAT1A.1296
&Z_S_TSTAR(ROW),Z_L_TSTAR(ROW),(Z_SOILT(ROW,LEVEL), AJS1F401.1545
& LEVEL=1,ST_LEVELS),Z_L_SFRU(ROW),Z_L_SBRU(ROW) AJS1F401.1546
ELSE ZONMAT1A.1299
WRITE(6,135)ROW,SPTS(ROW,1),LPTS(ROW,1),(Z_PMSL(ROW)*0.01), ZONMAT1A.1300
&(Z_PSTAR(ROW)*0.01),Z_SNOWD(ROW),Z_CANOPYW(ROW),Z_LSRN(ROW), ZONMAT1A.1301
&Z_CVRN(ROW),Z_TATMOS(ROW),(Z_QATMOS(ROW)*1000.),Z_SOILM(ROW), ZONMAT1A.1302
&Z_S_TSTAR(ROW),Z_L_TSTAR(ROW),(Z_SOILT(ROW,LEVEL), AJS1F401.1547
& LEVEL=1,ST_LEVELS) AJS1F401.1548
END IF ZONMAT1A.1305
IF (MOD(ROW,5).EQ.0) WRITE(6,'(1X)') ZONMAT1A.1306
END DO ZONMAT1A.1307
ZONMAT1A.1308
135 FORMAT(T2,I3,T6,I3,T10,I3,T15,F7.2,T23,F7.2,T31,F9.2,T41,F5.2,T47, ZONMAT1A.1309
&F5.3,T53,F5.3,T59,F6.2,T66,F5.3,T73,F6.2,T80,F6.2,T87,F6.2,T94, ZONMAT1A.1310
&F6.2,T101,F6.2,T108,F6.2,T115,F6.2,T122,2F5.2) ZONMAT1A.1311
WRITE(6,'(1X)') ZONMAT1A.1312
WRITE(6,'(A80/)') BANNER ZONMAT1A.1313
WRITE(6,215) ZONMAT1A.1314
215 FORMAT(T2,'ROW SEA LAND',T16,'S.HEAT (W/M2)',T34,'EVAP (MM/DAY)' ZONMAT1A.1315
&,T52,'PPTN (MM/DAY)',T70,'SNOW (MM/DAY)',T87,'AICE' ZONMAT1A.1316
&,T93,'HICE(M)',T101,'SUBLM(mm/day)',T115,'Cloud water whole') ZONMAT1A.1317
WRITE(6,220) ZONMAT1A.1318
220 FORMAT(T115,'atmos (g/kg*1000)') ZONMAT1A.1319
WRITE(6,225) ZONMAT1A.1320
225 FORMAT(T16,'ALL',T22,'SEA',T27,'LAND',T34,'ALL',T40,'SEA',T45, ZONMAT1A.1321
&'LAND',T52,'ALL',T57,'SEA',T63,'LAND',T70,'ALL',T76,'SEA',T81, ZONMAT1A.1322
&'LAND',T88,'SEA',T95,'SEA',T102,'SEA',T108,'LAND',T115,'LIQUID', ZONMAT1A.1323
& T123,'ICE'/) ZONMAT1A.1324
ZONMAT1A.1325
DO ROW=1,P_ROWS ZONMAT1A.1326
IF (ICLOUD.GT.1) THEN ARS2F304.161
WRITE(6,235)ROW,SPTS(ROW,1),LPTS(ROW,1),Z_SH(ROW),Z_S_SH(ROW), ZONMAT1A.1328
&Z_L_SH(ROW),(Z_EVAP(ROW)*86400.),(Z_S_EVAP(ROW)*86400.), ZONMAT1A.1329
&(Z_L_EVAP(ROW)*86400.),Z_PPTN(ROW), ZONMAT1A.1330
&Z_S_PPTN(ROW),Z_L_PPTN(ROW),Z_SNOW(ROW),Z_S_SNOW(ROW), ZONMAT1A.1331
&Z_L_SNOW(ROW),Z_S_AICE(ROW),Z_S_HICE(ROW),Z_S_SUBL(ROW), ZONMAT1A.1332
&Z_L_SUBL(ROW),1000000.*Z_CLLIQ(ROW),1000000.*Z_CLICE(ROW) ZONMAT1A.1333
ELSE ZONMAT1A.1334
WRITE(6,235)ROW,SPTS(ROW,1),LPTS(ROW,1),Z_SH(ROW),Z_S_SH(ROW), ZONMAT1A.1335
&Z_L_SH(ROW),(Z_EVAP(ROW)*86400.),(Z_S_EVAP(ROW)*86400.), ZONMAT1A.1336
&(Z_L_EVAP(ROW)*86400.),Z_PPTN(ROW), ZONMAT1A.1337
&Z_S_PPTN(ROW),Z_L_PPTN(ROW),Z_SNOW(ROW),Z_S_SNOW(ROW), ZONMAT1A.1338
&Z_L_SNOW(ROW),Z_S_AICE(ROW),Z_S_HICE(ROW),Z_S_SUBL(ROW), ZONMAT1A.1339
&Z_L_SUBL(ROW) ZONMAT1A.1340
END IF ZONMAT1A.1341
IF (MOD(ROW,5).EQ.0) WRITE(6,'(1X)') ZONMAT1A.1342
END DO ZONMAT1A.1343
235 FORMAT(T2,I3,T6,I3,T10,I3,T15,F5.1,T21,F5.1,T27,F5.1,T33,F5.3,T39, ZONMAT1A.1344
&F5.3,T45,F5.3,T51,F5.3,T57,F5.3,T63,F5.3,T69,F5.3,T75,F5.3,T81, ZONMAT1A.1345
&F5.3,T87,F5.3,T94,F5.3,T100,2(1X,F6.3),T115,2(F6.2,1X)) ZONMAT1A.1346
C ZONMAT1A.1347
IF ((IPRTWIND.EQ.1) .AND. (IRAD.EQ.1)) THEN ARS2F304.162
WRITE(6,'(1X)') ZONMAT1A.1349
WRITE(6,'(A80/)') BANNER ZONMAT1A.1350
WRITE(6,315) ZONMAT1A.1351
315 FORMAT(T2,'ROW SEA LAND UATMOS VATMOS',T35,'TAUX (N/M2)',T56, ZONMAT1A.1352
&'TAUY (N/M2)') ZONMAT1A.1353
WRITE(6,325) ZONMAT1A.1354
325 FORMAT(T18,'M/S',T25,'M/S',T31,'ALL',T38,'SEA',T45,'LAND',T52, ZONMAT1A.1355
&'ALL',T59,'SEA',T66,'LAND'/) ZONMAT1A.1356
ZONMAT1A.1357
DO ROW=1,U_ROWS ZONMAT1A.1358
WRITE(6,330)ROW,SPTS(ROW,2),LPTS(ROW,2),Z_UATMOS(ROW), ZONMAT1A.1359
&Z_VATMOS(ROW),Z_TAUX(ROW),Z_S_TAUX(ROW),Z_L_TAUX(ROW),Z_TAUY(ROW), ZONMAT1A.1360
&Z_S_TAUY(ROW),Z_L_TAUY(ROW) ZONMAT1A.1361
END DO ZONMAT1A.1362
330 FORMAT(T2,I3,T6,I3,T10,I3,T15,F6.2,T22,F6.2,T29,6F7.3) ZONMAT1A.1363
END IF ZONMAT1A.1364
C ZONMAT1A.1365
IF ((IPRTWIND.NE.1) .AND. (IRAD.GT.1)) THEN ARS2F304.163
WRITE(6,'(1X)') ZONMAT1A.1367
WRITE(6,'(A80/)') BANNER ZONMAT1A.1368
WRITE(6,331) ZONMAT1A.1369
331 FORMAT(T2,'ROW TOTAL SURFACE SURFACE SOLAR TOTAL TOA LW TOA', ZONMAT1A.1370
&' ALBEDO') ZONMAT1A.1371
WRITE(6,332) ZONMAT1A.1372
332 FORMAT(T7,2('Land',4X,'Sea',4X),2X,2('rad',7X)) ZONMAT1A.1373
ZONMAT1A.1374
DO ROW=1,U_ROWS ZONMAT1A.1375
WRITE(6,333)ROW,Z_L_SDTR(ROW),Z_S_SDTR(ROW),Z_L_SDSR(ROW), ZONMAT1A.1376
&Z_S_SDSR(ROW),Z_TDTR(ROW),Z_TOLR(ROW),Z_ALBEDO(ROW) AAS1F403.2
END DO ZONMAT1A.1378
333 FORMAT(T2,I3,1X,2(2(F6.1,1X),1X),1X,2(F7.2,2X),1X,F7.4) PXZONMAT.1
END IF ZONMAT1A.1380
ZONMAT1A.1381
IF ((IPRTWIND.EQ.1) .AND. (IRAD.GT.1)) THEN ARS2F304.164
WRITE(6,'(1X)') ZONMAT1A.1383
WRITE(6,'(A80/)') BANNER ZONMAT1A.1384
WRITE(6,335) ZONMAT1A.1385
335 FORMAT(T2,'ROW SEA LAND UATMOS VATMOS',T35,'TAUX (N/M2)',T56, ZONMAT1A.1386
&'TAUY (N/M2)',T72,'TOTAL SURFACE SURFACE SOLAR TOTAL TOA', ZONMAT1A.1387
&' LW TOA ALBEDO') ZONMAT1A.1388
WRITE(6,340) ZONMAT1A.1389
340 FORMAT(T18,'M/S',T25,'M/S',T31,'ALL',T38,'SEA',T45,'LAND',T52, ZONMAT1A.1390
&'ALL',T59,'SEA',T66,'LAND',T73,2('Land',4X,'Sea',4X), ZONMAT1A.1391
& 2X,2('rad',7X)) ZONMAT1A.1392
ZONMAT1A.1393
DO ROW=1,U_ROWS ZONMAT1A.1394
WRITE(6,345)ROW,SPTS(ROW,2),LPTS(ROW,2),Z_UATMOS(ROW), ZONMAT1A.1395
&Z_VATMOS(ROW),Z_TAUX(ROW),Z_S_TAUX(ROW),Z_L_TAUX(ROW),Z_TAUY(ROW), ZONMAT1A.1396
&Z_S_TAUY(ROW),Z_L_TAUY(ROW), ZONMAT1A.1397
&Z_L_SDTR(ROW),Z_S_SDTR(ROW),Z_L_SDSR(ROW), ZONMAT1A.1398
&Z_S_SDSR(ROW),Z_TDTR(ROW),Z_TOLR(ROW),Z_ALBEDO(ROW) AAS1F403.3
AAS1F403.4
AAS1F403.5
END DO ZONMAT1A.1400
345 FORMAT(T2,I3,T6,I3,T10,I3,T15,F6.2,T22,F6.2,T29,6F7.3, ZONMAT1A.1401
& T72,2(2(F6.1,1X),1X),1X,2(F7.2,2X),1X,F7.4) PXZONMAT.2
END IF ZONMAT1A.1403
C ZONMAT1A.1404
IF (IPRTWIND.EQ.1) THEN ZONMAT1A.1405
WRITE(6,'(1X)') ZONMAT1A.1406
WRITE(6,'(A80/)') BANNER ZONMAT1A.1407
WRITE(6,415) ZONMAT1A.1408
415 FORMAT(T2,'ROW U1',T14,'U2',T20,'U3',T26,'U4',T32,'U5',T38,'U6', ZONMAT1A.1409
&T44,'U7',T50,'U8',T56,'U9',T62,'U10',T68,'U11',T74,'U12',T80, ZONMAT1A.1410
&'U13',T86,'U14',T92,'U15',T98,'U16',T104,'U17',T110,'U18',T116, ZONMAT1A.1411
&'U19',T122,'U20') ZONMAT1A.1412
WRITE(6,425) ZONMAT1A.1413
425 FORMAT(T7,'M/S',T13,'M/S',T19,'M/S',T25,'M/S',T31,'M/S', ZONMAT1A.1414
&T37,'M/S',T43,'M/S',T49,'M/S',T55,'M/S',T61,'M/S',T67, ZONMAT1A.1415
&'M/S',T73,'M/S',T79,'M/S',T85,'M/S',T91,'M/S',T97, ZONMAT1A.1416
&'M/S',T103,'M/S',T109,'M/S',T115,'M/S',T121,'M/S'/) ZONMAT1A.1417
DO ROW=1,U_ROWS ZONMAT1A.1418
WRITE(6,435)ROW,(Z_U(ROW,LEVEL), LEVEL=1,P_LEVELS) ZONMAT1A.1419
IF (MOD(ROW,5).EQ.0) WRITE(6,'(1X)') ZONMAT1A.1420
END DO ZONMAT1A.1421
435 FORMAT(T1,I3,T5,F5.2,19F6.2) ZONMAT1A.1422
c ZONMAT1A.1423
WRITE(6,'(1X)') ZONMAT1A.1424
WRITE(6,'(A80/)') BANNER ZONMAT1A.1425
WRITE(6,515) ZONMAT1A.1426
515 FORMAT(T2,'ROW V1',T14,'V2',T20,'V3',T26,'V4',T32,'V5',T38,'V6', ZONMAT1A.1427
&T44,'V7',T50,'V8',T56,'V9',T62,'V10',T68,'V11',T74,'V12',T80, ZONMAT1A.1428
&'V13',T86,'V14',T92,'V15',T98,'V16',T104,'V17',T110,'V18',T116, ZONMAT1A.1429
&'V19',T122,'V20') ZONMAT1A.1430
WRITE(6,525) ZONMAT1A.1431
525 FORMAT(T7,'M/S',T13,'M/S',T19,'M/S',T25,'M/S',T31,'M/S', ZONMAT1A.1432
&T37,'M/S',T43,'M/S',T49,'M/S',T55,'M/S',T61,'M/S',T67, ZONMAT1A.1433
&'M/S',T73,'M/S',T79,'M/S',T85,'M/S',T91,'M/S',T97, ZONMAT1A.1434
&'M/S',T103,'M/S',T109,'M/S',T115,'M/S',T121,'M/S'/) ZONMAT1A.1435
ZONMAT1A.1436
DO ROW=1,U_ROWS ZONMAT1A.1437
WRITE(6,535)ROW,(Z_V(ROW,LEVEL), LEVEL=1,P_LEVELS) ZONMAT1A.1438
IF (MOD(ROW,5).EQ.0) WRITE(6,'(1X)') ZONMAT1A.1439
END DO ZONMAT1A.1440
535 FORMAT(T1,I3,T5,F5.2,19F6.2) ZONMAT1A.1441
ENDIF ZONMAT1A.1442
ZONMAT1A.1443
IF (IPRTKE.EQ.1) THEN ZONMAT1A.1444
WRITE(6,'(1X)') ZONMAT1A.1445
WRITE(6,'(A80/)') BANNER ZONMAT1A.1446
ZONMAT1A.1447
WRITE(6,540) ZONMAT1A.1448
540 FORMAT(T2,'ROW TKE1',T14,'TKE2',T20, 'TKE3',T26,'TKE4',T32, ZONMAT1A.1449
&'TKE5',T38,'TKE6', T44,'TKE7',T50,'TKE8',T56,'TKE9',T62, ZONMAT1A.1450
&'TKE10',T68,'TKE11',T74,'TKE12',T80, 'TKE13',T86,'TKE14',T92, ZONMAT1A.1451
&'TKE15',T98, 'TKE16',T104,'TKE17',T110,'TKE18',T116, ZONMAT1A.1452
&'TKE19',T122,'TKE20') ZONMAT1A.1453
ZONMAT1A.1454
WRITE(6,543) ZONMAT1A.1455
543 FORMAT(T7,'M2/S2',T13,'M2/S2',T19,'M2/S2',T25,'M2/S2',T31,'M2/S2', ZONMAT1A.1456
& T37,'M2/S2',T43,'M2/S2',T49,'M2/S2',T55,'M2/S2',T61,'M2/S2',T67, ZONMAT1A.1457
& 'M2/S2',T73,'M2/S2',T79,'M2/S2',T85,'M2/S2',T91,'M2/S2',T97, ZONMAT1A.1458
& 'M2/S2',T103,'M2/S2',T109,'M2/S2',T115,'M2/S2',T121,'M2/S2'/) ZONMAT1A.1459
ZONMAT1A.1460
DO ROW=1,U_ROWS ZONMAT1A.1461
WRITE(6,545) ROW,(Z_TKE_LEV(ROW,LEVEL)*0.5,LEVEL=1,P_LEVELS) ZONMAT1A.1462
IF(MOD(ROW,5).EQ.0) WRITE(6,'(1X)') ZONMAT1A.1463
END DO ZONMAT1A.1464
545 FORMAT(T1,I3,T7,F5.1,19F6.1) ZONMAT1A.1465
ZONMAT1A.1466
WRITE(6,'(1X)') ZONMAT1A.1467
WRITE(6,'(A80/)') BANNER ZONMAT1A.1468
ZONMAT1A.1469
WRITE(6,550) ZONMAT1A.1470
550 FORMAT(T2,'ROW ZKE1',T14,'ZKE2',T20, 'ZKE3',T26,'ZKE4',T32, ZONMAT1A.1471
&'ZKE5',T38,'ZKE6', T44,'ZKE7',T50,'ZKE8',T56,'ZKE9',T62, ZONMAT1A.1472
&'ZKE10',T68,'ZKE11',T74,'ZKE12',T80, 'ZKE13',T86,'ZKE14',T92, ZONMAT1A.1473
&'ZKE15',T98, 'ZKE16',T104,'ZKE17',T110,'ZKE18',T116, ZONMAT1A.1474
&'ZKE19',T122,'ZKE20') ZONMAT1A.1475
ZONMAT1A.1476
WRITE(6,543) ZONMAT1A.1477
ZONMAT1A.1478
DO ROW=1,U_ROWS ZONMAT1A.1479
WRITE(6,545) ROW,(Z_ZKE_LEV(ROW,LEVEL)*0.5,LEVEL=1,P_LEVELS) ZONMAT1A.1480
IF(MOD(ROW,5).EQ.0) WRITE(6,'(1X)') ZONMAT1A.1481
END DO ZONMAT1A.1482
ZONMAT1A.1483
WRITE(6,'(1X)') ZONMAT1A.1484
WRITE(6,'(A80/)') BANNER ZONMAT1A.1485
ZONMAT1A.1486
WRITE(6,560) ZONMAT1A.1487
560 FORMAT(T2,'ROW EKE1',T14,'EKE2',T20, 'EKE3',T26,'EKE4',T32, ZONMAT1A.1488
&'EKE5',T38,'EKE6', T44,'EKE7',T50,'EKE8',T56,'EKE9',T62, ZONMAT1A.1489
&'EKE10',T68,'EKE11',T74,'EKE12',T80, 'EKE13',T86,'EKE14',T92, ZONMAT1A.1490
&'EKE15',T98, 'EKE16',T104,'EKE17',T110,'EKE18',T116, ZONMAT1A.1491
&'EKE19',T122,'EKE20') ZONMAT1A.1492
ZONMAT1A.1493
WRITE(6,543) ZONMAT1A.1494
ZONMAT1A.1495
DO ROW=1,U_ROWS ZONMAT1A.1496
WRITE(6,545) ROW,(Z_EKE_LEV(ROW,LEVEL)*0.5,LEVEL=1,P_LEVELS) ZONMAT1A.1497
IF(MOD(ROW,5).EQ.0) WRITE(6,'(1X)') ZONMAT1A.1498
END DO ZONMAT1A.1499
ENDIF ZONMAT1A.1500
ZONMAT1A.1501
IF(IPRTTEMP.EQ.1) THEN ZONMAT1A.1502
WRITE(6,'(1X)') ZONMAT1A.1503
WRITE(6,'(A80/)') BANNER ZONMAT1A.1504
WRITE(6,615) ZONMAT1A.1505
615 FORMAT(T2,'ROW T1',T13,'T2',T19,'T3',T25,'T4',T31,'T5',T37,'T6', ZONMAT1A.1506
&T43,'T7',T49,'T8',T55,'T9',T61,'T10',T67,'T11',T73,'T12',T79, ZONMAT1A.1507
&'T13',T85,'T14',T91,'T15',T97,'T16',T103,'T17',T109,'T18',T115, ZONMAT1A.1508
&'T19',T121,'T20') ZONMAT1A.1509
WRITE(6,625) ZONMAT1A.1510
625 FORMAT(T6,'DEG.K',T12,'DEG.K',T18,'DEG.K',T24,'DEG.K',T30,'DEG.K', ZONMAT1A.1511
&T36,'DEG.K',T42,'DEG.K',T48,'DEG.K',T54,'DEG.K',T60,'DEG.K',T66, ZONMAT1A.1512
&'DEG.K',T72,'DEG.K',T78,'DEG.K',T84,'DEG.K',T90,'DEG.K',T96, ZONMAT1A.1513
&'DEG.K',T102,'DEG.K',T108,'DEG.K',T114,'DEG.K',T120,'DEG.K'/) ZONMAT1A.1514
DO ROW=1,P_ROWS ZONMAT1A.1515
WRITE(6,635)ROW,(Z_T(ROW,LEVEL), LEVEL=1,P_LEVELS) ZONMAT1A.1516
IF (MOD(ROW,5).EQ.0) WRITE(6,'(1X)') ZONMAT1A.1517
END DO ZONMAT1A.1518
635 FORMAT(T1,I3,T5,F5.1,19F6.1) ZONMAT1A.1519
ZONMAT1A.1520
IF(IPRTVAR.EQ.1) THEN ZONMAT1A.1521
WRITE(6,'(1X)') ZONMAT1A.1522
WRITE(6,'(A80/)') BANNER ZONMAT1A.1523
WRITE(6,640) ZONMAT1A.1524
640 FORMAT(T2,'ROW TVR1',T13,'TVR2',T19, 'TVR3',T25,'TVR4',T31, ZONMAT1A.1525
&'TVR5',T37,'TVR6', T43,'TVR7',T49,'TVR8',T55,'TVR9',T61, 'TVR10', ZONMAT1A.1526
&T67,'TVR11',T73,'TVR12',T79,'TVR13',T85,'TVR14',T91, 'TVR15',T97, ZONMAT1A.1527
& 'TVR16',T103,'TVR17',T109,'TVR18',T115, 'TVR19',T121,'TVR20') ZONMAT1A.1528
ZONMAT1A.1529
WRITE(6,643) ZONMAT1A.1530
643 FORMAT(T8,'K2',T14,'K2',T20,'K2',T26,'K2',T32,'K2', T38,'K2',T44, ZONMAT1A.1531
&'K2',T50,'K2',T56,'K2',T62,'K2',T68, 'K2',T74,'K2',T80,'K2', T86, ZONMAT1A.1532
&'K2',T92,'K2',T98, 'K2',T104,'K2',T110,'K2',T116,'K2',T122,'K2'/) ZONMAT1A.1533
ZONMAT1A.1534
DO ROW=1,P_ROWS ZONMAT1A.1535
WRITE(6,646) ROW,(Z_VAR_T_LEV(ROW,LEVEL),LEVEL=1,P_LEVELS) ZONMAT1A.1536
IF(MOD(ROW,5).EQ.0) WRITE(6,'(1X)') ZONMAT1A.1537
END DO ZONMAT1A.1538
646 FORMAT(T1,I3,T7,F5.0,7F6.0,12F6.1) ZONMAT1A.1539
ENDIF ZONMAT1A.1540
ENDIF ZONMAT1A.1541
ZONMAT1A.1542
IF(IPRTQ.EQ.1) THEN ZONMAT1A.1543
WRITE(6,'(1X)') ZONMAT1A.1544
WRITE(6,'(A80/)') BANNER ZONMAT1A.1545
WRITE(6,715) ZONMAT1A.1546
715 FORMAT(T2,'ROW Q1',T13,'Q2',T20,'Q3',T27,'Q4',T34,'Q5',T41,'Q6', ZONMAT1A.1547
&T47,'Q7',T53,'Q8',T59,'Q9',T65,'Q10',T71,'Q11',T77,'Q12',T83, ZONMAT1A.1548
&'Q13',T89,'Q14',T95,'Q15',T101,'Q16',T107,'Q17',T113,'Q18',T119, ZONMAT1A.1549
&'Q19',T125,'Q20') ZONMAT1A.1550
WRITE(6,725) ZONMAT1A.1551
725 FORMAT(T5,'GM/KG',T12,'GM/KG',T19,'GM/KG',T26,'GM/KG',T33,'GM/KG', ZONMAT1A.1552
&T40,'GM/KG',T46,'GM/KG',T52,'GM/KG',T58,'GM/KG',T64,'GM/KG',T70, ZONMAT1A.1553
&'GM/KG',T76,'GM/KG',T82,'GM/KG',T88,'GM/KG',T94,'GM/KG',T100, ZONMAT1A.1554
&'GM/KG',T106,'GM/KG',T112,'GM/KG',T118,'GM/KG',T124,'GM/KG'/) ZONMAT1A.1555
DO ROW=1,P_ROWS ZONMAT1A.1556
WRITE(6,735)ROW,(Z_Q(ROW,LEVEL)*1000., LEVEL=1,Q_LEVELS) ZONMAT1A.1557
IF (MOD(ROW,5).EQ.0) WRITE(6,'(1X)') ZONMAT1A.1558
END DO ZONMAT1A.1559
735 FORMAT(T1,I3,T5,F6.3,4F7.3,15F6.3) ZONMAT1A.1560
IF (IPRTVAR.EQ.1) THEN ZONMAT1A.1561
WRITE(6,'(1X)') ZONMAT1A.1562
WRITE(6,'(A80/)') BANNER ZONMAT1A.1563
ZONMAT1A.1564
WRITE(6,647) ZONMAT1A.1565
647 FORMAT(T2,'ROW QVR1',T13,'QVR2',T19, 'QVR3',T25,'QVR4',T31, ZONMAT1A.1566
&'QVR5',T37,'QVR6', T43,'QVR7',T49,'QVR8',T55,'QVR9',T60, 'QVR10', ZONMAT1A.1567
&T66,'QVR11',T72,'QVR12',T78,'QVR13',T84,'QVR14',T90, 'QVR15',T96, ZONMAT1A.1568
& 'QVR16',T102,'QVR17',T108,'QVR18',T114, 'QVR19',T120,'QVR20') ZONMAT1A.1569
ZONMAT1A.1570
WRITE(6,648) ZONMAT1A.1571
648 FORMAT(T6,'G2KG2',T12,'G2KG2',T18,'G2KG2',T24, 'G2KG2',T30, ZONMAT1A.1572
&'G2KG2', T36,'G2KG2',T42, 'G2KG2',T48,'G2KG2',T54,'G2KG2',T60, ZONMAT1A.1573
&'G2KG2', T66, 'G2KG2',T72,'G2KG2',T78,'G2KG2', T84, 'G2KG2',T90, ZONMAT1A.1574
&'G2KG2',T96, 'G2KG2',T102,'G2KG2', ZONMAT1A.1575
& T108,'G2KG2',T114,'G2KG2',T120,'G2KG2'/) ZONMAT1A.1576
ZONMAT1A.1577
DO ROW=1,P_ROWS ZONMAT1A.1578
WRITE(6,649)ROW,(Z_VAR_Q_LEV(ROW,LEVEL)*1.E6 ,LEVEL=1,Q_LEVELS) ZONMAT1A.1579
IF(MOD(ROW,5).EQ.0) WRITE(6,'(1X)') ZONMAT1A.1580
END DO ZONMAT1A.1581
649 FORMAT(T1,I3,T6,F5.1,5F6.1,4F6.2,10F6.3) ZONMAT1A.1582
ZONMAT1A.1583
ENDIF ZONMAT1A.1584
ENDIF ZONMAT1A.1585
CL ZONMAT1A.1586
CL 9.2 ** Quarter Global Means Output ** ZONMAT1A.1587
CL ZONMAT1A.1588
WRITE(6,'(1X)') ZONMAT1A.1589
WRITE(6,'(A80/)') BANNER ZONMAT1A.1590
WRITE(6,15) ZONMAT1A.1591
15 FORMAT(T40,'QUARTER GLOBAL MEAN STATISTICS'/) ZONMAT1A.1592
WRITE(6,116) ZONMAT1A.1593
116 FORMAT(T3,'ROWS',T12,'PMSL',T20,'PSTAR',T28,'SNOWD',T37,'CANOPYW', ZONMAT1A.1594
&T47,'SH (W/m2)',T62,'LSRN',T68,'CVRN', ZONMAT1A.1595
&T76,'SOILM',T82,'T*ALL',T89,'T*SEA',T96,'T*LAND',T103,'SOILT1', AJS1F401.1549
&T110,'SOILT2',T117,'SOILT3',T124,'SOILT4') AJS1F401.1550
WRITE(6,126) ZONMAT1A.1598
126 FORMAT(T4,'TO',T13,'MB',T21,'MB',T28,'MM',T38,'MM',T45,'All Land ZONMAT1A.1599
& Sea ',T62,'MM/',T68,'MM/',T76,'MM',T82,'DEG.K' ZONMAT1A.1600
&,T89,'DEG.K',T96,'DEG.K',T103,'DEG.K',T110,'DEG.K',T117,'DEG.K', AJS1F401.1551
& T124,'DEGK') AJS1F401.1552
WRITE(6,136) ZONMAT1A.1602
136 FORMAT(T62,'DAY',T68,'DAY') ZONMAT1A.1603
DO 310, QUART=1,4 ZONMAT1A.1604
WRITE(6,146)START_ROW(1,QUART),END_ROW(1,QUART), ZONMAT1A.1605
& (G_Q_PMSL(QUART)*0.01), ZONMAT1A.1606
&(G_Q_PSTAR(QUART)*0.01),G_Q_SNOWD(QUART),G_Q_CANOPYW(QUART), ZONMAT1A.1607
&G_Q_SH(QUART),G_Q_L_SH(QUART),G_Q_S_SH(QUART),G_Q_LSRN(QUART), ZONMAT1A.1608
&G_Q_CVRN(QUART),G_Q_SOILM(QUART),G_Q_TSTAR(QUART), ZONMAT1A.1609
&G_Q_S_TSTAR(QUART),G_Q_L_TSTAR(QUART), ZONMAT1A.1610
&(G_Q_SOILT(QUART,LEVEL),LEVEL=1,ST_LEVELS) AJS1F401.1553
310 CONTINUE ZONMAT1A.1612
146 FORMAT(T2,I3,':',T6,I3,T10,F7.2,T18,F7.2,T26,F9.3,T36,F7.3,T44, ZONMAT1A.1613
&F5.2,T50,F5.2,T56,F5.2,T62,F5.3,T68,F5.3,T76,F5.1,T82, ZONMAT1A.1614
&F6.2,T89,F6.2,T96,F6.2,T103,F6.2,T110,F6.2,T117,F6.2,T124,F6.2, AJS1F401.1554
&T131,F6.2) AJS1F401.1555
WRITE(6,117)(G_PMSL*0.01),(G_PSTAR*0.01),G_SNOWD,G_CANOPYW,G_SH, ZONMAT1A.1616
&G_L_SH,G_S_SH,G_LSRN,G_CVRN,G_SOILM,G_TSTAR, ZONMAT1A.1617
&G_S_TSTAR,G_L_TSTAR,(G_SOILT(LEVEL), LEVEL=1,ST_LEVELS) AJS1F401.1556
117 FORMAT(T3,'GLOBAL ',T10,F7.2,T18,F7.2,T26,F9.3,T36,F7.3,T44, ZONMAT1A.1619
&F5.2,T50,F5.2,T56,F5.2,T62,F5.3,T68,F5.3,T76,F5.1,T82, ZONMAT1A.1620
&F6.2,T89,F6.2,T96,F6.2,T103,F6.2,T110,F6.2,T117,F6.2,T124,F6.2) ZONMAT1A.1621
ZONMAT1A.1622
WRITE(6,'(1X)') ZONMAT1A.1623
CTITLE =' Rows Land Sea Evaporation(mm/day) Rainfall(mm ZONMAT1A.1624
&/day) Snowfall(mm/day) Sublimation(mm/day) Soil Evap' ZONMAT1A.1625
CHEAD1 =' to All Land Sea All Land ZONMAT1A.1626
&Sea All Land Sea All Land Sea (mm/day)' ZONMAT1A.1627
WRITE(6,'(1X)') ZONMAT1A.1628
WRITE(6,'(A130)') CTITLE ZONMAT1A.1629
WRITE(6,'(A130)') CHEAD1 ZONMAT1A.1630
DO 1310, QUART=1,4 ZONMAT1A.1631
WRITE(6,1146)START_ROW(1,QUART),END_ROW(1,QUART),QGLPTS(QUART,1), ZONMAT1A.1632
&QGSPTS(QUART,1),(G_Q_EVAP(QUART)*86400.), ZONMAT1A.1633
&(G_Q_L_EVAP(QUART)*86400.), ZONMAT1A.1634
&(G_Q_S_EVAP(QUART)*86400.),G_Q_PPTN(QUART),G_Q_L_PPTN(QUART), ZONMAT1A.1635
&G_Q_S_PPTN(QUART),G_Q_SNOW(QUART),G_Q_L_SNOW(QUART), ZONMAT1A.1636
&G_Q_S_SNOW(QUART),G_Q_SUBL(QUART),G_Q_L_SUBL(QUART), ZONMAT1A.1637
&G_Q_S_SUBL(QUART),G_Q_L_SOEV(QUART) ZONMAT1A.1638
1310 CONTINUE ZONMAT1A.1639
1146 FORMAT(T2,I3,':',T6,I3,2I6,T22,4(1X,3(1X,F5.3)),2X,F5.3) ZONMAT1A.1640
WRITE(6,1117) GLPTS(1),GSPTS(1),(G_EVAP*86400.),(G_L_EVAP*86400.), ZONMAT1A.1641
&(G_S_EVAP*86400.),G_PPTN,G_L_PPTN,G_S_PPTN,G_SNOW,G_L_SNOW, ZONMAT1A.1642
&G_S_SNOW,G_SUBL,G_L_SUBL,G_S_SUBL,G_L_SOEV ZONMAT1A.1643
1117 FORMAT(T2,'GLOBAL ',2I6,T22,4(1X,3(1X,F5.3)),2X,F5.3) ZONMAT1A.1644
ZONMAT1A.1645
IF(IPRTWIND.EQ.1) THEN ZONMAT1A.1646
WRITE(6,'(1X)') ZONMAT1A.1647
WRITE(6,216) ZONMAT1A.1648
216 FORMAT(T3,'ROWS',T10,'TATMOS QATMOS',T38,'UATMOS VATMOS',T57, ZONMAT1A.1649
&'TAUX (N/M2)',T78,'TAUY (N/M2)',T95,'TKE',T102,'ZKE',T109, ZONMAT1A.1650
&'EKE') ZONMAT1A.1651
WRITE(6,226) ZONMAT1A.1652
226 FORMAT(T10,'DEG.K',T17,'GM/KG',T40,'M/S',T47,'M/S',T53,'ALL',T60, ZONMAT1A.1653
&'SEA',T67,'LAND',T74,'ALL',T81,'SEA',T88,'LAND',T94,'M2/S2', ZONMAT1A.1654
&T101,'M2/S2',T108,'M2/S2'/) ZONMAT1A.1655
ZONMAT1A.1656
DO 320, QUART=1,4 ZONMAT1A.1657
WRITE(6,236)START_ROW(1,QUART),END_ROW(1,QUART),G_Q_TATMOS(QUART), ZONMAT1A.1658
&(G_Q_QATMOS(QUART)*1000.),START_ROW(2,QUART),END_ROW(2,QUART), ZONMAT1A.1659
&G_Q_UATMOS(QUART),G_Q_VATMOS(QUART),G_Q_TAUX(QUART), ZONMAT1A.1660
&G_Q_S_TAUX(QUART),G_Q_L_TAUX(QUART),G_Q_TAUY(QUART), ZONMAT1A.1661
&G_Q_S_TAUY(QUART),G_Q_L_TAUY(QUART), ZONMAT1A.1662
&G_Q_TKE(QUART)*0.5,G_Q_ZKE(QUART)*0.5,G_Q_EKE(QUART)*0.5 ZONMAT1A.1663
320 CONTINUE ZONMAT1A.1664
ZONMAT1A.1665
236 FORMAT(T2,I3,':',T6,I3,T10,F6.2,F6.2,T30,I3,':',T34,I3,T38,F6.2, ZONMAT1A.1666
&F6.2,6F7.3,3F7.2) ZONMAT1A.1667
WRITE(6,127)G_TATMOS,(G_QATMOS*1000.),G_UATMOS,G_VATMOS,G_TAUX, ZONMAT1A.1668
& G_S_TAUX,G_L_TAUX,G_TAUY,G_S_TAUY,G_L_TAUY, ZONMAT1A.1669
& G_TKE*0.5,G_ZKE*0.5,G_EKE*0.5 ZONMAT1A.1670
127 FORMAT(T3,'GLOBAL ',T10,F6.2,F6.2,T31,'GLOBAL',T38,F6.2,F6.2, ZONMAT1A.1671
&6F7.3,3F7.2) ZONMAT1A.1672
WRITE(6,'(1X)') ZONMAT1A.1673
WRITE (6,316) ZONMAT1A.1674
316 FORMAT(T3,'ROW',T11,'U1',T17,'U2',T23,'U3',T29,'U4',T35,'U5', ZONMAT1A.1675
&T41,'U6',T47,'U7',T53,'U8',T59,'U9',T65,'U10',T71,'U11',T77,'U12', ZONMAT1A.1676
&T83,'U13',T89,'U14',T95,'U15',T101,'U16',T107,'U17',T113,'U18', ZONMAT1A.1677
&T119,'U19',T125,'U20') ZONMAT1A.1678
WRITE(6,326) ZONMAT1A.1679
326 FORMAT(T11,'M/S',T17,'M/S',T23,'M/S',T29,'M/S',T35,'M/S', ZONMAT1A.1680
&T41,'M/S',T47,'M/S',T53,'M/S',T59,'M/S',T65,'M/S',T71, ZONMAT1A.1681
&'M/S',T77,'M/S',T83,'M/S',T89,'M/S',T95,'M/S',T101, ZONMAT1A.1682
&'M/S',T107,'M/S',T113,'M/S',T119,'M/S',T125,'M/S'/) ZONMAT1A.1683
DO QUART=1,4 ZONMAT1A.1684
WRITE(6,336)START_ROW(2,QUART),END_ROW(2,QUART), ZONMAT1A.1685
& (G_Q_U(QUART,LEVEL), LEVEL=1,P_LEVELS) ZONMAT1A.1686
END DO ZONMAT1A.1687
336 FORMAT(T2,I3,':',T6,I3,T10,F5.2,19F6.2) ZONMAT1A.1688
WRITE(6,137)(G_U(LEVEL), LEVEL=1,P_LEVELS) ZONMAT1A.1689
137 FORMAT(T3,'GLOBAL ',T10,F5.2,19F6.2) ZONMAT1A.1690
WRITE(6,'(1X)') ZONMAT1A.1691
WRITE (6,416) ZONMAT1A.1692
416 FORMAT(T3,'ROW',T11,'V1',T17,'V2',T23,'V3',T29,'V4',T35,'V5', ZONMAT1A.1693
&T41,'V6',T47,'V7',T53,'V8',T59,'V9',T65,'V10',T71,'V11',T77,'V12', ZONMAT1A.1694
&T83,'V13',T89,'V14',T95,'V15',T101,'V16',T107,'V17',T113,'V18', ZONMAT1A.1695
&T119,'V19',T125,'V20') ZONMAT1A.1696
WRITE(6,426) ZONMAT1A.1697
426 FORMAT(T11,'M/S',T17,'M/S',T23,'M/S',T29,'M/S',T35,'M/S', ZONMAT1A.1698
&T41,'M/S',T47,'M/S',T53,'M/S',T59,'M/S',T65,'M/S',T71, ZONMAT1A.1699
&'M/S',T77,'M/S',T83,'M/S',T89,'M/S',T95,'M/S',T101, ZONMAT1A.1700
&'M/S',T107,'M/S',T113,'M/S',T119,'M/S',T125,'M/S'/) ZONMAT1A.1701
DO QUART=1,4 ZONMAT1A.1702
WRITE(6,436)START_ROW(2,QUART),END_ROW(2,QUART), ZONMAT1A.1703
& (G_Q_V(QUART,LEVEL), LEVEL=1,P_LEVELS) ZONMAT1A.1704
END DO ZONMAT1A.1705
436 FORMAT(T2,I3,':',T6,I3,T10,F5.2,19F6.2) ZONMAT1A.1706
WRITE(6,147)(G_V(LEVEL), LEVEL=1,P_LEVELS) ZONMAT1A.1707
147 FORMAT(T3,'GLOBAL ',T10,F5.2,19F6.2) ZONMAT1A.1708
WRITE(6,1140) ZONMAT1A.1709
ENDIF ZONMAT1A.1710
IF(IPRTKE.EQ.1) THEN ZONMAT1A.1711
1140 FORMAT(T3,'KINETIC ENERGY IN M2/S2') ZONMAT1A.1712
WRITE(6,1141)(G_TKE_LEV(LEVEL)*0.5,LEVEL=1,P_LEVELS) ZONMAT1A.1713
1141 FORMAT(T3,'TKE',T10,F5.1,19F6.1) ZONMAT1A.1714
WRITE(6,1142)(G_ZKE_LEV(LEVEL)*0.5,LEVEL=1,P_LEVELS) ZONMAT1A.1715
1142 FORMAT(T3,'ZKE',T10,F5.1,19F6.1) ZONMAT1A.1716
WRITE(6,1143)(G_EKE_LEV(LEVEL)*0.5,LEVEL=1,P_LEVELS) ZONMAT1A.1717
1143 FORMAT(T3,'EKE',T10,F5.1,19F6.1) ZONMAT1A.1718
END IF ZONMAT1A.1719
ZONMAT1A.1720
IF(IPRTTEMP.EQ.1) THEN ZONMAT1A.1721
WRITE(6,'(1X)') ZONMAT1A.1722
WRITE (6,516) ZONMAT1A.1723
516 FORMAT(T3,'ROWS',T10,' T1 T2 T3 T4 T5 T6' ZONMAT1A.1724
&,T53,'T7 T8 T9 T10') ZONMAT1A.1725
WRITE(6,526) ZONMAT1A.1726
526 FORMAT(T10,'DEG.K DEG.K DEG.K DEG.K DEG.K DEG.K DEG.K', ZONMAT1A.1727
&T59,'DEG.K DEG.K DEG.K'/) ZONMAT1A.1728
DO QUART=1,4 ZONMAT1A.1729
WRITE(6,536)START_ROW(1,QUART),END_ROW(1,QUART), ZONMAT1A.1730
& (G_Q_T(QUART,LEVEL), LEVEL=1,10) ZONMAT1A.1731
END DO ZONMAT1A.1732
536 FORMAT(T2,I3,':',T6,I3,T9,10F7.2) ZONMAT1A.1733
WRITE(6,157)(G_T(LEVEL), LEVEL=1,10) ZONMAT1A.1734
157 FORMAT(T3,'GLOBAL ',T9,10F7.2) ZONMAT1A.1735
WRITE(6,1150) ZONMAT1A.1736
1150 FORMAT(T2,'TEMPERATURE VARIANCE IN DEG. K2') ZONMAT1A.1737
WRITE(6,1151) (G_VAR_T_LEV(LEVEL),LEVEL=1,10) ZONMAT1A.1738
1151 FORMAT(T9,10F7.2) ZONMAT1A.1739
WRITE(6,'(1X)') ZONMAT1A.1740
WRITE (6,546) ZONMAT1A.1741
546 FORMAT(T3,'ROWS',T10,' T11 T12 T13 T14 T15 T16' ZONMAT1A.1742
&,T53,'T17 T18 T19 T20') ZONMAT1A.1743
WRITE(6,556) ZONMAT1A.1744
556 FORMAT(T10,'DEG.K DEG.K DEG.K DEG.K DEG.K DEG.K DEG.K', ZONMAT1A.1745
&T59,'DEG.K DEG.K DEG.K'/) ZONMAT1A.1746
DO QUART=1,4 ZONMAT1A.1747
WRITE(6,566)START_ROW(1,QUART),END_ROW(1,QUART), ZONMAT1A.1748
& (G_Q_T(QUART,LEVEL), LEVEL=11,P_LEVELS) ZONMAT1A.1749
END DO ZONMAT1A.1750
ZONMAT1A.1751
566 FORMAT(T2,I3,':',T6,I3,T9,10F7.2) ZONMAT1A.1752
WRITE(6,167)(G_T(LEVEL), LEVEL=11,P_LEVELS) ZONMAT1A.1753
167 FORMAT(T3,'GLOBAL ',T9,10F7.2) ZONMAT1A.1754
WRITE(6,1160) ZONMAT1A.1755
1160 FORMAT(T3,'TEMPERATURE VARIANCE IN DEG. K2') ZONMAT1A.1756
WRITE(6,1161) (G_VAR_T_LEV(LEVEL),LEVEL=11,P_LEVELS) ZONMAT1A.1757
1161 FORMAT(T9,10F7.2) ZONMAT1A.1758
ENDIF ZONMAT1A.1759
ZONMAT1A.1760
IF(IPRTQ.EQ.1) THEN ZONMAT1A.1761
WRITE(6,'(1X)') ZONMAT1A.1762
WRITE (6,616) ZONMAT1A.1763
616 FORMAT(T3,'ROW',T11,'Q1',T18,'Q2',T25,'Q3',T32,'Q4',T39,'Q5', ZONMAT1A.1764
&T45,'Q6',T51,'Q7',T57,'Q8',T63,'Q9',T69,'Q10',T75,'Q11',T81,'Q12', ZONMAT1A.1765
&T87,'Q13',T93,'Q14',T99,'Q15',T105,'Q16',T111,'Q17',T117,'Q18', ZONMAT1A.1766
&T123,'Q19',T129,'Q20') ZONMAT1A.1767
WRITE(6,626) ZONMAT1A.1768
626 FORMAT(T10,'GM/KG',T17,'GM/KG',T24,'GM/KG',T31,'GM/KG',T38,'GM/KG' ZONMAT1A.1769
&,T44,'GM/KG',T50,'GM/KG',T56,'GM/KG',T62,'GM/KG',T68,'GM/KG',T74, ZONMAT1A.1770
&'GM/KG',T80,'GM/KG',T86,'GM/KG',T92,'GM/KG',T98,'GM/KG',T104, ZONMAT1A.1771
&'GM/KG',T110,'GM/KG',T116,'GM/KG',T122,'GM/KG',T128,'GM/KG'/) ZONMAT1A.1772
ZONMAT1A.1773
DO QUART=1,4 ZONMAT1A.1774
WRITE(6,636)START_ROW(1,QUART),END_ROW(1,QUART), ZONMAT1A.1775
& (G_Q_Q(QUART,LEVEL)*1000., LEVEL=1,Q_LEVELS) ZONMAT1A.1776
END DO ZONMAT1A.1777
636 FORMAT(T1,I3,':',T5,I3,T8,5F7.3,15F6.3) ZONMAT1A.1778
WRITE(6,177)(G_Q(LEVEL)*1000., LEVEL=1,Q_LEVELS) ZONMAT1A.1779
177 FORMAT(T2,'GLOBAL ',T8,5F7.3,15F6.3) ZONMAT1A.1780
WRITE(6,1170) ZONMAT1A.1781
1170 FORMAT(T3,'MOISTURE VARIANCE IN G2/KG2') ZONMAT1A.1782
WRITE(6,1171) (G_VAR_Q_LEV(LEVEL)*1.E6,LEVEL=1,Q_LEVELS) ZONMAT1A.1783
1171 FORMAT(T8,5F7.3,15F6.3) ZONMAT1A.1784
ENDIF ZONMAT1A.1785
ZONMAT1A.1786
WRITE(6,'(1X)') ZONMAT1A.1787
WRITE(6,1180) ZONMAT1A.1788
1180 FORMAT(T3,'GLOBAL MEANS') ZONMAT1A.1789
WRITE(6,1181) G_TKE*0.5, G_ZKE*0.5, G_EKE*0.5 ZONMAT1A.1790
1181 FORMAT(T5,'TKE',T10,F8.1,T19,'ZKE',T24,F8.1,T33,'EKE',F8.1, ZONMAT1A.1791
& '(M2 S-2)') ZONMAT1A.1792
WRITE(6,1182) G_VAR_T,G_VAR_Q*1.E6 ZONMAT1A.1793
1182 FORMAT(T5,'T_VARIANCE',T15,F8.1,T23,' K2 Q_VARIANCE',T38,F8.1, ZONMAT1A.1794
& T47,' G2/KG2') ZONMAT1A.1795
CL---------------------------------------------------------------------- ZONMAT1A.1796
C Extra variables added and only printed if available ZONMAT1A.1797
C ZONMAT1A.1798
IF (IPRTEXTRA.EQ.1) THEN ZONMAT1A.1799
IF(IHYDRO.GT.1) THEN ARS2F304.165
CTITLE ='Quarter-globe and Global means over land from hydrolog ZONMAT1A.1801
&y section' ZONMAT1A.1802
CHEAD1 =' Thro Fall Snowmelt Sur runoff Subsur Runof ZONMAT1A.1803
&f' ZONMAT1A.1804
CHEAD2 =' mm/day mm/day mm/day mm/day' ZONMAT1A.1805
WRITE(6,'(1X)') ZONMAT1A.1806
WRITE(6,'(A130)') CTITLE ZONMAT1A.1807
WRITE(6,'(A130)') CHEAD1 ZONMAT1A.1808
WRITE(6,'(A130)') CHEAD2 ZONMAT1A.1809
ZONMAT1A.1810
DO QUART=1,4 ZONMAT1A.1811
WRITE(6,2004)START_ROW(1,QUART),END_ROW(1,QUART), ZONMAT1A.1812
& G_Q_L_THRF(QUART),G_Q_L_SNML(QUART),G_Q_L_SFRU(QUART), ZONMAT1A.1813
& G_Q_L_SBRU(QUART) ZONMAT1A.1814
END DO ZONMAT1A.1815
2004 FORMAT(T1,I3,':',T5,I3,T8,4(3x,F8.3)) ZONMAT1A.1816
WRITE(6,2006) G_L_THRF,G_L_SNML,G_L_SFRU,G_L_SBRU ZONMAT1A.1817
2006 FORMAT(T2,'GLOBAL ',T8,4(3x,F8.3)) ZONMAT1A.1818
END IF ZONMAT1A.1819
ZONMAT1A.1820
IF(IRAD.GT.1) THEN ARS2F304.166
CTITLE =' Quarter-globe and Global means - Radiation in W/m2' ZONMAT1A.1822
CHEAD1 =' Total Surface Surface Solar Total ZONMAT1A.1823
& TOA LW TOA SW TOA Albedo' ZONMAT1A.1824
CHEAD2 =' All Land Sea All land Sea rad ZONMAT1A.1825
& rad rad' ZONMAT1A.1826
WRITE(6,'(1X)') ZONMAT1A.1827
WRITE(6,'(A130)') CTITLE ZONMAT1A.1828
WRITE(6,'(A130)') CHEAD1 ZONMAT1A.1829
WRITE(6,'(A130)') CHEAD2 ZONMAT1A.1830
DO I=1,4 ZONMAT1A.1831
WRITE(6,2014)START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1832
& G_Q_SDTR(I),G_Q_L_SDTR(I),G_Q_S_SDTR(I), ZONMAT1A.1833
& G_Q_SDSR(I),G_Q_L_SDSR(I),G_Q_S_SDSR(I), ZONMAT1A.1834
& G_Q_TDTR(I),G_Q_TOLR(I),G_Q_TOSW(I),G_Q_ALBEDO(I) ZONMAT1A.1835
END DO ZONMAT1A.1836
2014 FORMAT(T1,I3,':',T5,I3,T8,6(1X,F6.1),T52,3(F7.2,2X),F7.4) ZONMAT1A.1837
WRITE(6,2016) G_SDTR,G_L_SDTR,G_S_SDTR,G_SDSR,G_L_SDSR, ZONMAT1A.1838
& G_S_SDSR,G_TDTR,G_TOLR,G_TOSW,G_ALBEDO ZONMAT1A.1839
2016 FORMAT(T2,'GLOBAL ',T8,6(1X,F6.1),T52,3(F7.2,2X),F7.4) ZONMAT1A.1840
END IF ZONMAT1A.1841
ZONMAT1A.1842
IF (IFLUXL.GT.1) THEN ARS2F304.167
CTITLE =' Quarter-globe and Global means - Fluxes over land ZONMAT1A.1844
&and over whole atmosphere' ZONMAT1A.1845
CHEAD1 =' Soil Flux Snow Bud Can Bud Soil Bud Water(l ZONMAT1A.1846
&+s) Soil En Energy Flux (l+s)' ZONMAT1A.1847
CHEAD2 =' W/m2 mm/day mm/day mm/day mm/day ZONMAT1A.1848
& W/m2 W/m2' ZONMAT1A.1849
WRITE(6,'(1X)') ZONMAT1A.1850
WRITE(6,'(A130)') CTITLE ZONMAT1A.1851
WRITE(6,'(A130)') CHEAD1 ZONMAT1A.1852
WRITE(6,'(A130)') CHEAD2 ZONMAT1A.1853
DO I=1,4 ZONMAT1A.1854
WRITE(6,2024)START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1855
& G_Q_L_SOHF(I),G_Q_L_WFSS(I),G_Q_L_WFCA(I),G_Q_L_WFSO(I), ZONMAT1A.1856
& G_Q_WAFL(I),G_Q_L_ENFS(I),G_Q_ENFL(I) ZONMAT1A.1857
END DO ZONMAT1A.1858
2024 FORMAT(T1,I3,':',T5,I3,T8,7(F8.3,2x)) ZONMAT1A.1859
WRITE(6,2026) G_L_SOHF,G_L_WFSS,G_L_WFCA,G_L_WFSO,G_WAFL, ZONMAT1A.1860
& G_L_ENFS,G_ENFL ZONMAT1A.1861
2026 FORMAT(T2,'GLOBAL ',T8,7(F8.3,2x)) ZONMAT1A.1862
END IF ZONMAT1A.1863
C ZONMAT1A.1864
C Cloud water printouts ZONMAT1A.1865
C ZONMAT1A.1866
IF(ICLOUD.GT.1) THEN ARS2F304.168
CHEAD1 ='ROW 1 2 3 4 5 6 7 8 ZONMAT1A.1868
& 9 10 11 12 13 14 15 16 17 18 19 ZONMAT1A.1869
& 20' ZONMAT1A.1870
CHEAD2 =' to ' ZONMAT1A.1871
CTITLE =' Cloud liquid water - all levels (g/kg x1000)' ZONMAT1A.1872
WRITE(6,'(1X)') ZONMAT1A.1873
WRITE(6,'(A130)') CTITLE ZONMAT1A.1874
WRITE(6,'(A130)') CHEAD1 ZONMAT1A.1875
WRITE(6,'(A130)') CHEAD2 ZONMAT1A.1876
DO I=1,4 ZONMAT1A.1877
WRITE(6,3024)START_ROW(1,I),END_ROW(1,I),(1000000.*G_Q_CLLIQL(I, ZONMAT1A.1878
&LEVEL),LEVEL=1,Q_LEVELS) ZONMAT1A.1879
END DO ZONMAT1A.1880
WRITE(6,3026) (1000000.*G_CLLIQL(LEVEL),LEVEL=1,Q_LEVELS) ZONMAT1A.1881
CTITLE =' Cloud ice water - all levels (g/kg x1000)' ZONMAT1A.1882
WRITE(6,'(1X)') ZONMAT1A.1883
WRITE(6,'(A130)') CTITLE ZONMAT1A.1884
WRITE(6,'(A130)') CHEAD1 ZONMAT1A.1885
WRITE(6,'(A130)') CHEAD2 ZONMAT1A.1886
DO I=1,4 ZONMAT1A.1887
WRITE(6,3024)START_ROW(1,I),END_ROW(1,I),(1000000.*G_Q_CLICEL(I, ZONMAT1A.1888
&LEVEL),LEVEL=1,Q_LEVELS) ZONMAT1A.1889
END DO ZONMAT1A.1890
WRITE(6,3026) (1000000.*G_CLICEL(LEVEL),LEVEL=1,Q_LEVELS) ZONMAT1A.1891
CTITLE =' Cloud water whole atmosphere (g/kg x1000)' ZONMAT1A.1892
CHEAD1 =' ROW Liquid Ice' ZONMAT1A.1893
CHEAD2 =' to ' ZONMAT1A.1894
WRITE(6,'(1X)') ZONMAT1A.1895
WRITE(6,'(A130)') CTITLE ZONMAT1A.1896
WRITE(6,'(A130)') CHEAD1 ZONMAT1A.1897
WRITE(6,'(A130)') CHEAD2 ZONMAT1A.1898
DO I=1,4 ZONMAT1A.1899
WRITE(6,3030)START_ROW(1,I),END_ROW(1,I), ZONMAT1A.1900
& (1000000.*G_Q_CLLIQ(I)),(1000000.*G_Q_CLICE(I)) ZONMAT1A.1901
END DO ZONMAT1A.1902
WRITE(6,3032) (1000000.*G_CLLIQ),(1000000.*G_CLICE) ZONMAT1A.1903
3024 FORMAT(T1,I3,':',T5,I3,T8,20(F6.2)) ZONMAT1A.1904
3026 FORMAT(T2,'GLOBAL ',T8,20(F6.2)) ZONMAT1A.1905
3030 FORMAT(T1,I3,':',T5,I3,T8,2(3X,F6.2)) ZONMAT1A.1906
3032 FORMAT(T2,'GLOBAL ',T8,2(3X,F6.2)) ZONMAT1A.1907
ENDIF ZONMAT1A.1908
IF(ISICE.GT.1) THEN ARS2F304.169
CTITLE =' Quarter-globe means Sea-ice meaned over all sea po ZONMAT1A.1910
&ints' ZONMAT1A.1911
CHEAD1 =' Heat Flux Melt ht S Heat Sur Rad Sur SW ZONMAT1A.1912
&Energy Fl' ZONMAT1A.1913
CHEAD2 =' W/m2 W/m2 W/m2 W/m2 W/m2 ZONMAT1A.1914
& W/m2' ZONMAT1A.1915
WRITE(6,'(1X)') ZONMAT1A.1916
WRITE(6,'(A130)') CTITLE ZONMAT1A.1917
WRITE(6,'(A130)') CHEAD1 ZONMAT1A.1918
WRITE(6,'(A130)') CHEAD2 ZONMAT1A.1919
DO I=1,4 ZONMAT1A.1920
WRITE(6,3124)START_ROW(1,I),END_ROW(1,I),G_Q_SIHF(I), ZONMAT1A.1921
& G_Q_SIMH(I),G_Q_SISH(I),G_Q_SIST(I),G_Q_SISS(I),G_Q_SIEF(I) ZONMAT1A.1922
END DO ZONMAT1A.1923
WRITE(6,3126) G_SIHF,G_SIMH,G_SISH,G_SIST, ZONMAT1A.1924
& G_SISS,G_SIEF ZONMAT1A.1925
3124 FORMAT(T1,I3,':',T5,I3,T8,6(1X,F8.3)) ZONMAT1A.1926
3126 FORMAT(T2,'GLOBAL ',T8,6(1X,F8.3)) ZONMAT1A.1927
ENDIF ZONMAT1A.1928
ENDIF ZONMAT1A.1929
WRITE(6,'(1X)') ZONMAT1A.1930
CL---------------------------------------------------------------------- ZONMAT1A.1931
999 CONTINUE ZONMAT1A.1932
RETURN ZONMAT1A.1933
CL---------------------------------------------------------------------- ZONMAT1A.1934
END ZONMAT1A.1935
*ENDIF ZONMAT1A.1936