File:Gliese 12b example simulation locked ocean planet 1 r 1 1 1 1.png
Original file (1,040 × 704 pixels, file size: 435 KB, MIME type: image/png)
Captions
Summary
[edit]DescriptionGliese 12b example simulation locked ocean planet 1 r 1 1 1 1.png |
English: Gliese 12b example simulation as locked ocean planet |
Date | |
Source | Own work |
Author | Merikanto |
Exoplasim
https://github.com/alphaparrot/ExoPlaSim
Python3 exoplasim source code. Exoplasim compiled with pip
pip install git+https://github.com/alphaparrot/ExoPlaSim
Anaconda python 3.8 number of processors 1.
-
- exoplasim planet runner
- python 3, exoplasim
- basic run settings T21
-
- 29.10.2023 0000.0002a2
-
import math
import numpy as np
import exoplasim as exo
timestep=30.0
years=30
Sol=1361.5
- basic input params
startemp=3296 # assumption Earth-like
- Gliese 12 b
flux=1.63*1361.5
rotationperiod=12.76140
radius=0.958 ## rel to earth
mass=math.pow(radius, 3.7)
eccentricity=0
obliquity=0
fixedorbit=True
synchronous=True
- substellarlon=0
substellarlon=0
aquaplanet=True
desertplanet=False
vegetation=2 ## 0 no 1 static 2 dynamic
stormclim=False
aerosol=False
co2weathering=False
outgassing=0
evolveco2=False
maxsnow=-1 ## False ?
ozone=False
vegaccel=1
seaice=False
wetsoil=False
year=rotationperiod
glaciers1= {
"toggle": True,
"mindepth":2,
"initialh":-1
}
- landmap="Alderaan_surf_0129.sra"
- topomap="Alderaan_surf_0172.sra"
- landmap="sermak_surf_0129.sra"
- topomap="sermak_surf_0172.sra"
- landmap="five_surf_0129.sra"
- topomap="five_surf_0172.sra"
landmap=None
topomap=None
- relCO2=0.95
- relN2=1-relCO2
relCO2=280e-6
relO2=0.21
relN2=1-relCO2-relO2
- relCO2=0.95
- relO2=0.045
- relN2=1-relCO2-relO2
density=mass/np.power(radius,3)*5.519
vesc=np.sqrt(mass/radius)*11.186
geese=mass/(radius*radius)
geeseg=geese*9.80665
patm1=geese*mass*radius*radius ## estimation
patm2=math.pow(radius, 2.4)
gravity=geeseg
pressure=(patm1+patm2)/2
pN2=pressure*relN2
pCO2=pressure*relCO2
pO2=pressure*relO2
print(" insol ",round(flux,2) )
print(" mass ",round(mass,2) )
print(" radius ",round(radius,2) )
print(" density ",round(density,2) )
print(" vesc ",round(vesc,2) )
print(" gee_e ",round(geese,2) )
print(" Patm ",round(pressure,6) )
- quit(-1)
planeta = exo.Model(workdir="planeta_run",modelname="Planeta",ncpus=4,resolution="T21",outputtype=".nc")
planeta.configure(year=year, wetsoil=wetsoil,pO2=pO2, ozone=ozone, vegaccel=vegaccel, seaice=seaice, glaciers=glaciers1, maxsnow=maxsnow, evolveco2=evolveco2, outgassing=outgassing, co2weathering=co2weathering, vegetation=vegetation, stormclim=stormclim, aerosol=aerosol,landmap=landmap, topomap=topomap,startemp=startemp, flux=flux, eccentricity=eccentricity,obliquity=obliquity,fixedorbit=fixedorbit,synchronous=synchronous,substellarlon=substellarlon,rotationperiod=rotationperiod,radius=radius,gravity=gravity,aquaplanet=aquaplanet,desertplanet=desertplanet,pN2=pN2,pCO2=pCO2,timestep=timestep,snapshots=False,physicsfilter="gp|exp|sp")
print(" Run ...")
planeta.run(years=years,crashifbroken=True)
planeta.exportcfg()
planeta.run(years=years,crashifbroken=True)
planeta.finalize("Planeta",allyears=True,keeprestarts=True)
planeta.save()
Licensing
[edit]- You are free:
- to share – to copy, distribute and transmit the work
- to remix – to adapt the work
- Under the following conditions:
- attribution – You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.
File history
Click on a date/time to view the file as it appeared at that time.
Date/Time | Thumbnail | Dimensions | User | Comment | |
---|---|---|---|---|---|
current | 12:22, 27 May 2024 | 1,040 × 704 (435 KB) | Merikanto (talk | contribs) | Uploaded own work with UploadWizard |
You cannot overwrite this file.
File usage on Commons
There are no pages that use this file.