Planetary Code Collection

by Norbert Schorghofer


This program collection contains
* semi-implicit one-dimensional thermal model for planetary surfaces
* an explicit subsurface vapor diffusion and deposition model
* models for the equilibrium ice table on Mars
* fast method for subsurface-atmosphere vapor exchange for Mars
* fast (asynchronous) model for ice retreat on asteroids
* Monte-Carlo model for ballistic hops in the exosphere of the Moon and Ceres
* 3D model of shadowing, illumination, and scattering



==Mars Subsurface Ice Model (M-SIM)==

*Mars Thermal Model*

Mars/mars_thermal1d.f: (main program)
Mars/flux.f
Mars/marsorbit.f
Common/conductionQ.f
Common/conductionT.f
Common/tridag.for
Common/grids.f
Common/psvco2.f
Common/julday.for
Mars/input.par
Documentation: User Guide Part 1


*Vapor Diffusion Model*

Mars/exper_thermal1d.f: (main program)
Mars/vapordiffusioni.f
Mars/adsorption.f
Mars/exper.par
Documentation: User Guide Part 2
Documentation: Schorghofer, N. & Aharonson, O. (2005) J. Geophys. Res. 110, E05003, Appendices


*Equilibrium Ice Table*

Mars/mars_mapi.f: (main program)
Mars/mars_mapt2.f: (main program)
Mars/mars_mapiqh2v.f90: (main program)
Mars/jsub.f
Mars/jsubv.f90
Mars/marsorbit.f
Common/conductionQ.f
Common/conductionT.f
Common/tridag.for
Common/grids.f
Common/julday.for
Common/psv.f
Common/psvco2.f
Mars/mapgrid.dat
Documentation: Users Guide Section 3.1
Documentation: Schorghofer, N. & Aharonson, O. (2005) J. Geophys. Res. 110, E05003


*Fast Method for Subsurface Ice Dynamics*

Mars/stabgrow_fast.f90: (main program)
Mars/exper_fast.f90: (main program)
Mars/mars_fast.f90: (main program)
Mars/fast_modules.f90
Mars/fast_subs_univ.f90
Mars/fast_subs_exper.f90
Mars/fast_subs_mars.f90
Common/conductionQ.f
Common/conductionT.f
Common/tridag.for
Common/grids.f
Common/generalorbit.f
Common/psv.f
Common/psvco2.f
Common/derivs.f90
Mars/input_fast.par
Mars/lats.ph
Documentation: Schorghofer, N. (2010) Icarus 208, 598-607


==Other Models for Planetary Surfaces==

*Basic Thermal Model for Asteroids*

Asteroid/asteroid_thermal.f90: (main program)
Asteroid/oneasteroid_thermal1d.f90
Asteroid/insolonly.f90 
Common/flux_noatm.f90
Common/conductionQ.f
Common/tridag.for
Common/grids.f
Common/generalorbit.f
Documentation: User Guide Part 1


*Monte-Carlo Model for Surface-bounded Exosphere*

Exosphere/moon4.f90: (main program)
Exosphere/ceres_exo.f90: (main program)
Exosphere/montecarlo.f90
Exosphere/geogrid.f90
Exosphere/geogrid_D.f90
Common/subl_subs.f90
Common/gasdev.for
Common/ran2.for
Documentation: None
Documentation: User Guide Part 6


*Shadowing and Illumination*

Topo3D/shadows.f90: (main program)
Topo3D/fieldofviews.f90: (main program)
Topo3D/cratersQ_snapshot.f90: (main program)
Topo3D/cratersQ_moon.f90: (main program)
Topo3D/topos.f90
Topo3D/crater_modules.f90
Topo3D/crater_common.f90
Topo3D/shadow_subs.f90
Topo3D/fieldofview_subs.f90
Topo3D/model_subs.f90
Common/hpsort.for
Topo3D/topo40.xyz
Documentation: User Guide Part 5


*Asynchronous Model for Temperature, Impact Stirring, and Ice Loss on Asteroids*

Asteroid/asteroid_fast2.f90: (main program)
Asteroid/fast_modules_asteroid2.f90
Asteroid/fast_subs_asteroid2.f90 
Asteroid/impactstirring.f90
Common/flux_noatm.f90
Common/conductionQ.f
Common/tridag.for
Common/grids.f
Common/generalorbit.f
Common/psv.f
Common/derivs.f90
Common/ran2.f
Asteroid/lats.0
Documentation: Schorghofer, N. (2016) Icarus 276, 88-95.


*Irradiance Model for Terrestrial Analog*

EarthAnalogs/insol_earth.f90: (main program)
EarthAnalogs/insol_flat.f90: (simple main program)
EarthAnalogs/mk_atmosphere.f90
EarthAnalogs/sunpos.f90
Topo3D/topos.f90
Topo3D/crater_modules.f90 
Topo3D/crater_common.f90
Topo3D/model_subs.f90
Documentation: User Guide Sections 4.1 & 5.2-5.4



NOTE: Third party source code from Numerical Recipes is included in this distribution, but is covered by a separate copyright. These are files ending with .for.  Most of this code will not work without them. 



ACKNOWLEDGMENTS:

Mar 2016: Thanks to Lior Rubanenko for a bug report in cratersQ_*

2006: Troy Hudson discovered a grid-point offset in conductionT and conductionQ, which has been corrected.

2005: Thanks to Mischa Kreslavsky for providing correct formulas for energy balance on a slope.

* Many Thanks to Andy Vaught for developing an open-source Fortran 95 compiler (www.g95.org).  The early versions of this code were developed with this compiler.

2001: Samar Khatiwala invented an elegant implementation of the upper radiation boundary condition for the Crank-Nicolson method.

SUPPORT: This code development was supported by NASA, Caltech, and the University of Hawaii. Undoubtedly, some parts were written without support.