NCL处理数据

loftydongshi

请教大家给指点错误！附件中的数据spectral_irradiance_Lean_1610-2009_ann_c100405.nc一共有8个变量，具体文件说明如下

dimensions:

     ncl_scalar = 1

     time = 400  // unlimited

     ref_time_bound = 2

     wavelength = 3780

  variables:

     double tsi_ref ( ncl_scalar )

        time_op :        average

        bounds : ref_time_bound

        _FillValue :      -9999

        long_name :  average of tsi overref_time_bound inclusive

        units :     W/m^2

        average_op_ncl :   dim_avg overdimension: time

     double ref_time_bound ( ref_time_bound )

        units :     days since 0000-01-0100:00:00

        long_name :  reference timeinterval boundaries

        _FillValue :      -9999

     integer date ( time )

        format : YYYYMMDD

        _FillValue :      -9999

     double time ( time )

        calendar :       noleap

        axis :       T

        time_origin :  01-JAN-0000

        units :     days since 0000-01-0100:00:00

        _FillValue :      -9999

     double wavelength ( wavelength )

        units :     nm

        long_name :  Wavelength of bandcenter

        _FillValue :      -9999

     double band_width ( wavelength )

        units :     nm

        long_name :  Wavelength width ofband

        _FillValue :      -9999

     double ssi_ref ( wavelength )

        time_op :        average

        bounds : ref_time_bound

        _FillValue :      -9999

        long_name :  average of ssi overref_time_bound inclusive

        units :     mW/m^2/nm

        average_op_ncl :   dim_avg overdimension: time

     double tsi ( time )

        units :     W/m^2

        long_name :  Total Solar Irradianceat 1 a.u.

        _FillValue :      -9999

     double ssi ( time, wavelength )

        units :     mW/m^2/nm

        long_name :  Solar SpectralIrradiance at 1 a.u.

        _FillValue :      -9999

我主要是想对源数据ssi中的一部分进行放大，然后再将放大的这部分赋值给源数据，其他不变。

这是我的ncl脚本，处理之后我挨着变量看，数值都正确，但是把这个数据放到模式中，数据有错误，我想是不是格式有问题，请大家帮着看一下，十分的感谢。

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"

;---------------------

begin

;---------------------

in= addfile("/cmb/g5/dongshi/ssi/spectral_irradiance_Lean_1610-2009_ann_c100405_1.nc","w")

ssi= in->ssi

time=in->time

wavelength=in->wavelength

ntime=dimsizes(time)

nwavelength=dimsizes(wavelength)

printVarSummary(ssi)

cssi= new((/400,3780/),"double",ssi@_FillValue)

do ilt=0,ntime-1

  do ilw=0,nwavelength-1

    if((ilw.le.249).and.(ilw.gt.79)) then

    cssi(ilt,ilw)=ssi(ilt,ilw)*0.8

    else

    cssi(ilt,ilw)=ssi(ilt,ilw)

    end if

  end do

end do

copy_VarAtts(ssi,cssi)

printVarSummary(cssi)

in->ssi=cssi

end

longlivehj

“挨着变量看，数值都正确；放到模式中，数据有错误”，到底是正确还是错误啊？错误的详情？
程序看上去没啥问题，有两点看法供楼主参考：
1. 从nc文件的命名来看，楼主是做过备份了，这个好！但是，仍要小心对同一个文件重复操作，比如调试程序时反复运行程序。
2. 循环会大大降低程序执行效率。程序中的两层循环是没有必要的，直接指定数组下标进行修改。

loftydongshi

longlivehj 发表于 2014-9-5 12:00
“挨着变量看，数值都正确；放到模式中，数据有错误”，到底是正确还是错误啊？错误的详情？
程序看上去没 ...

谢谢版主。
运行是没有错误的，运行的数值也是正确的。
但是运行后的数据（spectral_irradiance_Lean_1610-2009_ann_c100405_1.nc）与原始数据（spectral_irradiance_Lean_1610-2009_ann_c100405.nc）除了数值之外，还有其他的问题。

昨天下午做了调整已经可以了
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"


;---------------------
begin
;---------------------

in= addfile("/cmb/g5/dongshi/ssi/spectral_irradiance_Lean_1610-2009_ann_c100405_20up.nc","w")

ssi= in->ssi

printVarSummary(ssi)

  do ilw=0,3779
    if((ilw .le. 249) .and. (ilw .gt. 79)) then
    ssi(:,ilw)=ssi(:,ilw)*1.2
    end if
  end do

in->ssi=ssi

end

loftydongshi

longlivehj 发表于 2014-9-5 12:00
“挨着变量看，数值都正确；放到模式中，数据有错误”，到底是正确还是错误啊？错误的详情？
程序看上去没 ...

是原来想复杂了，确实是循环嵌套了两层，这NCL直接罢工了。。

longlivehj

loftydongshi 发表于 2014-9-6 10:21
是原来想复杂了，确实是循环嵌套了两层，这NCL直接罢工了。。

不用循环，把下标范围指定一下！

loftydongshi

longlivehj 发表于 2014-9-6 11:18
不用循环，把下标范围指定一下！

我真是犯二了。。就是嘛。。。（：，70:249）不就得了。。。。转模式转的头都晕了。。。

走路睡觉的猫WD

正好学习一下

imaJUNation

{:eb502:}{:eb502:}{:eb502:}{:eb502:}{:eb502:}正好学习一下

四月雨之音

楼主我想问个问题，我在处理数据的时候需要把一个三维的（year,lat,lon)付给四维的（year,day,lat,lon)第九十一天的数值day=90,出现fatal: Assignment type mismatch, right hand side can't be coerced to type o left hand side三次。。怎么处理比较合适呢

loftydongshi

四月雨之音 发表于 2016-7-16 20:34
楼主我想问个问题，我在处理数据的时候需要把一个三维的（year,lat,lon)付给四维的（year,day,lat,lon)第九 ...

没大懂你的意思，第九十一天的数值？如果三维的数组付给四维数组时已经定义好单点某个时刻就不会出现这种问题了，可否把程序贴一下