[GRASS-SVN] r74488 - grass-addons/grass7/raster/r.green/r.green.gshp/r.green.gshp.theoretical

[svn_grass at osgeo.org]

Author: Giulia
Date: 2019-05-16 02:37:46 -0700 (Thu, 16 May 2019)
New Revision: 74488

Modified:
   grass-addons/grass7/raster/r.green/r.green.gshp/r.green.gshp.theoretical/r.green.gshp.theoretical.py
Log:
only positive values

Modified: grass-addons/grass7/raster/r.green/r.green.gshp/r.green.gshp.theoretical/r.green.gshp.theoretical.py
===================================================================
--- grass-addons/grass7/raster/r.green/r.green.gshp/r.green.gshp.theoretical/r.green.gshp.theoretical.py	2019-05-16 03:28:36 UTC (rev 74487)
+++ grass-addons/grass7/raster/r.green/r.green.gshp/r.green.gshp.theoretical/r.green.gshp.theoretical.py	2019-05-16 09:37:46 UTC (rev 74488)
@@ -222,6 +222,7 @@
 
 from grass.script import core as gcore
 from grass.script.utils import set_path
+from grass.script import mapcalc
 
 
 try:
@@ -329,6 +330,8 @@
     gpot.r_power(power, tc, ground_conductivity, ground_temperature,
                  fluid_limit_temperature, borehole_length, borehole_resistence,
                  gmax, execute=True, overwrite=OVER)
+    command = "{new} = if({old}<0, null(), {old})".format(old=power, new=power)
+    mapcalc(command, overwrite=True)
 
     energy = opts['energy']
     gpot.r_energy(energy, power, execute=True, overwrite=OVER)