integer igen_etan,igen_bdrag,igen_geoth
igen_geoth=0
& igen_geoth=iarr
if (igen_geoth.GT.0)
& call CTRL_MAP_GENARR2D(geothermalFlux,igen_geoth,myThid)