Index: trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/subslope_mola.F90
===================================================================
--- trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/subslope_mola.F90	(revision 3302)
+++ trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/subslope_mola.F90	(revision 3303)
@@ -60,5 +60,5 @@
 !-------------Building of theta_mola and psi_mola 
       ! Assume that the mola file is to ben found in "datadir"
-      CALL mola(trim(datadir)//"/",iim, &
+      CALL mola(trim(datadir)//"/",&
                 ierr,theta_mola,psi_mola,resol,iim_mola,jjm_mola)
 
@@ -111,5 +111,5 @@
 !==========================================================================================
 
- SUBROUTINE mola(dset,iim,ierr,theta_mola,psi_mola,resol,iim_mola,jjm_mola)
+ SUBROUTINE mola(dset,ierr,theta_mola,psi_mola,resol,iim_mola,jjm_mola)
 
 
@@ -124,11 +124,13 @@
 logical output_messages
       parameter (output_messages=.true.)
-      double precision resol ! MOLA pixel/degree resolution
-      integer jjm_mola, iim_mola ! # of longitude and latitude MOLA data values 
+      double precision resol
+!c      parameter(resol=16) ! MOLA pixel/degree resolution
+
+
+      integer jjm_mola, iim_mola   ! # of longitude and latitude MOLA data values 
 
 ! Arguments
 ! inputs
       character*(*) dset ! Path to MCD datafiles 
-      integer :: iim
 
 ! outputs
@@ -141,7 +143,11 @@
 
 ! Local variables
+
       real latitude  ! north latitude (degrees)m
       real longitude ! east longitude (degrees)
-      character*140 molafile ! MOLA datafile
+      character*140 molafile ! MOLA datafile 
+!      data molafile/'mola_32.2.nc'/
+!c      data molafile/'mola16.nc'/
+!c      data molafile/'mola32.nc'/
 !c      real invresol
 !c      parameter(invresol=1./resol)
@@ -162,10 +168,8 @@
       double precision lat,lon ! longitude and latitude, local values (in degrees)
       real topogrid(iim_mola,jjm_mola) ! altmola in 'real' version
-      character(2) :: iim_char
 
 
 !C 1.1. Open MOLA file
-         write(iim_char,'(i2)') iim
-         molafile = dset//'mola'//iim_char//'.nc'
+         molafile=dset//'mola64.nc'
          if (output_messages) then
            write(*,*)"Loading MOLA topography from file ", &
@@ -574,2 +578,3 @@
 
 END MODULE subslope_mola_mod
+