663573f911
Added adiabatic motion to convection
...
- Converted temperature to potential temperature for thermal advection
- Added optional argument for vertical smoothing in 3D FFT function
- Small changes to plotting
2020-09-16 20:57:03 +01:00
a124787a05
Smoothing added
...
Added FFT function to smooth velocity and temperature fields.
Also minor graphical improvements.
Added surface friction
Clamped velocities over pole
2020-09-09 21:00:57 +01:00
0cc169b000
Post stream fix
...
Minor plotting addition and improved efficiency of vertical derivative
2020-09-02 21:47:19 +01:00
b8a7f94e2c
Cythonised all the things!
...
Used cython to improve speed by factor of 10, fixed bug in advection, decamped functions to claude libraries
2020-09-02 21:26:31 +01:00
6135e0b849
Big update 26/8/20
...
Major fixes:
- Tweaked radiation scheme to fix representation of equations
- Added vertical motion assuming perturbation to hydrostatic equilibriu
- Implemented seasons / axial tilt
- Added stratospheric heating
Minor changes:
- Removed old unused code
- Added frictional forces throughout atmosphere
- Tweaked optical depth parameterisation
- Tweaked text output
- Removed vertical advection (for now)
2020-08-27 09:24:39 +01:00
23b65ae1a6
Fixed radiation scheme
...
- Fixed lowest level of radiation scheme
- Added plot functionality
- Removed vertical motion
- Added w to save file
2020-08-19 20:54:05 +01:00
3fd7ce1892
Added grey radiation scheme
...
Implemented 3D grey radiation and variable optical depth with latitude.
2020-08-12 20:53:41 +01:00
6e197187fc
Generalising height
...
Allowing for variable maximum height of the model, with initial conditions of density and temperature read from US standard atmosphere.
2020-08-10 13:55:54 +01:00
c29c0aa741
Added plotting functionality
...
Added 2D plots of temperature and velocity for each level and a zonal mean representation
2020-07-08 20:59:08 +01:00
4a5aa8abbf
3D atmosphere!
...
Added multiple layers to the atmosphere. These have adjacent radiation exchange but not cross-level exchange, and no vertical momentum calculated. Also tweaked the absorption in different levels of the atmosphere. Differential sections of the code have been generalised for 3D fields
2020-07-01 21:18:16 +01:00
72247c418d
Update toy_model.py
...
Stream update from 24/6/2020
- Fixed advection
- Fixed coriolis
- Added spinup
- Add time functionality
2020-06-24 21:08:16 +01:00
dd1cd70ee6
Update toy_model.py
...
Increased spin-up time, decreased resolution
2020-06-17 18:08:38 +01:00
d6f83e4558
Update toy_model.py
...
Increased heat capacity of all components, resulting in more realistic physics (i.e. no immediate blowup). Also implemented toggle for velocity calculations (now largely redundant). Various other fiddles.
2020-06-17 17:17:11 +01:00
df23337825
Update toy_model.py
...
Corrected geometric mistake with Coriolis, tidied up optional oceans, changed the order of operation such that all thermal processes take place before updating the pressure
2020-06-15 22:10:59 +01:00
bb6c00a2e5
Update toy_model.py
...
Fixed geometry of radiative forcing, testing problems with derivatives
2020-06-15 18:24:09 +01:00
59393eecbb
Add files via upload
...
First version of toy atmospheric model, currently without working advection!
2020-06-15 16:30:56 +01:00
Simon Clark