Opened 3 years ago
Closed 3 years ago
#2476 closed help (fixed)
convert hybrid pressure to pressure coordinate
Reported by: | ajd | Owned by: | willie |
---|---|---|---|
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 3 years ago by willie
- Description modified (diff)
- Keywords byte swapping added
- Owner changed from um_support to willie
- Status changed from new to accepted
comment:2 Changed 3 years ago by ajd
Hi Willie,
Many thanks, that worked!
Kind regards,
Andrea
comment:3 Changed 3 years ago by jeff
ppinfo -ieeebs pph
should also work on the original file
comment:4 Changed 3 years ago by willie
- Resolution set to fixed
- Status changed from accepted to closed
Note: See
TracTickets for help on using
tickets.
Hi Andrea,
It's byte swapped. If you copy it to ARCHER and then
The pph.out is in the right format.
Regards
Willie