jgr2014: Greenhouse effect and climate.

In brasmus/GHE: GHE

Greenhouse effect and climate.

Description

Functions for plotting scematics and carrying out some simple climate analysis plus a number of other functions.

p.hydrostatic the pressure at a given height under the assumption of hydrostatic balance. forcings Time series plot of solar forcing. gcr retrieves galactic cosmic ray data and sunspots the sunspot data. jgr2014 R-code for reproducing figures in the publication JGR2014. HITRAN function plotting IR spectral lines. solarsystem Emission temperature of the objectys in the solar system. Keelercurve time series of CO2. blackbody code to plot the Planck function - blackbody radiation co2 co2.temp plots scatteer plots between ln(CO2) and temperature.

Usage

p.hydrostatic(z,p0 = 1013.250,Temp = 288,g = 9.81,
              k = 1.38e-23,M = 0.027/6.022e+23)
sunspots(url="http://sidc.oma.be/DATA/monthssn.dat")
forcings()
gcr()
jgr2014()
HITRAN(param=2, plot = TRUE,xlim=c(1,20),ylim=c(0.1,100))
solarsystem()
Keelercurve(url="ftp://ftp.cmdl.noaa.gov/ccg/co2/trends/co2_mm_mlo.txt",
            urlbg="http://www.realclimate.org/images/raffinery.jpg")
blackbody(h=6.64e-34,k=1.38e-23,c=3.00e8,HITRAN = TRUE)
co2(url="ftp://ftp.cmdl.noaa.gov/ccg/co2/trends/co2_mm_mlo.txt",
    plot=TRUE,Temp=TRUE)
co2.temp()
vostoc()
T.e(R,A,S0=1367,rho=5.67e-8)
stand(x,ii=NULL)
ma.filt(x, n)
colorbar(yymm,fig=c(0.33,0.66,0.30,0.34),breaks,col)
prettygraph(t,y,col="black",tshadeoffset=0.1,yshadeoffset=0.1)
largesymbols(x,y,col="blue",cex=3)
IR(x,levels)
rot(x,angle=80)
sol(x,y,r=10)
piler(lw,x,y)
testrot()
gauss.filt(x, n)
elninoyrs(ylim=c(-10,10),yrs=NULL,col="grey90",lwd=7) 
prepareData(tempfil= "~/ferret_scripts/ncep-ta.nc") 
AtmOverturning() 
emissionHeigth(ncep=TRUE) 
Hulbert1931Tz(x,xb=3300,S=1.35e6,rho=5.67e-8,alpha=0.10) 
supportDiag() 
olr.corr()
GHE()
GHEschematic()
GeoChaos(X.0=c(-0.1,0.1,0),N=10000,R=0.1,nu=1,sigma=5)
LorenzSystem(X.0=c(-0.1,0.1,0),N=100000,sigma=10,beta=8/3, rho=28)
greenhouseeffect()
caldat(julian)
julday(mm, id, iyyy)

Arguments

z	height in m.
p0	mean pressure at mean sea-level
Temp	surface temperature
g	acceleration due to gravity in m/s^2
k	Boltzmann constant

.

M	mass ratio for air
url	URL of data source.
urlbg	URL of background graphics.
param	parameter
h	Planck constant in Js
c	speed of light in m/s
R	distance from the sun
A	planetary albedo
S0	'solar constant'
x	time series/horizontal coordinate or parameter for equation in Hulbert (1931)

.

y	time series/vertical coordinate

.

angle

angle in degrees

.

n	filter window width

.

levels

Vertical levels

.

t

time axis

.

lw

long-wave

.

HITRAN	TRUE: use the prvided HITRAN data
plot	TRUE: plot the results
fig	see par
yymm	year and month
breaks	see hist
tshadeoffset	graphics adjustmant
yshadeoffset	graphics adjustmant
cex	see par
r	radius
tempfil	File name for temperature data after processing with Ferret
ylim	see plot
lwd	see plot
yrs	years
col	see plot
ncep	FLAG - TRUE plot NCEP
xb	parameter for equation in Hulbert (1931)
S	The 'solar constant'
rho	Stefan-Boltzman's constant or parameters for the for Lorenz/GeoChaos system.
alpha	The albedo
X.0	initial conditions for Lorenz/GeoChaos system.
N	length of integration
sigma	parameters for the for Lorenz/GeoChaos system.
beta	parameters for the for Lorenz/GeoChaos system.
nu	parameters for the for Lorenz/GeoChaos system.
julian	julian day
mm	month
id	day
iyyy	year

Author(s)

R.E. Benestad

Examples

print("emissionHeigth()")
emissionHeigth()

print("AtmOverturning()")
AtmOverturning()

print("olr.corr()")
olr.corr()

# other plots...

print("jgr2014()")
jgr2014()

## Not run: 
print("co2()")
co2()

## End(Not run)
print("LorenzSystem()")
LorenzSystem()

print("GeoChaos()")
GeoChaos()

print("GHE()")
GHE()

print("greenhouseeffect()")
greenhouseeffect()

print("blackbody()")
blackbody()

print("solarsystem()")
solarsystem()

print("supportDiag()")
supportDiag() 

brasmus/GHE documentation built on Jan. 3, 2024, 6:03 p.m.