Philippe 04/04/2017: OpenACC: porting to OpenPOWER + optimisations

* use data zones instead of data create (to prevent crashes on OpenPOWER
  with PGI up to 17.3)
* use atomic and force loop parallelisation for COUNTJVxD_DEVICE subroutines
* force parallelisation of some loops that were not automatically parallelised
* added NOKRR7 define key to remove optional argument part of the code
  of OpenACC kernels (to prevent crashes on OpenPOWER)

src/MNH/ice_adjust.f90

@@ -290,9 +290,6 @@ LOGICAL,DIMENSION(SIZE(PEXNREF,1),SIZE(PEXNREF,2),SIZE(PEXNREF,3)) :: GTEMP
 !
 !$acc declare device_resident(ZEXNS,ZT,ZCPH,ZLV,ZLS,ZW1,ZW2,ZW3,ZW4,ZW5,ZW6,ZW7,ZCND,GTEMP)
 !
-#ifdef _OPENACC
-PRINT *,'OPENACC: ICE_ADJUST being implemented'
-#endif
 !-------------------------------------------------------------------------------
 !
 !*       1.     PRELIMINARIES

src/MNH/resolved_cloud.f90

@@ -1026,9 +1026,6 @@ CALL ABORT
    END IF
 !
   CASE ('ICE3')
-#ifdef _OPENACC
-PRINT *,'OPENACC: RESOLVED_CLOUD::ICE3 being implemented'
-#endif
 !
 !*       9.     MIXED-PHASE MICROPHYSICAL SCHEME (WITH 3 ICE SPECIES)
 !               -----------------------------------------------------