[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[IDV #OMI-310298]: Help Me please GDI _Index

This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.

  • Subject: [IDV #OMI-310298]: Help Me please GDI _Index
  • Date: Fri, 18 Aug 2023 14:53:56 -0600
> Hi,
> I used your software, how fix the following error?
> 
> def GDI_index(T950, T850, T700, T500, R950, R850, R700, R500, PSFC):
> # Cálculos para ThetaA, ThetaB y ThetaC
> ThetaA = T950 * (1000/950)**(2/7)
> ThetaB = 0.5 * (T850*(1000/850)**(2/7) + T700*(1000/700)**(2/7))
> ThetaC = T500 * (1000/950)**(2/7)
> 
> # Cálculos para EPTPa, EPTPb y EPTPc
> Ra = R950
> Rb = 0.5 * (R850 + R700)
> Rc = R500
> EPTPa = ThetaA * exp(((2.69*10**6)*Ra) / (1005.7*T850))
> EPTPb = ThetaB * exp(((2.69*10**6)*Rb) / (1005.7*T850)) - 10
> EPTPc = ThetaC * exp(((2.69*10**6)*Rc) / (1005.7*T850)) - 10
> 
> # Cálculos para ME y LE
> ME = EPTPc - 303
> LE = EPTPa - 303
> 
> # Cálculos para CBI
> if LE > 0:
> CBI = (6.5*10**-2) * LE * ME
> else:
> CBI = 0
> 
> # Cálculos para MWI
> if T500 - 263.15 > 0:
> MWI = -7 * (T500 - 263.15)
> else:
> MWI = 0
> 
> # Cálculos para II
> if (T950 - T700) + (EPTPb - EPTPa) > 0:
> II = 0
> else:
> II = 1.5 * (T950 - T700) + (EPTPb - EPTPa)
> 
> # Cálculos para TC
> TC = 18 - (9000 / PSFC - 500)
> 
> # Cálculo final del GDI
> GDI = CBI + MWI + II + TC
> return GDI
> 
> [cid:image001.png@01D9D056.D4DC6A00]
> Jesús L. Salazar G.
> Meteorólogo
> Vicepresidencia de Administración Recurso Hídrico (HI)
> Tel: +(507) 276-1729 - www.pancanal.com<http://www.pancanal.com/>
> e-mail: address@hidden<mailto:address@hidden>
> 
> 
> 
Please post the error message and provide some information of this formula.

Yuan

Ticket Details
===================
Ticket ID: OMI-310298
Department: Support IDV
Priority: Normal
Status: Open
===================
NOTE: All email exchanges with Unidata User Support are recorded in the Unidata 
inquiry tracking system and then made publicly available through the web.  If 
you do not want to have your interactions made available in this way, you must 
let us know in each email you send to us.