File:Earth like tidally locked ocean planet that orbits star teff 3500K S1 P1 1 1 1 1.png

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

import math import numpy as np import exoplasim as exo

timestep=30.0 years=40

Sol=1361.5

1. 1. basic input params

1. 1. startemp= 5778 # assumption Earth-like
   2. WARNING note M type star M5?
2. startemp=5778/2

1. 1. !!!!!

startemp=3500

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

flux=1367*1.0

rotationperiod=360 radius=1 ## rel to earth mass=1

eccentricity=0 obliquity=0

fixedorbit=True synchronous=True

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

1. landmap="Alderaan_surf_0129.sra"
2. topomap="Alderaan_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

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