f = np.where(h2 >= 5000) h2[f] = 5000 del f Jrho, Irho = rlon2.shape Mr2, Lr2 = rlon2.shape Lu2 = Lr2 - 1 Mu2 = Mr2 Lv2 = Lr2 Mv2 = Mr2 - 1 cosa = np.cos(angle) sina = np.sin(angle) del angle rmask2 = np.ma.masked_where(rmask2 == 0, rmask2) [Zsig, dZsig] = zlev(h2, run.theta_s, run.theta_b, run.tcline, run.klevels) ### Interpolating T, S to ROMS 3D S-COORD grid ############################### # Currently doing only TEMP and SALT. ZETA, UBAR, VBAR, U, V are zero arrays lN = run.klevels lt = np.size(run.time) ZETA = np.zeros([lt, Jrho, Irho]) UBAR = np.zeros([lt, Mu2, Lu2]) VBAR = np.zeros([lt, Mv2, Lv2]) U = np.zeros([lt, lN, Mu2, Lu2]) V = np.zeros([lt, lN, Mv2, Lv2]) z2 = np.zeros([N1, Jrho, Irho]) Zlev2 = np.zeros([N1, 1]) print ' \n' + '==> ' + ' INTERPOLATING TEMPERATURE ...\n' + ' '
# Modify the bathymetry f = np.where(h2 >= 5000) h2[f] = 5000; del f N = int(run.klevels) Jrho, Irho = rlon2.shape Mr2, Lr2 = rlon2.shape Lu2 = Lr2-1; Mu2 = Mr2 Lv2 = Lr2; Mv2 = Mr2-1 cosa = np.cos(angle); sina = np.sin(angle); del angle rmask2 = np.ma.masked_where(rmask2 == 0, rmask2) hu = griddata(rlon2.ravel(), rlat2.ravel(), h2.ravel(), ulon2, ulat2) hv = griddata(rlon2.ravel(), rlat2.ravel(), h2.ravel(), vlon2, vlat2) [Zsig,dZsig] = zlev(h2,run.theta_s,run.theta_b,run.tcline,run.klevels) [ZsigU,dZsigU] = zlev(hu,run.theta_s,run.theta_b,run.tcline,run.klevels) [ZsigV,dZsigV] = zlev(hv,run.theta_s,run.theta_b,run.tcline,run.klevels) ### Interpolating T, S to ROMS 3D S-COORD grid ############################### lN = run.klevels lt = np.size(run.time) ZETA = np.zeros([lt, Jrho, Irho]) UBAR = np.zeros([lt, Mu2, Lu2]) VBAR = np.zeros([lt, Mv2, Lv2]) TEMP = np.zeros([lt, N, Mv2, Lv2]) SALT = np.zeros([lt, N, Mv2, Lv2]) U = np.zeros([lt, N, Mu2, Lu2]) V = np.zeros([lt, N, Mv2, Lv2])
# DOING COMPUTATIONS TO INTERPOLATE THE FIELDS TO ROMS GRID ######### # Modify the bathymetry f = np.where(h2 >= 5000) h2[f] = 5000; del f Jrho, Irho = rlon2.shape Mr2, Lr2 = rlon2.shape Lu2 = Lr2-1; Mu2 = Mr2 Lv2 = Lr2; Mv2 = Mr2-1 cosa = np.cos(angle); sina = np.sin(angle); del angle rmask2 = np.ma.masked_where(rmask2 == 0, rmask2) [Zsig,dZsig] = zlev(h2,run.theta_s,run.theta_b,run.tcline,run.klevels) ### Interpolating T, S to ROMS 3D S-COORD grid ############################### # Currently doing only TEMP and SALT. ZETA, UBAR, VBAR, U, V are zero arrays lN = run.klevels lt = np.size(run.time) ZETA = np.zeros([lt, Jrho, Irho]) UBAR = np.zeros([lt, Mu2, Lu2]) VBAR = np.zeros([lt, Mv2, Lv2]) U = np.zeros([lt, lN, Mu2, Lu2]) V = np.zeros([lt, lN, Mv2, Lv2]) z2 = np.zeros([N1, Jrho, Irho]) Zlev2 = np.zeros([N1, 1])