Python ztosigmaの例

コード例 #1

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);

###

コード例 #2

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)

コード例 #3

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)

コード例 #4

ファイル: make_ini.py プロジェクト: rsoutelino/sandbox

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);

##############################################################################