Project

General

Profile

« Previous | Next » 

Revision cf724805

Added by Adam Wilson about 12 years ago

Resolved problem with temporary directory for Rscript call.

View differences:

climate/procedures/MOD06_L2_process.r
90 90
#RESAMPLING_TYPE =",ifelse(grepl("Flag|Mask|Quality",vars),"NN","CUBIC"),"
91 91
RESAMPLING_TYPE =NN
92 92
OUTPUT_PROJECTION_TYPE = SIN
93
 OUTPUT_PROJECTION_PARAMETERS = ( 6371007.181 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 )
93
OUTPUT_PROJECTION_PARAMETERS = ( 6371007.181 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 )
94 94
# projection parameters from http://landweb.nascom.nasa.gov/cgi-bin/QA_WWW/newPage.cgi?fileName=sn_gctp
95 95
ELLIPSOID_CODE = WGS84
96 96
OUTPUT_TYPE = HDFEOS
......
98 98
END
99 99

  
100 100
",sep="")
101
 
101

  
102 102
  ## if any remnants from previous runs remain, delete them
103 103
  if(length(list.files(tempdir(),pattern=basename(file)))>0)
104 104
    file.remove(list.files(tempdir(),pattern=basename(file),full=T))
......
269 269
##mod06(date,tile)
270 270

  
271 271
## run it for all dates
272
mclapply(notdone[1:100],mod06,tile)
272
mclapply(notdone,mod06,tile)
273 273

  
274 274

  
275 275
## quit R
climate/procedures/Pleiades.R
11 11

  
12 12
cat(paste("
13 13
#PBS -S /bin/bash
14
#PBS -l select=1:ncpus=4:mpiprocs=4:model=wes
14
#PBS -l select=64:ncpus=4:mpiprocs=4:model=wes
15 15
####old PBS -l select=64:ncpus=4:mpiprocs=4:model=wes
16 16
####### old: select=48:ncpus=8:mpiprocs=8:model=neh
17
#PBS -l walltime=0:10:00
17
#PBS -l walltime=10:00:00
18 18
#PBS -j oe
19 19
#PBS -m e
20 20
#PBS -V
21 21
####PBS -W group_list=s1007
22
#PBS -q devel
22
###PBS -q devel
23 23
###PBS -o log/log_^array_index^
24 24
#PBS -o log/log_DataCompile
25 25
#PBS -M adam.wilson@yale.edu
......
54 54
#system("cat ~/environmental-layers/climate/procedures/MOD06_L2_process.r")
55 55

  
56 56
## Submit it (and keep the pid)!
57
pid=system("qsub -q devel MOD06_process",intern=T); pid; pid=strsplit(pid,split="[.]")[[1]][1]
57
pid=system("qsub MOD06_process",intern=T); pid; pid=strsplit(pid,split="[.]")[[1]][1]
58 58

  
59 59
#system("qsub MOD06_process")
60 60

  
......
71 71

  
72 72
### copy the files back to Yale
73 73
system("scp 2_daily/* adamw@acrobates.eeb.yale.edu:/data/personal/adamw/projects/interp/")
74

  
75
system("scp  /tmp/Rtmp6I6tFn/MOD06_L2.A2000061.1615.051.2010273184629.hdf adamw@acrobates.eeb.yale.edu:/data/personal/adamw/projects/interp/")
76

  
77
list.files(" /tmp/Rtmp6I6tFn")

Also available in: Unified diff