alphaparrot / ExoPlaSim

Exoplanet Planet Simulator (PlaSim extended for different planet types (including tidally-locked) and evolution on geological timescales--glaciers and carbon cycle)
GNU General Public License v2.0
65 stars 12 forks source link

Example given in readme doesn't work on opensuse 15.5 #8

Open pkp24 opened 1 year ago

pkp24 commented 1 year ago

The example didn't work. This did work. `import exoplasim as exo

toi700d = exo.Model(workdir="toi700d_run",modelname="TOI-700d",ncpus=4,resolution="T21",outputtype=".npz") toi700d.configure(startemp=3480.0, flux=1167.0, # Stellar parameters eccentricity=0.,obliquity=0.,fixedorbit=True, # Orbital parameters synchronous=True,rotationperiod=37.426, # Rotation radius=1.19,gravity=11.9,aquaplanet=True, # Bulk properties pN2=1.47(1-360e-6),pCO2=1.47360e-6,ozone=False, # Atmosphere timestep=30.0,snapshots=720,physicsfilter="gp|exp|sp") # Model dynamics toi700d.exportcfg() toi700d.run(years=10,crashifbroken=True)`

alphaparrot commented 1 year ago

You're going to have to give a little more info on what you originally tried and what happened when it didn't work.

Specifically, please include the exact code that didn't run correctly, and the error messages you got when it crashed.

On Sun, Oct 22, 2023, 7:58 AM pkp24 @.***> wrote:

The example didn't work. This did work. `import exoplasim as exo

toi700d = exo.Model(workdir="toi700d_run",modelname="TOI-700d",ncpus=4,resolution="T21",outputtype=".npz") toi700d.configure(startemp=3480.0, flux=1167.0, # Stellar parameters eccentricity=0.,obliquity=0.,fixedorbit=True, # Orbital parameters synchronous=True,rotationperiod=37.426, # Rotation radius=1.19,gravity=11.9,aquaplanet=True, # Bulk properties pN2=1.47(1-360e-6),pCO2=1.47360e-6,ozone=False, # Atmosphere timestep=30.0,snapshots=720,physicsfilter="gp|exp|sp") # Model dynamics toi700d.exportcfg() toi700d.run(years=10,crashifbroken=True)`

— Reply to this email directly, view it on GitHub https://github.com/alphaparrot/ExoPlaSim/issues/8, or unsubscribe https://github.com/notifications/unsubscribe-auth/ADRQKQF7UA6SFBRLK73L2E3YAUC5JAVCNFSM6AAAAAA6K3CQVWVHI2DSMVQWIX3LMV43ASLTON2WKOZRHE2TKOBXGMYDAMQ . You are receiving this because you are subscribed to this thread.Message ID: @.***>

pkp24 commented 1 year ago

It seemed like it was failing as though it wasn't stable, and the fact that a preset tidally locked planet worked and the planet without any changes didn't seems to support that.

import exoplasim as exo mymodel = exo.Model(workdir="mymodel_testrun",modelname="mymodel",resolution="T21",layers=10,ncpus=4) mymodel.configure() mymodel.exportcfg() mymodel.run(years=100,crashifbroken=True) mymodel.finalize("mymodel_output")

What happened when I ran it: 7HKzCHg8.txt

u1145 commented 1 year ago

I just ran into the same issue when I tried running the example. The messages I got look very similar to what pkp24 got, with the warnings about numpy and things being out of bounds. I tried running it and now I only get the last part: output.txt