[Thread Prev][Thread Next][Index]

[ferret_users] calculation RMSD.



Hi folks,

I am trying to calculate the rmsd and bias using the following script but i am getting the following errors could you guys please help with the problem which i am facing.::::


yes? set mode metafie
yes?  use pottm2008.nc
yes? use pottmp.2008.nc
yes? set memory/size=500
 Cached data cleared from memory
yes? set reg/t=10-may-2008:31-may-2008/z=10
yes? let mtemp=temperature[d=1,l=@ave]
yes? let mpot=pottmp[d=2,l=@var]
yes? let sdtemp=temperature[d=1,l=@var]^0.5
yes? let sdpot=pottmp[d=2,l=@var]^0.5
yes? let dU=temperature[d=1]-pottmp[d=2,g=temperature[d=1]]
yes? let bias=dU[l=@ave]
yes? let rmsd=dU[l=@var]^0.5
yes? list mtemp,mpot,sdtemp,sdpot,bias,rmsd
 **ERROR: inconsistent sizes of data regions: X axis
          Variable 2 is I=001:360
          _expression_ is I=0001:2001
yes? go multi_view 3, 1, .8, .1, 0, .2, .1, .1
yes? plo/nolabel/vlim=-100:140/color=black/thick=3 temperature[d=1]
 **ERROR: dimensions improperly specified: specified data is not a line
           - its a 3D region: "plo/nolabel/vlim=-100:140/color=black/thick=3 temperature[d=1]"
yes? plo/nolabel/color=black/thick=3 temperature[d=1]
 **ERROR: dimensions improperly specified: specified data is not a line
           - its a 3D region: "plo/nolabel/color=black/thick=3 temperature[d=1]"
yes? plo/ov/nolabel/color=red/thick=3 pottmp[d=2]
 **ERROR: dimensions improperly specified: overlay is on a different axis
          "plo/ov/nolabel/color=red/thick=3 pottmp[d=2]"
yes?


the grids of the data are as follows:::

yes? use pottm2008.nc
yes? sh g temperature
    GRID GOZ1
 name       axis              # pts   start                end
 LONN2551_N551 LONGITUDE     2001mr   130W(-130)           30E(30)
 LAT959_2051 LATITUDE        1093 i   39.976S              39.976N
 LEV       DEPTH (m)           10 i-  0                    200
 TIME      TIME                32 r   08-MAY-2008 00:00    08-JUN-2008 00:00
 normal    E
 normal    F
yes? sh d
     currently SET data sets:
    1> ./pottm2008.nc  (default)
 name     title                             I         J         K         L         M         N
 TEMPERATURE
          potential temperature (degc)     1:2001    1:1093    1:10      1:32      ...       ...
 
yes? use pottm2008.nc
yes? sh g pottmp
 **ERROR: variable unknown or not in data set: POTTMP
yes? use pottmp.2008.nc
yes? sh g pottmp
    GRID GCN2
 name       axis              # pts   start                end
 LON       LONGITUDE          360mr   0.5E                 0.5W
 LAT       LATITUDE           418 r   74.5S                64.499N
 LEVEL     DEPTH (m)           40 i-  5                    4478
 TIME1     TIME                12 i   01-JAN-2008 00:00    01-DEC-2008 00:00
 normal    E
 normal    F
yes? sh d pottmp
 **ERROR: unknown data set:
yes? sh d
     currently SET data sets:
    1> ./pottm2008.nc
 name     title                             I         J         K         L         M         N
 TEMPERATURE
          potential temperature (degc)     1:2001    1:1093    1:10      1:32      ...       ...
 
    2> ./pottmp.2008.nc  (default)
 name     title                             I         J         K         L         M         N
 DATE                                      ...       ...       ...       1:12      ...       ...
 TIMEPLOT                                  ...       ...       ...       1:12      ...       ...
 POTTMP   Potential temperature            1:360     1:418     1:40      1:12      ...       ...
 
yes?


[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement