Opened 7 months ago

Closed 6 months ago

#2476 closed help (fixed)

convert hybrid pressure to pressure coordinate

Reported by: ajd Owned by: willie
Priority: normal Component: UM Tools
Keywords: byte swapping Cc:
Platform: Other UM Version: 4.5

Description (last modified by willie)

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


Change History (4)

comment:1 Changed 7 months ago by willie

  • Description modified (diff)
  • Keywords byte swapping added
  • Owner changed from um_support to willie
  • Status changed from new to accepted

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

comment:2 Changed 7 months ago by ajd

Hi Willie,

Many thanks, that worked!

Kind regards,
Andrea

comment:3 Changed 7 months ago by jeff

ppinfo -ieeebs pph

should also work on the original file

comment:4 Changed 6 months ago by willie

  • Resolution set to fixed
  • Status changed from accepted to closed
Note: See TracTickets for help on using tickets.