mpmd-Rscripts/write_mpmd_config_file.R at master · oliversus/mpmd-Rscripts · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
# R script to write config files for MPMD job on ECMWF supercomputer

# each MPMD job will be executed on an individual CPU on CRAY,
# thus each job needs an individual config file

# read command arguments
args = commandArgs( trailingOnly = T )

# associate local variables with appropriate entries in argument vector
year         = args[2]
month        = args[3]
prodtype     = args[4]
sensor       = args[5]
platform     = args[6]
jobID        = args[7]
ndays        = args[8]
cfg_dir      = args[9]
cfg_prefix   = args[10]
cfg_suffix   = args[11]
cfg_base     = args[12]
base_path    = args[13]
flist_l2_out = args[14]

# convert single digit months to double digits
month_mm = ifelse( nchar( month ) == 1, paste( "0", month, sep="" ), month )

# list L2 output folder matching filter criteria
L2_folder = list.files( base_path, pattern = paste( year, month_mm, "01_", sensor,
	            "_", platform, "_retrieval_*", sep="" ), ignore.case=T )

# if more than one matching output folder available, sort these as a function
# of creation time and select most recently created folder
if ( length( L2_folder ) > 1 ) {
   creation_time = file.info( paste( base_path, L2_folder, sep="/" ) )
   L2_folder = L2_folder[ order( creation_time$ctime, decreasing = T )[1] ]
}

# extract ID from L2 output folder
L2_id = paste( "ID", unlist( strsplit( L2_folder, "_ID" ) )[2], sep="" )

# create MPMD config files for each day
for (i in 1:ndays) {

    # convert single digit days to double digits
    i_dd = as.character(i)
    i_dd = ifelse(nchar(i_dd) == 1, paste("0", i_dd, sep=""), i_dd)

    # build config file name (consistent with write MPMD taskfile!)
    config_file = paste(cfg_dir, "/", cfg_prefix, cfg_base, i_dd, cfg_suffix, sep="")

    # open config file connection
    fileConn = file( config_file )

    # write data to config file
    writeLines( c(

	paste("YEAR=", year, sep=""),
	paste("MONTH=", month, sep=""),
	paste("DAY=", i, sep=""),
	paste("prodtype=", prodtype, sep=""),
	paste("sensor=", sensor, sep=""),
	paste("platform=", platform, sep=""),
	paste("filelist_level2_output=", flist_l2_out, sep=""),
	"local=F",
	"slon=0",
	"elon=18",
	"slat=42",
	"elat=53",
	"gridxloc=10",
	"gridyloc=10",
	"gridxl3=2",
	"gridyl3=2",
	"gridxl2b=10",
	"gridyl2b=10",
	paste("id=", L2_id, sep="")

	), fileConn)

    # close config file connection
    close(fileConn)

}