Python 为什么netCDF数据和ferret中的存储值不同
在Anaconda虚拟环境中 使用雪貂库Python 为什么netCDF数据和ferret中的存储值不同,python,netcdf,ferret,Python,Netcdf,Ferret,在Anaconda虚拟环境中 使用雪貂库 (FERRET) [sion@chukchi20 3.4.2]$ ferret NOAA/PMEL TMAP PyFerret v7.63 (optimized) Linux 4.15.0-1096-azure - 10/13/20 24-May-21 17:17 *** NOTE: Unable to create journal file ferret.jnl
(FERRET) [sion@chukchi20 3.4.2]$ ferret
NOAA/PMEL TMAP
PyFerret v7.63 (optimized)
Linux 4.15.0-1096-azure - 10/13/20
24-May-21 17:17
*** NOTE: Unable to create journal file ferret.jnl
**ERROR: invalid command: Unable to open journal file: ferret.jnl
yes? use soda3.4.2_mn_ocean_reg_1980.nc
*** NOTE: netCDF parent axis definition error
*** NOTE: Edges definition "st_edges_ocean" points to no existing axis
*** NOTE: Edge definitions ignored
*** NOTE: netCDF parent axis definition error
*** NOTE: Edges definition "sw_edges_ocean" points to no existing axis
*** NOTE: Edge definitions ignored
yes? show data
currently SET data sets:
1> ./soda3.4.2_mn_ocean_reg_1980.nc (default)
name title I J K L
TEMP Potential temperature 1:720 1:330 1:50 1:12
SALT Practical Salinity 1:720 1:330 1:50 1:12
WT dia-surface velocity T-points 1:720 1:330 1:50 1:12
SSH effective sea level (eta_t + pa 1:720 1:330 ... 1:12
MLT mixed layer depth determined by 1:720 1:330 ... 1:12
MLP Depth of potential density mixe 1:720 1:330 ... 1:12
MLS mixed layer depth determined by 1:720 1:330 ... 1:12
NET_HEATING
surface ocean heat flux coming 1:720 1:330 ... 1:12
PRHO potential density referenced to 1:720 1:330 1:50 1:12
U i-current 1:720 1:330 1:50 1:12
V j-current 1:720 1:330 1:50 1:12
TAUX i-directed wind stress forcing 1:720 1:330 ... 1:12
TAUY j-directed wind stress forcing 1:720 1:330 ... 1:12
从1980年到2019年,所有数据都被确认具有相同的形式。
因此,我创建了一个列表,使用下面的代码只对某些区域进行切片
(FERRET) [sion@chukchi20 ~]$ awk '{print "use \""$NF"\"; set region/x=100:300/y=-30:30/M=1; (save/file=Eq_Pac_TSUV.nc/append Temp,Salt,SSH,U,V); cancel data/all;"}' list > slicing.jnl
还有,在雪貂
(FERRET) [sion@chukchi20 ~]$ ferret
NOAA/PMEL TMAP
PyFerret v7.63 (optimized)
Linux 4.15.0-1096-azure - 10/13/20
24-May-21 17:06
yes? go slicing.jnl
use "/data00/MLearning/SJTuna/SODA/3.4.2/soda3.4.2_mn_ocean_reg_1980.nc"; set region/x=100:300/y=-30:30/M=1; (save/file=Eq_Pac_TSUV.nc/append Temp,Salt,SSH,U,V); cancel data/all;
*** NOTE: netCDF parent axis definition error
*** NOTE: Edges definition "st_edges_ocean" points to no existing axis
*** NOTE: Edge definitions ignored
*** NOTE: netCDF parent axis definition error
*** NOTE: Edges definition "sw_edges_ocean" points to no existing axis
*** NOTE: Edge definitions ignored
LISTing to file Eq_Pac_TSUV.nc
use "/data00/MLearning/SJTuna/SODA/3.4.2/soda3.4.2_mn_ocean_reg_1981.nc"; set region/x=100:300/y=-30:30/M=1; (save/file=Eq_Pac_TSUV.nc/append Temp,Salt,SSH,U,V); cancel data/all;
*** NOTE: netCDF parent axis definition error
*** NOTE: Edges definition "st_edges_ocean" points to no existing axis
*** NOTE: Edge definitions ignored
*** NOTE: netCDF parent axis definition error
*** NOTE: Edges definition "sw_edges_ocean" points to no existing axis
*** NOTE: Edge definitions ignored
LISTing to file Eq_Pac_TSUV.nc
... skip ...
use "/data00/MLearning/SJTuna/SODA/3.4.2/soda3.4.2_mn_ocean_reg_2018.nc"; set region/x=100:300/y=-30:30/M=1; (save/file=Eq_Pac_TSUV.nc/append Temp,Salt,SSH,U,V); cancel data/all;
*** NOTE: netCDF parent axis definition error
*** NOTE: Edges definition "st_edges_ocean" points to no existing axis
*** NOTE: Edge definitions ignored
*** NOTE: netCDF parent axis definition error
*** NOTE: Edges definition "sw_edges_ocean" points to no existing axis
*** NOTE: Edge definitions ignored
LISTing to file Eq_Pac_TSUV.nc
use "/data00/MLearning/SJTuna/SODA/3.4.2/soda3.4.2_mn_ocean_reg_2019.nc"; set region/x=100:300/y=-30:30/M=1; (save/file=Eq_Pac_TSUV.nc/append Temp,Salt,SSH,U,V); cancel data/all;
*** NOTE: netCDF parent axis definition error
*** NOTE: Edges definition "st_edges_ocean" points to no existing axis
*** NOTE: Edge definitions ignored
*** NOTE: netCDF parent axis definition error
*** NOTE: Edges definition "sw_edges_ocean" points to no existing axis
*** NOTE: Edge definitions ignored
LISTing to file Eq_Pac_TSUV.nc
最后,我检查了Eq_Pac_TSUV.nc文件的内容
yes? use Eq_Pac_TSUV.nc
yes? show data
currently SET data sets:
1> ./Eq_Pac_TSUV.nc (default)
name title I J K L
TEMP Potential temperature 1:400 1:120 1:50 1:499
SALT Practical Salinity 1:400 1:120 1:50 1:499
SSH effective sea level (eta_t + pa 1:400 1:120 ... 1:499
U i-current 1:401 1:121 1:50 1:499
V j-current 1:401 1:121 1:50 1:499
以U,V的名义。一、 J,有一个不同的范围。
为什么范围是1:401,1:121?
我所期望的是1:400,1:121
如果你知道答案,请帮助我