convert_float2double_mitgcm.py

#!/usr/bin/env python

import numpy as np
import struct

def write_new_data(data,fileout):
    fid   = open(fileout, "wb")
    flatdata = data.flatten()
    for kk in np.arange(len(flatdata)):
        tmp = struct.pack('>d',flatdata[kk])
        fid.write(tmp)
    fid.close()
    return None

def convert_to_real8(filein, fileout):
    # read the original data in real4
    data = np.fromfile(filein,dtype='>f')
    # write file in real8
    write_new_data(data,fileout)
    return None

# usage
# convert_to_real8(indir+filein,outdir+fileout)