Changeset 2274 in lmdz_wrf for trunk/tools/module_ForDiagnostics.f90

Timestamp:

Dec 28, 2018, 3:35:45 PM (6 years ago)

Author:

lfita

Message:

Adding:

• `WRFbnds': grid cell bounds from XLONG_U, XLAT_U, XLONG_V, XLAT_V as intersection of their related parallels and meridians

File:

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

trunk/tools/module_ForDiagnostics.f90

@@ -16 +16 @@
 ! compute_cape_afwa4D: Subroutine to use WRF phys/module_diag_afwa.F `buyoancy' subroutine to compute 
 !   CAPE, CIN, ZLFC, PLFC, LI
+! compute_cellbnds: Subroutine to compute cellboundaries using wind-staggered lon, lats as 
+!   intersection of their related parallels and meridians
 ! compute_cllmh4D3: Computation of low, medium and high cloudiness from a 4D CLDFRA and pressure being 
 !   3rd dimension the z-dim 
@@ -1051 +1053 @@
   END SUBROUTINE compute_range_faces
 
+  SUBROUTINE compute_cellbnds(dx, dy, sdx, sdy, ulon, ulat, vlon, vlat, xbnds, ybnds)
+! Subroutine to compute cellboundaries using wind-staggered lon, lats as intersection of their related
+!   parallels and meridians
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: dx, dy, sdx, sdy
+    REAL(r_k), DIMENSION(sdx, dy), INTENT(in)            :: ulon, ulat
+    REAL(r_k), DIMENSION(dx, sdy), INTENT(in)            :: vlon, vlat
+    REAL(r_k), DIMENSION(dx, dy, 4), INTENT(out)         :: xbnds, ybnds
+
+! Local
+    INTEGER                                              :: i,j,iv
+    INTEGER                                              :: ix,ex,iy,ey
+    CHARACTER(len=2), DIMENSION(4)                       :: Svertex
+    INTEGER, DIMENSION(4,2,2,2)                          :: indices
+    REAL(r_k), DIMENSION(2)                              :: ptintsct
+    REAL(r_k), DIMENSION(2,2)                            :: merid, paral
+    LOGICAL                                              :: intsct
+
+!!!!!!! Variables
+! dx, dy: un-staggered dimensions
+! sdx, sdy: staggered dimensions
+! ulon, ulat: x-wind staggered longitudes and latitudes
+! vlon, vlat: y-wind staggered longitudes and latitudes
+! xbnds, ybnds: x and y cell boundaries
+
+    fname = 'compute_cellbnds'
+
+    ! Indices to use indices[SW/NW/NE/SE, m/p, x/y, i/e]
+    Svertex = (/ 'SW', 'NW', 'NE', 'SE' /)
+
+    ! SW
+    indices(1,1,1,1) = 0
+    indices(1,1,1,2) = 0
+    indices(1,1,2,1) = -1
+    indices(1,1,2,2) = 0
+    indices(1,2,1,1) = -1
+    indices(1,2,1,2) = 0
+    indices(1,2,2,1) = -1
+    indices(1,2,2,2) = -1
+    ! NW
+    indices(2,1,1,1) = 0
+    indices(2,1,1,2) = 0
+    indices(2,1,2,1) = 0
+    indices(2,1,2,2) = 1
+    indices(2,2,1,1) = -1
+    indices(2,2,1,2) = 0
+    indices(2,2,2,1) = 1
+    indices(2,2,2,2) = 1
+    ! NE
+    indices(3,1,1,1) = 1
+    indices(3,1,1,2) = 1
+    indices(3,1,2,1) = 0
+    indices(3,1,2,2) = 1
+    indices(3,2,1,1) = 0
+    indices(3,2,1,2) = 1
+    indices(3,2,2,1) = 1
+    indices(3,2,2,2) = 1
+    ! SE
+    indices(4,1,1,1) = 1
+    indices(4,1,1,2) = 1
+    indices(4,1,2,1) = -1
+    indices(4,1,2,2) = 0
+    indices(4,2,1,1) = 0
+    indices(4,2,1,2) = 1
+    indices(4,2,2,1) = -1
+    indices(4,2,2,2) = -1
+
+    DO i=1,dx
+      DO j=1,dy
+        DO iv=1,4
+
+          ix = MAX(i+indices(iv,1,1,1),1)
+          !ex = MIN(i+indices(iv,1,1,2),dx)
+          ex = i+indices(iv,1,1,2)
+          iy = MAX(j+indices(iv,1,2,1),1)
+          ey = MIN(j+indices(iv,1,2,2),dy)
+ 
+          merid(1,1) = ulon(ix,iy)
+          merid(1,2) = ulat(ix,iy)
+          merid(2,1) = ulon(ex,ey)
+          merid(2,2) = ulat(ex,ey)
+
+          ix = MAX(i+indices(iv,2,1,1),1)
+          ex = MIN(i+indices(iv,2,1,2),dx)
+          iy = MAX(j+indices(iv,2,2,1),1)
+          !ey = MIN(i+indices(iv,2,2,2),dy)
+          ey = j+indices(iv,2,2,2)
+          paral(1,1) = vlon(ix,iy)
+          paral(1,2) = vlat(ix,iy)
+          paral(2,1) = vlon(ex,ey)
+          paral(2,2) = vlat(ex,ey)
+
+          CALL intersection_2Dlines(merid, paral, intsct, ptintsct)
+          IF (.NOT.intsct) THEN
+            msg = 'not interection found for ' // Svertex(iv) // ' vertex'
+            CALL ErrMsg(msg, fname, -1)
+          END IF
+          xbnds(i,j,iv) = ptintsct(1)
+          ybnds(i,j,iv) = ptintsct(2)
+        END DO
+      END DO
+    END DO
+
+  END SUBROUTINE
+
   SUBROUTINE compute_Koeppen_Geiger_climates(dx, dy, dt, pr, tas, climatesI, climatesS, climlegend)
   ! Subroutine to compute the Koeppen-Geiger climates after: