File:Gliese 12b locked lot of land little sea example simulation 1 r 1 1 1 1.png

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

import math import numpy as np import exoplasim as exo

timestep=30.0 years=30

Sol=1367

1. 1. basic input params

startemp= 3296 # 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

flux=1.63*1367 rotationperiod=12.76140 radius=0.96 ## rel to earth

1. mass=0.85

mass=math.pow(radius, 3.7)

eccentricity=0 obliquity=0

fixedorbit=True synchronous=True substellarlon=180

1. substellarlon=0

aquaplanet=False 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 }

1. landmap="Alderaan_surf_0129.sra"
2. topomap="Alderaan_surf_0172.sra"

landmap="g3_surf_0129.sra" topomap="g3_surf_0172.sra"

1. landmap=None
2. 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. patm=geese*mass*radius*radius ## estimation

patm1=geese*mass*radius*radius ## estimation patm2=math.pow(radius, 2.4)

patm=(patm1+patm2)/2

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="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()