Esclonamento de máquinas de cogeração utilzando programação inteira mista

As centrais termoelétricas convencionais convertem apenas parte do combustível consumido na produção de energia elétrica, sendo que outra parte resulta em perdas sob a forma de calor. Neste sentido, surgiram as unidades de cogeração, ou Combined Heat and Power (CHP), que permitem reaproveitar a energia dissipada sob a forma de energia térmica e disponibilizá-la, em conjunto com a energia elétrica gerada, para consumo doméstico ou industrial, tornando-as mais eficientes que as unidades convencionais Os custos de produção de energia elétrica e de calor das unidades CHP são representados por uma função não-linear e apresentam uma região de operação admissível que pode ser convexa ou não-convexa, dependendo das caraterísticas de cada unidade. Por estas razões, a modelação de unidades CHP no âmbito do escalonamento de geradores elétricos (na literatura inglesa Unit Commitment Problem (UCP)) tem especial relevância para as empresas que possuem, também, este tipo de unidades. Estas empresas têm como objetivo definir, entre as unidades CHP e as unidades que apenas geram energia elétrica ou calor, quais devem ser ligadas e os respetivos níveis de produção para satisfazer a procura de energia elétrica e de calor a um custo mínimo. Neste documento são propostos dois modelos de programação inteira mista para o UCP com inclusão de unidades de cogeração: um modelo não-linear que inclui a função real de custo de produção das unidades CHP e um modelo que propõe uma linearização da referida função baseada na combinação convexa de um número pré-definido de pontos extremos. Em ambos os modelos a região de operação admissível não-convexa é modelada através da divisão desta àrea em duas àreas convexas distintas. Testes computacionais efetuados com ambos os modelos para várias instâncias permitiram verificar a eficiência do modelo linear proposto. Este modelo permitiu obter as soluções ótimas do modelo não-linear com tempos computationais significativamente menores. Para além disso, ambos os modelos foram testados com e sem a inclusão de restrições de tomada e deslastre de carga, permitindo concluir que este tipo de restrições aumenta a complexidade do problema sendo que o tempo computacional exigido para a resolução do mesmo cresce significativamente.

Spray Pyrolysed Zinc Oxide Thin Films : Effects of Doping and Ion Beam Irradiation

In recent years scientists have made rapid and significant advances in the field of semiconductor physics. One of the most important fields of current interest in materials science is the fundamental aspects and applications of conducting transparent oxide thin films (TCO). The characteristic properties of such coatings are low electrical resistivity and high transparency in the visible region. The first semitransparent and electrically conducting CdO film was reported as early as in 1907 [1]. Though early work on these films was performed out of purely scientific interest, substantial technological advances in such films were made after 1940. The technological interest in the study of transparent semiconducting films was generated mainly due to the potential applications of these materials both in industry and research. Such films demonstrated their utility as transparent electrical heaters for windscreens in the aircraft industry. However, during the last decade, these conducting transparent films have been widely used in a variety of other applications such as gas sensors [2], solar cells [3], heat reflectors [4], light emitting devices [5] and laser damage resistant coatings in high power laser technology [6]. Just a few materials dominate the current TCO industry and the two dominant markets for TCO’s are in architectural applications and flat panel displays. The architectural use of TCO is for energy efficient windows. Fluorine doped tin oxide (FTO), deposited using a pyrolysis process is the TCO usually finds maximum application. SnO2 also finds application ad coatings for windows, which are efficient in preventing radiative heat loss, due to low emissivity (0.16). Pyrolitic tin oxide is used in PV modules, touch screens and plasma displays. However indium tin oxide (ITO) is mostly used in the majority of flat panel display (FPD) applications. In FPDs, the basic function of ITO is as transparent electrodes. The volume of FPD’s produced, and hence the volume of ITO coatings produced, continues to grow rapidly. But the current increase in the cost of indium and the scarcity of this material created the difficulty in obtaining low cost TCOs. Hence search for alternative TCO materials has been a topic of active research for the last few decades. This resulted in the development of binary materials like ZnO, SnO2, CdO and ternary materials like II Zn2SnO4, CdSb2O6:Y, ZnSO3, GaInO3 etc. The use of multicomponent oxide materials makes it possible to have TCO films suitable for specialized applications because by altering their chemical compositions, one can control the electrical, optical, chemical and physical properties. But the advantages of using binary materials are the easiness to control the chemical compositions and depositions conditions. Recently, there were reports claiming the deposition of CdO:In films with a resistivity of the order of 10-5 ohm cm for flat panel displays and solar cells. However they find limited use because of Cd-Toxicity. In this regard, ZnO films developed in 1980s, are very useful as these use Zn, an abundant, inexpensive and nontoxic material. Resistivity of this material is still not very low, but can be reduced through doping with group-III elements like In, Al or Ga or with F [6]. Hence there is a great interest in ZnO as an alternative of ITO. In the present study, we prepared and characterized transparent and conducting ZnO thin films, using a cost effective technique viz Chemical Spray Pyrolysis (CSP). This technique is also suitable for large area film deposition. It involves spraying a solution, (usually aqueous) containing soluble salts of the constituents of the desired compound, onto a heated substrate.

AN INVESTIGATION ON AIR-SEA INTERACTION IN THE INDIAN OCEAN DURING CONTRASTING MONSOONS

The study is undertaken with an objective to investigate the linkage between air-sea fluxes in the Indian Ocean and monsoon forcing. Since the monsoon activity is linked to fluxes, the variability of surface marine meteorological fields under the variable monsoon conditions is also studied. The very objective of the present study is to document various sea surface parameters of the Indian Ocean and to examine the anomalies found in them. Hence it is attempted to relate the anomaly to the variability of monsoon over India, highlighting the occasion of contrasting monsoon periods. The analysis of anomalies of surface meteorological fields such as SST, wind speed and direction, sea level pressure and cloud cover for contrasting monsoons are also studied. During good monsoon years, the pressure anomalies are negative indicating a fall in SLP during pre-monsoon and monsoon months. The interaction of the marine atmosphere with tropical Indian Ocean and its influence on ISMR continue to be an area of active research.

Development of a Technique for Estimating Rain Parameters from Rain Noise Measurements

The measurement of global precipitation is of great importance in climate modeling since the release of latent heat associated with tropical convection is one of the pricipal driving mechanisms of atmospheric circulation.Knowledge of the larger-scale precipitation field also has important potential applications in the generation of initial conditions for numerical weather prediction models Knowledge of the relationship between rainfall intensity and kinetic energy, and its variations in time and space is important for erosion prediction. Vegetation on earth also greatly depends on the total amount of rainfall as well as the drop size distribution (DSD) in rainfall.While methods using visible,infrared, and microwave radiometer data have been shown to yield useful estimates of precipitation, validation of these products for the open ocean has been hampered by the limited amount of surface rainfall measurements available for accurate assessement, especially for the tropical oceans.Surface rain fall measurements(often called the ground truth)are carried out by rain gauges working on various principles like weighing type,tipping bucket,capacitive type and so on.The acoustic technique is yet another promising method of rain parameter measurement that has many advantages. The basic principle of acoustic method is that the droplets falling in water produce underwater sound with distinct features, using which the rainfall parameters can be computed. The acoustic technique can also be used for developing a low cost and accurate device for automatic measurement of rainfall rate and kinetic energy of rain.especially suitable for telemetry applications. This technique can also be utilized to develop a low cost Disdrometer that finds application in rainfall analysis as well as in calibration of nozzles and sprinklers. This thesis is divided into the following 7 chapters, which describes the methodology adopted, the results obtained and the conclusions arrived at.

Sea Surface Temperature-Convection Relationship in Tropical Oceans with Special Emphasis to Intraseasonal Variability of Indian Summer Monsoon

The SST convection relation over tropical ocean and its impact on the South Asian monsoon is the first part of this thesis. Understanding the complicated relation between SST and convection is important for better prediction of the variability of the Indian monsoon in subseasonal, seasonal, interannual, and longer time scales. Improved global data sets from satellite scatterometer observations of SST, precipitation and refined reanalysis of global wind fields have made it possible to do a comprehensive study of the SST convection relation. Interaction of the monsoon and Indian ocean has been discussed. A coupled feedback process between SST and the Active-Break cycle of the Asian summer monsoon is a central theme of the thesis. The relation between SST and convection is very important in the field of numerical modeling of tropical rainfall. It is well known that models generally do very well simulating rainfall in areas of tropical convergence zones but are found unable to do satisfactory simulation in the monsoon areas. Thus in this study we critically examined the different mechanisms of generation of deep convection over these two distinct regions.The study reported in chapter 3 has shown that SST - convection relation over the warm pool regions of Indian and west Pacific oceans (monsoon areas) is in such a way that convection increases with SST in the SST range 26-29 C and for SST higher than 29-30 C convection decreases with increase of SST (it is called Waliser type). It is found that convection is induced in areas with SST gradients in the warm pool areas of Indian and west Pacific oceans. Once deep convection is initiated in the south of the warmest region of warm pool, the deep tropospheric heating by the latent heat released in the convective clouds produces strong low level wind fields (Low level Jet - LLJ) on the equatorward side of the warm pool and both the convection and wind are found to grow through a positive feedback process. Thus SST through its gradient acts only as an initiator of convection. The central region of the warm pool has very small SST gradients and large values of convection are associated with the cyclonic vorticity of the LLJ in the atmospheric boundary layer. The conditionally unstable atmosphere in the tropics is favorable for the production of deep convective clouds.

Development of solar air heaters & thermal energy storage system for drying applications in food processing industries

In the present work, the author has designed and developed all types of solar air heaters called porous and nonporous collectors. The developed solar air heaters were subjected to different air mass flow rates in order to standardize the flow per unit area of the collector. Much attention was given to investigate the performance of the solar air heaters fitted with baffles. The output obtained from the experiments on pilot models, helped the installation of solar air heating system for industrial drying applications also. Apart from these, various types of solar dryers, for small and medium scale drying applications, were also built up. The feasibility of ‘latent heat thermal energy storage system’ based on Phase Change Material was also undertaken. The application of solar greenhouse for drying industrial effluent was analyzed in the present study and a solar greenhouse was developed. The effectiveness of Computational Fluid Dynamics (CFD) in the field of solar air heaters was also analyzed. The thesis is divided into eight chapters.

Assessment of Energy Conservation in Indian Cement Industry and Forecasting of Co2 Emissions

Cement industry ranks 2nd in energy consumption among the industries in India. It is one of the major emitter of CO2, due to combustion of fossil fuel and calcination process. As the huge amount of CO2 emissions cause severe environment problems, the efficient and effective utilization of energy is a major concern in Indian cement industry. The main objective of the research work is to assess the energy cosumption and energy conservation of the Indian cement industry and to predict future trends in cement production and reduction of CO2 emissions. In order to achieve this objective, a detailed energy and exergy analysis of a typical cement plant in Kerala was carried out. The data on fuel usage, electricity consumption, amount of clinker and cement production were also collected from a few selected cement industries in India for the period 2001 - 2010 and the CO2 emissions were estimated. A complete decomposition method was used for the analysis of change in CO2 emissions during the period 2001 - 2010 by categorising the cement industries according to the specific thermal energy consumption. A basic forecasting model for the cement production trend was developed by using the system dynamic approach and the model was validated with the data collected from the selected cement industries. The cement production and CO2 emissions from the industries were also predicted with the base year as 2010. The sensitivity analysis of the forecasting model was conducted and found satisfactory. The model was then modified for the total cement production in India to predict the cement production and CO2 emissions for the next 21 years under three different scenarios. The parmeters that influence CO2 emissions like population and GDP growth rate, demand of cement and its production, clinker consumption and energy utilization are incorporated in these scenarios. The existing growth rate of the population and cement production in the year 2010 were used in the baseline scenario. In the scenario-1 (S1) the growth rate of population was assumed to be gradually decreasing and finally reach zero by the year 2030, while in scenario-2 (S2) a faster decline in the growth rate was assumed such that zero growth rate is achieved in the year 2020. The mitigation strategiesfor the reduction of CO2 emissions from the cement production were identified and analyzed in the energy management scenarioThe energy and exergy analysis of the raw mill of the cement plant revealed that the exergy utilization was worse than energy utilization. The energy analysis of the kiln system showed that around 38% of heat energy is wasted through exhaust gases of the preheater and cooler of the kiln sysetm. This could be recovered by the waste heat recovery system. A secondary insulation shell was also recommended for the kiln in the plant in order to prevent heat loss and enhance the efficiency of the plant. The decomposition analysis of the change in CO2 emissions during 2001- 2010 showed that the activity effect was the main factor for CO2 emissions for the cement industries since it is directly dependent on economic growth of the country. The forecasting model showed that 15.22% and 29.44% of CO2 emissions reduction can be achieved by the year 2030 in scenario- (S1) and scenario-2 (S2) respectively. In analysing the energy management scenario, it was assumed that 25% of electrical energy supply to the cement plants is replaced by renewable energy. The analysis revealed that the recovery of waste heat and the use of renewable energy could lead to decline in CO2 emissions 7.1% for baseline scenario, 10.9 % in scenario-1 (S1) and 11.16% in scenario-2 (S2) in 2030. The combined scenario considering population stabilization by the year 2020, 25% of contribution from renewable energy sources of the cement industry and 38% thermal energy from the waste heat streams shows that CO2 emissions from Indian cement industry could be reduced by nearly 37% in the year 2030. This would reduce a substantial level of greenhouse gas load to the environment. The cement industry will remain one of the critical sectors for India to meet its CO2 emissions reduction target. India’s cement production will continue to grow in the near future due to its GDP growth. The control of population, improvement in plant efficiency and use of renewable energy are the important options for the mitigation of CO2 emissions from Indian cement industries

Proporción de pacientes que presentan hipotermia peri operatoria durante un reemplazo articular de cadera, rodilla u hombro

Introducción: La hipotermia perioperatoria se ha documentado como factor de riesgo para el aumento de la morbimortalidad de los pacientes aumentando morbilidad miocárdica, riesgo de infección, pérdidas sanguíneas y tiempo de hospitalización. La aplicación de anestésicos toma relevancia ya que causa la pérdida de control central de la temperatura. Nuestro objetivo con este estudio fue describir la proporción de casos de hipotermia en la población sometida a un reemplazo articular durante un periodo de cuatro meses. Materiales y métodos: Se realizó un estudio de cohorte prospectivo. La población a estudio fueron los pacientes que fueron sometidos a un reemplazo total de cadera, rodilla u hombro. Se registró la temperatura central en el momento previo a la inducción anestésica, 30, 60 y 90 minutos después, al finalizar el procedimiento y al ingresar a recuperación. Se reportó el porcentaje de pacientes con hipotermia en cada tiempo. Resultados: Se analizaron en total 88 pacientes, el 55,7% fue llevado a cirugía de cadera, 39,7% de rodilla y 4,5% de hombro. El tipo de anestesia más utilizado fue general y la duración promedio de anestesia fue 164 minutos. La medición de la temperatura central se realizó en nasofaringe, esófago o tímpano. La proporción de pacientes que presentaron hipotermia en la inducción fue 21,6%, a 30 minutos 83%, a 60 minutos 73,9%, a 90 minutos 68,2%, al finalizar 59,1% y en recuperación 58%. Se realizó una prueba Chi cuadrado comparando las proporciones entre la inducción y los cinco periodos posteriores, se encontró que la proporción de pacientes con hipotermia en los cinco tiempos posteriores tuvo una diferencia estadísticamente significativa (p=0,00) comparada con la proporción de pacientes con hipotermia durante la inducción. Conclusión: En los pacientes sometidos a un reemplazo articular la hipotermia fue una condición prevalente posterior a la aplicación de los anestésicos sistémicos. Los dispositivos de calentamiento intraoperatorio usados actualmente son insuficientes para evitar la hipotermia, lo que indica concordancia con la literatura en cuanto a las recomendaciones de calentamiento perioperatorio, con énfasis en el precalentamiento, para prevenir la caída significativa de la temperatura y la morbimortalidad asociada.

African easterly wave variability and its relationship to Atlantic tropical cyclone activity

Automatic tracking of vorticity centers in European Centre for Medium-Range Weather Forecasts analyses has been used to develop a 20-yr climatology of African easterly wave activity. The tracking statistics at 600 and 850 mb confirm the complicated easterly wave structures present over the African continent. The rainy zone equatorward of 15 degreesN is dominated by 600-mb activity, and the much drier Saharan region poleward of 15 degreesN is more dominated by 850-mb activity. Over the Atlantic Ocean there is just one storm track with the 600- and 850-mb wave activity collocated. Based on growth/decay and genesis statistics, it appears that the 850-mb waves poleward of 15 degreesN over land generally do not get involved with the equatorward storm track over the ocean. Instead, there appears to be significant development of 850-mb activity at the West African coast in the rainy zone around (10 degreesN, 10 degreesW), which, it is proposed, is associated with latent heat release. Based on the tracking statistics, it has been shown that there is marked interannual variability in African easterly wave (AEW) activity. It is especially marked at the 850-mb level at the West African coast between about 10 degrees and 15 degreesN, where the coefficient of variation is 0.29. For the period between 1985 and 1998, a notable positive correlation is seen between this AEW activity and Atlantic tropical cyclone activity. This correlation is particularly strong for the postreanalysis period between 1994 and 1998. This result suggests that Atlantic tropical cyclone activity may be influenced by the number of AEWs leaving the West African coast, which have significant low-level amplitudes, and not simply by the total number of AEWs.

Life-cycle simulations of shallow frontal waves and the impact of deformation strain

The life-cycle of shallow frontal waves and the impact of deformation strain on their development is investigated using the idealised version of the Met Office non-hydrostatic Unified Model which includes the same physics and dynamics as the operational forecast model. Frontal wave development occurs in two stages; first, a deformation strain is applied to a front and a positive potential vorticity (PV) strip forms, generated by latent heat release in the frontal updraft; second, as the deformation strain is reduced the PV strip breaks up into individual anomalies. The circulations associated with the PV anomalies cause shallow frontal waves to form. The structure of the simulated frontal waves is consistent with the conceptual model of a frontal cyclone. Deeper frontal waves are simulated if the stability of the atmosphere is reduced. Deformation strain rates of different strengths are applied to the PV strip to determine whether a deformation strain threshold exists above which frontal wave development is suppressed. An objective method of frontal wave activity is defined and frontal wave development was found to be suppressed by deformation strain rates $\ge 0.4\times10^{-5}\mbox{s}^{-1}$. This value compares well with observed deformation strain rate thresholds and the analytical solution for the minimum deformation strain rate needed to suppress barotropic frontal wave development. The deformation strain rate threshold is dependent on the strength of the PV strip with strong PV strips able to overcome stronger deformation strain rates (leading to frontal wave development) than weaker PV strips.

An objective climatology of the dynamical forcing of polar lows in the Nordic seas

A new objective climatology of polar lows in the Nordic (Norwegian and Barents) seas has been derived from a database of diagnostics of objectively identified cyclones spanning the period January 2000 to April 2004. There are two distinct parts to this study: the development of the objective climatology and a characterization of the dynamical forcing of the polar lows identified. Polar lows are an intense subset of polar mesocyclones. Polar mesocyclones are distinguished from other cyclones in the database as those that occur in cold air outbreaks over the open ocean. The difference between the wet-bulb potential temperature at 700 hPa and the sea surface temperature (SST) is found to be an effective discriminator between the atmospheric conditions associated with polar lows and other cyclones in the Nordic seas. A verification study shows that the objective identification method is reliable in the Nordic seas region. After demonstrating success at identifying polar lows using the above method, the dynamical forcing of the polar lows in the Nordic seas is characterized. Diagnostics of the ratio of mid-level vertical motion attributable to quasi-geostrophic forcing from upper and lower levels (U/L ratio) are used to determine the prevalence of a recently proposed category of extratropical cyclogenesis, type C, for which latent heat release is crucial to development. Thirty-one percent of the objectively identified polar low events (36 from 115) exceeded the U/L ratio of 4.0, previously identified as a threshold for type C cyclones. There is a contrast between polar lows to the north and south of the Nordic seas. In the southern Norwegian Sea, the population of polar low events is dominated by type C cyclones. These possess strong convection and weak low-level baroclinicity. Over the Barents and northern Norwegian seas, the well-known cyclogenesis types A and B dominate. These possess stronger low-level baroclinicity and weaker convection.

Understanding African easterly waves: a moist singular vector approach

Moist singular vectors (MSV) have been applied successfully to predicting mid-latitude storms growing in association with latent heat of condensation. Tropical cyclone sensitivity has also been assessed. Extending this approach to more general tropical weather systems here, MSVs are evaluated for understanding and predicting African easterly waves, given the importance of moist processes in their development. First results, without initial moisture perturbations, suggest MSVs may be used advantageously. Perturbations bear similar structural and energy profiles to previous idealised non-linear studies and observations. Strong sensitivities prevail in the metrics and trajectories chosen, and benefits of initial moisture perturbations should be appraised. Copyright © 2009 Royal Meteorological Society

Climate entropy budget of the HadCM3 atmosphere-ocean general circulation model and of FAMOUS, its low resolution version

The entropy budget is calculated of the coupled atmosphere–ocean general circulation model HadCM3. Estimates of the different entropy sources and sinks of the climate system are obtained directly from the diabatic heating terms, and an approximate estimate of the planetary entropy production is also provided. The rate of material entropy production of the climate system is found to be ∼50 mW m−2 K−1, a value intermediate in the range 30–70 mW m−2 K−1 previously reported from different models. The largest part of this is due to sensible and latent heat transport (∼38 mW m−2 K−1). Another 13 mW m−2 K−1 is due to dissipation of kinetic energy in the atmosphere by friction and Reynolds stresses. Numerical entropy production in the atmosphere dynamical core is found to be about 0.7 mW m−2 K−1. The material entropy production within the ocean due to turbulent mixing is ∼1 mW m−2 K−1, a very small contribution to the material entropy production of the climate system. The rate of change of entropy of the model climate system is about 1 mW m−2 K−1 or less, which is comparable with the typical size of the fluctuations of the entropy sources due to interannual variability, and a more accurate closure of the budget than achieved by previous analyses. Results are similar for FAMOUS, which has a lower spatial resolution but similar formulation to HadCM3, while more substantial differences are found with respect to other models, suggesting that the formulation of the model has an important influence on the climate entropy budget. Since this is the first diagnosis of the entropy budget in a climate model of the type and complexity used for projection of twenty-first century climate change, it would be valuable if similar analyses were carried out for other such models.

Illustration of micro-scale advection using grid-pattern mini-lysimeters

Twenty-five small soil-filled perspex boxes arranged in a square, with dwarf sunflowers growing in them, were used to study micro-scale advection. Hydrological heterogeneity was introduced by applying two different amounts of irrigation water (low-irrigation, L, versus high-irrigation, H). The nine central boxes (4 H, 4 L and I bare box) were precision weighing lysimeters, yielding diurnal measurements of evaporation. After the onset of soil water stress, a large difference in latent heat flux (up to 4-fold) was observed between the lysimeters of the H and L treatments, mainly caused by large differences between H and L stomatal conductance values. This resulted in micro-advection, causing H soil-sunflower systems to evaporate well above equilibrium latent heat flux. The occurrence of micro-advective enhancement was reflected in large values of the Priestley-Taylor constant (often larger than 2.0) and generally negative values of sensible heat flux for the H treatment. (c) 2005 Elsevier B.V. All rights reserved.

Atmosphere feedbacks during ENSO in a coupled GCM with a modified atmospheric convection scheme

The too diverse representation of ENSO in a coupled GCM limits one’s ability to describe future change of its properties. Several studies pointed to the key role of atmosphere feedbacks in contributing to this diversity. These feedbacks are analyzed here in two simulations of a coupled GCM that differ only by the parameterization of deep atmospheric convection and the associated clouds. Using the Kerry–Emanuel (KE) scheme in the L’Institut Pierre-Simon Laplace Coupled Model, version 4 (IPSL CM4; KE simulation), ENSO has about the right amplitude, whereas it is almost suppressed when using the Tiedke (TI) scheme. Quantifying both the dynamical Bjerknes feedback and the heat flux feedback in KE, TI, and the corresponding Atmospheric Model Intercomparison Project (AMIP) atmosphere-only simulations, it is shown that the suppression of ENSO in TI is due to a doubling of the damping via heat flux feedback. Because the Bjerknes positive feedback is weak in both simulations, the KE simulation exhibits the right ENSO amplitude owing to an error compensation between a too weak heat flux feedback and a too weak Bjerknes feedback. In TI, the heat flux feedback strength is closer to estimates from observations and reanalysis, leading to ENSO suppression. The shortwave heat flux and, to a lesser extent, the latent heat flux feedbacks are the dominant contributors to the change between TI and KE. The shortwave heat flux feedback differences are traced back to a modified distribution of the large-scale regimes of deep convection (negative feedback) and subsidence (positive feedback) in the east Pacific. These are further associated with the model systematic errors. It is argued that a systematic and detailed evaluation of atmosphere feedbacks during ENSO is a necessary step to fully understand its simulation in coupled GCMs.