[GRASS-dev] python script problem...

Massimiliano Cannata massimiliano.cannata at supsi.ch
Tue Oct 16 11:18:42 EDT 2007


Hi all,
I'm new to python and I'm trying to build up a script...
My question is:
After I setted up an option of type double, How can I get it in a python 
variable?

h = float(os.getenv('GIS_OPT_h')) it doesn't work, it require a string 
or numeric attribute !?!?

how to do?

thanks,
Maxi
------------------------
Here is my script:
-----------------------
#!/usr/bin/python
# -*- coding:utf-8 -*-
import math

#%Module
#%  description: g.parser test script
#%  keywords: keyword1, keyword2
#%End
#%option
#%  key: h
#%  type: double
#%  description: Profondità del mare in prossimità della costa h
#%  required : yes
#%end
#%option
#%  key: ho
#%  type: double
#%  description: Profondità del mare in corrispondenza dell'epicentro
#%  required : yes
#%end
#%option
#%  key: Ho
#%  type: double
#%  description: Altezza d'onda nel punto di generazione dello tsunami
#%  required : yes
#%end

import sys
import os

def main():

#### add your code here ####
    print ""
    print "Value of GIS_OPT_h: %s" % os.getenv("GIS_OPT_h")
    print "Value of GIS_OPT_ho:  %s" % os.getenv("GIS_OPT_ho")
    print "Value of GIS_OPT_Ho:  %s" % os.getenv("GIS_OPT_Ho")

#### it works and print the user input ####

h = float(os.getenv('GIS_OPT_h'))
ho = os.getenv('GIS_OPT_ho')
Ho = os.getenv('GIS_OPT_Ho')

#### NOT WORKING ! i guess i need something similar to atof in C ANSI ####

Co=math.sqrt(9.8 * h)
print "celerità dell'onda=",Co

H= Ho * ((ho/h)**(0.25))
print "Altezza d'onda a riva=",H

U=Co * math.sqrt((1+H/h))
print "Velocità d'onda a riva=",U

L= 4 * 2.415 *h /(math.sqrt(3*H/h))
print "Lunghezza di un'onda solitaria=",L

k1=U**2 * L/(19.6)
print "k1=",k1

k2= U**2 * H/(39.2)
print "k2=",k2

#### end of your code ####
if __name__ == "__main__":
    if ( len(sys.argv) <= 1 or sys.argv[1] != "@ARGS_PARSED@" ):
        os.execvp("g.parser", [sys.argv[0]] + sys.argv)
    else:
    main();





-- 
-----------------------------------------------------
Dr. Massimiliano Cannata
Environmental & Geomatic Engineer

Institute of Earth Scences - SUPSI
Trevano, C.P. 72, CH-6952 Canobbio, SWITZERLAND

Tel:    +41 (0)58 / 666 62 14  
Fax:    +41 (0)58 / 666 62 09
E-mail: massimiliano.cannata at supsi.ch

Web: 
http://www.ist.supsi.ch
http://istgis.ist.supsi.ch:8001/geomatica/
-------------------------------------------------------





More information about the grass-dev mailing list