Python 为什么netCDF数据和ferret中的存储值不同

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

在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
 **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

如果你知道答案,请帮助我