The HARRY_READ_ME.txt file

Part 21

21. A little experimentation goes a short way..

I tried using the 'stn' option of anomdtb.for. Not completely sure what
it's supposed to do, but no matter as it didn't work:

crua6[/cru/cruts/rerun1/data/cruts/rerun_pre] ./anomdtb

> ***** AnomDTB: converts .dtb to anom .txt for gridding *****

> Enter the suffix of the variable required:
.pre
> Will calculate percentage anomalies.
> Select the .cts or .dtb file to load:
pre.0312031600H.dtb

> Specify the start,end of the normals period:
1961,1990
> Specify the missing percentage permitted:
25
> Data required for a normal: 23
> Specify the no. of stdevs at which to reject data:
5
> Select outputs (1=.cts,2=.ann,3=.txt,4=.stn):
4
> Check for duplicate stns after anomalising? (0=no,>0=km range)
8
> Select the .stn file to save:
pre.fromanomdtb.stn
> Enter the correlation decay distance:
450
> Submit a grim that contains the appropriate grid.
> Enter the grim filepath:
cru_ts_2_10.1961-1970.pre

> Grid dimensions and domain size: 720 360 67420
> Select the first,last years AD to save:
1901,2002
> Operating...

> NORMALS MEAN percent STDEV percent
> .dtb 2635548 29.6
> .cts 4711327 52.8 7325296 82.2
> PROCESS DECISION percent %of-chk
> no lat/lon 20761 0.2 0.2
> no normal 1585342 17.8 17.8
> out-of-range 20249 0.2 0.3
> duplicated 317035 3.6 4.3
> accepted 6972308 78.2
> Calculating station coverages...
> ##### WithinRange: Alloc: DataB #####
forrtl: severe (174): SIGSEGV, segmentation fault occurred
crua6[/cru/cruts/rerun1/data/cruts/rerun_pre]

..knowing how long it takes to debug this suite - the experiment
endeth here. The option (like all the anomdtb options) is totally
undocumented so we'll never know what we lost.


Go on to part 22, back to index or Email search