--- trunk/libf/dyn3d/start_init_phys_m.f90	2010/04/06 17:52:58	32
+++ trunk/libf/dyn3d/start_init_phys_m.f90	2011/03/24 11:52:41	42
@@ -5,7 +5,7 @@
 
   IMPLICIT NONE
 
-  REAL, ALLOCATABLE, SAVE, DIMENSION(:, :):: qsol_2d
+  REAL, ALLOCATABLE, SAVE:: qsol_2d(:, :)
 
 CONTAINS
 
@@ -21,7 +21,7 @@
 
     REAL, intent(out):: tsol_2d(:, :)
 
-    !  LOCAL
+    ! Variables local to the procedure:
 
     INTEGER fid_phys, iml_phys, jml_phys
     REAL, ALLOCATABLE, DIMENSION(:, :):: lon_phys, lat_phys
@@ -31,9 +31,6 @@
     INTEGER:: itau(1)
     INTEGER::  llm_tmp, ttm_tmp
 
-    CHARACTER(len=120) physfname
-    LOGICAL:: check=.TRUE.
-
     REAL, ALLOCATABLE:: lon_rad(:), lat_rad(:)
     REAL, ALLOCATABLE:: lon_ini(:), lat_ini(:)
     REAL, ALLOCATABLE:: var_ana(:, :)
@@ -43,17 +40,14 @@
 
     print *, "Call sequence information: start_init_phys"
     if (any(shape(tsol_2d) /= (/iim + 1, jjm + 1/))) stop "start_init_phys"
-    physfname = 'ECPHY.nc'
-    IF ( check ) print *, 'Opening the surface analysis'
-    CALL flininfo(physfname, iml_phys, jml_phys, llm_tmp, ttm_tmp, fid_phys)
+    CALL flininfo('ECPHY.nc', iml_phys, jml_phys, llm_tmp, ttm_tmp, fid_phys)
 
     ALLOCATE(lat_phys(iml_phys, jml_phys))
     ALLOCATE(lon_phys(iml_phys, jml_phys))
     ALLOCATE(levphys_ini(llm_tmp))
 
-    CALL flinopen_nozoom(iml_phys, jml_phys,  &
-         llm_tmp, lon_phys, lat_phys, levphys_ini, ttm_tmp,  &
-         itau, date, dt, fid_phys)
+    CALL flinopen_nozoom(iml_phys, jml_phys, llm_tmp, lon_phys, lat_phys, &
+         levphys_ini, ttm_tmp, itau, date, dt, fid_phys)
 
     DEALLOCATE(levphys_ini)
 
@@ -79,9 +73,8 @@
        lat_ini(:) = lat_phys(1, :) 
     ENDIF
 
-    !   We get the two standard varibales
-    !   Surface temperature
-    ! 'ST'            : Surface temperature
+    ! We get the two standard variables
+    ! 'ST': surface temperature
     CALL flinget(fid_phys, 'ST', iml_phys, jml_phys,  &
          llm_tmp, ttm_tmp, 1, 1, var_ana)
     CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, var_ana)