Comments (7)
@mcthreems - A quick clarification first - I'm a little confused as to what circumstance this will apply. Presumably we're thinking of Titan lakes, so that if the bucket depth is locally shallower than the prescribed mld then we adopt a lower value, with the lowest possible value being the land mld? My Earth perspective was confused at first, because over Earth's oceans you're imagining an essentially infinite bucket depth, so there's never any way it would dry out over ocean.
In terms of why this is happening, have you tried setting the land mld instead of the bucket depth as a first try? Does that still give you bad temperatures? That will hopefully tell you that it's that part of your code that is causing the problems.
Another quick one to check would be to try outputting the mld straight after this calculation, before any temperature changes are done. Would be good to check that the mld is a reasonable number (is it perhaps a strangely large number by accident?)
In terms of other possibilities, the model does tend to get a little upset when there are large contrasts in mld as this means large contrasts in temperatures -> high winds -> CFL errors. Might be worth trying lots of smaller timesteps to see if this is possibly the issue.
I guess you might have tried some of these things already. Hopefully there's something there that'll be useful.
from isca.
Also, is it failing on the first timestep, or is it failing after a few? I can imagine this kind of setup leading to a feedback loop of very rapid bucket emptying and evaporation if it's running in a dry atmosphere for a few timesteps.
from isca.
Have you tried initialising with a really deep and full bucket? That might help it run for a few timesteps where you know this condition won't be satisfied.
from isca.
Regarding your first question, we are running simulations with a pseudo-aquaplanet version of Earth in which we vary the initial bucket depth. The pseudo is because we have the model treat all grid cells as land cells, even though we can have up to 1000m of water in them. The way the model currently handles this is just applying the land prefactor for every cell. We wanted to avoid that land prefactor until a cell is either dry or would have a calculated heat capacity lower than the dry value, hence the system I described.
from isca.
After trying some of your suggestions, I finally realized the issue was a simple mistake on my part. When I was implementing the "where" functions with the bucket depth array, I assumed the dimensions would match the order of the output dimensions. As it turns out, they are reversed from the output dimensions, and thus the time dimension is the 3rd rather than the 1st. Once I switched the "size" function to the 3rd dimension, the model ran fine (up until later issues with crashing due to what looks like cold temperatures).
from isca.
Regarding those later issues, I'll paste in the error message here just for reference. If this constitutes a separate issues post we can close this one and I'll repost these on a new one.
2020-01-16 23:06:25,950 - isca - DEBUG - PE, location, Textr(curr,future): 3: 43 1 4 100.326 99.942
2020-01-16 23:06:25,950 - isca - DEBUG - PE, location, Uextr(curr,future): 3: 43 1 4 -151.763 -152.305
2020-01-16 23:06:25,951 - isca - DEBUG -
2020-01-16 23:06:25,951 - isca - DEBUG - FATAL from PE 3: spectral_dynamics: temperatures out of valid range
from isca.
We figured out that the heat capacity calculation we were editing is in the init function, meaning it isn't done each time step. We've changed it to be located in the function proper such that it is done each time step. I'll close this for now.
from isca.
Related Issues (20)
- Tidying Up Docs HOT 1
- Compile Errors with GCC version 10> HOT 2
- Compiling issue HOT 2
- Main front-page Readme should be updated with up-to-date installation instructions HOT 2
- gfortran float to integer type conversion problems HOT 2
- x86_64-conda-linux-gnu-cc can't find liblto_plugin.so on Ubuntu HOT 3
- Min/Max Time Averaging
- Compile issue, error with make HOT 2
- sh 2.0.0 breaks Isca installation HOT 3
- utime no longer a class in cftime
- Compiling ISCA on Mac laptop HOT 4
- Be careful with using the frierson_solar_rad option with Socrates.
- Compiling ISCA for the first time, issue with helpers.py HOT 10
- Is it possible to use Isca model to simulate atmoshperic responses to a diabatic heating source with a given seasonal background state?
- Should topography be blocking out portions of the atmosphere? HOT 4
- Compile Errors when running test (held_suarez_test_case.py) HOT 5
- Typo
- Socrates Test Case run failure
- Unable to achieve an Earth-like tsurf gradient using SOCRATES and and frierson_solar_rad = True HOT 2
- Segmentation fault caused by unset `tau` in "dry_convection_nml"
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from isca.