File:Gliese 12b example simulation locked desert planet 1 r 1 1 1 1.png

https://github.com/alphaparrot/ExoPlaSim

insol    2219.24
mass     0.85
radius   0.96
density  5.36
vesc     10.56
gee_e    0.93
Patm     0.815051

2. 1. exoplasim planet runner
   2. python 3, exoplasim
4. basic run settings T21
5. 1. 29.10.2023 0000.0002a2

import math import numpy as np import exoplasim as exo

timestep=30.0 years=30

Sol=1361.5

1. 1. basic input params

startemp=3296 # assumption Earth-like

1. 1. 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

1. substellarlon=0

substellarlon=0 aquaplanet=False desertplanet=True 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 }

1. landmap="Alderaan_surf_0129.sra"
2. topomap="Alderaan_surf_0172.sra"
3. landmap="sermak_surf_0129.sra"
4. topomap="sermak_surf_0172.sra"
5. landmap="five_surf_0129.sra"
6. topomap="five_surf_0172.sra"

landmap=None topomap=None

1. relCO2=0.95
2. relN2=1-relCO2

relCO2=280e-6 relO2=0.21 relN2=1-relCO2-relO2

1. relCO2=0.95
2. relO2=0.045

1. 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) )

1. quit(-1)

planeta = exo.Model(workdir="planeta_run",modelname="Planeta",ncpus=1,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()