Energias de ligação do estado fundamental de doadores rasos em super-redes são consideradas teoricamente com o auxílio de um procedimento variacional que leva em conta a mistura do contínuo de estados da minibanda à qual o mesmo está associado. Os cálculos são realizados para um grande número de parâmetros de super-rede e qualquer posição da impureza na mesma. É mostrado que a dependência da energia de ligação com os vários parâmetros envolvidos pode ser completamente explicada em termos de um modelo simples unidimensional (tight-binding) onde a largura da respectiva minibanda de condução e a energia de ligação no Limite de poço isolado são os únicos parâmetros relevantes. A extrema concordância quantitativa entre as energias de ligação derivadas deste modelo e as obtidas pelo método variacional mais rigoroso vem enfatizar o papel fundamental desempenhado pela largura de minibanda com o único parâmetro relevante ao se levar em conta os efeitos de tunelamento existentes nas super-redes. / A variational procedure which takes into account the mixing of a continuum of subband states has been used to investigate the binding energies of shallow donors in superlattices. The calculations where performed for a wide range of superlattices parameters and impurity positions. It is shown that the dependence of the binding energy upon the various superlattice parameters can be completely explained in terms of a simple onedimensional tight-binding model where the bandwidth of the respective conduction subband and the binding energy in the isolated quantum well are the only relevant parameters. The quantitative overall agreement between the binding energies derived from this model and those found variationally is excellent and emphasizes the fundamental role played by the bandwidth as the only relevant parameter accounting for the tunneling effects.

CoordenaÃÃo de AperfeÃoamento de Pessoal de NÃvel Superior / O dimensionamento de fundaÃÃes superficiais requer definiÃÃo de forma clara e objetiva, para o solo em estudo, da capacidade de carga e a previsÃo dos recalques a que se submete a estrutura, baseando-se nas propriedades mecÃnicas obtidas por ensaios de laboratÃrio e campo. O ensaio de placa à uma das formas mais seguras para se obter essas informaÃÃes, a partir de carregamento direto, que produz, em escala reduzida, o comportamento da futura fundaÃÃo. Os ensaios realizados nesta pesquisa foram executados em trÃs etapas: caracterizaÃÃo geotÃcnica (granulometria, determinaÃÃo da massa especÃfica, limite de consistÃncia, compactaÃÃo proctor normal e Ãndice de suporte califÃrnia (ISC)); ensaios especiais de laboratÃrio (adensamento e cisalhamento); e ensaios de campo (provas de carga). O objetivo deste trabalho à avaliar a partir da realizaÃÃo de provas de carga diretas com solo na umidade natural e inundado, se previsÃes de potencial de colapso e expansÃo, realizadas a partir de ensaios edomÃtricos simples e duplos sÃo capazes de dar indicaÃÃes concordantes. Foi observado que resultados de ensaios edometricos duplos, realizados no solo silto argiloso da cidade de IcÃ, proporcionaram estimativas mais concordantes para as situaÃÃes de colapso do solo. AlÃm disso, a tÃcnica da compactaÃÃo de solos problemÃticos se mostrou adequada para minorar, ou mesmo anular, problemas relacionados ao colapso e a expansÃo do solo estudado. / The design of shallow foundations requires defining clearly and objectively to the soil under study, the load capacity and the prediction of settlements that undergoes structure, based on the mechanical properties obtained from laboratory testing and field. The plaque assay is one of the safest ways to get this information from direct loading, which produces, in reduced scale, the behavior of the Future Foundation. The tests performed in this study were performed in three steps: geotechnical characterization (particle size, density determination, consistency limits, standard Proctor compaction and California bearing ratio (ISC)​​); Special laboratory tests (density and shear); and field trials (load tests). The objective of this study is to evaluate the evidence from conducting direct load with soil and flooded in natural moisture, if forecasts of potential collapse and expansion, made ​​from single and double edomÃtricos tests are able to give consistent evidence. It was observed that results of double edometricos trials, conducted in silty clay soil of the city Ico, provided more consistent estimates for situations of soil collapse. Furthermore, the technique of compression problematic soils was adequate to mitigate or even nullify, problems related to the collapse and expansion of the studied soil.

Introdução: Extubação é a retirada da via aérea artificial. Desmame é o processo de transição da ventilação artificial para a respiração espontânea. A falha na extubação exerce um impacto adverso na clínica dos pacientes, além disso, a reintubação acrescenta dias no CTI (Centro de Terapia Intensiva) e na permanência hospitalar. O Índice de respiração rápida e superficial (IRRS) é um parâmetro fidedigno nas referências sobre o tema, devendo ser incluído em protocolos de desmame. Objetivos: Avaliar a capacidade do IRRS como preditor de sucesso de extubação da ventilação mecânica invasiva (VMI), em uma população geral de pacientes críticos e também relacionar este índice a uma população de pacientes subdivididos em diferentes comorbidades, tais como: Doentes neurológicos, cardíacos, pulmonares e outras. Material e Métodos: Foram incluídos neste trabalho pacientes que necessitaram de ventilação mecânica invasiva por um período maior que 48 horas e que foram considerados aptos para extubação pela equipe médica assistente. Foram excluídos pacientes que não fossem extubados em até 6 horas após ventilação espontânea em tubo T, paciente com traqueostomia e com doença neuromuscular periférica. Após, o paciente era colocado em ventilação espontânea com tubo T , e em seguida era registrada a relação f/Vt (freqüência respiratória/volume corrente) sendo esta calculada no 1º (IRRS 1) e 30o (IRRS 2) minuto do teste, e calculado o Delta do IRRS. Os dados foram expressos em freqüência e percentual, média e desvio padrão e mediana (percentil 25-percentil 75), com nível de significância p<0,05 e curvas ROC (Receiver-operating Characteristic) para as variáveis do IRRS. Resultados: De um total de 504 pacientes extubados, 403 (80,0%) obtiveram sucesso no processo de desmame da VMI e 101 (20,0%) falharam. Quando comparamos o grupo sucesso com o grupo insucesso, encontramos os seguintes resultados: Média de idade 56 versus 61 (p<0,025) , IRRS 1 56,5 versus 68,6 (p=0,004), IRRS 2 59,8 versus 81,8 (p=< 0,001), Delta do IRRS 3,4 versus 13,2 (p=0,011), dias de CTI 12 versus 16 (p=0,001), dias de hospitalização 29,5 versus 33 (p=0,185), óbito 11% versus 29% (p< 0,001 ). Em relação às curvas ROC encontramos os seguintes valores para as áreas sobre a curva IRRS 1- 0,592 (IC 95% 0,527 – 0,656) p=0,005; IRRS 2- 0,619 (IC 95% 0,552 – 0,686) p=< 0,001; Delta do IRRS- 0,567 (IC 95% 0,499 – 0,635) p=0,039. Não encontramos diferenças significativas quando analisamos os dados dividindo os pacientes por subgrupos de comorbidades. Conclusão: A utilização do IRRS como preditor de sucesso, pode colaborar para a identificação dos pacientes que terão sucesso no processo de extubação da ventilação mecânica invasiva. Este índice não deve ser usado isoladamente, mas deve ser incluído em protocolos de desmame. Não encontramos diferenças significativas no comportamento do IRRS nas diferentes comorbidades avaliadas neste estudo. / Introduction: Extubation is the removal of artificial airway. Weaning is the process of transition from artificial to spontaneous breathing. Extubation failure has an adverse impact on patients' clinical evolution; furthermore, re-intubation increases ICU (Intensive Care Unit) and hospitalization time. Rapid Shallow Breathing Index (RSBI) is a reliable parameter in references on the theme and should be included in weaning protocols. Objectives: To assess the capacity of RSBI as a predictor of successful extubation from invasive mechanical ventilation (IMV) in a general population of critical patients, and relate that rate to a population of patients subdivided into different co-morbidities such as neurological, cardiac, pulmonary patients and others. Material and Methods: This study includes patients who required invasive mechanical ventilation for a period of more than 48 hours and were deemed fit for extubation by the assistant medical team. It excluded patients who were not extubated within 6 hours after spontaneous ventilation with a T-piece, patients with tracheostomy and those with peripheral neuromuscular disease. After extubation, the patient was placed in spontaneous ventilation with T-piece, and then the f/Vt (respiratory rate/Tidal volume) ratio was recorded, calculated at the 1st (RSBI 1) and 30th (RSBI 2) minute of the trial and calculated Delta RSBI. The data are expressed in frequency and percentage, average, standard deviation, and median (25th percentile-75th percentile), with a significance level of p<0.05 and ROC (Receiver-operating Characteristic) curves for the variables of RSBI. Results: From a total of 504 patients extubated, 403 (80.0%) achieved successful extubation from IMV and 101 (20.0%) failed. When comparing the success group and the failure group, the following results are found: average age 56 versus 61 (p<0.025); RSBI 1 56.5 versus 68.6 (p=0.004); RSBI 2 59.8 versus 81.8 (p=< 0.001); Delta RSBI 3.4 versus 13.2 (p=0.011); days in ICU 12 versus 16 (p=0.001); days of hospitalization 29.5 versus 33 (p=0.185); death 11% versus 29% (p< 0.001 ). In relation to ROC curves we find the following values for the areas under the curve RSBI 1- 0,592 (IC 95% 0,527 – 0,656) p=0,005; RSBI 2- 0,619 (IC 95% 0,552 – 0,686) p=< 0,001; Delta RSBI- 0,567 (IC 95% 0,499 – 0,635) p=0,039. No significant differences were found when the data was analyzed dividing the patients into co-morbidity subgroups. Conclusion: The use of the RSBI as a predictor of success, can help identify patients who will succeed in extubation from invasive mechanical ventilation. This index should not be used alone, but should be included in the weaning protocols. No significant differences were found in RSBI behavior for different co-morbidities assessed in this study.

Este estudo foi realizado com base em coletas sazonais no período de um ano (2010-2011), em três pontos na Lagoa do Peixe e um ponto na Lagoa do Ruivo, lagoa à montante, com objetivo de avaliar a existência de padrões espaciais e/ou sazonais na estrutura fitoplanctônica. Este ambiente lagunar subtropical está incluso no Parque Nacional da Lagoa do Peixe (sul do Brasil) e é conectado ao Oceano Atlântico por uma barra de areia localizada na porção central da laguna. A salinidade variou amplamente de condições limnéticas à euhalinas, sendo observada em média, salinidade mais baixa no inverno, e mais elevada no verão, exceto na Lagoa do Ruivo no verão, que mostrou tendência oposta. Em relação à espacialidade, a salinidade apresentou tendência à um gradiente crescente dos pontos extremos da laguna em direção a barra, apresentando menor variação nesta (mesohalino à euhalino). A comunidade fitoplanctônica esteve formada por sete classes com 136 táxons genéricos à infra-específicos. A densidade variou de 653 a 114.829 ind.mL-1, sendo os menores valores registrados próximos a barra. Os grupos abundantes em relação a densidade foram as cianobactérias e diatomáceas. A biomassa variou de 0,3 a 24,1 mg.L-1, sendo o menor valor registrado no inverno e o maior no outono, ambos próximo a barra. A elevada biomassa no outono ocorreu devido o predomínio de euglenofíceas. Espacialmente, a Lagoa do Ruivo (menor salinidade) apresentou maior riqueza. O grupo das diatomáceas apresentou a maior riqueza total. Análises de ordenação (ACP, RDA e ACC) para verificação de tendências sazonais e/ou espaciais foram realizadas, sugerindo existência de tendências a padrões, tanto sazonais, quando espaciais na estrutura e distribuição fitoplanctônica do sistema lagunar da Lagoa do Peixe, onde a salinidade foi a variável que contribui para estes gradientes. / This study was based on seasonal sampling during one year (2010-2011) at three points in Peixe Lagoon and one point in Ruivo Lagoon, located upstream, to evaluate the spatial and seasonal patterns of phytoplankton communities. These subtropical lagoon environments are included in the National Park of Peixe Lagoon (southern Brazil) and are connected to the Atlantic Ocean by a sand bar located at the central portion of the lagoon. Salinity ranged widely from limnetic to euhaline conditions, on average from lower salinity in winter to higher in summer being observed, except at Ruivo Lagoon in summer, which showed an opposite trend. Regarding spatiality, salinity tended to an increasing gradient from extreme points of the lagoon toward the bar, with less variation at this points (mesohaline to euhaline). The phytoplankton taxonomic community was composed of seven classes with 136 generic to infraspecific taxa. Density ranged from 653 to 114,829 ind.mL-1, with the lowest values recorded near the bar. The most abundant groups in density were cyanobacteria and diatoms. Biomass ranged from 0.3 to 24.1 mg.L-1, with the lowest value recorded in winter and higher, in the fall, both near the bar. The high biomass in fall was due to the predominance of euglenoids. Spatially, the Ruivo Lagoon (lower salinity) showed greater richness. Diatoms groups showed the highest total richness. Ordination analysis (PCA, RDA and CCA) were performed, suggesting the existence seasonal and spatial patterns in relation to the structure and distribution in phytoplankton of this lagoon system, where salinity was the variable that most contributed to these gradients.

Este trabalho aborda a solução numérica das equações lineares de águas rasas. O método dos elementos finitos e utilizado para a discretização espacial das equações que modelam o problema, e para a discretização temporal, o esquema semi-implícito de Crank-Nicolson é empregado. Além de alguns conceitos comuns quando se trabalha com escoamentos geofísicos, são descritas também a formulação das equações de águas rasas, sua linearização e uma solução analítica para um caso onde o parâmetro de Coriolis é nulo. A escolha adequada de pares de elementos finitos é a principal dificuldade quando se trabalha com esse método para a resolução da equação de águas rasas. Assim, é discutido o uso de quatro pares de elementos finitos e técnicas de estabilização para contornar o surgimento de modos espúrios na solução discreta. Os resultados numéricos são realizados com auxílio do software FreeFem++, onde se pode notar a capacidade dos pares de elementos de reproduzirem o escoamento, através da solução discreta, além das propriedades de conservação de massa e energia de cada discretização. / This work is about the numerical solution of the linear shallow water equations. The finite element method is used for spatial discretization of the equations that model the problem and for the time discretization the semi-implicit Crank-Nicolson scheme is used. Besides the concepts related to geophysical flows, the formulation of the shallow water equations, their linearization and an analytical solution for a case where the Coriolis parameter is zero are also described. The appropriate choice of a pair of finite elements is the main difficulty when working with this method for solving the shallow water equations. The use of four pairs of finite elements and stabilization techniques to circumvent the appearance of spurious modes in the discrete solution are discussed. The numerical results are performed using the software FreeFem++, where one can notice the ability of the elements to represent the discrete solution and mass and energy conservation of each discretization.

Made available in DSpace on 2014-10-09T12:32:48Z (GMT). No. of bitstreams: 0 / Made available in DSpace on 2014-10-09T14:08:43Z (GMT). No. of bitstreams: 1 03530.pdf: 1830989 bytes, checksum: b0ff295a4249c851927486f820ebe929 (MD5) / Tese (Doutoramento) / IPEN/T / Instituto de Pesquisas Energeticas e Nucleares - IPEN/CNEN-SP

nÃo hà / Identificar se as comunidades estÃo estruturadas de forma aleatÃria ou com base em interaÃÃes biolÃgicas trata-se de um dos aspectos mais importantes em ecologia . Neste sentido a identificaÃÃo dos eventos envolvidos em sucessÃo de espÃcies pode contribuir de forma significativa. As lagoas temporÃrias do semi-Ãrido brasileiro sÃo ambientes Ãmidos temporÃrios, com caracterÃsticas peculiares que permitem uma rÃpida sucessÃo de espÃcies, uma vez que, apresentam variaÃÃes ambientais, fase Ãmida e seca, responsÃveis pela seleÃÃo de espÃcies de acordo com suas formas de vida. Estes ambientes apresentam comunidades, dentre as quais se destacam as macrÃfitas aquÃticas, adaptadas aos perÃodos sazonais de inundaÃÃo. Assim, teve-se como objetivo descrever as variaÃÃes na estrutura de uma comunidade de macrÃfitas aquÃticas, em uma lagoa temporÃria, durante seu perÃodo de inundaÃÃo, relacionando as alteraÃÃes na dominÃncia e composiÃÃo de espÃcies com variÃveis ambientais. Como tambÃm, levantar consideraÃÃes sobre a co-ocorrÃncia de espÃcies atravÃs do c-score. As amostragens de dados foram realizadas quinzenalmente atravÃs do mÃtodo dos pontos a fim de se obter dados quanto a estratificaÃÃo, composiÃÃo e dominÃncia de espÃcies na comunidade , totalizando13 amostragens. AlÃm dos dados da vegetaÃÃo foram tambÃm obtidos dados de pH, condutividade, turbidez, alcalinidade, nitrogÃnio total, nitrato, fosfato total, ortofosfato, oxigÃnio dissolvido e demanda bioquÃmica de oxigÃnio da Ãgua. Estes foram relacionados com os dados de dominÃncia de espÃcies mostrando que as variÃveis ambientais afetam de formas diferentes as espÃcies estudadas. O c-score foi calculado tanto para cada coleta individualmente como para todo o perÃodo de inundaÃÃo. Nos estÃgios do perÃodo de inundaÃÃo predominaram as interaÃÃes biolÃgicas na comunidade e considerando todo o perÃodo o c-score mostra uma ocupaÃÃo das espÃcies de acordo com um padrÃo de sucessÃo / Identify whether communities are structured randomly or based on bio logical interactions is one of the most important aspect in ecology. In this sense the identification of events involved in the succession of species can significantly contribute. The temporary ponds in the semi-arid environments are temporary wet, with distinctive features that allow a rapid species succession, as they present environmental variations, wet and dry phase, responsible for the selection of species according to their ways of life. These environments present communities, which emphasize the aquatic macrophytes, adapted to seasonal flooding periods. So, it came up to describe the variations in the structure of a aquatic macrophytes community in a temporary pond, during the flooding period, relating the changes in dominance and species composition to environmental variables. As well, raising considerations of co-occurrence of species through the c-score. Data samples were made fortnightly by the method of points in order to obtain data on the stratification, composition of dominant species in the community, totaling 13 samples. In addition to the vegetation data were also obtained data for pH, conductivity, turbidity, alkalinity, total nitrogen, nitrate, total phosphate, orthophosphate, dissolved oxygen and biochemical oxygen demand of water. These were related to the data of dominance of species showing that environmental variables affect different forms of the species studied. The c-score was both calculated for each sample individually and for the duration of flooding. For the each sample separate ly prevailed biological interactions in the community and considering the whole period the c-score shows an occupation of the species according to a pattern of succession

Neste trabalho de pesquisa foi realizado um estudo do método de deconvolução visando melhor adequação à situação encontrada na escala de investigação rasa para a estimativa da forma de onda e da resposta impulsiva da Terra. Procedimentos determinísticos e estatísticos (métodos cegos) foram avaliados sobre dados sintéticos e dados reais adquiridos com fontes de impacto e com uma fonte pseudo-aleatória. / In this research work a study of the method of deconvolution was conducted in order to improve the adequacy to the shallow subsurface scale of investigation to the estimate of the seismic wavelet and of the earth impulse response. Deterministic and statistical (blind) procedures were evaluated over synthetic and real data acquired with impact sources and a pseudo-random source.

Orientador: Jose Alexandre Diniz / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Eletrica e de Computação / Made available in DSpace on 2018-08-21T20:16:28Z (GMT). No. of bitstreams: 1 Biasotto_Cleber_D.pdf: 3795677 bytes, checksum: e4e2355f73ec8fb2ad7677cdf17f5daf (MD5) Previous issue date: 2012 / Resumo: Atualmente para atender à necessidade de fabricação dos sensores, dispositivos eletrônicos e circuitos integrados com dimensões micro e nanométricas, novos processo de custos e de thermal budgets reduzidos são necessários. Este trabalho apresenta o desenvolvimento de alguns destes novos processos alternativos para aplicação nesta fabricação. O trabalho está dividido em quatro partes: a primeira parte apresenta a obtenção e a caracterização de filmes isolantes de nitreto de silício para aplicação em microsensores, tais como o sensor de pressão. Estes filmes foram obtidos sobre substratos de Si em baixa temperatura (20°C) utilizando-se um reator de plasma do tipo ECRCVD (Electron Cyclotron Resonance - Chemical Vapor Deposition). Normalmente reatores do tipo Low Pressure Chemical Vapor Deposition (LPCVD) ou Plasma Enhanced - CVD (PECVD) em temperaturas maiores que 600ºC e 250ºC, respectivamente, são utilizados para essa aplicação. A caracterização dos plasmas ECR, que foram usados para as deposições dos nitretos, e a fabricação de membranas suspensas com estes filmes são apresentadas. A segunda parte apresenta a fabricação e a caracterização de diodos p+-n fabricados em camada de SiGe crescida por LPCVD sobre substrato de Si. Processo este alternativo em substituição aos executados em reatores epitaxiais de alto custo. Na terceira parte deste trabalho, é apresentado o desenvolvimento de processos em baixa temperatura para aplicação em diodos e tecnologia MOS (Metal-Oxide-Semiconductor). São apresentadas a fabricação e a caracterização elétrica dos capacitores MOS, utilizando as tecnologias ALD (Atomic Layer Deposition) e ICP (Inductively Coupled Plasma) para a obtenção em baixa temperatura dos dielétricos high-k de Al2O3 e SiON de porta MOS, respectivamente. Na quarta parte, são apresentadas também, a fabricação e a caracterização elétrica de diodos n+p utilizando a tecnologia de recozimento a laser. O desenvolvimento de capacitores MOS e diodos possibilitou a fabricação (usando processos em baixa temperatura (?400oC)) de n- e p-MISFETs (Metal- Insulator-Semiconductor Field Effect Transistors), como também a fabricação de um protótipo de transistor MOS de alta velocidade baseado em silício germânio chamado D-DotFET (Disposable Dot Field Effect Transistor). Os processos alternativos desenvolvidos nesta tese apresentam um enorme potencial para aplicação nas próximas gerações de dispositivos CMOS (Complementary Metal Oxide Semiconductor) de dimensões sub-22 nm / Abstract: Nowadays, to attend the needs of the fabrication of sensors, electronic devices and integrated circuits with dimensions of micro and nanometrics, new processes of reduced costs and thermal budgets are needed. This work presents the development of some of these alternative processes for this fabrication. This work is divided in four parts: the first part presents the synthesis and characterization of insulating films of silicon nitride for application in microsensors, such as pressure sensors. These films were deposited on Si substrates at low temperature (20°C) using an ECR-CVD (Electron Cyclotron Resonance - Chemical Vapor Deposition) plasma reactor. Normally, Low Pressure Chemical Vapor Deposition (LPCVD) or Plasma Enhanced CVD - (PECVD) reactors are used for this application with high temperature process higher than 600oC and 250oC, respectively. The characterization of ECR plasmas, which were used to get the silicon nitrides, and the fabrication of suspended membranes based on these nitrides are presented. The second part presents the fabrication and the characterization of p+-n silicon germanium (SiGe) diodes fabricated on SiGe layers, which were grown by LPCVD on Si substrate. The grown of SiGe layers by LPCVD is an alternative process to replace the high cost of epitaxial reactors. In the third part of this work is presented the development of low-temperature processes for application in diodes and MOS (Metal-Oxide-Semiconductor) technology. The fabrication at low temperature and electrical characterization of MOS capacitors, using technologies as: ALD (Atomic Layer Deposition) and ICP (Inductively Coupled Plasma) to get the Al2O3 and SiON high-k gate dielectrics of MOS capacitors are presented, respectively. In the four part, the fabrication and electrical characterization of n+-p diodes using the process of laser annealing are presented as well. The development of MOS capacitors and diodes have become feasible the fabrication (using processes at low temperature (? 400oC)) of n- and p-MISFETs (Metal- Insulator -Semiconductor Field Effect Transistors) and also the fabrication of a high speed MOS transistor prototype based on silicon germanium named D-DotFET (Disposable Dot Field Effect Transistor). In conclusion, the alternative processes developed in this thesis have shown to be a huge potential for application in next generations of CMOS (Complementary Metal Oxide Semiconductor) devices with sub- 22 nm dimensions / Doutorado / Eletrônica, Microeletrônica e Optoeletrônica / Doutor em Engenharia Elétrica

Ice sheets and glaciers constitute an enormous water storage, currently corresponding to a potential sea level rise of almost 70 meters if all ice was to melt completely. The ice sheets are dynamic components of the global climate system and numerical modeling is a useful tool that can help us understand and predict how the ice sheets develop. The most accurate model available for ice sheets is given by the Stokes equations, but to solve them for a real ice sheet on a relevant time scale would be way too computationally costly. Instead approximations of the Stokes equations are used such as the Shallow Ice Approximation (SIA). The SIA is valid for areas where the aspect ratio, the ice thickness divided by the horizontal extent of the ice, is small. In this project equations for temperature and surface evolution were implemented in a Matlab version of SIA. The model already had algorithms implemented for computation of stresses, velocities and pressures for an ice sheet with fixed geometry and temperature. Implementation of temperature and free surface equations also made the problem time-dependent. The result was evaluated by solving a simple test problem and comparing the solution to a full Stokes solution obtained with the code ElmerIce. The SIA solution was closer to the Stokes solution when the aspect ratio ε and slope α were decreased simultaneously such that ε=arctanα, but a similar improvement was also obtained when only the slope was decreased. The differences between the two solutions were satisfyingly small for both temperature, surface location and velocities for an aspect ratio of ε= 7.8 10−4 and ε=arctanα.