#####################################################################
# SCRIPT START ######################################################
# Basic Settings:
filenamestr = '_frc.nc'
filetypestr = 'ROMS Forcing file'
# READING PREVIOUSLY BUILT RELEVANT FILES: ###########################
# metadata ascii file
# OA-created netcdf initial T, S file
# grid netcdf file
print ' \n' + '==> ' + ' READING ASCII METADATA FILE ...\n' + ' '
run = run_setup('../phd9_run.setup')
print ' \n' + '==> ' + ' READING ECMWF NETCDF FILE ...\n' + ' '
datafile = netCDF4.Dataset(run.datadir + run.frc_filename)
# assigning some variables from data file
lon = datafile.variables['longitude'][:] - 360
lat = datafile.variables['latitude'][:]
u10 = datafile.variables['u10']
v10 = datafile.variables['v10']
lon, lat = np.meshgrid(lon, lat)
u = datafile.variables['u10'][:] * u10.scale_factor + u10.add_offset
v = datafile.variables['v10'][:] * v10.scale_factor + v10.add_offset
u = u.squeeze(); v = v.squeeze()
#####################################################################
# SCRIPT START ######################################################
# Basic Settings:
filenamestr = '_bry.nc'
filetypestr = 'ROMS Boundary Conditions file'
# READING PREVIOUSLY BUILT RELEVANT FILES: ###########################
# metadata ascii file
# OA-created netcdf initial T, S file
# grid netcdf file
print ' \n' + '==> ' + ' READING ASCII METADATA FILE ...\n' + ' '
run = run_setup('../phd16_run.setup')
print ' \n' + '==> ' + ' READING FEATURE MODEL FIELD ...\n' + ' '
datafile = sp.loadmat(run.datadir + run.ini_filename)
# assigning some variables from data file
Zlev = datafile['z'][:].ravel(); Zlev = np.abs(Zlev); Zlev = -Zlev
N1 = Zlev.size
lon = datafile['lon'][:]
lat = datafile['lat'][:]
temp = datafile['temp'][:]
salt = datafile['salt'][:]
u = datafile['u'][:]
v = datafile['v'][:]
ubar = datafile['ubar'][:]
vbar = datafile['vbar'][:]
dx = x - x0
dx = abs(dx)
fn = np.where( dx == dx.min() )
return fn[0][0]
# PYTHON SCRIPT START ######################################################
# Input parameters
expt = 'phd18'
zlev = 2000
print ' \n' + '==> ' + ' READING ASCII METADATA FILE ...\n' + ' '
roms = run_setup(expt + '_run.setup')
print ' \n' + '==> ' + ' READING GRID NETCDF FILE ...\n' + ' '
grd = netCDF4.Dataset(roms.rundir + roms.run_name + '_grd.nc')
# assigning some variables from grid file
lon = grd.variables['lon_rho'][:]
lat = grd.variables['lat_rho'][:]
latu = grd.variables['lat_u'][:]
latv = grd.variables['lat_v'][:]
lonu = grd.variables['lon_u'][:]
lonv = grd.variables['lon_v'][:]
h = grd.variables['h'][:]
print ' \n' + '==> ' + ' READING CHOSEN ROMS OUTPUT NETCDF FILE ...\n' + ' '
avg = netCDF4.Dataset(roms.rundir + roms.run_name + '_avg.nc')
dz = np.diff(z, axis=0)
aux = dz[0, :]
aux.shape = (1, np.size(aux))
dz = np.concatenate((dz, aux), axis=0)
transp = np.abs(dx) * np.abs(dz) * v
transp = transp.sum()
transp = transp / 1e6
return transp
# PYTHON SCRIPT START ######################################################
print ' \n' + '==> ' + ' READING ASCII METADATA FILE ...\n' + ' '
roms = run_setup(expt + '_run.setup')
print ' \n' + '==> ' + ' READING GRID NETCDF FILE ...\n' + ' '
grdfile = netCDF4.Dataset(roms.rundir + roms.run_name + '_grd.nc')
# assigning some variables from grid file
lon = grdfile.variables['lon_rho'][:]
lat = grdfile.variables['lat_rho'][:]
latu = grdfile.variables['lat_u'][:]
latv = grdfile.variables['lat_v'][:]
lonu = grdfile.variables['lon_u'][:]
lonv = grdfile.variables['lon_v'][:]
h = grdfile.variables['h'][:]
print ' \n' + '==> ' + ' READING CHOSEN ROMS OUTPUT NETCDF FILE ...\n' + ' '
outfile = netCDF4.Dataset(roms.rundir + roms.run_name + '_' + filetype + '.nc')
vst = 0.03 # contour interval
niN = -12 # latitude for NBUC section
niC = -17 # latitude for RC Eddy
niS = -19 # longitude for Abrolhos Eddy
XlimN = (-38, -35) # x axis limits for northern section
XlimC = (-39, -35) # x axis limits for southern section
XlimS = (-39, -35) # x axis limits for eastern section
Zlim = (-1500, 0) # z axis limits
print ' \n' + '==> ' + ' READING ASCII METADATA FILE ...\n' + ' '
pathname = "/home/rsoutelino/myroms/phd_run/"
roms = run_setup(pathname + expt + '_run.setup')
print ' \n' + '==> ' + ' READING GRID NETCDF FILE ...\n' + ' '
grdfile = netCDF4.Dataset(roms.rundir + roms.run_name + '_grd.nc')
# assigning some variables from grid file
lon = grdfile.variables['lon_rho'][:]
lat = grdfile.variables['lat_rho'][:]
h = grdfile.variables['h'][:]
grdfile = netCDF4.Dataset(roms.rundir + 'phd13_grd.nc')
h2 = grdfile.variables['h'][:]
fig1 = plt.figure(facecolor='w', figsize=(16,10))
m = Basemap(projection='merc',
llcrnrlat = lat.min(), urcrnrlat = lat.max(),
# SCRIPT START ######################################################
# Basic Settings:
bry_filenamestr = '_bry.nc'
bry_filetypestr = 'ROMS Boundary Forcing file'
# READING PREVIOUSLY BUILT RELEVANT FILES: ###########################
# metadata ascii file
# grid netcdf file
print ' \n' + '==> ' + ' READING ASCII METADATA FILE ...\n' + ' '
Run = run_setup('piloto.setup')
print ' \n' + '==> ' + ' READING GRID NETCDF FILE ...\n' + ' '
grdfile = netCDF4.Dataset(Run.datadir + Run.run_name + '_grd.nc')
# WRITING THE NETCDF FILE ####################################################
# Based on "bry_unlimit.cdl" NETCDF sample structure
# some computings regarding boundary limits:
print ' \n' + '==> ' + ' WRITING NETCDF Bry FILE ...\n' + ' '
ncfile = netCDF4.Dataset(run.datadir + run.run_name + bry_filenamestr, mode='w',
clobber='true', format='NETCDF3_CLASSIC')
from roms_setup import hanning_smoother_coef2d, FX, FY
# SCRIPT START ######################################################
# Basic Settings:
filenamestr = '_grd.nc'
filetypestr = 'ROMS Grid file'
# READING PREVIOUSLY BUILT RELEVANT FILES: ###########################
# metadata ascii file
# OA-created netcdf initial T, S file
# grid netcdf file
print ' \n' + '==> ' + ' READING ASCII METADATA FILE ...\n' + ' '
run = run_setup('../remo_relax.setup')
grid = nc.Dataset('/data_0/REMO/ROMS_tsheld/mesoscale_grd.nc')
lon = grid.variables['lon_rho'][:]
lat = grid.variables['lat_rho'][:]
coast = sp.loadmat('costa.mat')
xc = coast['ncst'][:,0]
yc = coast['ncst'][:,1]
dl = run.resolution
print('Pick the four corners of the grid starting\n from lower left in clockwise direction')
plt.figure()
