*&^%=USER ALIAS FILE @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@ INTERNATIONAL DESKS ALIAS.USR @@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@ EQUATION-PROGRAMMING FILE FOR @@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@ WINGRIDDS V3.2, LAST UPDATED @@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@ ON 23MAR2018 BY JOSE M GALVEZ @@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@ TERRAIN MASK - START @@@@@@@@@@@@@@@@@@@@@@@@@@20151216 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@ The terrain masks block calculations in elevations that are higher @@ than a threshold inside the grid, to avoid plotting ficticious @@ values whenever low-level information is interpolated into regions @@ that are inside the grid terrain. @@-------------------------------------------------------------------- @@ MASK FOR ELEVATIONS HIGHER THAN 750 METERS ABOVE SEA LEVEL MSKT=SADC 1 ZPOS SSBC 1 SMLC -1 ZPOS SSBC 750 HGHT SFC @@ MASK 975 HPA MS97=SADC 1 ZPOS SSBC 1 ZNEG SSBC 975 PSFC @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@ TERRAIN MASK (MAR2018) - END @@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@ MOISTURE CONVERGENCE - START @@@@@@@@@@@@@@ JOSE_M_GALVEZ 20180323 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@ Computation of the convergence of the flux of mixing ratio (COMA). @@ This relates to the low-level convergence of moisture, and it is @@ very useful to assess the triggering of convection by low-level @@ converging features whenever moisture is available. Note that low- @@ level convergence can also occur in regions with limited moisture, @@ which is ignored by this calculation. COMA is the convergence of @@ the flux of mixing ratio for any level. COLL is the average for @@ the 1000-850 hPa layer, particularly useful in the tropics during @@ the dry or transition season. For the rainy season or deeper-layer @@ moisture convergence, COL2 can be used, as it considers 700 hPa. @@-------------------------------------------------------------------- COMA=DVRG FLUX MIXR WIND COLL=SDVC 3 SSUM COMA 925 SSUM COMA 850 COMA 1000 COL2=SAVR SAVR COMA 700 COMA 850 SAVR COMA 925 COMA 1000 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@ MOISTURE CONVERGENCE - END @@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ DENSITY CALC START @@@@@@@@@@@@@@@@@@@@@@@@@@ JOSE_M_GALVEZ 20160615 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ This simple calculation estimates the atmospheric density @ considering the equation of state and moisture content. @--------------------------------------------------------------------- @@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 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@ The GDI (GALVEZ-DAVISON INDEX) @@@@@@@ JOSE_M_GALVEZ 20140710 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ %--------------------------------------------------------- % 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 %------------------------------------------------------------------- % GDI WITHOUT TERRAIN CORRECTION % Note that the terrain correction CCOO might require adjustment % upon model or analysis used. For example, as of 2018, adjustment % for the ECMWF and WRF model data seems to require more subtle % (less ample) values of CCOO than those used for the GFS data. % Resolution also matters. Note that the goal of CCOO is only for % forecaster visualization purposes. %------------------------------------------------------------------- GDIN=SSUM TTII IIII %------------------------------------------------------------------- % GDI (WITH TERRAIN CORRECTION) %------------------------------------------------------------------- GDIF=SSUM CCOO SSUM TTII IIII @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@ GDI (GALVEZ-DAVISON INDEX) - END @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@ TRADE WIND INVERSION HEIGHT/STRENGTH - START @@@@@@@@20180315 @@@@@@@@@@@@@@@@@@@ MAR.13.2017 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@ This set of equations estimate 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 what convection will do if it reaches @@ top of the moist layer. Example: If the inversion is weak and dry @@ air is present over it, dry air entrainment into the marine layer @@ might play a role in decreasing potential precipitation amounts. If @@ the inversion is strong, vertical mixing across it might be limited @@ allowing for a dense stratocumulus cloud deck. Shallow moist layers @@ are related to less precipitation than deep ones. @@-------------------------------------------------------------------- @@ 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 @@ Estimates the maximum pressure where the gentle lapse rate @@ can be found. This implies that the equations look for the @@ lowest inversion, ignoring those above. @@-------------------------------------------------------------------- ILEV=SMAX IC70 SMAX IC75 SMAX IC80 SMAX IC85 SMAX IC90 IC10 @@-------------------------------------------------------------------- @@ Part 3. MID-LEVEL WINDS @@ To evaluate the location and structure of mid ridges and @@ troughs, which relate to the structure of the trade wind @@ inversion. @@-------------------------------------------------------------------- WI78=VAVR WIND 700 WIND 600 WI56=VAVR WIND 500 WIND 550 WIML=VAVR WW76 WW54 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@ TRADE WIND INVERSION HEIGHT/STRENGTH - END @@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@ HAIL MACRO EQUATIONS (GR01) @@@@@@@@@@@@@@@@@@ NESTOR_SANTAYANA AND @@ MACRO PARA GRANIZO @@@@@@@@@@@@@@@@@@@@@@@@ JOSE_M_GALVEZ 20151216 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ These equations define variables that are used by the Wingridds @ macro GR01.CMD, designed to evaluate the potential for hail in @ South America. The macro was first developed by Nestor Santayana @ (Uruguay Weather Service, 2012-2013), and adjusted by J.M.Galvez @ on 2014. It helps to assess the potential for the occurrence of @ hail in subtropical (20-40 degrees) latitudes in South America. @--------------------------------------------------------------------- @ MOISTURE CONVERGENCE @--------------------------------------------------------------------- COML=COLL MSKC=SNEG SSBC 1 ZNEG SADC 1 ZPOS SADC 20 COML DATA @--------------------------------------------------------------------- @ VERTICAL VELOCITY MASK @--------------------------------------------------------------------- VVTS=SMTH SMLC 1+4 SDVC 4 SSUM VVEL 700 SSUM VVEL 600 SSUM VVEL 500 VVEL 400 MVTS=SNEG SSBC 1 ZNEG SADC 1 ZPOS SADC 10 VVTS MVT2=SNEG SSBC 1 ZNEG SADC 1 ZPOS SADC 25 VVTS @--------------------------------------------------------------------- @ PRECIPITABLE WATER MASK @--------------------------------------------------------------------- MSKP=SADC 1 ZPOS SSBC 1 ZNEG SSBC 20 SMLC 10 PWAT @--------------------------------------------------------------------- @ LIFTED INDEX MASK @--------------------------------------------------------------------- MSKL=SADC 1 SNEG ZNEG SADC 1 ZPOS SADC 3 SMTH SMTH LIFT @--------------------------------------------------------------------- @ FINAL MASK INCLUDING LIFTED INDEX @--------------------------------------------------------------------- MSKF=SMLT MSKP SMLT MSKT SMLT MVTS MSKL MSK2=SMLT MSKP SMLT MSKT SMLT MVT2 MSKL @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ HAIL INDEX - END @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@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 400 WIND 300 WWUU=VAVR WW43 WW22 DV22=SAVR DVRG WIND 200 DVRG WIND 250 DV43=SAVR DVRG WIND 400 DVRG WIND 300 DDUU=SMLC 6000 SMLC 1000 SAVR DV22 DV43 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@ WINDS AND DIVERGENCE AVERAGED IN THE VERTICAL - END @@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@ ENGE MACRO EQUATIONS: ICING POTENTIAL @@@@@@@@@@@@@@@@@@@@@@@@@@@ JOSE_M_GALVEZ @@ ECUACIONES DEL MACRO ENGE: POTENCIAL DE ENGELAMIENTO @@@@@@@@@@@@@@@@@ 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. These equations represent an alternative approach @ to evaluate the potential for icing. The associated macro (ICIN.CMD) @ 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 the 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 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@ BELOW THIS LINE LIE EQUATIONS THAT REQUIRE DOCUMENTATION AND THAT @@@ ARE USED WITH LESS FREQUENCY AT THE WPC INTERNATIONAL DESKS @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@ 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 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @----------------------------------------------------------------- @ 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 @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@