Changeset 4225 for dynamico_lmdz/simple_physics/phyparam

Timestamp:

Jan 10, 2020, 3:07:23 PM (5 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

Files:

• Makefile (modified) (3 diffs)
• physics/logging.F90 (modified) (2 diffs)
• physics/turbulence.F90 (modified) (3 diffs)

dynamico_lmdz/simple_physics/phyparam/Makefile

@@ -1 @@
-# Example (bash) :
-# clear ; make clean ; F90=ifort FFLAGS=-fast time make -j
+# Examples (bash) :
+# clear ; make clean ; 
+# F90=ifort FFLAGS=-fast time make -j
+# F90=gfortran F90FLAGS=-ffree-line-length-none gmake 
+# known to compile with gfortran version 7, ifort version 17
+
+# if XCodeML F_Front is installed, set FFRONT=F_Front to parse code into XML
+FFRONT ?= true
 
 VPATH = obj
@@ -9 +15 @@
 
 obj :
-        @mkdir obj include lib
+        @mkdir obj include lib xml
 
 clean :
-        rm -f obj/* include/* lib/*
+        rm -f obj/* include/* lib/* xml/*
 
 %.so : $(OBJECTS)
@@ -28 +34 @@
 obj/phyparam_mod.o     : MAIN solar.o radiative_sw.o radiative_lw.o comgeomfi.o writefield_mod.o
 
+SRC = $(basename $(notdir $@))
 obj/%.o: physics/%.F90
-        $(F90) -I include -fPIC -c $(F90FLAGS) $< -o $@
-        @mv $(basename $(notdir $@)).mod include
+        $(F90) -I include -fPIC $(F90FLAGS) -c $<
+        @mv $(SRC).mod include
+        @mv $(SRC).o obj
+        @cpp -DXCODEML $< > xml/$(SRC).f90
+        @$(FFRONT) -M xml xml/$(SRC).f90 > xml/$(SRC).xml

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)