7f56b3c96d
Stable; runs without crashing when one filter is enabled; crashes if both are enabled.
Adiv Paradise2020-09-17 17:38:53 -04:00
c4366f3657
Added tunings to better-capture extant storms and avoid generating too much output
Adiv Paradise2020-09-15 14:52:09 -04:00
353115faf3
Storm diagnostic computations turned off by default; bugfixes; slightly expanded namelist--by default use all available trigger metrics (most restrictive) for writing high-cadence output.
Adiv Paradise2020-08-22 20:01:34 -04:00
7998c9c902
Everything runs now and produces numbers that looks halfway reasonable. Some of the magnitudes seem suspect; perhaps a unit conversion problem. Should cross-check against PlaSim reanalysis.
Adiv Paradise2020-08-19 20:40:39 -04:00
64efc0100b
Routines written; attached to PlaSim mainframe. Still needs to be tested and debugged.
Adiv Paradise2020-08-18 11:29:45 -04:00
5bae6ba8ed
Initial translation of Daniel Gilford's Python version of Kerry Emanuel's potential intensity code to PlaSim-suitable Fortran, framework for hurricanemod.f90
Adiv Paradise2020-08-12 22:23:39 -04:00
04583e421b
Added compile options for T85, T127, and T170
Adiv Paradise2020-07-17 18:49:13 -04:00
43c3d99aa8
Fixed problem with snow albedo on land not being set properly
Adiv Paradise2020-07-07 16:53:47 -04:00
b4386bbf78
Fix bug with necham not being MPI broadcast
Adiv Paradise2020-06-17 17:38:47 -04:00
2e98bf4b9e
Added high-cadence output mode (USE EXTREME CAUTION; POTENTIAL FOR MACHINE-BREAKING IO AND HUGE OUTPUT FILES IF MISUSED), and attempted bugfixes for a problem where prescribed albedos wouldn't be respected if a stellar temperature was set.
Adiv Paradise2020-05-05 15:28:08 -04:00
852de87949
Updated rayleigh scattering so that effective optical depth correctly reflects not only the change in weighting between the two shortwave bands, but also the wavelength-dependence of the scattering cross-section.
Adiv Paradise2020-04-30 14:22:03 -04:00
05614259eb
Added switch to (now default) compute the broadband albedo for each surface type and use that for both bands. This *should* result in 5772 K behavior matching PlaSim default behavior exactly. Also retuned spectral reflectivities to produce defaults at 5772 K.
Adiv Paradise2020-03-01 15:02:49 -06:00
7e688a6f36
Adjusted spectral calculation to exclude fluxes at wavelengths less than 316 nm; this is to mirror UV absorption by upper atmosphere, and so that spectral partitioning at 5772 K matches PlaSim's defaults
Adiv Paradise2020-02-28 17:56:45 -06:00
8009579a35
Orography smoothing will not progress if Mars mode is enabled.
Adiv Paradise2020-02-27 17:11:41 -06:00
674e5d3a9e
Bugfix to address a dependency the substellar longitude had on rotation rate for tidally-locked planets.
Adiv Paradise2020-01-22 14:17:53 -05:00
f3d44aaf6b
Fixed units error on custom spectrum import for stellar spectra.
Adiv Paradise2019-12-09 18:28:35 -05:00
5514da2339
Fixed bug where changing the rotation rate would not change the diurnal cycle period, resulting in incomplete diurnal cycles and asymmetric insolation, if a slowdown factor wasn't also included. Now specifying rotation rate is sufficient.
Adiv Paradise2019-11-22 16:08:46 -05:00
d84b1a461f
Merge branch 'hybridvertical' into starcolor
Adiv Paradise2019-11-19 17:30:51 -05:00
1aff29ec30
Adjusted how rayleigh drag works with stratosphere, and added option for stratospheric newtonian cooling
Adiv Paradise2019-11-19 17:17:55 -05:00
9afca5a352
Attempts at bugfixing; may need to tweak units/scaling on input spectra
Adiv Paradise2019-11-05 14:33:37 -05:00
747a3725aa
Added namelist options and code to use custom input spectra (untested)
Adiv Paradise2019-11-01 12:07:18 -04:00
ea4e36abda
Changed how spectra are interpolated to fit PlaSim hi-res wavelength bins
Adiv Paradise2019-10-31 23:22:23 -04:00
ecf914733b
Added python utility for PlaSim-friendly stellar spectra from PHOENIX model spectra, plus samples
Adiv Paradise2019-10-31 18:30:51 -04:00
d461842137
Changes to no-ice mode, and added an array-reading utility to utilities.f90 (prepping for the future)
Adiv Paradise2019-10-31 16:57:17 -04:00
58383f69eb
Attempting to have nice=0 turn into no sea ice....
Adiv Paradise2019-10-17 18:38:55 -04:00
c38de13002
Bugfixes. Spectral version works now, but seems to be more robust with double precision, and prone to crashes with single precision.
Adiv Paradise2019-09-25 20:07:45 -04:00
1cd90a87f4
Implemented spectral albedos, using weighted combinations of JPL Ecostress spectra.
Adiv Paradise2019-09-23 14:10:50 -04:00
9f69f6e970
Fixed bug where surface temperature in the time-averaged output file was erroneously low on the first write. This was a purely diagnostic bug--no impact on physics.
Adiv Paradise2019-08-14 17:22:17 -04:00
7992cfe14c
Added namelist switch for pressure broadening (default on, so last commit's change by default will have no effect).
Adiv Paradise2019-06-27 16:16:47 -04:00
134ef1de7a
Changed effective GHG amount calculation in radmod, so that rather than use (p/pref), it uses (sigma/sigma_ref). This removes a surface pressure dependence from GHG calculations.
Adiv Paradise2019-06-27 12:40:08 -04:00
b285ae5f8c
Added batch job creation utility for HPC environments. Currently PBS/Torque and Slurm are supported. To configure the system, buildhpc.py should be run.
Adiv Paradise2019-04-12 18:11:16 -04:00
8f9ee06550
Reverting that fractional sea ice change from ages ago; it might be causing more problems than it solved.
Adiv Paradise2019-03-26 01:45:29 -04:00
e534acdd7b
Added switch for 20-layer Rayleigh drag at top of atmosphere
Adiv Paradise2019-03-03 15:18:20 -05:00
07f15ae523
Fixes to spin-up timestepping in plasim.f90. Code runs now.
Adiv Paradise2019-02-22 21:17:11 -05:00
8ff264a9c9
First full attempt at calendar overhaul. Currently broken.
Adiv Paradise2019-02-22 17:42:30 -05:00
084a6f788b
Starting to implement new calendar system to permit non-standard day lengths
Adiv Paradise2019-02-21 20:29:54 -05:00
b08ce767f4
Added g0/g factor to Rayleigh scattering to de-couple surface gravity, added cloud critical humidity controls, and added flux tendencies and convergence to output.
Adiv Paradise2019-02-13 10:03:40 -05:00
3ecccb087d
In original bottom-scattering setup, Rayleigh scattering now scales with total atmospheric pressure.
Adiv Paradise2018-11-15 11:56:33 -05:00
94ada2f0f7
Fixed burn7 postprocessor to allow layer fluxes to be output with NLEV+1 layers (at layer interfaces). Should also work for any other interface variables!
Adiv Paradise2018-11-03 18:59:09 -04:00
48080b80c5
Added layer-by-layer radiation output to snapshots, and updated burn7 build process. Burn7 depends on netcdf C++ bindings, which are no longer included in standard distributions. This model will now ship with patched netcdf-c++ and build scripts. On initial install, run build_init.sh; run compile.sh if making subsequent changes to burn7.cpp.
Adiv Paradise2018-11-03 16:41:48 -04:00
3d7c6f2aa8
Added neqsig=3 option, which increases level resolution near surface as well as near tropopause. It's like the pseudologarithmic neqsig=2 option, but with more emphasis on the 'pseudo'. Use with 20 levels and mpstep=5.0. Also fixed some bugs.
Adiv Paradise2018-10-07 15:59:13 -05:00
28dde679e5
Adjusted default top reference pressure to 75 hPa. Also note that a reduced timestep is necessary in many cases with logarithmic spacing and high surface pressures.
Adiv Paradise2018-10-02 16:29:56 -04:00
6749c41ca0
Bugfixes in glacier module--snow depth changes getting stored as rates rather than amounts would have given too-small changes in snow depth
Adiv Paradise2018-10-02 15:55:34 -04:00
131229fe1e
Added option to space vertical levels pseudo-logarithmically, anchored to a user-specified upper reference pressure (default 150 hPa). The upper two layers are at slightly higher pressures than if the spacing were truly logarithmic, which is necessary to ensure half-levels and full-levels can be appropriately spaced. Logarithmic spacing is turned on with neqsig=2 in plasim_namelist, and setting ptop in plasim_namelist (in Pa).
Adiv Paradise2018-10-02 15:33:28 -04:00
d4d343e5e5
Added bypass so nglacier doesn't have to be on for oroini to get called.
Adiv Paradise2018-08-19 17:02:06 -04:00
57190427a8
Added OROSCALE to landmod_namelist--apparently it had been removed from surfmod_namelist at some point?
Adiv Paradise2018-08-13 17:45:57 -04:00
3b5589255f
Added desert planet mode. Also made the soil heat capacity visible to landmod_namelist.
Adiv Paradise2018-08-09 12:51:32 -04:00
5bd4716a0d
Add low-IO flag and nstpw parameter, which override the model's output cadence. With nlowio=1 and nstpw>0, the code will only write output every nstpw timesteps, and will compute time-averages for all variables to avoid loss of fidelity.
Adiv Paradise2018-07-14 15:34:57 -04:00
Adiv Paradise