convert hybrid pressure to pressure coordinate

[ajd]

Hi,

I'm trying to convert model output from FAMOUS (which I think is at UM4.5) from hybrid pressure levels to pressure levels.

NCL has a function to do this provided with the A and B coefficients.

Following this thread, http://cms.ncas.ac.uk/ticket/338, it appears that these can be found in the header. However, when I attempt to print the header info using ppinfo, I get the following error:

 Field number            1
 Field code =****** (    16777216 )
 No processing on field (           0 )
 Date range of field ********:**** to ********:0000
 Day number range of field *********** to ***********
 Time indicator =    201326592
 Forecast period (hours) =    103715586
 Level code =    134217728
 Value of vertical coordinate =     0.0000
 Cross section grid (    16777216 )
 Global field (           0 )
 Number of rows =    620756992
 Number of grid points in each row =    805306368
 Total length of data =   -268042240
 Unpacked Data (           0 )
 Missing data indicator =   4.6206416E-41
 MKS scaling factor =   4.6006030E-41
 -1995243520   201326592    16777216   100663296           0  1258356736
    16777216   201326592    16777216           0           0  1258356736
   201326592   103715586  -268042240    16777216           0   620756992
   805306368           0           0    33554432    16777216           0
           0   134217728           0   436673554           0  -268042240
   167772160    33554432  -201261056           0           0           0
           0  1459879936    16777216  1091043840           0  1245642752
           0           0    16777216
  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00
  7.1300868E-41  8.9776989E-41  0.0000000E+00  0.0000000E+00  0.0000000E+00
  6.4664319E-41  0.0000000E+00  0.0000000E+00  6.8251643E-41  5.7666234E-41
  8.6364827E-41  8.6185461E-41  4.6206416E-41  4.6006030E-41
forrtl: severe (39): error during read, unit 11, file /net/plane/export/plane/data-19/robin/robinia/um12/um_archive/SHa20/pph
Image              PC                Routine            Line        Source             
ppinfo             0000000000454CDA  Unknown               Unknown  Unknown

Do you know what could be going wrong? The data is on local jasmin in Reading.

Many thanks,
Andrea

[willie]

Hi Andrea,

It's byte swapped. If you copy it to ARCHER and then

$UMDIR/bin/bigend -32 pph pph.out

The pph.out is in the right format.

Regards
Willie

[ajd]

Hi Willie,

Many thanks, that worked!

Kind regards,
Andrea

[jeff]

ppinfo -ieeebs pph

should also work on the original file