Displaying 7 results from an estimated 7 matches for "tstep".
Did you mean:
step
2011 Sep 07
1
Error: in routine alloca() there is a stack overflow: thread 0, max 535822282KB, used 0KB, request 24B
...get.var.ncdf(nc.ori,varid="MEOH")
### Here is the output associated with the print ("fininshed open file") and
print(nc) portion of the script:
[1] "fininshed open file"
[1] "file /users/me/latlong_spec2_WOCONV_150s_2004187.ncf has 6 dimensions:"
[1] "TSTEP Size: 2"
[1] "DATE-TIME Size: 2"
[1] "LAY Size: 1"
[1] "VAR Size: 30"
[1] "ROW Size: 1200"
[1] "COL Size: 2496"
[1] "------------------------"
[1] "file
/users/huxxx249/MEGANv2.04/Output/latlong_spec2_WOCONV_150s_2004...
2006 Sep 22
4
Creating Movies with R
...ma_pos)*exp(-1/2*(x-expect_position)^2/sigma_pos)
}
#### Now the physical parameters
tau<-0.1
beta<-1/tau
St<-tau ### since I am in dimensionless units and tau is already in
units of 1/|alpha|
D=2e-2
q<-2*beta^2*D
############### Now the grid in space and time
time<-5 # time extent
tsteps<-501 # time steps
newtime<-seq(0,time,len=tsteps)
#### Now the things specific for the dynamics along x
lam1<- -beta/2*(1+sqrt(1+4*St))
lam2<- -beta/2*(1-sqrt(1+4*St))
xmin<- -0.5
xmax<-0.5
x0<-0.1
vx0<-x0
nx<-101 ## grid intervals along x
newx<-seq(xmin,xmax,len=nx)...
2010 Apr 26
1
finite difference scheme for 2D differential equations
...on using a code like the
following. But I want to write a similar code for,say, a two dimensional
diffusion equation. Any kind of help/advice is highly appreciated.
Here is how I did for the 1D equation (in this case for 1D diffusion)
Equation: dc/dt = D*d^2c/dx^2
dx = 10
dt = 0.1
D = 15
n = 20
tstep = 200
C = 50
dt = 0.1
dx = 10
conc<-matrix(0,tstep,n)
conc[1,1:10] = C
for (i in 2:tstep) {
for (j in 1:n){
if (j==1){
conc[i,j] = conc[i-1,j] - dt*(D/dx^2)*(conc[i-1,j] - conc[i-1,j+1])
}
conc[i,j] = conc[i,j]
if(j >1 && j < n){
conc[i,j] = conc[i-1,j] + dt*(D/dx^2)* (conc[i-1,(...
2011 Oct 31
2
3 Overlayed simple plots
Dear all,
I am plotting 3 plots into the same x and y axis.
I want the first one to be painted red with a continuous line
The second one green with a continuous line
and the third one blue with a continuous line
plot(max_power(data),ylim=c(-120,-20))
par(new=T)
plot(min_power(data),ylim=c(-120,-20))
par(new=T)
plot(mean_power(data),ylim=c(-120,-20))
par(new=F)
Is it also a
2018 May 29
1
Difficulty in writing R code for one pool dynamic model
...ide the value of pool could be (23.8, 24.8, 25.0 or 24.7). Thevalue I got does not correspond to any of the choices given above. ?I don't know where I got it wrong.?
Below is the Ro code I tried to estimate thepool size of A at time 30.
#Step 1: Specify time?
Tstart <- 1
Tstop <- 30
Tstep <- 1
#Step 2: Set variables corresponding to pools
inA<-5
A <- inA
t <- Tstart
#Step 3: Set up rate parameters
A<-5
kinA <- 5
koutA <- 0.2
#Step 4: Write arrays that track pool size overtime
val <- array(1, dim=c(1, length(seq(Tstart,Tstop, by=Tstep))))
val
#St...
2008 Jul 09
0
"Rotated Lat-Lon" projection in mapproj
...sep="
"),intern=TRUE))
## Projected coordinates
coords_lamb<-mapproject(lon,lat,projection="lambert",parameters=c(60.,30.))
xs<-t(matrix(coords_lamb$x,ncol=NCC,nrow=NRC))[,1]
ys<-t(matrix(coords_lamb$y,ncol=NCC,nrow=NRC))[1,]}
## Plot map
filled.contour(xs,ys,Tsd[,,tstep],levels=fillevs,col=cpal,plot.axes={
map('world',projection="lambert",parameter=c(60.,30.),add=TRUE,col="grey");
contour(xs,ys,Tm[,,tstep]-273.15,levels=c(-20:20)*3,add=TRUE);
axis(1,labels=F,tick=F);axis(2,labels=F,tick=F)},main=title) #,asp=1.)
Appreciate your help...
2016 Oct 12
4
[test-suite] making polybench/symm succeed with "-Ofast" and "-ffp-contract=on"
On Wed, Oct 12, 2016 at 10:53 AM, Hal Finkel <hfinkel at anl.gov> wrote:
> I don't think that Clang/LLVM uses it by default on x86_64. If you're using -Ofast, however, that would explain it. I recommend looking at -O3 vs -O0 and make sure those are the same. -Ofast enables -ffast-math, which can legitimately cause differences.
>
The following tests pass at "-O3" and