SPI_Example/FSTD_InterpTime.tcl at main · kuohsienchang/SPI_Example · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
#!/bin/sh
# the next line restarts using tclsh \
exec $SPI_PATH/tclsh "$0" "$@"
#============================================================================
# Environnement Canada
# Centre Meteorologique Canadien
# 2100 Trans-Canadienne
# Dorval, Quebec
#
# Projet     : Exemple de scripts.
# Fichier    : FSTD_InterpTime.tcl
# Creation   : Mai 2000 - J.P. Gauthier - CMC/CMOE
# Description: Interpoler des chamsp dans le temps
#
# Parametres :
#
# Retour:
#
# Remarques  :
#
#============================================================================

package require TclData
#package require TclGeoEER
package require Logger

Log::Start [info script] 0.1

file delete $env(CI_DATA_OUT)/FSTD_InterpTime.fstd

#----- Ouvrir les fichiers d'entree (1) sortie (2)
fstdfile open 1 read  $env(CI_DATA_IN)/2005102612_012
fstdfile open 2 write $env(CI_DATA_OUT)/FSTD_InterpTime.fstd

#----- Lire les champs aux temps T(0) et T(1)
fstdfield read T0 1 -1 "" 12000 18 -1  "" ES
fstdfield read T1 1 -1 "" 12000 24 -1  "" ES

#----- Recuperer la date de validitee au temps T(0)
set stamp0 [fstdfield define T0 -DATEV]

#----- Interpoler a toute les demi-heure sur 3 heures
for { set hour 0.1 } { $hour < 3.0 } { set hour [expr $hour+0.1] } {

   set stamp [fstdstamp incr $stamp0 $hour]
   puts stderr $hour=$stamp
   fstdfield timeinterp TX T0 T1 $stamp
   fstdfield write TX 2 -32 False
}

fstdfile close 1 2

Log::End