##################################################################### # 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()