Changeset 4225 for dynamico_lmdz/simple_physics/phyparam/physics

Timestamp:

Jan 10, 2020, 3:07:23 PM (6 years ago)

Author:

dubos

Message:

simple_physics : minor changes to compile with gfortran and parse code using XCodeML F_Front

Location:

dynamico_lmdz/simple_physics/phyparam/physics

Files:

• logging.F90 (modified) (2 diffs)
• turbulence.F90 (modified) (3 diffs)

dynamico_lmdz/simple_physics/phyparam/physics/logging.F90

@@ -24 +24 @@
   END INTERFACE
 
-  ! This module provides a default implementation of flush_plugin but the top-level driver is welcome to override it.  
+#ifdef XCODEML
+  ! XCodeML cannot parse procedure pointers                                                                                                                                             
+
+#define default_flush_plugin flush_plugin
+#define default_log_gridpoint log_gridpoint_plugin
+
+#else
+
+  ! This module provides a default implementation of flush_plugin but the top-level driver is welcome to override it.                                                                    
   PROCEDURE(plugin), POINTER :: flush_plugin => default_flush_plugin
 
-  ! The top-level driver MUST provide an implementation for log_gridpoint_plugin
+  ! The top-level driver MUST provide an implementation for log_gridpoint_plugin                                                                                                         
   PROCEDURE(plugin_log_gridpoint), POINTER :: log_gridpoint_plugin => NULL()
+
+#endif
 
   INTEGER, PARAMETER :: linesize=10000, logging_bufsize=100
@@ -69 +79 @@
   END SUBROUTINE log_gridpoint
 
+  SUBROUTINE default_log_gridpoint(index, line)
+    INTEGER, INTENT(IN) :: index ! index of gridpoint                                                                                                                                    
+    CHARACTER(*), INTENT(OUT) :: line
+    line=''
+  END SUBROUTINE default_log_gridpoint
+
 END MODULE logging

dynamico_lmdz/simple_physics/phyparam/physics/turbulence.F90

@@ -248 +248 @@
     IF(lwrite) THEN
        ig=ngrid/2+1
-       WRITELOG(*,*) 'Pression (mbar) ,altitude (km),u,v,theta, rho dz'
+       WRITELOG(*,*) 'Pression (mbar) altitude (km),u,v,theta, rho dz'
        DO ilay=1,nlay
           WRITELOG(*,*) .01*pplay(ig,ilay),.001*pzlay(ig,ilay), &
                pu(ig,ilay),pv(ig,ilay),ph(ig,ilay),za(ig,ilay)
        ENDDO
-       WRITELOG(*,*) 'Pression (mbar) ,altitude (km),zb'
+       WRITELOG(*,*) 'Pression (mbar) altitude (km),zb'
        DO ilev=1,nlay
           WRITELOG(*,*) .01*pplev(ig,ilev),.001*pzlev(ig,ilev), &
@@ -291 +291 @@
        WRITELOG(*,*) 'coefficients Cd pour v et h'
        WRITELOG(*,*) zcdv(ngrid/2+1),zcdh(ngrid/2+1)
-       WRITELOG(*,*) ,'coefficients K pour v et h'
+       WRITELOG(*,*) 'coefficients K pour v et h'
        DO ilev=1,nlay
           WRITELOG(*,*) zkv(ngrid/2+1,ilev),zkh(ngrid/2+1,ilev)
@@ -429 +429 @@
     IF(lwrite) THEN
        WRITELOG(*,*) 
-       WRITELOG(*,*) ,'Diagnostique de la diffusion verticale'
-       WRITELOG(*,*) ,'h avant et apres diffusion verticale'
-       WRITELOG(*,*) ,ptsrf(ngrid/2+1),ztsrf2(ngrid/2+1)
+       WRITELOG(*,*) 'Diagnostique de la diffusion verticale'
+       WRITELOG(*,*) 'h avant et apres diffusion verticale'
+       WRITELOG(*,*) ptsrf(ngrid/2+1),ztsrf2(ngrid/2+1)
        DO  ilev=1,nlay
           WRITELOG(*,*) ph(ngrid/2+1,ilev),zh(ngrid/2+1,ilev)