Changeset 4236 for dynamico_lmdz/simple_physics/phyparam/physics/phyparam_mod.F90

Timestamp:

Jan 22, 2020, 9:46:46 PM (5 years ago)

Author:

dubos

Message:

simple_physics : some Python bindings

File:

• dynamico_lmdz/simple_physics/phyparam/physics/phyparam_mod.F90 (modified) (5 diffs)

dynamico_lmdz/simple_physics/phyparam/physics/phyparam_mod.F90

@@ -34 +34 @@
        &            pplev,pplay,pphi,           &
        &            pu,pv,pt,                   &
-       &            pdu,pdv,pdt,pdpsrf)
+       &            pdu,pdv,pdt,pdpsrf) BIND(C, name='phyparam_phyparam')
     USE phys_const, ONLY : g, rcp, r, unjours
     USE surface,    ONLY : soil
@@ -45 +45 @@
     !   20 parameters GCM for planetary atmospheres.
     !   It includes:
-    !   raditive transfer (long and shortwave) for CO2 and dust.
+    !   radiative transfer (long and shortwave) for CO2 and dust.
     !   vertical turbulent mixing
     !   convective adjsutment
+    !   heat diffusion in the soil
     !
     !   author: Frederic Hourdin 15 / 10 /93
     !=======================================================================
 
-    INTEGER, INTENT(IN) ::      &
+    INTEGER, INTENT(IN), VALUE :: &
          ngrid,                 & ! Size of the horizontal grid.
          nlayer                   ! Number of vertical layers.
-    LOGICAL, INTENT(IN) ::      &
+    LOGICAL, INTENT(IN), VALUE  :: &
          firstcall,             & ! True at the first call
          lastcall                 ! True at the last call
-    REAL, INTENT(IN)    ::      &
+    REAL, INTENT(IN), VALUE     ::      &
          rjourvrai,             & ! Number of days counted from the North. Spring equinox
          gmtime,                & ! time of the day in seconds
-         ptimestep,             & ! timestep (s)
+         ptimestep                ! timestep (s)
+    REAL, INTENT(IN) :: &
          pplev(ngrid,nlayer+1), & ! Pressure at interfaces between layers (pa)
          pplay(ngrid,nlayer),   & ! Pressure at the middle of the layers (Pa)
@@ -431 +433 @@
   END SUBROUTINE radiative_tendencies
 
-  SUBROUTINE alloc(ngrid, nlayer)
+  SUBROUTINE alloc(ngrid, nlayer) BIND(C, name='phyparam_alloc')
+    !$cython header void phyparam_alloc(int, int);
+    !$cython wrapper def alloc(ngrid, nlayer) : phy.phyparam_alloc(ngrid, nlayer)
     USE astronomy, ONLY : iniorbit
     USE surface, ONLY : zc,zd
-    INTEGER, INTENT(IN) :: ngrid, nlayer
+    INTEGER, INTENT(IN), VALUE :: ngrid, nlayer
     LOGICAL, PARAMETER :: firstcall=.TRUE.
     ! allocate arrays for internal state
@@ -448 +452 @@
   END SUBROUTINE alloc
 
-  SUBROUTINE precompute
+  SUBROUTINE precompute() BIND(C, name='phyparam_precompute')
+    !$cython header void phyparam_precompute();
+    !$cython wrapper def precompute() : phy.phyparam_precompute()
     USE surface
     ! precompute time-independent arrays
@@ -458 +464 @@
   END SUBROUTINE precompute
 
-  SUBROUTINE coldstart(ngrid, ptimestep)
+  SUBROUTINE coldstart(ngrid, ptimestep) BIND(C, name='phyparam_coldstart')
+    !$cython header void phyparam_coldstart(int, double);
+    !$cython wrapper def coldstart (ngrid, timestep): phy.phyparam_coldstart(ngrid, timestep)
     ! create internal state to start a run without a restart file
     USE surface, ONLY : soil
-    INTEGER, INTENT(IN) :: ngrid
-    REAL, INTENT(iN)    :: ptimestep
+    INTEGER, INTENT(IN), VALUE :: ngrid
+    REAL, INTENT(IN),    VALUE :: ptimestep
     tsurf(:)   = tsoil_init
     tsoil(:,:) = tsoil_init