File indexing completed on 2018-03-02 18:37:18 UTC
view on githubraw file Latest commit add29e06 on 2018-01-31 20:35:05 UTC
a176fcc394 Andr*0001 #include "AIM_OPTIONS.h"
0002
0003
0004
0005
0006 SUBROUTINE AIM_DIAGNOSTICS_INIT( myThid )
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016 IMPLICIT NONE
0017
0018
0019 #include "EEPARAMS.h"
0020
0021
0022
0023
0024 INTEGER myThid
0025
0026
0027 #ifdef ALLOW_DIAGNOSTICS
0028
0029
0030
ee828acb9e Jean*0031
a176fcc394 Andr*0032
0033
0034
0035
0036
0037 INTEGER diagNum
ee828acb9e Jean*0038 INTEGER diagMate
a176fcc394 Andr*0039 CHARACTER*8 diagName
0040 CHARACTER*16 diagCode
0041 CHARACTER*16 diagUnits
0042 CHARACTER*(80) diagTitle
0043
0044
02b51191e4 Jean*0045
a176fcc394 Andr*0046 diagName = 'DIABT '
02b51191e4 Jean*0047 diagTitle =
d2b04ae38b Jean*0048 & 'Pot. Temp. Tendency (Mass-Weighted) from Diabatic Processes'
520b5fccfc Jean*0049 diagUnits = 'K/s '
5644fea420 Jean*0050 diagCode = 'SMr MR '
ee828acb9e Jean*0051 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0052 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0053
0054 diagName = 'DIABQ '
02b51191e4 Jean*0055 diagTitle =
d2b04ae38b Jean*0056 & 'Spec.Humid. Tendency (Mass-Weighted) from Diabatic Processes'
520b5fccfc Jean*0057 diagUnits = 'g/kg/s '
5644fea420 Jean*0058 diagCode = 'SMr MR '
ee828acb9e Jean*0059 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0060 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0061
520b5fccfc Jean*0062 diagName = 'RADSW '
02b51191e4 Jean*0063 diagTitle =
520b5fccfc Jean*0064 & 'Temperature Tendency due to Shortwave Radiation (TT_RSW)'
0065 diagUnits = 'K/s '
a176fcc394 Andr*0066 diagCode = 'SM ML '
ee828acb9e Jean*0067 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0068 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0069
520b5fccfc Jean*0070 diagName = 'RADLW '
02b51191e4 Jean*0071 diagTitle =
520b5fccfc Jean*0072 & 'Temperature Tendency due to Longwave Radiation (TT_RLW)'
0073 diagUnits = 'K/s '
0074 diagCode = 'SM ML '
ee828acb9e Jean*0075 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0076 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0077
0078 diagName = 'DTCONV '
520b5fccfc Jean*0079 diagTitle = 'Temperature Tendency due to Convection (TT_CNV)'
a176fcc394 Andr*0080 diagUnits = 'K/s '
5644fea420 Jean*0081 diagCode = 'SM ML '
ee828acb9e Jean*0082 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0083 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0084
0085 diagName = 'TURBT '
02b51191e4 Jean*0086 diagTitle =
520b5fccfc Jean*0087 & 'Temperature Tendency due to Turbulence in PBL (TT_PBL)'
0088 diagUnits = 'K/s '
a176fcc394 Andr*0089 diagCode = 'SM ML '
ee828acb9e Jean*0090 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0091 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0092
0093 diagName = 'DTLS '
02b51191e4 Jean*0094 diagTitle =
520b5fccfc Jean*0095 & 'Temperature Tendency due to Large-scale condens. (TT_LSC)'
0096 diagUnits = 'K/s '
a176fcc394 Andr*0097 diagCode = 'SM ML '
ee828acb9e Jean*0098 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0099 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0100
0101 diagName = 'DQCONV '
520b5fccfc Jean*0102 diagTitle = 'Spec. Humidity Tendency due to Convection (QT_CNV)'
a176fcc394 Andr*0103 diagUnits = 'g/kg/s '
0104 diagCode = 'SM MR '
ee828acb9e Jean*0105 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0106 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0107
0108 diagName = 'TURBQ '
02b51191e4 Jean*0109 diagTitle =
520b5fccfc Jean*0110 & 'Spec. Humidity Tendency due to Turbulence in PBL (QT_PBL)'
0111 diagUnits = 'g/kg/s '
a176fcc394 Andr*0112 diagCode = 'SM ML '
ee828acb9e Jean*0113 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0114 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0115
0116 diagName = 'DQLS '
02b51191e4 Jean*0117 diagTitle =
520b5fccfc Jean*0118 & 'Spec. Humidity Tendency due to Large-Scale Condens. (QT_LSC)'
0119 diagUnits = 'g/kg/s '
a176fcc394 Andr*0120 diagCode = 'SM ML '
ee828acb9e Jean*0121 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0122 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0123
0124 diagName = 'RELHUM '
0125 diagTitle = 'Relative Humidity (fraction)'
520b5fccfc Jean*0126 diagUnits = '0-1 '
5644fea420 Jean*0127 diagCode = 'SMR MR '
ee828acb9e Jean*0128 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0129 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0130
520b5fccfc Jean*0131 diagName = 'TSR '
0132 diagTitle = 'Top-of-atm. net Shortwave Radiation (+=dw)'
0133 diagUnits = 'W/m^2 '
a176fcc394 Andr*0134 diagCode = 'SM P U1 '
ee828acb9e Jean*0135 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0136 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0137
0138 diagName = 'OLR '
520b5fccfc Jean*0139 diagTitle = 'Outgoing Longwave Radiation (+=up)'
0140 diagUnits = 'W/m^2 '
0141 diagCode = 'SM P U1 '
ee828acb9e Jean*0142 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0143 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0144
0145 diagName = 'RADSWG '
520b5fccfc Jean*0146 diagTitle = 'Net Shortwave Radiation at the Ground (+=dw)'
0147 diagUnits = 'W/m^2 '
a176fcc394 Andr*0148 diagCode = 'SM P L1 '
ee828acb9e Jean*0149 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0150 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0151
0152 diagName = 'RADLWG '
520b5fccfc Jean*0153 diagTitle = 'Net Longwave Radiation at the Ground (+=up)'
0154 diagUnits = 'W/m^2 '
a176fcc394 Andr*0155 diagCode = 'SM L1 '
ee828acb9e Jean*0156 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0157 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0158
0159 diagName = 'HFLUX '
520b5fccfc Jean*0160 diagTitle = 'Sensible Heat Flux (+=up)'
0161 diagUnits = 'W/m^2 '
a176fcc394 Andr*0162 diagCode = 'SM L1 '
ee828acb9e Jean*0163 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0164 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0165
0166 diagName = 'EVAP '
520b5fccfc Jean*0167 diagTitle = 'Surface Evaporation (g/m2/s)'
0168 diagUnits = 'g/m^2/s '
a176fcc394 Andr*0169 diagCode = 'SM L1 '
ee828acb9e Jean*0170 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0171 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0172
0173 diagName = 'PRECON '
520b5fccfc Jean*0174 diagTitle = 'Convective Precipitation (g/m2/s)'
0175 diagUnits = 'g/m^2/s '
a176fcc394 Andr*0176 diagCode = 'SM P L1 '
ee828acb9e Jean*0177 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0178 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0179
0180 diagName = 'PRECLS '
520b5fccfc Jean*0181 diagTitle = 'Large Scale Precipitation (g/m2/s)'
a176fcc394 Andr*0182 diagUnits = 'g/m^2/s '
0183 diagCode = 'SM M1 '
ee828acb9e Jean*0184 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0185 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0186
bd70561229 Davi*0187 #ifdef ALLOW_CLOUD_3D
0188 diagName = 'CLDCLW '
0189 diagTitle = 'Local Cloud Fraction (0-1)'
0190 diagUnits = '0-1 '
ee828acb9e Jean*0191 diagCode = 'SM P ML '
0192 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0193 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
bd70561229 Davi*0194 #endif
0195
ee828acb9e Jean*0196 diagName = 'CLDFRC '
0197 diagTitle = 'Total Cloud Fraction (0-1)'
0198 diagUnits = '0-1 '
0199 diagCode = 'SM P M1 '
0200 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0201 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0202
02b51191e4 Jean*0203 diagName = 'CLDPRS '
520b5fccfc Jean*0204 diagTitle = 'Cloud Top Pressure (normalized)'
0205 diagUnits = '0-1 '
ee828acb9e Jean*0206 diagCode = 'SM PC M1 '
0207 diagMate = diagNum
0208 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0209 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
a176fcc394 Andr*0210
0211 diagName = 'CLDMAS '
520b5fccfc Jean*0212 diagTitle = 'Cloud-base Mass Flux (kg/m^2/s)'
0213 diagUnits = 'kg/m^2/s '
893181d3ee Jean*0214 diagCode = 'SM P L1 '
ee828acb9e Jean*0215 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0216 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0217
520b5fccfc Jean*0218 diagName = 'DRAG '
0219 diagTitle = 'Surface Drag Coefficient (kg/m^2/s)'
0220 diagUnits = 'kg/m^2/s '
893181d3ee Jean*0221 diagCode = 'SM P L1 '
ee828acb9e Jean*0222 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0223 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0224
0225 diagName = 'WINDS '
02b51191e4 Jean*0226 diagTitle = 'Surface Wind Speed (m/s)'
520b5fccfc Jean*0227 diagUnits = 'm/s '
0228 diagCode = 'SM P L1 '
ee828acb9e Jean*0229 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0230 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0231
0232 diagName = 'TS '
520b5fccfc Jean*0233 diagTitle = 'near Surface Air Temperature (K)'
02b51191e4 Jean*0234 diagUnits = 'K '
a176fcc394 Andr*0235 diagCode = 'SM L1 '
ee828acb9e Jean*0236 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0237 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0238
0239 diagName = 'QS '
520b5fccfc Jean*0240 diagTitle = 'near Surface Specific Humidity (g/kg)'
0241 diagUnits = 'g/kg '
a176fcc394 Andr*0242 diagCode = 'SM P L1 '
ee828acb9e Jean*0243 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0244 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0245
0246 diagName = 'ENPREC '
520b5fccfc Jean*0247 diagTitle ='Energy flux associated with precip. (snow, rain Temp)'
0248 diagUnits = 'W/m^2 '
a176fcc394 Andr*0249 diagCode = 'SM M1 '
ee828acb9e Jean*0250 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0251 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0252
02b51191e4 Jean*0253 diagName = 'ALBVISDF'
520b5fccfc Jean*0254 diagTitle = 'Surface Albedo (Visible band) (0-1)'
a176fcc394 Andr*0255 diagUnits = '0-1 '
520b5fccfc Jean*0256 diagCode = 'SM P L1 '
ee828acb9e Jean*0257 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0258 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
520b5fccfc Jean*0259
0260 diagName = 'DWNLWG '
0261 diagTitle =
0262 & 'Downward Component of Longwave Flux at the Ground (+=dw)'
0263 diagUnits = 'W/m^2 '
0264 diagCode = 'SM P L1 '
ee828acb9e Jean*0265 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0266 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0267
7f98c35e47 Davi*0268 diagName = 'UPSWG '
0269 diagTitle =
0270 & 'Upward Shortwave Flux at the Ground (+=up)'
0271 diagUnits = 'W/m^2 '
0272 diagCode = 'SM P L1 '
0273 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0274 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0275
a176fcc394 Andr*0276 diagName = 'SWCLR '
520b5fccfc Jean*0277 diagTitle = 'Clear Sky Temp. Tendency due to Shortwave Radiation'
0278 diagUnits = 'K/s '
0279 diagCode = 'SM ML '
ee828acb9e Jean*0280 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0281 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0282
0283 diagName = 'LWCLR '
520b5fccfc Jean*0284 diagTitle = 'Clear Sky Temp. Tendency due to Longwave Radiation'
0285 diagUnits = 'K/s '
a176fcc394 Andr*0286 diagCode = 'SM ML '
ee828acb9e Jean*0287 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0288 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0289
520b5fccfc Jean*0290 diagName = 'TSRCLR '
0291 diagTitle = 'Clear Sky Top-of-atm. net Shortwave Radiation (+=dw)'
0292 diagUnits = 'W/m^2 '
a176fcc394 Andr*0293 diagCode = 'SM P U1 '
ee828acb9e Jean*0294 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0295 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0296
0297 diagName = 'OLRCLR '
520b5fccfc Jean*0298 diagTitle = 'Clear Sky Outgoing Longwave Radiation (+=up)'
0299 diagUnits = 'W/m^2 '
0300 diagCode = 'SM P U1 '
ee828acb9e Jean*0301 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0302 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0303
0304 diagName = 'SWGCLR '
02b51191e4 Jean*0305 diagTitle =
520b5fccfc Jean*0306 & 'Clear Sky Net Shortwave Radiation at the Ground (+=dw)'
0307 diagUnits = 'W/m^2 '
a176fcc394 Andr*0308 diagCode = 'SM P L1 '
ee828acb9e Jean*0309 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0310 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0311
0312 diagName = 'LWGCLR '
520b5fccfc Jean*0313 diagTitle =
0314 & 'Clear Sky Net Longwave Radiation at the Ground (+=up)'
0315 diagUnits = 'W/m^2 '
a176fcc394 Andr*0316 diagCode = 'SM L1 '
ee828acb9e Jean*0317 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0318 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0319
7f98c35e47 Davi*0320 diagName = 'UPSWGclr'
0321 diagTitle =
0322 & 'Upward Shortwave Flux at the Ground (+=up)'
0323 diagUnits = 'W/m^2 '
0324 diagCode = 'SM P L1 '
0325 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0326 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0327
a176fcc394 Andr*0328 diagName = 'UFLUX '
02b51191e4 Jean*0329 diagTitle = 'Zonal Wind Surface Stress (N/m^2)'
520b5fccfc Jean*0330 diagUnits = 'N/m^2 '
ee828acb9e Jean*0331 diagCode = 'UU L1 '
0332 diagMate = diagNum + 2
0333 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0334 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
a176fcc394 Andr*0335
0336 diagName = 'VFLUX '
02b51191e4 Jean*0337 diagTitle = 'Meridional Wind Surface Stress (N/m^2)'
520b5fccfc Jean*0338 diagUnits = 'N/m^2 '
ee828acb9e Jean*0339 diagCode = 'VV L1 '
0340 diagMate = diagNum
0341 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0342 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
1bc0e5d60a Davi*0343
0344 diagName = 'dKE_Ustr'
ee828acb9e Jean*0345 diagTitle =
1bc0e5d60a Davi*0346 & 'KE sink due to Zonal Surf. Stress and Strat. drag (<0)'
0347 diagUnits = 'W/m^2 '
ee828acb9e Jean*0348 diagCode = 'SU L1 '
0349 diagMate = diagNum + 2
0350 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0351 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
1bc0e5d60a Davi*0352
0353 diagName = 'dKE_Vstr'
ee828acb9e Jean*0354 diagTitle =
0355 & 'KE sink due to Merid. Surf. Stress and Strat. drag (<0)'
1bc0e5d60a Davi*0356 diagUnits = 'W/m^2 '
ee828acb9e Jean*0357 diagCode = 'SV L1 '
0358 diagMate = diagNum
0359 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0360 I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
a176fcc394 Andr*0361
0362 diagName = 'DTSIMPL '
02b51191e4 Jean*0363 diagTitle = 'Surf. Temp Change after 1 implicit time step'
520b5fccfc Jean*0364 diagUnits = 'K '
0365 diagCode = 'SM P L1 '
ee828acb9e Jean*0366 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0367 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a176fcc394 Andr*0368
4c105e462e Davi*0369 #ifdef ALLOW_AIM_CO2
0370 diagName = 'aim_pCO2'
0371 diagTitle = 'Atmospheric CO2'
0372 diagUnits = 'ppm '
0373 diagCode = 'SM P L1 '
0374 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0375 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
0376 #endif
0377
0378
7f98c35e47 Davi*0379 diagName = 'FSOL '
0380 diagTitle = 'Incoming Solar radiation'
0381 diagUnits = 'W/m^2 '
0382 diagCode = 'SM P L1 '
0383 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
0384 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
4c105e462e Davi*0385
0386
a176fcc394 Andr*0387 #endif /* ALLOW_DIAGNOSTICS */
0388
0389 RETURN
0390 END