File:Earth like 90 degrees obliquity desert planet july tsa 1 1 1 1b.png

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

import math import numpy as np import exoplasim as exo

timestep=30.0 years=200

Sol=1361.5

1. 1. basic input params

startemp=5778 # assumption Earth-like

1. 1. Trappist 1 d

1. flux=1.12*1367
2. rotationperiod=4.05
3. radius=0.78 ## rel to earth
4. mass=0.39
5. Trappist 1e
6. insol 884.98
7. mass 0.69
8. radius 0.92
9. density 4.89
10. vesc 9.69
11. gee_e 0.82
12. Patm 0.647368

flux=1361.5

pressure=1

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

1. relCO2=280e-6
2. relO2=0.21
3. relN2=1-relCO2-relO2

1. relCO2=0.95 ## lots of co2
2. relO2=0.001 ##
3. relN2=1-relCO2-relO2

relCO2=280e-6

relO2=0.21 relN2=1-relCO2-relO2

radius=1 ## rel to earth mass=1

rotationperiod=1 ##

1. eccentricity=0.0167

eccentricity=0.0

1. obliquity=23.44

obliquity=90

aquaplanet=False desertplanet=False fixedorbit=True synchronous=False

1. substellarlon=0

substellarlon=0

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=360

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

density=mass/np.power(radius,3)*5.519 vesc=np.sqrt(mass/radius)*11.186 geese=mass/(radius*radius) geeseg=geese*9.80665 gravity=geeseg

if(pressure==0): patm1=geese*mass*radius*radius ## estimation patm2=math.pow(radius, 2.4) 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=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()