Changeset 1052 for codes/icosagcm/devel/src/diagnostics/compute_diagnostics.f90

Timestamp:

09/23/20 16:00:58 (4 years ago)

Author:

dubos

Message:

devel : diagnose divergence and vorticity

File:

• codes/icosagcm/devel/src/diagnostics/compute_diagnostics.f90 (modified) (2 diffs)

codes/icosagcm/devel/src/diagnostics/compute_diagnostics.f90

@@ -27 +27 @@
     END SUBROUTINE comp_temperature
 
+    SUBROUTINE comp_vorticity(ue, vort)
+      IMPORT
+      REAL(rstd),INTENT(IN)  :: ue(:,:)
+      REAL(rstd),INTENT(OUT) :: vort(:,:)
+    END SUBROUTINE comp_vorticity
+
+    SUBROUTINE comp_divergence(ue, div)
+      IMPORT
+      REAL(rstd),INTENT(IN)  :: ue(:,:)
+      REAL(rstd),INTENT(OUT) :: div(:,:)
+    END SUBROUTINE comp_divergence
+
     SUBROUTINE comp_hydro_press(rhodz, theta_rhodz, ps, p)
       IMPORT
@@ -48 +60 @@
   PROCEDURE(comp_pression),    POINTER, SAVE :: compute_pression => NULL(), compute_pression_mid => NULL()
   PROCEDURE(comp_temperature), POINTER, SAVE :: compute_temperature => NULL()
+  PROCEDURE(comp_vorticity),   POINTER, SAVE :: compute_vorticity => NULL()
+  PROCEDURE(comp_divergence),  POINTER, SAVE :: compute_divergence => NULL()
   PROCEDURE(comp_hydro_press), POINTER, SAVE :: compute_hydrostatic_pressure => NULL()
   PROCEDURE(comp_vert_interp), POINTER, SAVE :: compute_vertical_interp => NULL()
 
-  PUBLIC :: compute_rhodz, compute_pression, compute_pression_mid, compute_temperature, compute_hydrostatic_pressure, &
-       compute_vertical_interp
+  PUBLIC :: compute_rhodz, compute_pression, compute_pression_mid, compute_temperature, &
+       compute_vorticity, compute_divergence, compute_hydrostatic_pressure, compute_vertical_interp
 
 END MODULE compute_diagnostics_mod