return(raster_name)

}

  #This functions mosaics tiles/files give a list of files. 

  #There are four options to mosaic:   use_sine_weights,use_edge,use_linear_weights, unweighted

  #Sine weights fits sine fuctions across rows and column producing elliptical/spherical patterns from center

  #             here e.g. dailyTmax and date!!

  #11)out_dir: output directory, default is NULL

  #12)algorithm: use R or python function

  #

  #OUTPUT:

  # Object is produced with 3 components:

    ###### PART 4: compute the weighted mean with the mosaic function #####

    if(algorithm=="python"){

      #The file to do the merge is /nobackupp6/aguzman4/climateLayers/sharedCode/gdal_merge_sum.py. Sample call below.

                                                    module_name="gdal_merge_sum.py",

                                                    input_file=filename_list_mosaics_prod_weights_m,

                                                    out_mosaic_name=out_mosaic_name_prod_weights_m)

      cmd_str2 <- mosaic_prod_weights_obj$cmd_str

      #write out python command used for mosaicing

      cmd_mosaic_logfile <- file.path(out_dir,paste("cmd_mosaic_",out_suffix,".txt",sep=""))

      writeLines(cmd_str1,con=cmd_mosaic_logfile) #weights files to mosaic 

      #writeLines(cmd_str2,con=file.path(out_dir,paste("cmd_mosaic_",out_suffix,".txt",sep=""))) #weights files to mosaic 

      cat(cmd_str2, file=cmd_mosaic_logfile, append=TRUE, sep = "\n")

    }

    if(algorithm=="R"){