利用Python读取netCDF文件

#!/usr/bin/python

from netCDF4 import Dataset

Prf_file_path  = '/home3/jmzhang/ex_nc_read/FY3C_GNOSX_GBAL_L2_ATP_MLT_NUL_20170518_0448_G21XX_MS.NC'
out_raw_data   = '/home3/jmzhang/ex_nc_read/outfile.txt'

def GetPrfVar(nc_file_name): 

    nc_file       = Dataset(nc_file_name, 'r')
    lat           = nc_file.lat
    lon           = nc_file.lon
    MSL_alt       = nc_file.variables['MSL_alt'][:]
    Temp          = nc_file.variables['Temp'][:]
    Pres          = nc_file.variables['Pres'][:]

    return (lat, lon, MSL_alt, Temp, Pres)

(lat, lon, MSL_alt, Temp, Pres) = GetPrfVar(Prf_file_path)

file_put = open(out_raw_data, 'a')
for k in range(len(MSL_alt)):
                
    file_put.write('{:+.10f}\t\t{:+.10f}\t\t{:+.10f}\n'.format(MSL_alt[k], Temp[k], Pres[k]))

file_put.close()