[GRASS5] Re: r.sun improvements

Jaro Hofierka hofierka at geomodel.sk
Thu Apr 10 05:20:02 EDT 2003


Paul, thanks for the fix. It works fine for me.

Jaro

Paul Kelly wrote:

>Hello Jaro and Hamish
>
>On Thu, 10 Apr 2003, Jaro Hofierka wrote:
>
>  
>
>>However, there is another problem with memory associated with proj
>>library. If you have latitude defined via region settings (i.e. you
>>don't define lat and latin options) then memory usage is unbelievably
>>high. I guess it is some memory allocation/freeing problem in  proj
>>library. I have used proj library functions to get latitude/longitute
>>for each ratser cell using examples from grass programmers manual. I
>>    
>>
>
>I didn't realise r.sun calls the proj functions so often. Maybe it needs
>to free the memory used for the key/value pairs read from the PROJ_INFO
>file. Does the patch below help? I don't use r.sun so can't really test
>it. When I have time I mean to update the example in the programmer's manual.
>
>Paul
>
>Index: main.c
>===================================================================
>RCS file: /grassrepository/grass/src/raster/r.sun/main.c,v
>retrieving revision 1.11
>diff -u -r1.11 main.c
>--- main.c	7 Apr 2003 15:15:20 -0000	1.11
>+++ main.c	10 Apr 2003 08:28:22 -0000
>@@ -1463,6 +1463,11 @@
> 	        if (pj_get_kv(&oproj, out_proj_info, out_unit_info) < 0)
> 	            G_fatal_error("Unable to set up lat/long projection parameters");
>
>+                G_free_key_value( in_proj_info );
>+                G_free_key_value( in_unit_info );
>+                G_free_key_value( out_proj_info );
>+                G_free_key_value( out_unit_info );
>+
> 	        if(pj_do_proj(&longitude,&latitude,&iproj,&oproj) < 0)
> 	        {
> 	                fprintf(stderr,"Error in pj_do_proj\n");
>
>_______________________________________________
>grass5 mailing list
>grass5 at grass.itc.it
>http://grass.itc.it/mailman/listinfo/grass5
>
>  
>





More information about the grass-dev mailing list