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
|
# Negatively/Positively Perturbed Forecast
transient pNumber=identificationNumber;
set identificationOfOriginatingGeneratingSubCentre = 98 ;
set localDefinitionNumber = 1 ;
set marsClass="od";
set marsStream="kwbc";
#print "marsType=[marsType:l]";
if ( marsType == 2 || marsType == 3 ) {
set marsType=11;
}
# Control Forecast
if ( marsType == 1 ) {
set marsType="cf";
}
set experimentVersionNumber="0001";
set numberOfForecastsInEnsemble=11;
set perturbationNumber=pNumber;
# Geopotential height
if ( ( level == 500 || level == 1000 ) && indicatorOfParameter == 7 ) {
print "found indicatorOfParameter=[indicatorOfParameter] level=[level] date=[date]";
set gribTablesVersionNo = 128;
set indicatorOfParameter = 129 ;
print "indicatorOfParameter=[indicatorOfParameter] level=[level] date=[date]";
write "out.grib";
print;
}
# U wind
if ( level == 250 && indicatorOfParameter == 33 ) {
print "found indicatorOfParameter=[indicatorOfParameter] level=[level] date=[date]";
set gribTablesVersionNo = 128;
set indicatorOfParameter = 131 ;
print "indicatorOfParameter=[indicatorOfParameter] level=[level] date=[date]";
write "out.grib";
print;
}
# V wind
if ( level == 250 && indicatorOfParameter == 34 ) {
print "found indicatorOfParameter=[indicatorOfParameter] level=[level] date=[date]";
set gribTablesVersionNo = 128;
set indicatorOfParameter = 132;
print "indicatorOfParameter=[indicatorOfParameter] level=[level] date=[date]";
write "out.grib";
print;
}
# Relative humidity
if ( level == 700 && indicatorOfParameter == 52 ) {
print "found indicatorOfParameter=[indicatorOfParameter] level=[level] date=[date]";
set gribTablesVersionNo = 128;
set indicatorOfParameter = 157 ;
print "indicatorOfParameter=[indicatorOfParameter] level=[level] date=[date]";
write "out.grib";
print;
}
# Temperature
if ( level == 850 && indicatorOfParameter == 11 ) {
print "found indicatorOfParameter=[indicatorOfParameter] level=[level] date=[date]";
set gribTablesVersionNo = 128;
set indicatorOfParameter = 130;
print "indicatorOfParameter=[indicatorOfParameter] level=[level] date=[date]";
write "out.grib";
print;
}
|
