source: branches/2015/dev_r5102_CNRS11_TestCase/NEMOGCM/TOOLS/TEST_CASES/CREATE_FILES/make_bathy.py @ 5346

#!/opt/local/bin/python

import os,sys
from netCDF4 import Dataset as netcdf
import numpy as np

resname='r025'

# input file
fcoord='coordinates_'+str(resname)+'.nc'

# output file
fbathy='bathy_'+str(resname)+'.nc'

# Reading coordinates file
nccoord=netcdf(fcoord,'r')
nav_lon=nccoord.variables['glamt']
nav_lat=nccoord.variables['gphit']
LON1= nav_lon.shape[1]
LAT1= nav_lon.shape[0]

print '   creating bathymetry file  '+fbathy+'...'

# Creating bathy netcdf file
nc=netcdf(fbathy,'w')
nc.createDimension('y',LAT1)
nc.createDimension('x',LON1)

cdflon=nc.createVariable('nav_lon','f',('y','x'))
cdflat=nc.createVariable('nav_lat','f',('y','x'))
cdfbathy=nc.createVariable('Bathymetry','f',('y','x'))

cdflon[:,:]=nav_lon[:,:]
cdflat[:,:]=nav_lat[:,:]
cdfbathy[:,:]=5000

nc.close()
nccoord.close()