*&^%=USER ALIAS FILE @███████████████████████████████████████████████████████████████████ @███████████████████████████████████████████████████████████████████ @ ALIAS.USR USED AT THE WPC INTERNATIONAL DESKS @ This version of the ALIAS is exclusively adapted to run the @ operational Wingridds Macro website generated in March 2019 @ at https://www.wpc.ncep.noaa.gov/international/wng/ @ Last updated: 22 December 2022 @ By: Jose Manuel Galvez @███████████████████████████████████████████████████████████████████ @███████████████████████████████████████████████████████████████████ @ FRONTAL IDENTIFICATION WITH ABS MOISTURE PARAMETERS WBT1092=SAVR WBTC 1000 WBTC 925 MIX1092=SAVR MIXR 1000 MIXR 925 DWP1092=SAVR DWPT 1000 DWPT 925 @█████████████████████████████████████████████████████████ @█ ONDAS TROPICALES (MACRO OTRP) ████████████████ 20221222 @█████████████████████████████████████████████████████████ %▒▒▒▒▒ VON1: Viento promedio 850-700-600 hPa VON1=VDVC 3 VSUM BKNT 850 VSUM BKNT 700 BKNT 600 %▒▒▒▒▒ VRO1: Vort. Rel. promedio 850-700-600-500 hPa por 10^6. VRO1=SMLC 1000 SMLC 1000 SDVC 4 SSUM RVRT WIND 500 SSUM RVRT WIND 850 SSUM RVRT WIND 700 RVRT WIND 600 %▒▒▒▒▒ OTRCOMA: Moisture flux convg. 925-850-700 hPa OTRCOMA=SMLC 100 SMLC 1000 SMLC 1000 SDVC 3 SSUM COMA 925 SSUM COMA 850 COMA 700 %▒▒▒▒▒ OTRRELH: Rel. Humidity averaged over 850-600 hPa x 10 OTRRELH=SMLC 10 SDVC 3 SSUM RELH 600 SSUM RELH 850 RELH 700 @█████████████████████████████████████████████████████████ @████████████████████████████████████████████████████████████████ @█ TROPW: WPC INTL DESKS TROPICAL WAVE TRACKING ALGORITHM ██ @█ Developed by Clarke and Galvez, based on the Tropical Wave ██ @█ Tracking algorithm developed after a meeting with Andrew ██ @█ Levine from TAFB (NHC) on September 2023 ▒▒▒ 20240708██ @████████████████████████████████████████████████████████████████ %▒▒▒▒▒ Potential vorticity averaged over the 315 and 320 isentrope %SMLT LATMULT is to multiply times -1 for the South Hemisphere TWPV2=SMLT LATMULT SDVC 2 SSUM VPOT I320 VPOT I315 %▒▒▒▒▒ Curvature vorticity CURVOR=SSUM SMLC 0.9 SAVR WSHD 700 WSHD 850 SMLC 1.1 SAVR RVRT WIND 700 RVRT WIND 850 %▒▒▒▒▒ Relative vort (VOREL), shear vort (VOSHR) and curvature vort (VOCUR) all times 10^-6 %VOREL=SMLC 1000 SMLC 1000 SAVR RVRT WIND 700 RVRT WIND 850 %VOSHR=SMLC 1000 SMLC 1000 SAVR WSHD 700 WSHD 850 VOREL=SMLC 1000 SMLC 1000 SDVC 3 SSUM RVRT WIND 600 SSUM RVRT WIND 700 RVRT WIND 850 VOSHR=SMLC 1000 SMLC 1000 SDVC 3 SSUM WSHD 600 SSUM WSHD 700 WSHD 850 VOCUR=SSUM VOSHR VOREL %▒▒▒▒▒ Average Winds 850-700 VAVR78=VAVR BKNT 700 BKNT 850 %▒▒▒▒▒ Latitude correction factor for Pot Vorticity LATMULT=SADC 1 ZPOS SADC -2 ZNEG ALAT %▒▒▒▒▒ Curvature vorticity corrected for latitude VOCURLA=SMLT LATMULT VOCUR @████████████████████████████████████████████████ END OF TROPW █ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@ TRADE WIND INVERSION HEIGHT/STRENGTH - START @@@@@@@@20180315 @@@@@@@@@@@@@@@@@@@ MAR.13.2017 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@ This macro aids with the discrimination of the height and strength @@ of the trade wind inversion. The objective is to help with the @@ assessment of the depth of the moist layer, and with the type of @@ convective regime that is expected. @@-------------------------------------------------------------------- @@Part 1. INVERSION LEVEL DISCRIMINATORS BY COMPUTING @@ DT/DZ IN 50-HPA LAYER @@-------------------------------------------------------------------- IT70=ZNEG SADC 2.5 SDIF TEMP 650 TEMP 700 IC70=SMLC 675 SADC 1 ZPOS SSBC 1 IT70 IT75=ZNEG SADC 2.5 SDIF TEMP 700 TEMP 750 IC75=SMLC 725 SADC 1 ZPOS SSBC 1 IT75 IT80=ZNEG SADC 2.5 SDIF TEMP 750 TEMP 800 IC80=SMLC 780 SADC 1 ZPOS SSBC 1 IT80 IT85=ZNEG SADC 2.5 SDIF TEMP 800 TEMP 850 IC85=SMLC 840 SADC 1 ZPOS SSBC 1 IT85 IT90=ZNEG SADC 2.5 SDIF TEMP 850 TEMP 900 IC90=SMLC 895 SADC 1 ZPOS SSBC 1 IT90 IT95=ZNEG SADC 2.5 SDIF TEMP 900 TEMP 950 IC95=SMLC 945 SADC 1 ZPOS SSBC 1 IT95 IT10=ZNEG SADC 2.5 SDIF TEMP 925 TEMP 975 IC10=SMLC 975 SADC 1 ZPOS SSBC 1 IT10 @@-------------------------------------------------------------------- @@ Part 2. INVERSION LEVEL @@-------------------------------------------------------------------- @ILEV=SMAX IC70 SMAX IC75 SMAX IC80 SMAX IC85 SMAX IC90 IC10 @WINDS WI78=VAVR WIND 700 WIND 600 WI56=VAVR WIND 500 WIND 550 WIML=VAVR WW76 WW54 @@Part 1. INVERSION LEVEL DISCRIMINATORS BY COMPUTING @@ DT/DZ IN 50-HPA LAYER @@-------------------------------------------------------------------- IDISC=3 NIT70=ZNEG SADC IDISC SDIF TEMP 650 TEMP 700 NIC70=SMLC 675 SADC 1 ZPOS SSBC 1 IT70 NIT75=ZNEG SADC IDISC SDIF TEMP 700 TEMP 750 NIC75=SMLC 725 SADC 1 ZPOS SSBC 1 IT75 NIT80=ZNEG SADC IDISC SDIF TEMP 750 TEMP 800 NIC80=SMLC 780 SADC 1 ZPOS SSBC 1 IT80 NIT85=ZNEG SADC IDISC SDIF TEMP 800 TEMP 850 NIC85=SMLC 840 SADC 1 ZPOS SSBC 1 IT85 NIT90=ZNEG SADC IDISC SDIF TEMP 850 TEMP 900 NIC90=SMLC 895 SADC 1 ZPOS SSBC 1 IT90 NIT95=ZNEG SADC IDISC SDIF TEMP 900 TEMP 950 NIC95=SMLC 945 SADC 1 ZPOS SSBC 1 IT95 NIT10=ZNEG SADC IDISC SDIF TEMP 925 TEMP 975 NIC10=SMLC 975 SADC 1 ZPOS SSBC 1 IT10 @@-------------------------------------------------------------------- @@ Part 2. INVERSION LEVEL @@-------------------------------------------------------------------- ILEV=SMAX NIC70 SMAX NIC75 SMAX NIC80 SMAX NIC85 SMAX NIC90 NIC10 @@-------------------------------------------------------------------- @@ Part 3. RELATIVE HUMIDITY MID-LEVELS @@-------------------------------------------------------------------- RELMID=SMTH SMTH SMTH SMTH SMTH SDVC 3 SSUM RELH 70 SSUM RELH 600 RELH 500 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@ TRADE WIND INVERSION HEIGHT/STRENGTH - END @@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@ CENTRAL ANDEAN PRECIPITATION ACCUMULATION INDEX @@@@@@@@@@ @@@@@@@@ CAPI @@@@@@@@@@ @@@@@@@@ 24 February 2020 @@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@BEGIN %%% CAPI is developed in February 2020 to enhance the forecasting of %%% precipitation accumulation in the central Andes. An index CAPRE, %%% the main calculation, considers 7 factors that are summed %%% algebraically: %%% 1.CAMIXR - Mixing Ratio 600-400 hPa. %%% 2.CARELH - Relative humidity 700-300 hPa. %%% 3.CAPWAT - Precipitable water correction. %%% 4.CALIFT - Instability. %%% 5.CACOMA - Moisture convergence 700-500 hPa. %%% 6.CADDUU - Upper Divergence 300-200 hPa. %%% 7.CASHER - Shear 500-300 hPa. %%% Once the index is calculated, the final index includes a GDI %%% correction since the GDI works better for lower elevations of %%% the Andes. The balance of CAPI is ~2/3 CAPRE and 1/3 GDI. %%% %%%================================================================ %%% 1.MIXING RATIO (CAMIXR) - Moisture Variable %%% Role: Moisture availability for precipitation. %%% Layer: 600-400 hPa layer seems to adjust best. %%% Smoothing: No. %%% Capping: No. %%%================================================================ CAMIXS=SSUM MIXR 400 SSUM MIXR 500 SSUM MIXR 600 SSUM MIXR 550 MIXR 700 CAMIXR=SMLC 6+3 SDVC 5 CAMIXS %%% %%%================================================================ %%% 2.RELATIVE HUMIDITY (CARELH) - Moisture Variable %%% Role: Relates to environment ready to saturate rapidly, but %%% can also be a proxy for model-resolved convection. Low %%% values are proxys for dry air entrainment and in some %%% instances elevated inversions. %%% Layer: 700-300 hPa layer seems to adjust best. %%% Smoothing: No %%% Capping: Yes. At a maximum of 120 in CAPI units %%%================================================================ CARHIN=SSBC 65 SDVC 5 SSUM RELH 300 SSUM RELH 400 SSUM RELH 500 SSUM RELH 600 RELH 700 CARELH=SADC 120 ZPOS SSBC 120 SMLC 4 CARHIN %%% %%%================================================================ %%% 3.MOISTURE: PRECIPITABLE WATER (CAPWAT) %%% Role: Corrects for excessive values over elevated terrain in %%% the Altiplano, and low values over the mountain slopes. %%% Defines a terrain mask first (CAMSKT) to correct in %%% model terrain >1100 mASL and avoid producing huge %%% values in areas adjacent to the Andes. %%% Smoothing: 1 (Not more to preserve fine structures) %%% Capping: Yes. %%%================================================================ CAMSKT=SADC 1 ZPOS SSBC 1 ZNEG SSBC 1100 HGHT SFC CAPWAT=SMTH SMLT CAMSKT SSBC 260 SMLC 100 PWAT %%% %%%================================================================ %%% 4.INSTABILITY: LIFTED (CALIFT) %%% Role: Instability for convection. %%% Smoothing: No. %%% Capping: Yes, at -4 because values are generally recurrently %%% excessive in higher parts of the altiplano. %%%================================================================ CALIFT=SMLC 21 SNEG SSBC 4 ZNEG SADC 4 LIFT %%% %%%================================================================ %%% 5.TRIGGER: MOISTURE CONVERGENCE (CACOMA) %%% Role: Surface-based moist ascent/convection trigger. %%% Layer: 700-500 hPa layer seems to adjust best. %%% Smoothing: 1 (Not more to preserve fine structures) %%% Capping: No %%%================================================================ CACOMA=SMTH SNEG SMLC 6+8 SDVC 5 SSUM COMA 650 SSUM COMA 700 SSUM COMA 500 SSUM COMA 550 COMA 600 %%% %%%================================================================ %%% 6.TRIGGER: UPPER DIVERGENCE (CADDUU) %%% Role: Enhance ascent in the column/proxy for model convection %%% Layer: 200-300 hPa layer seems to adjust best. %%% Smoothing: 3 (Noisy field and conv. param. dependent) %%% Capping: Yes. To a max of 120. This is to avoid excessive %%% values during convective feedback situations. %%%================================================================ CADDU1=SDVC 3 SSUM DVRG WIND 200 SSUM DVRG WIND 250 DVRG WIND 300 CADDU2=SMTH SMTH SMTH SMLC 5+6 CADDU1 CADDUU=SADC 120 ZPOS SSBC 120 CADDU2 %%% %%%================================================================ %%% 7.SHEAR (CASHER) %%% Role: Vertical structure of storm cells. Large shear can %%% lead to some severity and locally heavy rain, but %%% generally decreases rainfall totals over a wider area. %%% Layer: 500-300 hPa layer seems to adjust best. %%% Smoothing: 1 (Not more to preserve fine structures) %%% Capping: Yes. %%%================================================================ CASHE1=MAGN VDIF BKNT 300 BKNT 500 CASHE2=SMLC 5 SADC 10 SMTH SNEG CASHE1 CASHER=SSBC 100 ZNEG SADC 100 CASHE2 %%% %%%================================================================ %%% CAPI: FINAL CALCULATION OF THE INDEX (Feb.24.2020) %%% CAPI1: SUM OF ALL THE FACTORS %%% CAPREF: Final Index (2/3 CAPI, 1/3 GDI) %%%================================================================ CAPI1=SSUM CAPWAT SSUM CALIFT SSUM CADDUU SSUM CACOMA SSUM CASHER SSUM CAMIXR CARELH CAPREF=SSUM SMLC 0.45 GDIF SMLC 0.95 CAPI1 %%% %%%================================================================ %%% WINDSPEED FOR CAPI %%% Average of the 600-300 hPa layer. Informs about the origin of %%% the air mass and the movement of the storm cells. %%%================================================================ CAWIND=VDVC 5 VSUM BKNT 300 VSUM BKNT 400 VSUM BKNT 500 VSUM BKNT 550 BKNT 600 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@ CENTRAL ANDEAN PRECIPITATION ACCUMULATION INDEX @@@@@@@@@@ @@@@@@@@ CAPI @@@@@@@@@@ @@@@@@@@ 24 February 2020 @@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@END @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@ FRONT ANALYSIS @@@@@@@@@@ @@@@@@@@ (FRONT) @@@@@@@@@@ @@@@@@@@ 2 March 2020 @@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@BEGIN %%% (1) MASKS FOR TERRAIN ELEVATION %%% Masks are first developed for terrain of different elevations %%% -F_ML allows plots for terrain under 1000 mASL %%% -F_MH allows plots for terrain over 1000 mASL F_ML=SNEG ZPOS SSBC 1 ZNEG SSBC 1000 HGHT SFC F_MH=SADC 1 ZPOS SSBC 1 ZNEG SSBC 1000 HGHT SFC %%%================================================================ %%% (2) THICKNESS COMPONENT (F_THC) - Thermal Aspect %%% -Layer Low-levels: 900-1000 hPa. %%% -Layer Upper-Levels: 925-700 hPa. %%% -Thickness are multiplied to enhance gradients. %%% -Constants are calculated empirically to match the %%% field structure using observations and analyses. %%% -Higher terrain thickness are smoothed twice to %%% reduce 'noise' associated with surface roughness. %%%================================================================ %LOW LEVELS F_THL1=SSBC 750 SDIF HGHT 900 HGHT 1000 F_THL2=SMLC 1-2 SMLT F_THL1 F_THL1 F_THL3=SMLT F_ML F_THL2 %ELEVATED TERRAIN F_THH1=SSBC 1870 SDIF HGHT 700 HGHT 925 F_THH2=SMLC 1-3 SMLT F_THH1 F_THH1 F_THH3=SMLT F_MH SMTH SMTH F_THH2 %COMBINATION F_THC=SSUM F_THL3 F_THH3 %%%================================================================ %%% (3) DEWPOINT COMPONENT (F_DWP) - Moisture Aspect %%% -Layer Low-levels: 950-1000 hPa. %%% -Layer Upper-Levels: 925-850 hPa. %%% -Dewpoints are multiplied to enhance gradients. %%% -Constants are calculated empirically to match the %%% field structure using observations and analyses. %%% -Higher terrain dewpoints are smoothed twice to %%% reduce 'noise' associated with surface roughness. %%%================================================================ F_DWL1=SSBC 650 SMLC 1-2 SMLT DWPK 1000 DWPK 950 F_DWL2=SMLT F_ML F_DWL1 F_DWH1=SSBC 580 SMLC 1-2 SMLT DWPK 850 DWPK 925 F_DWH2=SMLT F_MH SMTH SMTH SMTH F_DWH1 %COMBINATION F_DWP=SMLC 2 SSUM F_DWL2 F_DWH2 %%%================================================================ %%% (4) AIRMASS FIELD 'ALPHA' (F_ALPP) %%% -Multiplies the dewpoint and and thickness components %%% (F_THC and F_DWP) and caps them at 1100 (dimensionless) %%% to limit gradients in the tropics that associate with %%% orographic features and low dewpoints generated by %%% processes that are not frontal. %%%================================================================ F_ALP=SMLC 1-2 SMLT F_THC F_DWP F_ALPP=SADC 1100 ZPOS SSBC 1100 F_ALP %%%================================================================ %%% (5) ADDITIONAL MASKS %%% -F_THC: Mask for large values of thickness that are %%% to limit identification of fronts in the deep tropics. %%% All the values larger than 290 are masked out when %%% calculating the gradient (F_GAL) later on. %%% -F_OCO: Orographic correction mask to maskout the %%% effects of gradients associated with the terrain. %%%================================================================ F_THCM=SNEG ZPOS SSBC 1 ZNEG SSBC 290 F_THC F_OCO=SMLC 5+4 MGRD HGHT SFC %%%================================================================ %%% (6) MAGNITUDE OF THE GRADIENT (F_GAL) %%% -F_GAL: Calculates the magnitude of the gradient %%% to limit identification of fronts in the deep tropics. %%% -F_OCO: Orographic correction mask to maskout the %%% effects of gradients associated with the terrain. %%%================================================================ F_GAL=SMLT F_THCM SMLC 9+5 MGRD SMTH F_ALPP %%%================================================================ %%% (7) RELATIVE HUMIDITY ADDITION %%% -F_RELH: Average of relative humidity in the lower %%% troposphere. %%% -F_GRRH: Gradient of relative humidity. %%%================================================================ F_RELH=SDVC 4 SSUM RELH 950 SSUM RELH 900 SSUM RELH 800 RELH 700 F_GRRH=SDVC 2.0 SQRT SMLC 8+4 MGRD F_RELH %%%================================================================ %%% (8) FRONT GRADIENT MAGNITUDE FIELD (F_GRAD) %%% -F_GALM: Mask to maskout the gradients that are too %%% low (F_GAL<30) and highlight larger gradients. %%% -F_GRAD: Adds front gradient to relative humidity %%% gradient, substracts the orography-related gradients, %%% and multiplies by the mask F_GALM to hgihlight large %%% gradients only. %%%================================================================ F_GALM=ZNEG SADC 1 ZPOS SSBC 30 F_GAL F_GRAD=SMTH SMLT F_GALM SDIF SSUM F_GAL F_RELH F_OCO @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@ FRONTS PART2 - HIGH LATITUDES @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ %%% New field without elevated terrain component F_ALPG=SMLC 1-2 SMLT SMLC 2 F_THL2 F_DWL1 F_ALPPG=SADC 1100 ZPOS SSBC 1100 F_ALPG %%% Latitude enhancement F_LAT=SADC 1 SDVC 35 ZNEG SSBC 30 ABSV ALAT F_GALG=SMLT F_THCM SMLC 9+5 MGRD SMTH F_ALPPG F_GRADG=SMTH SMLC 0.9 SMTH SMLT F_LAT SMLT F_GALM SSUM F_GALG F_RELH @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@ FRONT ANALYSIS @@@@@@@@@@ @@@@@@@@ (FRONT) @@@@@@@@@@ @@@@@@@@ 2 March 2020 @@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@END @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@ POTENTIAL FOR SEVERE WEATHER AND HAIL @@@@@@@@@@ @@@@@@@@ GR02T @@@@@@@@@@ @@@@@@@@ 06 January 2020 @@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@BEGIN %%% %%% This algorithm has 2 parts. The generation of a mask for strong %%% convection (GMKTT) and population of this mask with enhancers %%% of severity. %%% %%%================================================================ %% PART 1. ENVIRONMENTS SUITABLE FOR STRONG DEEP MOIST CONVECTION %% This is assumes that the following 5 processes are present: %% 1a) Sufficient moisture in the column [PWAT > 20mm] %% 1b) Sufficient deep-layer instability [LI <0C] %% 1c) Dynamically-induced ascent [OMEGA_600-300 < 0Pa/s] %% 1d) Sufficient saturation in mid-levels [RELH_700-500 > 60%] %% 1e) Sufficient mid-level instability/cold air [T_600 <+2C] %%%================================================================ %% 1a) PRECIPITABLE WATER MASK: PWAT > 20mm %% Assumes that at least 20mm of precipitable water %% are required to sustain deep moist convection. GMKPW=SADC 1 ZPOS SSBC 1 ZNEG SSBC 20 SMLC 10 PWAT %%-------------------------------------------------------------- %% 1b) LIFTED INDEX MASK: LIFT < 0C %% Negative Lifted: Instability favorable for deep %% convection, although not necesarilly severe. %%2021Oct GMKLI=ZNEG SADC 1 SNEG ZNEG LIFT GMKLI=ZNEG SADC 1 SNEG ZNEG SSBC 1 LIFT %%-------------------------------------------------------------- %% 1c) OMEGA MASK: OMEGA_600-300 < 0 Pa/s (negative=ascent) %% Looks for negative omegas, but not necesarrily the value. %% Dynamically-induced ascent favors vertical growth of %% convection. Weight is given to the mid-upper troposphere %% due to the important role of upper jets and other upper %% features that enhance upper divergence and ascent in %% the region. Smoothing is applied since OMGA is a noisy %% field in NWP and could be off by tenths of km due to the %% mesoscalar/microscalar nature of ascending updrafts. GM_OMGS=SMTH SSUM OMGA 500 SSUM OMGA 600 SSUM OMGA 400 OMGA 300 %%% Changed from SMLC 1+3 to SMLC 1+1 on 12 Dec 2020 - Windows 10 update. %%2021Oct GM_OMGA=SMLC 1+1 SDVC 4 GM_OMGS GM_OMGA=SMLC 10 SDVC 4 GM_OMGS GMKOM=SNEG ZPOS SMTH SMTH SMTH GM_OMGA %%-------------------------------------------------------------- %% (1d) RELATIVE HUMIDITY MASK: RELH_700-500 > 60% %% Saturation favors moist updrafts which are important %% for the growth of water droplets and hail stones. GRHML=SMTH SDVC 4 SSUM RELH 850 SSUM RELH 700 SSUM RELH 600 RELH 500 %%2021Oct GMKRH=ZNEG SSBC 60 GRHML GMKRH=ZNEG SSBC 50 GRHML %%-------------------------------------------------------------- %% (1e) MID-LEVEL TEMPERATURE MASK: T_600 <+2C %% Indicates that air is cool enough in the mid-troposphere %% to sustain strong convection and hail stones that could %% reach the ground. GMKTE=SNEG ZPOS SSBC 2 TEMP 600 CIN1 GMKTE1=SADC 1 ZPOS SSBC 1 GMKTE %%------------------------------------------------------------------ %% MASK FOR THE POTENTIAL OF STRONG DEEP MOIST CONVECTION %% The mask (variable GMKTT) contains 1 where the five conditions %% calculated above are present. This means that the environment %% is suitable for strong deep moist convection. Otherwise it %% contains zeros. Note that applies a terrain mask MSKT. %%------------------------------------------------------------------ GMKTTA=SMLT GMKOM SMLT GMKPW SMLT GMKTE SMLT GMKLI GMKRH GMKTT=SMLC 10 SADC .1 ZPOS SSBC .1 SMLT MSKT GMKTTA %% %%%================================================================ %% PART 2. ENHANCERS OF SEVERITY AND THE POTENTIAL FOR HAIL %% Values are added to the mask GMKTT when any of the %% following conditions are present: %% 2a) Enhanced 700-500 hPa Lapse Rate [LR700-500>15C] %% 2b) Cool mid-level temperatures. [T_500hPa<-8C] %% 2c) Negative Lifted Index %% 2d) Dynamically-induced ascent %% 2e) Upper divergence enhancement %% 2f) Enhanced moisture convergence [COMA_950-700<-0.5] %% 2g) Enhanced shear %%%================================================================ %% 2a) 700-500 Lapse Rate. Larger values enhance violent ascent %% in the region of hail formation, favoring hail growth. GM_LR75=SMTH SDIF TEMP 700 TEMP 500 GM_MKLR=ZNEG SMLC 7 SSBC 15 GM_LR75 MSKMLKR=SADC 1 ZPOS SADC -1 GM_MKLR %Before 2020Mar11 GM_MKLR=ZNEG SMLC 4 SSBC 15 GM_LR75 %Before 2021Oct GM_MKLR=ZNEG SMLC 7 SSBC 15 GM_LR75 %%-------------------------------------------------------------- %% (2b) Temperature 500hPa. Cooler values enhance instability and %% limit melting of hail stone allowing for larger stones. %%2021Oct GM_T500=SMLC -2 SADC 8 TEMP 500 GM_T500=SMLC -2 SADC 7 TEMP 500 %%-------------------------------------------------------------- %% (2c) Deep-layer Instability: LIFTED < -3 %GM_LI=SMLC 4 ZNEG SNEG LIFT %%2021Oct GM_LI=SMLC 6 ZNEG SNEG LIFT GM_LI=SMLC 7 ZNEG SNEG SSBC 1 LIFT %%-------------------------------------------------------------- %% (2d) Stronger dynamically-induced ascent (Omega < ) GMKOM2=SMLC 3 SNEG ZPOS SADC 2 SMTH GM_OMGA %%2021Oct GMKOM2=SNEG ZPOS SADC 4 SMTH SMTH GM_OMGA %%-------------------------------------------------------------- %% (2f) Upper Divergence Enhancer %The commented line was valid until 2020Mar11 - Added more weight %to upper divergence on March 12, based on California Severe %thunderstorms event and on comments from Nestor Santayana %G_DDUU=SMLC 4+4 SMTH SMTH SMTH SSUM DVRG WIND 300 SSUM DVRG WIND 250 DVRG WIND 200 %%2021Oct G_DDUU=SMLC 1+5 SMTH SMTH SMTH SSUM DVRG WIND 300 SSUM DVRG WIND 250 DVRG WIND 200 G_DDUU=SMLC 4000 SMLC 100 SMTH SMTH SDVC 3 SSUM DVRG WIND 300 SSUM DVRG WIND 250 DVRG WIND 200 %%2021Oct GM_DDU=ZNEG SSBC 1 G_DDUU GM_DDU=ZNEG G_DDUU %%-------------------------------------------------------------- %% (2c) Moisture convergence. GM_COMA=SMLC 1000 SMLC 1000 SMLC 100 SDVC 5 SSUM COMA 700 SSUM COMA 800 SSUM COMA 850 SSUM COMA 900 COMA 950 %%2021Oct GM_MKCO=ZNEG SMLC -0.5 SMTH GM_COMA GM_MKCO=ZNEG SMLC -3 SMTH GM_COMA %%-------------------------------------------------------------- %% (2e) Shear/cizalla. %%2021Oct GM_CIZA=SAVR MAGN VDIF BKNT 500 BKNT 1000 MAGN VDIF BKNT 700 BKNT 1000 %%2021Oct GM_MKCZ=SMLC .1 ZNEG SMTH SSBC 40 GM_CIZA GM_CIZA=SAVR MAGN VDIF BKNT 500 BKNT 1000 MAGN VDIF BKNT 700 BKNT 1000 GM_MKCZ=SMLC 3 ZNEG SSBC 25 GM_CIZA %% %%================================================================= %% PART 3. POTENTIAL FOR SEVERITY (GRF1) %%================================================================= GRF1A=SSUM GM_LI SSUM GMKOM2 SSUM GM_MKCO SSUM GM_MKLR GM_T500 GRF1=SMLT GMKTT SSUM GM_DDU SSUM GM_MKCZ GRF1A %%================================================================= %% PART 4. WIND IN THE HAIL FALLING REGION (GRF1) %%================================================================= % WIND IN THE HAIL FALLING REGION %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% GM_HWND=VDVC 3 VSUM BKNT 800 VSUM BKNT 700 BKNT 600 GM_ADGR=SMLC 2+3 ADVT GRF1 GM_HWND GRF2=ZNEG SSUM GRF1 GM_ADGR @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@ MACRO DE GRANIZO - GR03 (09 DECEMBER 2019) - END @@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@ TEST FOR A HAIL MACRO THAT CONSIDERS SHEAR AND TEMPERATURES - BEGIN @@@@ 06 DECEMBER 2019 - @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ PROB=VDIF BKNT 500 BKNT 1000 CIZALLA06=MAGN VDIF BKNT 500 BKNT 1000 CIZALLA03=MAGN VDIF BKNT 700 BKNT 1000 @@ Mascara cizalla > 20 y > 40 GR_MKCIZ20=SADC 1 ZPOS SSBC 1 ZNEG SADC -20 SMTH CIZALLA03 GR_MKCIZ40=SADC 1 ZPOS SSBC 1 ZNEG SADC -40 SMTH CIZALLA03 @@ Mascara cizalla total GR_MKCIZTOTAL=SSUM GR_MKCIZ20 GR_MKCIZ40 @@ Mascara PWAT GR_MKPW20=SADC 1 ZPOS SSBC 1 ZNEG SADC -20 SMLC 10 PWAT GR_MKPW40=SADC 1 ZPOS SSBC 1 ZNEG SADC -40 SMLC 10 PWAT @@ Mascara cizalla total GR_MKPWTOTAL=SSUM GR_MKPW20 GR_MKPW40 @@ Mascara TEMPERATURA @----TEMP 600<2C GR_MKT600=ZNEG SADC 1 ZPOS SNEG SADC -2 TEMP 600 GR_MKT500=ZNEG SADC 1 ZPOS SNEG SADC 6 TEMP 500 GR_MKT5002=ZNEG SADC 1 ZPOS SNEG SADC 9 TEMP 500 GR_MKTFN=SSUM GR_MKT600 SSUM GR_MKT500 GR_MKT5002 @@ Mascara LAPSE RATE GR_MKLR16=ZNEG SADC 1 ZPOS SADC -16 SDIF TEMP 700 TEMP 500 GR_MKLR19=ZNEG SADC 1 ZPOS SADC -19 SDIF TEMP 700 TEMP 500 GR_MKLRFN=SSUM GR_MKLR19 GR_MKLR16 @---------------------------MASKOUTS @@ Mascara OMEGA GR_OMGA=SMTH SMLC 1+4 SDVC 4 SSUM OMGA 700 SSUM OMGA 500 SSUM OMGA 400 OMGA 300 GR_MKOMGA=SNEG SSBC 1 ZNEG SADC 1 ZPOS SADC 5 GR_OMGA @@ Mascara humedad relativa troposfera media GR_HUMETM=SDVC 3 SSUM RELH 700 SSUM RELH 600 RELH 500 GR_MKHTM=ZNEG SADC 1 ZPOS SADC -40 GR_HUMETM @@ Mascara LIFTED GR_MKLF-1=SNEG SSBC 1 ZNEG SADC 1 ZPOS SADC 1 LIFT GR_MKLF-2=SNEG SSBC 1 ZNEG SADC 1 ZPOS SADC 2 LIFT GR_MKLF-4=SNEG SSBC 1 ZNEG SADC 1 ZPOS SADC 4 LIFT GR_MKLF-6=SNEG SSBC 1 ZNEG SADC 1 ZPOS SADC 6 LIFT GR_MKLFFN=SSUM GR_MKLF-2 SSUM GR_MKLF-4 GR_MKLF-6 @@ Multiplier Masks GR_MASKOUT=SMLT GR_MKOMGA SMLT GR_MKLF-1 GR_MKHTM @@ Mascara Final GR_MKFN=SMLT GR_MASKOUT SSUM GR_MKLRFN SSUM GR_MKTFN SSUM GR_MKLFFN SSUM GR_MKPWTOTAL GR_MKCIZTOTAL @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@ TEST FOR A HAIL MACRO THAT CONSIDERS SHEAR AND TEMPERATURES - END @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@ VARIABLES FOR HAIL MACROS (GR02 and GR03) @@@@@@@@@@@@@@ @@@@@@@@@ VARIABLES PARA MACROS DE GRANIZO (GR02 y GR03) @@@@@@@@@@@@@@ @@@@@@@@@ Galvez & Santayana, 12-Sep-2019 @@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ %%%%% These are the variables defined to be able to operate the macros %%%%% GR02 and GR03. These are designed to detect the potential for hail %%%%% in South America. The algorithm was initially developed by Nestor %%%%% Santayana (INUMET) in 2012-13, and adjusted by Jose Manuel Galvez %%%%% in 2014. The macro used from 2014 through June 2019 was GR01. But %%%%% due to changes on GFS variable VVEL, the macro GR01 does not work %%%%% anylonger. The recommended macros are GR02 and GR03. %%%%% %%%%%==================================================================== %%%%% TERRAIN MASK (MSKT) %%%%% Blocks areas where the surface is higher than 750 mASL. %%%%% -MSKT: Contains ones in gridpoints where HGHT SFC < 750 mASL. %%%%% Otherwise MSKT contains zeros. %%%%%==================================================================== MSKT=SADC 1 ZPOS SSBC 1 SMLC -1 ZPOS SSBC 750 HGHT SFC %%%%% %%%%%==================================================================== %%%%% PRECIPITABLE WATER MASK (MSKP) %%%%% Assumes sufficient moisture available for hail when PWAT>20mm. %%%%% -MSKP: Contains ones in gridpoints where PWAT > 20mm. %%%%% Otherwise MSKP contains zeros. %%%%%==================================================================== MSKP=SADC 1 ZPOS SSBC 1 ZNEG SSBC 20 SMLC 10 PWAT %%%%% %%%%%==================================================================== %%%%% MID-TROPOSPHERIC TEMPERATURE MASK (MSKN) %%%%% Assumes that 500hPa temperatures < -7C are sufficiently cold %%%%% to sustain instability and cold air favorable for hail. %%%%% -MSKN: Contains ones in gridpoints where Temp_500 < -7C. %%%%% Otherwise MSKN contains zeros. %%%%%==================================================================== MSKN=ZNEG SADC 1 SNEG ZNEG SADC 7 TEMP 500 GMT1=ZNEG SADC 1 SNEG ZNEG SADC 6.5 TEMP 500 %%%%% %%%%%==================================================================== %%%%% 700-500 LAPSE RATE (LMTP) AND LAPSE RATE MASK (MSKR) %%%%% Assumes that lapse rates >16C in this layer favor updrafts that %%%%% are sufficiently strong to favor hail growth. %%%%% -LMTP: 500-700 hPa Lapse Rate or Temperature Difference. %%%%% -MSKR: Contains ones in gridpoints where LMTP > 16C. %%%%% Otherwise MSKR contains zeros. %%%%%==================================================================== LMTP=SDIF TEMP 700 TEMP 500 MSKR=SADC 1 ZPOS SSBC 1 ZNEG SADC -16 LMTP GMR1=SADC 1 ZPOS SSBC 1 ZNEG SSBC 15.5 LMTP %%%%% %%%%%==================================================================== %%%%% AVERAGED VERTICAL VELOCITY (VVTS) %%%%% Vertical velocity in pascals/s, omega (OMGA), is averaged %%%%% over the 300-700hPa layer and stored in the variable VVTS. %%%%% This is just a simple scalar average. It is multiplied by %%%%% 10^4 SMLC 1+4) and smoothed once (SMTH) for 1DEG GFS. The %%%%% variable OMGA was the variable VVEL before 12 June 2019. %%%%%==================================================================== SOMG=SSUM OMGA 700 SSUM OMGA 500 SSUM OMGA 400 OMGA 300 VVTS=SMTH SMLC 10000 SDVC 4 SOMG %%%%% %%%%% %%%%%==================================================================== %%%%% LIFTED INDEX MASKS (GML1, 2 AND 3) %%%%% They contain ones where LIFT < -0.5, -3 and -6. They contain %%%%% zeros otherwise. The more negative LIFT, the higher chance for %%%%% severity. %%%%%==================================================================== GML1=SADC 1 SNEG ZNEG SADC 1 ZPOS SADC -1 LIFT GML2=SADC 1 SNEG ZNEG SADC 1 ZPOS SADC 3 LIFT GML3=SADC 1 SNEG ZNEG SADC 1 ZPOS SADC 6 LIFT %%%%%==================================================================== %%%%% VERTICAL VELOCITY MASKS (MVT1, MVT2 and MVT3) %%%%% The following are masks of vertical velocity. %%%%% -MVT1 contains ones in regions where VVTS > 10^-4 Pa/s %%%%% -MVT2 contains ones in regions where VVTS > 20^-4 Pa/s %%%%% -MVT2 contains ones in regions where VVTS > 30^-4 Pa/s %%%%% -Otherwise they contain zeros. %%%%%==================================================================== MVT0=SNEG SSBC 1 ZNEG SADC 1 ZPOS SADC 5 VVTS MVT1=SNEG SSBC 1 ZNEG SADC 1 ZPOS SADC 10 VVTS MVT2=SNEG SSBC 1 ZNEG SADC 1 ZPOS SADC 40 VVTS MVT3=SNEG SSBC 1 ZNEG SADC 1 ZPOS SADC 70 VVTS %%%%% %%%%%==================================================================== %%%%% FINAL MASKS FOR GR02 %%%%% They contains ones when the following intercept: %%%%% -MSKT: Terrain below 750 mASL. %%%%% -MSKP: Precipitable Water > 20mm. %%%%% -GMF1: Lifted<-0.5C (GML1) and Omegas<-10^-3 Pa/s (MVT1). %%%%% Used to allow the contouring of moisture convergence, %%%%% 500hPa mixing ratios and the plotting of green boxes. %%%%% Indicates areas prone for strong thunderstorms. %%%%% -GMF2: Lifted<-3C (GML2) and Omegas<-20^-3 Pa/s (MVT3). %%%%% Used to plot the red contour that shows a risk for strong %%%%% updrafts. Implies increased risk for severity and hail. %%%%% -GMF3: Lifted<-6C (GML3) and Omegas<-20^-3 Pa/s (MVT3). %%%%% Used to plot the fucsia contour that shows a high risk for %%%%% strong updrafts and associated enhanced potential for hail. %%%%% -GMF4: Lifted<-3C (GML3) and Omegas<-15^-3 Pa/s (MVT2). %%%%% Used to allow the plotting of colored boxes inside the %%%%% red colored contours (moderate risk for hail). %%%%% -GMF5: Lifted<-6C (GML3) and Omegas<-15^-3 Pa/s (MVT2). %%%%% Used to allow the plotting of colored boxes inside the %%%%% fuscia colored contours (high risk for hail). %%%%%==================================================================== GMF1=SMLT MSKMLKR SMLT GMKTE1 SMLT MSKT SMLT MSKP SMLT MVT1 GML1 GMF2=SMLT MSKMLKR SMLT GMKTE1 SMLT MSKT SMLT MSKP SMLT MVT3 GML2 GMF3=SMLT MSKMLKR SMLT GMKTE1 SMLT MSKT SMLT MSKP SMLT MVT3 GML3 GMF4=SMLT MSKMLKR SMLT GMKTE1 SADC 1 ZPOS SSBC 1 SMLC 10 SMLT MSKT SMLT MSKP SMLT MVT2 GML2 GMF5=SMLT MSKMLKR SMLT GMKTE1 SADC 1 ZPOS SSBC 1 SMLC 10 SMLT MSKT SMLT MSKP SMLT MVT2 GML3 %%%%% %%%%%==================================================================== %%%%% FINAL MASKS FOR GR03 %%%%% They produce ones when the following intercept: %%%%% -MSKT: Terrain below 750 mASL. %%%%% -MSKP: Precipitable Water > 20mm. %%%%% -MSKN: Temp_500<-7C (Unstable & sufficiently cold) %%%%% -MSKR: Lapse Rate_700-500>15C (Unstable=strong updrafts) %%%%% -MVT1, MVT2 and MVT3: Omegas (Ascent) %%%%% -MSKL: Lifted Index<-1.5C (Unstable) %%%%% -Otherwise they contain zeros. %%%%%==================================================================== @@Lapse Rate Mask G3R1=SADC 1 ZPOS SSBC 1 ZNEG SADC -15 LMTP G3M1=SMLT MSKN SMLT G3R1 SMLT MSKT SMLT MSKP SMLT MVT1 GML2 G3M2=SMLT MSKN SMLT G3R1 SMLT MSKT SMLT MSKP SMLT MVT1 GML3 G3M3=SMLT MSKN SMLT G3R1 SMLT MSKT SMLT MSKP SMLT MVT3 GML3 %%%%% %%%%%==================================================================== %%%%% MOISTURE CONVERGENCE IN 950-700 hPa LAYER (MCNV) %%%%% Is not used to generate a mask. But is plotted to describes the %%%%% potential for ascending moist parcels, which feed low-level %%%%% moisture into developing deep convection. Can favor an increase %%%%% in mid-level moisture that enhances hail formation. Moisture %%%%% convergence are negative values of COMA and MCNV. Units of both %%%%% COMA and MCNV: kg/kg * m^-2 * s^-1 * 10^-8 %%%%% -COMA=Divergence of the Flux of Mixing Ratio %%%%% -MCNV=COMA Averaged in the 950-700 hPa layer. %%%%%==================================================================== COMA=DVRG FLUX MIXR WIND COM5=SSUM COMA 950 SSUM COMA 900 SSUM COMA 850 SSUM COMA 800 COMA 700 MCNV=SMLC 1+8 SDVC 5 COM5 %%%%% %%%%%==================================================================== %%%%% UPPER DIVERGENCE IN THE 200-400 hPa LAYER (DDU3) %%%%% Does not generate a mask, but is plotted by GR02 and GR03 as it %%%%% indicates areas where strong dynamically-induced ascent might %%%%% be present. %%%%%==================================================================== %%%%% DDU1: Is the wind divergence averaged over the 200-400 hPa layer %%%%% multiplied by 10^7. It is smoothed to remove noise and %%%%% preserve the most robust signal. Units: 10^7 s^-1. %%%%%-------------------------------------------------------------------- DV22=SAVR DVRG WIND 200 DVRG WIND 250 DV43=SAVR DVRG WIND 400 DVRG WIND 300 DDUU=SMLC 1+7 SAVR DV22 DV43 DDU1=SMTH SMTH SMTH SMTH SMTH DDUU %%%%%-------------------------------------------------------------------- %%%%% Large positive and large negative values are truncated for %%%%% contouring purposes. %%%%%-------------------------------------------------------------------- DDU2=SSBC 050 ZNEG SADC 050 DDU1 DDU3=SADC 320 ZPOS SSBC 320 DDU2 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@ CODE FOR HAIL MACRO VARIABLES (GR02,GR03) ENDS @@@@@@@@@@@@ @@@@@@@@@@@ FIN DEL CODIGO PARA LOS MACROS (GR02,GR03) @@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@Rodrigo Mathias - Potencial de Chuvas fortes para Rio de Janeiro - START @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@Mask for relative humidity >90% at 900 hPa MKRH=ZNEG SADC 1 ZPOS SSBC 90 RELH 900 @@Mask for 10m winds >5m/s MKWN=ZNEG SADC 1 ZPOS SSBC 5 YCMP WIND 10M @@Sum of both masks. Risk arises when MKCJ>=2 MKCJ=SSUM MKRH MKWN @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@Rodrigo Mathias - Potencial de Chuvas fortes para Rio de Janeiro - END @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@ MOISTURE CONVERGENCE @@@@@@@@@@@@@@@@@@@@@@@@@@@ JMGALVEZ 20190528 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@BEGIN %% Computes the divergence of the flux of mixing ratio. Relates to %% the triggering of convection by the convergence of water vapor. %% COLL is the average for the 1000-850 hPa Layer %%-------------------------------------------------------------------- COMA=DVRG FLUX MIXR WIND COLL=SDVC 3 SSUM COMA 850 SSUM COMA 925 COMA 1000 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@ MOISTURE CONVERGENCE @@@@@@@@@@@@@@@@@@@@@@@@@@@ JMGALVEZ 20190528 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@END @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@ ENHANCED GDI @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ JMGALVEZ 20190528 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@BEGIN %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 1.BASE CALCULATIONS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%-------------------------------------------------------------------- %% 1.1.OROGRAPHIC CORRECTION FACTOR %% It is used to tamper values that are ficticiously high over %% elevated terrain. Its a fraction: 0.5 at ~3200m. It is used %% to limit excessive values of moisture convergence generated %% over mountainous terrain (Part 2.2). %%-------------------------------------------------------------------- OCO1=SNEG SSBC 1 SMLC 2 SMLC 1-4 HGHT SFC %%-------------------------------------------------------------------- %% 1.2.MASK GDI>25 %% Masks GDI values that are lower than 25. This allows to later %% select regions that have the potential for deep convection by %% masking regions where convection is generally shallow. With %% the exception of moisture convergence adjustments, the rest of %% the adjustments are applied to large (>25) GDI values only. %%-------------------------------------------------------------------- MGD4=SADC 1 ZPOS SSBC 1 ZNEG SSBC 25 GDIF %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 2.MIXING RATIO FLUX CONVERGENCE CONTRIBUTION (SMCF and SMCH) %% Mixing ratio flux convergence in the 1000-700hPa layer is an %% important factor for tropical convection. It plays a role in the %% triggering of convection but also by contributing to the amount %% of moisture available for convection %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%-------------------------------------------------------------------- %% 2.1.MIXING RATIO FLUX CONVERGENCE IN THE 950-700 HPA LAYER %% Simple average in the 950-700 layer. %%-------------------------------------------------------------------- MCLX=SSUM COMA 950 SSUM COMA 900 SSUM COMA 850 SSUM COMA 800 COMA 700 MCL1=SDVC 5 MCLX %%-------------------------------------------------------------------- %% 2.2.SCALING AND CORRECTION FOR OROGRAPHY %%-------------------------------------------------------------------- MCL2=SMLT OCO1 SMLC 1+9 MCL1 %%-------------------------------------------------------------------- %% 2.3.SERIES TRUNCATION %% Truncates the field when excessive values are present (MCL3). %%-------------------------------------------------------------------- MCL3=SADC 220 ZPOS SSBC 440 ZNEG SADC 220 MCL2 %%-------------------------------------------------------------------- %% 2.4.SMOOTHING (ONLY FOR PLOTTING, NOT FOR FURTHER CALCULATIONS) %%-------------------------------------------------------------------- MCL4=SMTH MCL3 MCL5=SMTH SMTH SMTH SMTH SMTH MCL4 %%-------------------------------------------------------------------- %% 2.5.ADJUSTMENT FACTOR FOR THE GDI %% Conversion to GDI value ranges. %% +Was 0.04 before August 14, 2019 %%-------------------------------------------------------------------- SMCF=SMLC 0.08 SNEG SMTH SMTH MCL3 %%-------------------------------------------------------------------- %% 2.6.ADJUSTMENT FACTOR FOR THE GDI - HIRES DOMAIN %% The weight provided for the high resolution domain is lower %% than the weight used for 1 degree domains, given that %% resolution has a strong impact on the values of convergence. %% +Was 0.02 before August 14, 2019 %%-------------------------------------------------------------------- SMCH=SMLC 0.04 SNEG SMTH SMTH MCL3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 3.REDUCTION BY DRY AIR IN THE MID-UPPER TROPOSPHERE (AJRH) %% Tropical convection is sensitive to dry air entrainment. %% Since the GDI considers up to 500 hPa, this adjustment factor %% uses 300-400 hPa relative humidity. A reduction factor is %% computed for regions where the GDI is higher than 25 (MGD4) %% whenever the average RH in the 400-300 layer is lower than 70% %% (SSBC 70 SAVR RELH 400 RELH 300). The scaling constant is 0.2. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% AJRH=SMLT MGD4 SMLC 0.2 ZPOS SSBC 70 SAVR RELH 400 RELH 300 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 4.REDUCTION BY MID-UPPER CONVERGENT PATTERN (MUC2 and MUC4) %% Convergence in the upper troposphere is used to limit GDI %% values. In the deep tropics, the levels of 250-150 hPa work %% best due to a high tropopause. Elsewhere, the 300-200 hPa %% layer tends to work well. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%-------------------------------------------------------------------- %% MUC2: Reduction under convergent pattern in the 300-200 hPa layer %%-------------------------------------------------------------------- MUC1=SMLC 8+5 SSUM DVRG WIND 200 SSUM DVRG WIND 300 DVRG WIND 250 MUC2=SMTH SMTH SMLT MGD4 ZPOS MUC1 MUCH=SMLC 0.25 smth SMTH SMTH MUC2 %%-------------------------------------------------------------------- %% MUC4: Reduction under convergent pattern in the 250-150 hPa layer %%-------------------------------------------------------------------- MUC3=SMLC 8+5 SSUM DVRG WIND 200 SSUM DVRG WIND 150 DVRG WIND 250 MUC4=SMTH SMTH SMLT MGD4 ZPOS MUC3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 5. ENHANCEMENT BY PRECIPITABLE WATER > 30MM (PWCF) % 20210505 - Added SMLC 0.1 to PWAT, because the units changed in the new GFS file. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% PWCF=SMLC 2 ZNEG SSBC 3 SMLC 0.1 PWAT %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 6. ENHANCED GDI %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%-------------------------------------------------------------------- %% GDIS (Standard adjustment) %% -MUC2=Convergence in the 300-200hPa layer %% -AJRH=Unsaturated air in the 300-400 hPa layer %% -SMCF=Moisture divergence/convergence in the 950-700hPa layer %% -PWCF=Precipitable Water>30mm %%-------------------------------------------------------------------- GDIS=SSUM MUC2 SSUM AJRH SSUM SMCF SSUM PWCF GDIF %%-------------------------------------------------------------------- %% GDIT (Deep tropics variant) %% Considers upper convergent pattern in the %% 150-250 layer (MUC4) versus 300-200 hPa (MUC2). %%-------------------------------------------------------------------- GDID=SSUM MUC4 SSUM AJRH SSUM SMCF SSUM PWCF GDIF %%-------------------------------------------------------------------- %% GDIH (High resolution domain variant) %% -SMCH instead of SMCF: Uses a lower weight for the moisture %% conv. contribution which is highly sensitive to resolution. %% -MUCH instead of MUC2: Decreases the contribution of upper %% convergence which is also highly sensitive to resolution. %%-------------------------------------------------------------------- GDIH=SSUM MUCH SSUM AJRH SSUM SMCH SSUM PWCF GDIF %%-------------------------------------------------------------------- %% GDIT (Standard minus upper convergent adjustment) %% -This is used for the GDI advection mask %% -AJRH=Unsaturated air in the 300-400 hPa layer %% -SMCF=Moisture divergence/convergence in the 950-700hPa layer %% -PWCF=Precipitable Water>30mm %%-------------------------------------------------------------------- GDIT=SSUM AJRH SSUM SMCF SSUM PWCF GDIF %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % 7. GDI ADVECTION %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%-------------------------------------------------------------------- %% GDI MASK FOR VALUES>30 %%-------------------------------------------------------------------- GDM1=ZNEG SADC -30 GDIT GDM2=SMTH SMTH SMTH GDM1 %%-------------------------------------------------------------------- %% STEERING FLOW %%-------------------------------------------------------------------- SFUU=VDVC 4 VSUM BKNT 400 VSUM BKNT 300 VSUM BKNT 250 BKNT 200 SFMM=VDVC 4 VSUM BKNT 850 VSUM BKNT 700 VSUM BKNT 600 BKNT 500 SFL1=VAVR SFUU SFMM SFU2=VDVC 4 VSUM BKNT 600 VSUM BKNT 500 VSUM BKNT 400 BKNT 300 SFM2=VDVC 4 VSUM BKNT 850 VSUM BKNT 750 VSUM BKNT 700 BKNT 650 SFL2=VAVR SFU2 SFM2 STFL=VDVC 4 VSUM BKNT 850 VSUM BKNT 700 VSUM BKNT 500 BKNT 300 %%-------------------------------------------------------------------- %% GDI ADVECTION BY THE MEAN-LAYER FLOW %% ADGD: Uses 850-200 hPa flow. %% ADG2: Uses 850-300 hPa flow. Gives more weight to the mid- %% levels of the troposphere (for tropical islands, where %% convection is not as violent as over land, but the %% volume of water in the low-mid troposphere is so large %% that increases the relative weight of the mid-levels %% on the steering of convection) %%-------------------------------------------------------------------- ADGD=SMLC 1+5 ADVT GDM1 SFL1 ADG2=SMLC 1+5 ADVT GDM2 SFL2 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@ ENHANCED GDI @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ JMGALVEZ 20190528 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@BEGIN @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@ STRATUS FORECASTING @@@@@@@@@@@@@@@@@@@@@@@@@@@@ JMGALVEZ 20190530 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@BEGIN % % METHOD % % 1.DETERMINE THE STRENGTH OF THE INVERSION % A first approach is a simple difference between the 850 and 1000 % hPa temperatures. To establish stability, the lapse rate must be % larger than the moist adiabatic lapse rate. For temperatures in % the 10-20C range, the 1000-850 hPa lapse rates range within % -Dry adiabatic: -14 to -12 C % -Moist adiabatic: -6 to -8 C % Thus a lapse rate gentler than -5C will be considered as stable. % % 1.1.GENERATE A MASK TO SHOW ONLY LAPSE RATES GENTLER THAN -5C %===================================================================== MKST=SADC 1 ZPOS SSBC 1 ZNEG SADC 5 SDIF TEMP 850 TEMP 1000 % % 1.2.GENERATE A MASK TO SHOW ONLY LAPSE RATES GENTLER THAN -5C %===================================================================== M2ST=ZNEG SADC 1 ZPOS SDIF TEMP 850 TEMP 1000 %===================================================================== % % 2.DETERMINE IS THERE IS SATURATION AND CLOUDS % Clouds are favored in saturated (RH close to 100%) environments. % Relative humidity (RELH) can be used to assess the degree of % saturation and cloud water mixing ratio (CWMR) for the liquid % water content in clouds. If either or both values are high, the % chance for overcast increases. The peak of saturation and cloud % formation occurs near the base of the inversion. Since this % oscillates between the surface and 850 hPa, a way to determine % if the base of the inversion has clouds is to find the maximum % relative humidity and cloud water mixing ratios under 850 hPa. % % 2.1.MAXIMUM RELATIVE HUMIDITY UNDER 850 HPA %===================================================================== XRHX=SMAX RELH 900 SMAX RELH 925 SMAX RELH 975 RELH 950 XRHM=SMAX RELH 1000 SMAX RELH 850 XRHX %===================================================================== % % 2.2.MAXIMUM CLOUD WATER MIXING RATIO UNDER 850 HPA %===================================================================== XCWX=SMAX CWMR 900 SMAX CWMR 925 SMAX CWMR 975 CWMR 950 XCWM=SMLC 1+5 SMAX CWMR 1000 SMAX CWMR 850 XCWX %===================================================================== % % 2.3.MASKS FOR RELATIVE HUMIDITY %===================================================================== MKR1=ZNEG SADC 1 ZPOS SSBC 88 XRHM % % 2.4.MASKS FOR CLOUD WATER %===================================================================== MKC1=SADC 1 ZPOS SSBC 1 ZNEG XCWM MKC2=ZNEG SADC 1 ZPOS SADC -20 XCWM %===================================================================== % % 3.LARGE SCALE SUBSIDENCE % Relates to the strength of the inversion. The more large scale % subsidence (positive VVEL) relates to a strengthening of the % inversion, less vertical mixing, and favors cloud maintenance or % increase. The more ascent relates to a destabilization of the % inversion, increase in the vertical mixing and thus potential % for cloud dissipation. % % 3.1.AVERAGE VERTICAL VELOCITY %===================================================================== VVRX=SMLC 1+3 SSUM VVEL 950 SSUM VVEL 925 VVEL 900 MKV1=SADC 1 ZPOS SSBC 1 ZNEG SSBC 10 SMTH VVRX % % 4.MOISTURE CONVERGENCE/DIVERGENCE % Favors cloud formation/dissipation % % 4.1. %===================================================================== MCRX=SMLC 1+8 SSUM COMA 1000 SSUM COMA 975 COMA 950 MKM1=SSBC 1 ZNEG SADC 1 ZPOS MCRX MKM2=SADC 1 ZPOS SSBC 1 ZNEG SSBC 10 MCRX %===================================================================== % % % 5.SHEAR IN THE BOUNDARY LAYER (975-900) % Enhances vertical mixing and cloud dissipation. % % 5.1.MAGNITUDE OF THE 900-975 HPA SHEAR %===================================================================== MSHX=MAGN VDIF BKNT 900 BKNT 975 MKS1=SNEG SADC 3 ZPOS SSBC 3 ZNEG SSBC 12 MSHX %===================================================================== % % % 6.WIND SPEED NEAR THE SURFACE - EXTRACTION OF OCEAN MOISTURE % Motivation: Analyses and forecasts of the vertical and horizontal % structure of relative humidity in the cold current % marine boundary layer are limited. Thus, relative % humidity alone cannot be used as a predictor or as % as a diagnosis tool for cloud formation/dissipation. % % Hypothesis: Near-surface relative humidity can be better % represented by considering the flux of water vapor % from the ocean surface. Given weak water vapor % gradients in the marine layer, the role of the speed % of the wind near the surface of the ocean should be % the main predictor. The strength of the inversion % will also matter: a stronger inversion will trap % moisture and speed-up the saturation process. % % 6.1.WIND SPEED NEAR THE SURFACE (WSNS) %===================================================================== %MASK WIND SPEED MWS1=SADC 1 ZPOS SADC -1 ZNEG SSBC 10 WSPK %ADVECTION OF WIND SPEED AWSP=SMLC 1+5 ADVT WSPK WIND 1000 %MASK WIND ADVECTION MWA1=SADC 1 ZPOS SSBC 1 ZNEG AWSP %EFFECTS OF WIND SPEED MWD1=SSUM MWA1 MWS1 MWD2=SMLC 3 SSBC 0.5 MWD1 %WIND DIVERGENCE AT 950 WDV1=SMTH SMLC 1+6 DVRG WIND 950 %WIND DIVERGENCE MASK, WHEN VALUES > 4 WDF1=SADC 2 ZPOS SADC -2 ZNEG SSBC 4 WDV1 %ADVECTION OF WIND DIVERGENCE WDA1=SADC 2 ZPOS SADC -2 ZNEG SADC -2 SMLC 1+4 ADVT WDV1 WIND 950 % 950 WIND DIVERGENCE INDEX WDI1=SNEG SSUM WDF1 WDA1 %===================================================================== % % % 7.STRATUS INDEX % Considers a simple sum of subindices % a)MKC1: Max cloud water values > 10^-5 kg/kg % b)MKC2: Max cloud water values > 20x10^-5 kg/kg % c)MKR1: Max relative Humidity values > 90% %===================================================================== %MKSP=SMLT MKST SSUM MKV1 SSUM MKR1 SSUM MKC1 MKC2 %MKSP=SMLT MKST SSUM MKM1 SSUM MKV1 SSUM MKR1 SSUM MKC1 MKC2 MKSP=SMLT MKST SSUM WDI1 SSUM MWD1 SSUM MKS1 SSUM MKM2 SSUM MKM1 SSUM MKV1 SSUM MKR1 SSUM MKC1 MKC2 %MKSP=SMLT MKST SSUM MKM2 SSUM MKM1 SSUM MKV1 SSUM MKR1 SSUM MKC1 MKC2 %%OTHER EXPLORATORY CALCULATIONS %ISENTROPIC ISN1=SMLC 1+4 SSUM ADVT PRES WIND I297 ADVT PRES WIND I298 ISN2=ISN1 %% MASK TERRAIN %===================================================================== MTR1=ZNEG SADC 1 SNEG ZNEG SSBC 500 HGHT SFC %%AVERAGE CLOUD WATER MIXING RATIOS %===================================================================== CWM1=SDVC 4 SSUM CWMR 975 SSUM CWMR 950 SSUM CWMR 925 CWMR 900 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@ MOISTURE CONVERGENCE @@@@@@@@@@@@@@@@@@@@@@@@@@@ JMGALVEZ 20190528 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@END @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@HAIL MACRO - GR02 - DOES NOT WORK @@@@@@@@@@@@@@@@@@@@@@@@@20190410 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @------MOISTURE CONVERGENCE @COML=SMTH SMLC 1+8 SDVC 4 SSUM COMA 1000 SSUM COMA 925 SSUM COMA 850 COMA 700 @MMC1=SNEG SSBC 1 ZNEG SADC 1 ZPOS COML @MMC2=SNEG SSBC 1 ZNEG SADC 1 ZPOS SADC 10 COML @------VVEL MASK @VVTS=SMTH SMLC 1+4 SDVC 4 SSUM VVEL 700 SSUM VVEL 600 SSUM VVEL 500 VVEL 400 @MVV1=SNEG SSBC 1 ZNEG SADC 1 ZPOS SADC 50 VVTS @MVV2=SNEG SSBC 1 ZNEG SADC 1 ZPOS SADC 200 VVTS @------SHEAR MASK @SHER=MAGN VDIF BKNT 500 BKNT 1000 @MSH1=SADC 1 ZPOS SADC -1 ZNEG SADC -20 SHER @MSH2=SADC 1 ZPOS SADC -1 ZNEG SADC -25 SHER @------500 TEMP MASK @MSKT=SNEG SADC -1 ZNEG SADC 1 ZPOS SADC 9 TEMP 500 @------PRECIPITABLE WATER MASK @MSKP=SADC 1 ZPOS SSBC 1 ZNEG SSBC 20 SMLC 10 PWAT @------MASK LIFTED @MSL1=SADC 1 SNEG ZNEG SADC 1 ZPOS SADC 1.5 SMTH SMTH LIFT @MSL2=SADC 1 SNEG ZNEG SADC 1 ZPOS SADC 3 SMTH SMTH LIFT @======================================================= @ FINAL MASKS @ MSK2 is narrow, and depicts regions of hail potential. @ MSK1 is broader, to plot variables that setup the @ environment for hail. @======================================================= MSK1=SMTH SMLT MSL1 SMLT MSKP SMLT MSKT SMLT MSH1 SMLT MVV1 MMC1 MSK2=SMLT MSL2 SMLT MSKP SMLT MSKT SMLT MSH2 SMLT MVV2 MMC2 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @HAIL MACRO - END @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ SHEAR VECTOR AND LAYER-AVERAGED FLOW - BEGIN @@@@@@@@@@JMG -20190401 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ SHUU=VDIF BKNT 250 BKNT 500 SHMM=VDIF BKNT 500 BKNT 850 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ SHEAR VECTOR AND LAYER-AVERAGED FLOW - END @@@@@@@@@@@@JMG -20190401 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@ TERRAIN MASK - START @@@@@@@@@@@@@@@@@@@@@@@@@@20151216 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@ Masks calculations in elevations that are higher than ~750 mASL in @@ the grid, to avoid plotting ficticious values whenever low-level @@ information is interpolated into regions that are under the terrain @@ in reality. @@-------------------------------------------------------------------- MSKT=SADC 1 ZPOS SSBC 1 SMLC -1 ZPOS SSBC 750 HGHT SFC @@ MASK 700 HPA MS70=SADC 1 ZPOS SSBC 1 SMLC -1 ZPOS SSBC 3000 HGHT SFC @@ MASK 975 HPA MS97=SADC 1 ZPOS SSBC 1 ZNEG SSBC 975 PSFC @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@ TERRAIN MASK (MAR2018) - END @@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ DENSITY CALC BEGIN @@@@@@@@@@@@@@@@@@@@@@@@@@ JOSE_M_GALVEZ 20160615 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@PRESSURE IN PASCALS PREP=SMLC 100 PRES @@TEMPERATURE IN KELVIN TEMK=SADC 0.15 SADC 273 TEMP @@NUMERATOR NUM1=SMLT PREP SADC 1 MIXR @@DENOMINATOR DEN1=SMLC 287 SMLT TEMK SADC 1 SMLC 1.61 MIXR @@DENSITY DENS=SMLT NUM1 INVS DEN1 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ DENSITY CALC END @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@ GDI (GALVEZ-DAVISON INDEX) @@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@ Updated July 10 2014 @@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ %--------------------------------------------------------- % LAYER AVERAGES OF TEMPERATURE AND MIXING RATIOS %--------------------------------------------------------- POTC=THTA 500 MXRC=MIXR 500 POTB=SAVR THTA 700 THTA 850 MXRB=SAVR MIXR 700 MIXR 850 POTA=THTA 925 MXRA=MIXR 925 %--------------------------------------------------------- % EPT (EQUIVALENT POTENTIAL TEMP) PROXY FACTORS %--------------------------------------------------------- %%% PROXY FOR TEMPERATURE LCL IN KELVIN (TEMP 850HPA) LCLK=SADC 0.15 SADC 273 TEMP 850 %%% DENOMINATOR FOR EXPONENTIAL EXDE=SMLC 1006 LCLK %%% EXPONENTIALS FOR LAYERS C, B, A EXPC=EXPP INVS SDVD EXDE SMLC 1000 SMLC 2690 MXRC EXPB=EXPP INVS SDVD EXDE SMLC 1000 SMLC 2690 MXRB EXPA=EXPP INVS SDVD EXDE SMLC 1000 SMLC 2690 MXRA %%% EPT PROXYS (EQUIVALENT POTENTIAL TEMPERATURES) EPTC=SADC -10 SMLT POTC EXPC EPTB=SADC -10 SMLT POTB EXPB EPTA=SMLT POTA EXPA %--------------------------------------------------------- % TI (THERMAL INDEX) %--------------------------------------------------------- %%% MC (MID-TROPOSPHERIC EPT CORE, BASED ON LAYER C EPT) MMCC=SSBC 303 EPTC %%% E (ENHANCEMENT BASED ON LAYER A EPT) EEEE=ZNEG SMLC .065 SSBC 303 EPTA %%% C (CORE FACTOR) CCFF=SMLT MMCC EEEE %%% MT (MID-LEVEL TEMPERATURE THRESHOLD) MTTT=ZNEG SADC 10 TEMP 500 %%% MW (MID-LEVEL WARMING FACTOR) MWWW=SNEG SMLC 7.0 MTTT %%% TI (THERMAL INDEX) TTII=SSUM MWWW CCFF %--------------------------------------------------------- % II (INVERSION INDEX) %--------------------------------------------------------- %%% SSSS (STABILITY FACTOR) SSSS=SMLC 1.5 SDIF TEMP 925 TEMP 700 %%% DDDD (DRYING FACTOR) DDDD=SMLC 1.5 SDIF EPTB EPTA %%% II (INVERSION INDEX) IIII=SNEG ZNEG SNEG SSUM SSSS DDDD %--------------------------------------------------------- % CORRECTION FOR VISUALIZATION OVER MOUNTAINS %--------------------------------------------------------- %CCOO=SNEG SSBC 18 INVS SDVC 9000 SSBC 500 PSFC CCOO=SNEG SSBC 18 INVS SDVC 6000 SSBC 500 PSFC %--------------------------------------------------------- % GDI %--------------------------------------------------------- GDIN=SSUM TTII IIII GDIF=SSUM CCOO SSUM TTII IIII @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@ GDI (GALVEZ-DAVISON INDEX) - END @@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@ WINDS AND DIVERG. AVERAGED IN THE VERTICAL - START @@@@@20140314 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@ LOW-LEVEL WINDS USING 5 LEVELS WLL5=VDVC 5 VSUM WIND 750 VSUM VSUM WIND 800 WIND 850 VSUM WIND 950 WIND 900 @@@@@@@@@@@@@ LOWER WINDS WW11=VAVR WIND 975 WIND 1000 WW99=VAVR WIND 950 WIND 925 WWLL=VAVR WW11 WW99 DDLL=SMLC 1000 SMLC 1000 DVRG WIND WWLL @@@@@@@@@@@@@ LOW/MID WINDS WW19=VAVR WIND 1000 WIND 950 WW88=VAVR WIND 900 WIND 850 WWML=VAVR WW19 WW88 DV19=SAVR DVRG WIND 1000 DVRG WIND 950 DV88=SAVR DVRG WIND 900 DVRG WIND 850 DDML=SMLC 1000 SMLC 1000 SAVR DV19 DV88 @@@@@@@@@@@@@ MID WINDS WW76=VAVR WIND 700 WIND 600 WW54=VAVR WIND 500 WIND 450 WWMM=VAVR WW76 WW54 @@@@@@@@@@@@@ MID/UPPER WINDS WW34=VAVR WIND 300 WIND 400 WW56=VAVR WIND 600 WIND 500 WWMU=VAVR WW56 WW34 DV34=SAVR DVRG WIND 300 DVRG WIND 400 DV56=SAVR DVRG WIND 600 DVRG WIND 500 DDMU=SMLC 6000 SMLC 1000 SAVR DV56 DV34 @@@@@@@@@@@ UPPER WINDS WW22=VAVR WIND 200 WIND 250 WW43=VAVR WIND 350 WIND 300 WWUU=VAVR WW43 WW22 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@ WINDS AND DIVERGENCE AVERAGED IN THE VERTICAL - END @@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ MACRO ENGE [Engelamiento/Icing] @ Last Update: JMGalvez 20160729 @ ========================================== @ Depicts the potential for icing using @ (a) temperatures between 0C and -20C @ (b) wet bulb temperatures between 0 and -20C @ (c) negative omegas (ascent) @ (d) cloud water mixing ratios @ (e) temperatures between -2 and -12C for enhanced icing risk @ Calculates data for four layers [700, 600, 500 and 400 hPa] @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@ ENGE CALCULATIONS BEGIN @@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @MAXIMUM NEGATIVE OMEGAS IN EACH THREE-LEVEL LAYER [mb/s] MXO7=SMIN VVEL 750 SMIN VVEL 700 VVEL 650 MXO6=SMIN VVEL 650 SMIN VVEL 600 VVEL 550 MXO5=SMIN VVEL 550 SMIN VVEL 500 VVEL 450 MXO4=SMIN VVEL 450 SMIN VVEL 400 VVEL 350 @BINARY MASKS TO DEPICT NEGATIVE OMEGAS (ASCENT) AS 1 MMO7=SADC 1 ZPOS SADC -1 SNEG ZPOS SMLC 1+4 MXO7 MMO6=SADC 1 ZPOS SADC -1 SNEG ZPOS SMLC 1+4 MXO6 MMO5=SADC 1 ZPOS SADC -1 SNEG ZPOS SMLC 1+4 MXO5 MMO4=SADC 1 ZPOS SADC -1 SNEG ZPOS SMLC 1+4 MXO4 @BINARY MASKS TO DEPICT NEGATIVE OMEGAS (ASCENT) LARGER THAN 3microbar/s M2O7=SNEG SADC -1 ZNEG SADC 1 ZPOS SMTH SADC 30 SMLC 1+4 MXO7 M2O6=SNEG SADC -1 ZNEG SADC 1 ZPOS SMTH SADC 30 SMLC 1+4 MXO6 M2O5=SNEG SADC -1 ZNEG SADC 1 ZPOS SMTH SADC 30 SMLC 1+4 MXO5 M2O4=SNEG SADC -1 ZNEG SADC 1 ZPOS SMTH SADC 30 SMLC 1+4 MXO4 @AVERAGE CLOUD WATER MIXING RATIOS IN g/m3 CWM7=SMLC 1000 SMAX CWMR 750 SMAX CWMR 700 CWMR 650 CWM6=SMLC 1000 SMAX CWMR 650 SMAX CWMR 600 CWMR 550 CWM5=SMLC 1000 SMAX CWMR 550 SMAX CWMR 500 CWMR 450 CWM4=SMLC 1000 SMAX CWMR 450 SMAX CWMR 400 CWMR 350 @MASKS TO DEPICT CLOUD WATER MIXING RATIOS > 1g/m3 MCW7=SADC 1 ZPOS SADC -1 ZNEG SADC -1 SMLC 100 CWM7 MCW6=SADC 1 ZPOS SADC -1 ZNEG SADC -1 SMLC 100 CWM6 MCW5=SADC 1 ZPOS SADC -1 ZNEG SADC -1 SMLC 100 CWM5 MCW4=SADC 1 ZPOS SADC -1 ZNEG SADC -1 SMLC 100 CWM4 @AVERAGE TEMPERATURES IN C TEM7=SDVC 3 SSUM TEMP 750 SSUM TEMP 700 TEMP 650 TEM6=SDVC 3 SSUM TEMP 650 SSUM TEMP 600 TEMP 550 TEM5=SDVC 3 SSUM TEMP 550 SSUM TEMP 500 TEMP 450 TEM4=SDVC 3 SSUM TEMP 450 SSUM TEMP 400 TEMP 350 @MASKS TO DEPICT TEMPS BETWEEN 0 AND -20C MTE7=SADC 1 SNEG ZNEG SADC 1 SMLT ZNEG SADC 20 TEM7 ZPOS TEM7 MTE6=SADC 1 SNEG ZNEG SADC 1 SMLT ZNEG SADC 20 TEM6 ZPOS TEM6 MTE5=SADC 1 SNEG ZNEG SADC 1 SMLT ZNEG SADC 20 TEM5 ZPOS TEM5 MTE4=SADC 1 SNEG ZNEG SADC 1 SMLT ZNEG SADC 20 TEM4 ZPOS TEM4 @MASKS TO DEPICT TEMPS BETWEEN -2 AND -12C MTT7=SADC 1 SNEG ZNEG SADC 1 SMLT ZNEG SADC 12 TEM7 ZPOS SADC 2 TEM7 MTT6=SADC 1 SNEG ZNEG SADC 1 SMLT ZNEG SADC 12 TEM6 ZPOS SADC 2 TEM6 MTT5=SADC 1 SNEG ZNEG SADC 1 SMLT ZNEG SADC 12 TEM5 ZPOS SADC 2 TEM5 MTT4=SADC 1 SNEG ZNEG SADC 1 SMLT ZNEG SADC 12 TEM4 ZPOS SADC 2 TEM4 @AVERAGE WET BULB TEMPERATURES IN C WBT7=SDVC 3 SSUM WBTC 750 SSUM WBTC 700 WBTC 650 WBT6=SDVC 3 SSUM WBTC 650 SSUM WBTC 600 WBTC 550 WBT5=SDVC 3 SSUM WBTC 550 SSUM WBTC 500 WBTC 450 WBT4=SDVC 3 SSUM WBTC 450 SSUM WBTC 400 WBTC 350 @MASKS TO DEPICT WET BULB TEMPS BETWEEN 0 AND -20C MWB7=SADC 1 SNEG ZNEG SADC 1 SMLT ZNEG SADC 20 WBT7 ZPOS WBT7 MWB6=SADC 1 SNEG ZNEG SADC 1 SMLT ZNEG SADC 20 WBT6 ZPOS WBT6 MWB5=SADC 1 SNEG ZNEG SADC 1 SMLT ZNEG SADC 20 WBT5 ZPOS WBT5 MWB4=SADC 1 SNEG ZNEG SADC 1 SMLT ZNEG SADC 20 WBT4 ZPOS WBT4 @--------------------------------------------------------------------- @MASK CWMR, TEMP AND WET BULB TEMP @Gives 5 when the five variables fall within thresholds. @Gives twice as much weight to omegas higher than -1 microbar/s (ascent) @--------------------------------------------------------------------- MIC7=SSUM MTT7 SSUM SMLC 2 M2O7 SSUM MWB7 SSUM MTE7 MCW7 MIC6=SSUM MTT6 SSUM SMLC 2 M2O6 SSUM MWB6 SSUM MTE6 MCW6 MIC5=SSUM MTT5 SSUM SMLC 2 M2O5 SSUM MWB5 SSUM MTE5 MCW5 MIC4=SSUM MTT4 SSUM SMLC 2 M2O4 SSUM MWB4 SSUM MTE4 MCW4 @--------------------------------------------------------------------- @ADJUSTED CLOUD WATER MIXING RATIOS IN g/m3 - EMPIRICAL ADJUSTMENT @Adds 0.4 to CWMR if temperature is between -2 and -12 @Adds 0.4 if omegas are larger than 3 microbars/s @--------------------------------------------------------------------- CWF7=SMLT M2O7 SMLT MTE7 SSUM SMLC 0.4 M2O7 SSUM SMLC 0.4 MTT7 CWM7 CWF6=SMLT M2O6 SMLT MTE6 SSUM SMLC 0.4 M2O6 SSUM SMLC 0.4 MTT6 CWM6 CWF5=SMLT M2O5 SMLT MTE5 SSUM SMLC 0.4 M2O5 SSUM SMLC 0.4 MTT5 CWM5 CWF4=SMLT M2O4 SMLT MTE4 SSUM SMLC 0.4 M2O4 SSUM SMLC 0.4 MTT4 CWM4 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@ ENGE CALCULATIONS END @@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ MACRO ICIN - BEGIN @@@@@@@@@@@@@@@@@@@@@@@@@@ JOSE_M_GALVEZ 20160620 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ ICIN=Icing @ Depicts areas with the potential for icing based upon @ relative humidity > 70%, wet builb temperatures of 0 tot -20C, @ and negative omegas. Considers severe icing when the magnitude @ of negative omegas is larger than 6microbar/s @--------------------------------------------------------------------- @MASK FOR RELH GRTN 70 MK01=SADC 1 ZPOS SADC -1 ZNEG SADC -70 RELH @MASK FOR OMEGAS MK02=SADC 1 ZPOS SADC -1 ZNEG SNEG SMLC 1+5 SMTH VVEL MK03=SADC 1 SNEG ZNEG SADC 1 ZPOS SADC 600 SMLC 1+5 SMTH VVEL @MASK FOR WET BULB TEMPS OF 0 TO -20C MK04=ZNEG SADC 20 WBTC MK05=ZPOS WBTC MK06=SADC 1 SNEG ZNEG SADC 1 SMLT MK05 MK04 @MASK ICING MKIC=SMLT MK06 SMLT MK01 MK02 MKIE=SMLT MK06 SMLT MK01 MK03 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ MACRO ICIN - END @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@ COSTA ECUADOR/NORTE PERU - START @@@@@@@@@@ JOSE_M_GALVEZ 20170315 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ OROG=SNEG ZNEG SNEG SSBC 700 PSFC IRM1=SMLC 1+4 SDVC 4 SSUM MIXR 700 SSUM MIXR 850 SSUM MIXR 950 MIXR 900 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@ COSTA ECUADOR/NORTE PERU - END @@@@@@@@@@@@ JOSE_M_GALVEZ 20170315 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@ INDICE ANDES CENTRALES - RICARDO DURAN - JOSE GALVEZ - START @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ V85A=VDVC 4 VSUM VSUM WIND 500 WIND 600 VSUM WIND 700 WIND 850 V85B=SMLC 5 SNEG XCMP V85A R85A=SDVC 4 SSUM SSUM RELH 500 RELH 600 SSUM RELH 700 RELH 850 R85B=ZNEG SSBC 60 R85A DALT=SMTH SMLC 2+5 DVRG WIND 200 CFRM=SMTH SMLC 1+8 SNEG DVRG FLUX MIXR WIND 600 IND1=SSUM CFRM SSUM DALT SSUM R85B V85B @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@ INDICE ANDES CENTRALES - RICARDO DURAN - JOSE GALVEZ - END @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ PARAGUAY SEVERE WEATHER - START @@@@@@@@@@@@@ JOSE_M_GALVEZ 20150315 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@Shear SH03=MAGN VDIF WIND 700 WIND 1000 SH06=MAGN VDIF WIND 450 WIND 1000 SHC3=SADC 12 ZNEG SSBC 12 SH03 SHC6=SADC 20 ZNEG SSBC 20 SH06 SHSU=SSUM SHC6 SHC3 @@Low-level convergence CV89=SMLC 1+6 SAVR DVRG WIND 850 DVRG WIND 925 @@Vvel VVNB=SMLC 1+3 VVEL 925 VVBM=SMLC 3+2 SSUM VVEL 700 SSUM VVEL 925 VVEL 850 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@PARAGUAY SEVERE WEATHER - END @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@STORM RELATIVE HELICITY 0-3km - START - might have problems @@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ CVEC=VDVC 4 VSUM WIND 700 VSUM WIND 800 VSUM WIND 900 WIND 975 @@3KM (700) UWL7=XCMP VDIF WIND 700 CVEC VSH7=SDIF YCMP WIND 750 YCMP WIND 650 VWL7=YCMP VDIF WIND 700 CVEC USH7=SDIF XCMP WIND 750 XCMP WIND 650 HEL7=SSUM SMLT VWL7 USH7 SNEG SMLT UWL7 VSH7 @@2KM (800) UWL8=XCMP VDIF WIND 800 CVEC VSH8=SDIF YCMP WIND 850 YCMP WIND 750 VWL8=YCMP VDIF WIND 800 CVEC USH8=SDIF XCMP WIND 850 XCMP WIND 750 HEL8=SSUM SMLT VWL8 USH8 SNEG SMLT UWL8 VSH8 @@1.5KM (900) UWL9=XCMP VDIF WIND 900 CVEC VSH9=SDIF YCMP WIND 950 YCMP WIND 850 VWL9=YCMP VDIF WIND 900 CVEC USH9=SDIF XCMP WIND 950 XCMP WIND 850 HEL9=SSUM SMLT VWL9 USH9 SNEG SMLT UWL9 VSH9 @@SFC UWL1=XCMP VDIF WIND 975 CVEC VSH1=SDIF YCMP WIND 1000 YCMP WIND 950 VWL1=YCMP VDIF WIND 975 CVEC USH1=SDIF XCMP WIND 1000 XCMP WIND 950 HEL1=SSUM SMLT VWL1 USH1 SNEG SMLT UWL1 VSH1 HELI=SNEG SSUM HEL1 SSUM HEL9 SSUM HEL7 HEL8 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @STORM RELATIVE HELICITY - END @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ GDI-CHILE (GALVEZ-DAVISON INDEX FPR CHILE) @@@@@@@@@@@@@@@@20160415 @ Use with care, still requires validation as of Mar 2018@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ %%% EPT PROXYS (EQUIVALENT POTENTIAL TEMPERATURES) POT1=SDVC 3 SSUM THTA 925 SSUM THTA 850 THTA 700 MXR1=SDVC 3 SSUM MIXR 925 SSUM MIXR 850 MIXR 700 %%% EPT PROXYS (EQUIVALENT POTENTIAL TEMPERATURES) EX1A=EXPP INVS SDVD EXDE SMLC 1000 SMLC 2690 MXR1 EP1A=SMLT POT1 EX1A %%% E (ENHANCEMENT BASED ON LAYER A EPT) EEE1=ZNEG SMLC .350 SSBC 303 EP1A %%% C (EPT CORE INDEX) CCF1=SMLT MMCC EEE1 %%% MT (MID-LEVEL TEMPERATURE THRESHOLD) MTT1=ZNEG SADC 15 TEMP 500 %%% MW (MID-LEVEL WARMING FACTOR) MWW1=SNEG SMLC 20 MTT1 %%% TI (THERMAL INDEX) TTI1=SSUM MWW1 CCF1 %%% GDI ADJUSTED FOR CHILE GDIC=SSUM CCOO TTI1 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@ GDI (GALVEZ-DAVISON INDEX - CHILE) - END @@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@OTHER CALCULATIONS - START @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ CMBB=SMLC 1+8 SDVC 0004 SSUM COMA 800 SSUM COMA 850 SSUM COMA 900 COMA 975 @----------------------------------------------------------------- @ K-INDEX @----------------------------------------------------------------- KINX=SSUM DWPT 850 SSUM SDIF TEMP 850 TEMP 500 SDIF DWPT 700 TEMP 700 @----------------------------------------------------------------- @ TOTAL-TOTALS @----------------------------------------------------------------- TOTI=SSUM SDIF TEMP 850 TEMP 500 SDIF DWPT 850 TEMP 500 @----------------------------------------------------------------- @ AVERAGED MID-LAYER VVEL AND RELH @----------------------------------------------------------------- RELM=SAVR SAVR RELH 700 RELH 600 SAVR RELH 500 RELH 400 VVLM=SAVR SAVR VVEL 700 VVEL 600 SAVR VVEL 500 VVEL 400 @----------------------------------------------------------------- @ GDI AVERAGED OVER TIME - May not run due to limitations on the @ number of commands that can be used in Wingridds V3.2 @----------------------------------------------------------------- GDIP=SDVC 5 SSUM GDIF F12 SSUM GDIF F18 SSUM GDIF F24 SSUM GDIF F30 GDIF F36 @----------------------------------------------------------------- @ MID-LEVEL VORTICITY (500-600 HPA) @----------------------------------------------------------------- VT56=SMLC 1+6 SAVR RVRT WIND 500 RVRT WIND 600 @----------------------------------------------------------------- @ VERTICAL VELOCITY @----------------------------------------------------------------- VV54=SAVR VVEL 500 VVEL 400 VV67=SAVR VVEL 700 VVEL 600 VVML=SMLC 1000 SAVR VV54 VV67 @----------------------------------------------------------------- @ DEFORMATION @----------------------------------------------------------------- DEFO=SQRT SSUM SMLT WSHD WSHD SMLT WSTD WSTD @----------------------------------------------------------------- @ MOISTURE CONVERGENCE @----------------------------------------------------------------- CM18=SMLC 1+8 SDVC 0003 SSUM COMA 925 SSUM COMA 850 COMA 975 vv18=SDVC 0003 SSUM VVEL 600 SSUM VVEL 700 VVEL 850 @@------------------------------------------------------------- @ POTENTIAL VORTICITY (VORTICIDAD POTENCIAL) @@------------------------------------------------------------- VPOT=SMLC -1+6 SMLT WVRT LEVL SDVD THTA PRES LDIF LY-1 @@------------------------------------------------------------- @ POTENTIAL VORTICITY ADVECTION @ (ADVECCION DE VORTICIDAD POTENCIAL) @@------------------------------------------------------------- AVPT=SMLC 1+3 ADVT VPOT WIND @@------------------------------------------------------------- @ POTENTIAL VORTICITY AVERAGES OVER MANY DAYS @ (PROMEDIOS DE VORT POTENCIAL SOBRE VARIOS DIAS) @@------------------------------------------------------------- VPT1=SDIF VPOT F24 VPOT F00 VPT2=SDIF VPOT F48 VPOT F24 VPT3=SDIF VPOT F72 VPOT F48 VPT4=SDIF VPOT F96 VPOT F72 VPTT=SDVC 3 SSUM VPT2 SSUM VPT3 VPT1 @----------------------------------------------------------------- @ THERMAL ADVECTION AVERAGES OVER SEVERAL DAYS @ (PROMEDIOS DE ADVECCION TERMAL SOBRE VARIOS DIAS) @----------------------------------------------------------------- ATM1=SAVR ADVT TEMP WIND B015 F00 ADVT TEMP WIND B015 F12 ATM2=SAVR ADVT TEMP WIND B015 F24 ADVT TEMP WIND B015 F36 ATM3=SAVR ADVT TEMP WIND B015 F48 ADVT TEMP WIND B015 F60 ATM4=SAVR ADVT TEMP WIND B015 F72 ADVT TEMP WIND B015 F84 ATMT=SDVC 4 SSUM ATM1 SSUM ATM2 SSUM ATM3 ATM4 @----------------------------------------------------------------- @ OTHER VARIABLES @----------------------------------------------------------------- AMDV=DPOS DVRG FLUX PRES LDIF AGEO LAVE GMDV=DPOS DVRG FLUX PRES LDIF GEOS LAVE TMDV=DPOS DVRG FLUX PRES LDIF WIND LAVE THKX=RSLV SDIF HGHT 500 SMLC 7.5 SSBC 1000 PMSL SVLV PVRT=C5-6 SMLT RLTN -1 MGTN -49 VPDF SMLT WVRT LEVL SDVD THTA PRES LDIF LY-1 VPDF=SDIF PRES LVL0 PRES 0000 VENT=WIND VNTO=WIND KEPS=LAST KEPV=LAST THXP=PBOX&THKX CI60&PMSL CIN4 PBX1=BOX1 LSTN 0.25 GRTN 0.01 PBX2=BOX2 LSTN 0.50 GRTN 0.25 PBX3=BOX3 LSTN 1.00 GRTN 0.50 PBX4=BOX4 LSTN 2.00 GRTN 1.00 PBX5=BOX5 GRTN 2.00 PBOX=CLRG PBX5 LAST&CLRF PBX4 LAST&CLRF PBX3 LAST&CLRE PBX2 LAST&CLRE PBX1 TPCI TMPL=SDIF TEMP SMT9 TEMP LSTX=LAST& INRI=SDVD SMLT THTA SMLT SCPY MAGN WIND LDIF SMLT HGHT LDIF THTA LDIF LY-1 RIN1=SDVD[SDVD(THTA LDIF,HGHT LDIF)THTA LAVE]LY-1 RIN2=SMLT SCPY SDVD MAGN WIND LDIF HGHT LDIF LY-1 RIG2=SMLT SCPY SDVD MAGN GEOS LDIF HGHT LDIF LY-1 IGRI=SDVD SMLT THTA SMLT SCPY MAGN GEOS LDIF SMLT HGHT LDIF THTA LDIF LY-1 QVC1=SMLC -1 SSUM SMLT DSDX XCMP GEOS DSDX TEMP SMLT DSDX YCMP GEOS DSDY TEMP QVC2=SMLC -1 SSUM SMLT DSDY XCMP GEOS DSDX TEMP SMLT DSDY YCMP GEOS DSDY TEMP MSFC=SSUM NORM GEOS SMLT FFFF DIST MSKS=SSUM RGTN 0000 MLTN 0000 SDIF PSFC PRES BOXS=BOX9 LT00 MGTN 0000 SDIF PSFC PRES CLRF AVRD=ADVT RVRT WIND WIND VNM1=SMTH SMLC 1+6 SAVR RVRT WIND 500 RVRT WIND 600 WNM1=VAVR BKNT 500 BKNT 600 GUS2=VDIF WIND 10M VMLC 0.70 WIND 925 D850=SMTH SMLC 1+6 SAVR DVRG WIND 850 DVRG WIND 925 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ OTHER CALCULATIONS - END @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@