Changeset 2641 in lmdz_wrf

Timestamp:

Jun 28, 2019, 9:38:35 PM (5 years ago)

Author:

lfita

Message:

Adding:

• `curl2D_1o': Subroutine to compute the first order curl of a 2D vectorial field
• `divergence2D_1o': Subroutine to compute the first order divergence of a 2D vectorial field
• `gradient2D_1o': Subroutine to compute the first order gradient of a 2D field

File:

• trunk/tools/module_scientific.f90 (modified) (2 diffs)

trunk/tools/module_scientific.f90

@@ -18 +18 @@
 ! coincident_polygon: Subroutine to provide the intersection polygon between two polygons
 ! crossingpoints_polys: Subroutine to provide the crossing points between two polygons
+! curl2D_1o: Subroutine to compute the first order curl of a 2D vectorial field
 ! distanceRK: Function to provide the distance between two points
+! divergence2D_1o: Subroutine to compute the first order divergence of a 2D vectorial field
 ! FindMinimumR_K*: Function returns the location of the minimum in the section between Start and End.
 ! fill3DI_2Dvec: Subroutine to fill a 3D integer matrix from a series of indices from a 2D matrix
 ! fill3DR_2Dvec: Subroutine to fill a 3D float matrix from a series of indices from a 2D matrix
+! gradient2D_1o: Subroutine to compute the first order gradient of a 2D field
 ! grid_within_polygon: Subroutine to determine which grid cells from a matrix lay inside a polygon
 ! grid_spacepercen: Subroutine to compute the space-percentages of a series of grid cells (B) which lay inside another 
@@ -7602 +7605 @@
   END SUBROUTINE crossingpoints_polys
 
+  SUBROUTINE gradient2D_1o(dx, dy, var, dsx, dsy, grad)
+  ! Subroutine to compute the first order gradient of a 2D field
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: dx, dy
+    REAL(r_k), DIMENSION(dx,dy), INTENT(in)              :: var, dsx, dsy
+    REAL(r_k), DIMENSION(dx,dy,2), INTENT(out)           :: grad
+
+! Local 
+    INTEGER                                              :: i, j
+
+!!!!!!! Variables
+! dx, dy: shape of the 2D field
+! var: variable
+! dsx, dsy: matrices of distances betweeen grid points along each axis
+! grad: gradient
+
+    fname = 'gradient2D_1o'
+
+    grad = zeroRK
+
+    DO i=2, dx-1
+      DO j=2, dy-1
+        grad(i,j,:) = (/ (var(i+1,j)-var(i,j))/dsx(i,j), (var(i,j+1)-var(i,j))/dsy(i,j) /)
+      END DO
+    END DO
+
+  END SUBROUTINE gradient2D_1o
+
+  SUBROUTINE divergence2D_1o(dx, dy, dsx, dsy, xvar, yvar, diver)
+  ! Subroutine to compute the first order divergence of a 2D vectorial field
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: dx, dy
+    REAL(r_k), DIMENSION(dx,dy), INTENT(in)              :: xvar, yvar, dsx, dsy
+    REAL(r_k), DIMENSION(dx,dy), INTENT(out)             :: diver
+
+! Local 
+    INTEGER                                              :: i, j
+
+!!!!!!! Variables
+! dx, dy: shape of the 2D field
+! xvar, yvar: x-y components of vectorial variable
+! dsx, dsy: matrices of distances betweeen grid points along each axis
+! diver: divergence
+
+    fname = 'divergence2D_1o'
+
+    diver = zeroRK
+
+    DO i=2, dx-1
+      DO j=2, dy-1
+        diver(i,j) = (xvar(i+1,j)-xvar(i,j))/dsx(i,j) + (yvar(i,j+1)-yvar(i,j))/dsy(i,j)
+      END DO
+    END DO
+
+  END SUBROUTINE divergence2D_1o
+
+  SUBROUTINE curl2D_1o(dx, dy, dsx, dsy, xvar, yvar, curl)
+  ! Subroutine to compute the first order curl of a 2D vectorial field
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: dx, dy
+    REAL(r_k), DIMENSION(dx,dy), INTENT(in)              :: xvar, yvar, dsx, dsy
+    REAL(r_k), DIMENSION(dx,dy), INTENT(out)             :: curl
+
+! Local 
+    INTEGER                                              :: i, j
+
+!!!!!!! Variables
+! dx, dy: shape of the 2D field
+! xvar, yvar: x-y components of vectorial variable
+! dsx, dsy: matrices of distances betweeen grid points along each axis
+! curl: curl
+
+    fname = 'curl2D_1o'
+
+    curl = zeroRK
+
+    DO i=1, dx-1
+      DO j=1, dy-1
+        curl(i,j) = (yvar(i+1,j)-yvar(i,j))/dsx(i,j) - (xvar(i,j+1)-xvar(i,j))/dsy(i,j)
+      END DO
+    END DO
+
+  END SUBROUTINE curl2D_1o
+
 END MODULE module_scientific