diff --git a/mean_rdr.f90 b/mean_rdr.f90
new file mode 100644
index 0000000..65a3cbc
--- /dev/null
+++ b/mean_rdr.f90
@@ -0,0 +1,24 @@
+      REAL(r8) FUNCTION mean_rad(k,rad_xy)
+      use rcfl
+      use grd_str
+      
+      IMPLICIT NONE
+      INTEGER(i4):: j,i,count
+      INTEGER(i4),INTENT(IN):: k
+      REAL(r8),INTENT(IN):: rad_xy(grd%im,grd%jm,1)
+      REAL(r8):: meanxy
+     
+      print*, k
+      count=0
+      meanxy=0
+      do j=1,grd%jm
+        do i=1,grd%im
+          if(rad_xy(i,j,1)>0.0001) then
+            count=count+1
+            meanxy=meanxy+rad_xy(i,j,1)
+          end if
+        enddo
+      enddo
+      mean_rad=meanxy/count
+      return     
+      END FUNCTION mean_rad
diff --git a/rdrcorr.f90 b/rdrcorr.f90
new file mode 100644
index 0000000..eb13c03
--- /dev/null
+++ b/rdrcorr.f90
@@ -0,0 +1,78 @@
+subroutine rdrcorr
+
+!---------------------------------------------------------------------------
+!                                                                          !
+!    Copyright 2015 Anna teruzzi, OGS Trieste                              !
+!                                                                          !
+!    This file is part of OceanVar.                                          !
+!                                                                          !
+!    OceanVar is free software: you can redistribute it and/or modify.     !
+!    it under the terms of the GNU General Public License as published by  !
+!    the Free Software Foundation, either version 3 of the License, or     !
+!    (at your option) any later version.                                   !
+!                                                                          !
+!    OceanVar is distributed in the hope that it will be useful,           !
+!    but WITHOUT ANY WARRANTY; without even the implied warranty of        !
+!    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         !
+!    GNU General Public License for more details.                          !
+!                                                                          !
+!    You should have received a copy of the GNU General Public License     !
+!    along with OceanVar.  If not, see <http://www.gnu.org/licenses/>.       !
+!                                                                          !
+!--------------------------------------------------------------------------- 
+
+
+  use set_knd
+  use drv_str
+  use netcdf
+  use grd_str
+  use cns_str
+  use filenames
+  use rcfl
+
+  implicit none
+
+  INTEGER(i4)                    :: stat, ncid, idvar,imr,jmr,kmr
+
+!write(*,*)trim(RCORR_FILE)
+    stat = nf90_open(trim(RCORR_FILE), NF90_NOWRITE, ncid)
+    if (stat /= nf90_noerr) call netcdf_err(stat)
+
+! Get dimensions 
+      stat = nf90_inq_dimid (ncid, 'longitude', idvar)
+    if (stat /= nf90_noerr) call netcdf_err(stat)
+      stat = nf90_inquire_dimension (ncid, idvar, len = imr)
+    if (stat /= nf90_noerr) call netcdf_err(stat)
+      stat = nf90_inq_dimid (ncid, 'latitude', idvar)
+    if (stat /= nf90_noerr) call netcdf_err(stat)
+      stat = nf90_inquire_dimension (ncid, idvar, len = jmr)
+    if (stat /= nf90_noerr) call netcdf_err(stat)
+!Laura add depth
+      stat = nf90_inq_dimid (ncid, 'depth', idvar)
+    if (stat /= nf90_noerr) call netcdf_err(stat)
+      stat = nf90_inquire_dimension (ncid, idvar, len = kmr)
+    if (stat /= nf90_noerr) call netcdf_err(stat)
+
+! Check on dimensions
+    if ((imr .ne. grd%im).OR.(jmr.ne.grd%jm)) then
+       write(drv%dia,*)'Error: dimensions of rcorr different from grid ones'
+       call f_exit(24)
+    endif
+
+
+!  Allocate rcorr arrays
+     ALLOCATE ( rcf%Lxyz(grd%im,grd%jm,kmr))
+
+       stat = nf90_inq_varid (ncid, 'radius', idvar)
+    if (stat /= nf90_noerr) call netcdf_err(stat)
+      stat = nf90_get_var (ncid,idvar,rcf%Lxyz)
+    if (stat /= nf90_noerr) call netcdf_err(stat)
+
+!laura from km to meter
+    where (rcf%Lxyz<=0.0001) 
+       rcf%Lxyz=rcf%Lxyz/1000
+    end where
+     rcf%Lxyz= rcf%Lxyz*1000  !from km to meter
+
+end subroutine rdrcorr
+