SALT = np.zeros([lt, N, Mv2, Lv2])
U = np.zeros([lt, N, Mu2, Lu2])
V = np.zeros([lt, N, Mv2, Lv2])
z2 = np.zeros([N1, Jrho, Irho])
Zlev2 = np.zeros([N1, 1])
print ' \n' + '==> ' + ' INTERPOLATING TEMPERATURE ...\n' + ' '
for k in np.arange(0, N1, 1):
print 'TEMP: Z Level = ' + str(-1*Zlev[k]) + ' m'
z1 = np.squeeze(temp[k,:,:])
z2[N1-k-1,:,:] = griddata(lon.ravel(),lat.ravel(),z1.ravel(),rlon2,rlat2)
Zlev2[N1-k-1] = Zlev[k]
print ' \n' + '==> ' + ' INTERPOLATING TEMP FROM Z --> S COORD ...\n' + ' '
TEMP = ztosigma(z2,Zsig,Zlev2); del z1, z2
###
z2 = np.zeros([N1, Jrho, Irho])
print ' \n' + '==> ' + ' INTERPOLATING SALINITY ...\n' + ' '
for k in np.arange(0, N1, 1):
print 'SALT: Z Level = ' + str(-1*Zlev[k]) + ' m'
z1 = np.squeeze(salt[k,:,:])
z2[N1-k-1,:,:] = griddata(lon.ravel(),lat.ravel(),z1.ravel(),rlon2,rlat2)
Zlev2[N1-k-1] = Zlev[k]
print ' \n' + '==> ' + ' INTERPOLATING SALT FROM Z --> S COORD ...\n' + ' '
SALT = ztosigma(z2,Zsig,Zlev2);
###
U = np.zeros([lt, N, Mu2, Lu2])
V = np.zeros([lt, N, Mv2, Lv2])
z2 = np.zeros([N1, Jrho, Irho])
Zlev2 = np.zeros([N1, 1])
print ' \n' + '==> ' + ' INTERPOLATING TEMPERATURE ...\n' + ' '
for k in np.arange(0, N1, 1):
print 'TEMP: Z Level = ' + str(-1 * Zlev[k]) + ' m'
z1 = np.squeeze(temp[k, :, :])
z2[N1 - k - 1, :, :] = griddata(lon.ravel(), lat.ravel(), z1.ravel(),
rlon2, rlat2)
Zlev2[N1 - k - 1] = Zlev[k]
print ' \n' + '==> ' + ' INTERPOLATING TEMP FROM Z --> S COORD ...\n' + ' '
TEMP = ztosigma(z2, Zsig, Zlev2)
del z1, z2
###
z2 = np.zeros([N1, Jrho, Irho])
print ' \n' + '==> ' + ' INTERPOLATING SALINITY ...\n' + ' '
for k in np.arange(0, N1, 1):
print 'SALT: Z Level = ' + str(-1 * Zlev[k]) + ' m'
z1 = np.squeeze(salt[k, :, :])
z2[N1 - k - 1, :, :] = griddata(lon.ravel(), lat.ravel(), z1.ravel(),
rlon2, rlat2)
Zlev2[N1 - k - 1] = Zlev[k]
print ' \n' + '==> ' + ' INTERPOLATING SALT FROM Z --> S COORD ...\n' + ' '
SALT = ztosigma(z2, Zsig, Zlev2)
V = np.zeros([lt, lN, Mv2, Lv2])
z2 = np.zeros([N1, Jrho, Irho])
Zlev2 = np.zeros([N1, 1])
print ' \n' + '==> ' + ' INTERPOLATING TEMPERATURE ...\n' + ' '
for k in np.arange(0, N1, 1):
print 'TEMP: Z Level = ' + str(-1 * Zlev[k]) + ' m'
z1 = np.squeeze(temp[:, :, k])
z2[N1 - k - 1, :, :] = griddata(lon.ravel(), lat.ravel(), z1.ravel(),
rlon2, rlat2)
Zlev2[N1 - k - 1] = Zlev[k]
print ' \n' + '==> ' + ' INTERPOLATING TEMP FROM Z --> S COORD ...\n' + ' '
TEMP = ztosigma(z2, Zsig, Zlev2)
del z1, z2
###
z2 = np.zeros([N1, Jrho, Irho])
print ' \n' + '==> ' + ' INTERPOLATING SALINITY ...\n' + ' '
for k in np.arange(0, N1, 1):
print 'SALT: Z Level = ' + str(-1 * Zlev[k]) + ' m'
z1 = np.squeeze(salt[:, :, k])
z2[N1 - k - 1, :, :] = griddata(lon.ravel(), lat.ravel(), z1.ravel(),
rlon2, rlat2)
Zlev2[N1 - k - 1] = Zlev[k]
print ' \n' + '==> ' + ' INTERPOLATING SALT FROM Z --> S COORD ...\n' + ' '
SALT = ztosigma(z2, Zsig, Zlev2)
V = np.zeros([lt, lN, Mv2, Lv2])
z2 = np.zeros([N1, Jrho, Irho])
Zlev2 = np.zeros([N1, 1])
print ' \n' + '==> ' + ' INTERPOLATING TEMPERATURE ...\n' + ' '
for k in np.arange(0, N1, 1):
print 'TEMP: Z Level = ' + str(-1*Zlev[k]) + ' m'
z1 = np.squeeze(temp[:,:,k])
z2[N1-k-1,:,:] = griddata(lon.ravel(),lat.ravel(),z1.ravel(),rlon2,rlat2)
Zlev2[N1-k-1] = Zlev[k]
print ' \n' + '==> ' + ' INTERPOLATING TEMP FROM Z --> S COORD ...\n' + ' '
TEMP = ztosigma(z2,Zsig,Zlev2); del z1, z2
###
z2 = np.zeros([N1, Jrho, Irho])
print ' \n' + '==> ' + ' INTERPOLATING SALINITY ...\n' + ' '
for k in np.arange(0, N1, 1):
print 'SALT: Z Level = ' + str(-1*Zlev[k]) + ' m'
z1 = np.squeeze(salt[:,:,k])
z2[N1-k-1,:,:] = griddata(lon.ravel(),lat.ravel(),z1.ravel(),rlon2,rlat2)
Zlev2[N1-k-1] = Zlev[k]
print ' \n' + '==> ' + ' INTERPOLATING SALT FROM Z --> S COORD ...\n' + ' '
SALT = ztosigma(z2,Zsig,Zlev2);
##############################################################################