Philippe 09/11/2017: mppdb: bug correction + minor improvements

src/LIB/SURCOUCHE/src/mode_mppdb.f90

@@ -73,7 +73,7 @@ CONTAINS
 
     !NMNH_COMM_WORLD = MPI_COMM_WORLD
 
-    ! If already initialized , no think to do
+    ! If already initialized, nothing to do
     IF (MPPDB_INITIALIZED) RETURN
     !
     MPPDB_INITIALIZED = .TRUE.

src/MNH/prep_nest_pgd.f90

@@ -171,8 +171,6 @@ TYPE(TPTR2FILE),DIMENSION(:),ALLOCATABLE,TARGET :: TZFILENESTPGD ! Output PGD fi
 !
 CALL MPPDB_INIT()
 !
-CALL MPPDB_INIT()
-!
 CALL VERSION
 CPROGRAM='NESPGD'
 !

src/MNH/prep_pgd.f90

@@ -153,7 +153,6 @@ CPROGRAM='PGD   '
 !
 TZFILE => NULL()
 !
-CALL MPPDB_INIT()
 !*    1.      Set default names and parallelized I/O
 !             --------------------------------------
 !

src/MNH/retrieve2_nest_infon.f90

@@ -261,9 +261,6 @@ IF(.NOT.LCARTESIAN) THEN
                  XXHAT(JPHEXT+1),XYHAT(JPHEXT+1), &
                  ZLAT2,ZLON2)
 
-  !20131024 MPPDB CHECK
-  CALL MPPDB_CHECK2D(ZPGDLAT1,"retrieve2_nest_info:ZPGDLAT1",PRECISION)
-  CALL MPPDB_CHECK2D(ZPGDLON1,"retrieve2_nest_info:ZPGDLON1",PRECISION)
 !
 !*      1.2   latitude and longitude of all local flux points (model1)
 !             --------------------------------------------------
@@ -273,6 +270,8 @@ IF(.NOT.LCARTESIAN) THEN
   CALL SM_LATLON(XPGDLATOR,XPGDLONOR,                                 &
                  SPREAD(XPGDXHAT(:),2,IJU),SPREAD(XPGDYHAT(:),1,IIU), &
                  ZPGDLAT1(:,:),ZPGDLON1(:,:))
+  CALL MPPDB_CHECK2D(ZPGDLAT1,"retrieve2_nest_info:ZPGDLAT1",PRECISION)
+  CALL MPPDB_CHECK2D(ZPGDLON1,"retrieve2_nest_info:ZPGDLON1",PRECISION)
 ENDIF  
 !
 !*      1.3   KXOR_C_ll, KYOR_C_ll - origin (global) of son model in father grid