# Temperature
if ( level == 850 && indicatorOfParameter == 11 ) {
print "found indicatorOfParameter=[indicatorOfParameter] level=[level] date=[date]";
transient oldtype = type ;
set identificationOfOriginatingGeneratingSubCentre=98;
set gribTablesVersionNo = 128;
set indicatorOfParameter = 130;
set localDefinitionNumber=1;
set marsClass="od";
set marsStream="kwbc";
# Negatively/Positively Perturbed Forecast
if ( oldtype == 2 || oldtype == 3 ) {
set marsType="pf";
set experimentVersionNumber="4001";
}
# Control Forecast
if ( oldtype == 1 ) {
set marsType="cf";
set experimentVersionNumber="0001";
}
set numberOfForecastsInEnsemble=11;
write;
print "indicatorOfParameter=[indicatorOfParameter] level=[level] date=[date]";
print;
}