File:Faint young sun earth temperature july S0p7 P 1p0bar co2 95percent 1 1 1 1.png
Original file (1,488 × 992 pixels, file size: 469 KB, MIME type: image/png)
Captions
Summary
[edit]DescriptionFaint young sun earth temperature july S0p7 P 1p0bar co2 95percent 1 1 1 1.png |
English: Faint young sun and Earth temperature of July if S=0.7 sol P=1.0bar CO2=99 per cent |
Date | |
Source | Own work |
Author | Merikanto |
Python3 source code
-
- exoplasim planet runner
- python 3, exoplasim
- basic run settings T21
-
- 17.10.2023 0000.0002a
-
import math
import numpy as np
import exoplasim as exo
timestep=30.0
years=40
Sol=1361.5
- basic input params
- startemp= 5778 # assumption Earth-like
- WARNING note M type star M5?
- startemp=5778
- !!!!!
startemp=3500
- Trappist 1 d
- flux=1.12*1367
- rotationperiod=4.05
- radius=0.78 ## rel to earth
- mass=0.39
flux=1367*1.0
rotationperiod=360
radius=1 ## rel to earth
mass=1
eccentricity=0
obliquity=0
fixedorbit=True
synchronous=True
- substellarlon=180
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=True
vegaccel=1
seaice=True
wetsoil=True
year=rotationperiod
glaciers1= {
"toggle": True,
"mindepth":2,
"initialh":-1
}
- landmap="Alderaan_surf_0129.sra"
- topomap="Alderaan_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
- patm=geese*mass*radius*radius ## estimation
patm=1
gravity=geeseg
pressure=patm
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(patm,6) )
planeta = exo.Model(workdir="planet_run",modelname="Planet",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, takes long time ...")
planeta.run(years=years,crashifbroken=True)
planeta.exportcfg()
planeta.run(years=years,crashifbroken=True)
planeta.finalize("Planet",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.
- share alike – If you remix, transform, or build upon the material, you must distribute your contributions under the same or compatible license as the original.
File history
Click on a date/time to view the file as it appeared at that time.
Date/Time | Thumbnail | Dimensions | User | Comment | |
---|---|---|---|---|---|
current | 11:26, 29 October 2023 | 1,488 × 992 (469 KB) | Merikanto (talk | contribs) | Update | |
11:25, 29 October 2023 | 1,488 × 992 (469 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.