Changeset 2642 in lmdz_wrf

Timestamp:

Jun 28, 2019, 9:40:39 PM (5 years ago)

Author:

lfita

Message:

Adding:

• `compute_front_R04d3': Subroutine to compute presence of a front following Rodrigues et al.(2004)
• `var_front_R04': Subroutine to compute presence of a front following Rodrigues et al.(2004)

Location:

trunk/tools

Files:

• module_ForDiagnostics.f90 (modified) (2 diffs)
• module_ForDiagnosticsVars.f90 (modified) (2 diffs)

trunk/tools/module_ForDiagnostics.f90

@@ -32 +32 @@
 ! compute_fog_RUC: Computation of fog and visibility following RUC method Smirnova, (2000)
 ! compute_fog_FRAML50: fog and visibility following Gultepe and Milbrandt, (2010)
+! compute_front_R04d3: Subroutine to compute presence of a front following Rodrigues et al.(2004)
 ! compute_massvertint1D: Subroutine to vertically integrate a 1D variable in eta vertical coordinates
 ! compute_psl_ecmwf: Compute sea level pressure using ECMWF method following Mats Hamrud and Philippe Courtier [Pa]
@@ -1549 +1550 @@
   END SUBROUTINE compute_tws_RK4
 
+  SUBROUTINE compute_front_R04d3(d1, d2, d3, tas, uas, vas, dsx, dsy, front)
+! Subroutine to compute presence of a front following Rodrigues et al.(2004), Rev. Bras.  Geofis. 22, 
+!    135-151
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: d1, d2, d3
+    REAL(r_k), DIMENSION(d1,d2,d3), INTENT(in)           :: tas, uas, vas, dsx, dsy
+    INTEGER, DIMENSION(d1,d2,d3), INTENT(out)            :: front
+ 
+! Local
+
+!!!!!!! Variables
+! tas: 2-m air temperature [K]
+! uas: 10-m eastward wind [ms-1]
+! vas: 10-m northward wind [ms-1]
+! dxs: grid spacing betweeen grid points along x-axis [m]
+! dys: grid spacing betweeen grid points along y-axis [m]
+
+    fname = 'compute_front_R04d3'
+
+    CALL var_front_R04(d1, d2, d3, tas, uas, vas, dsx, dsy, front)
+
+    RETURN
+
+  END SUBROUTINE compute_front_R04d3
+
 END MODULE module_ForDiagnostics

trunk/tools/module_ForDiagnosticsVars.f90

@@ -34 +34 @@
 ! var_fog_RUC: fog and visibility following RUC method Smirnova, (2000)
 ! var_fog_FRAML50: fog and visibility following Gultepe and Milbrandt, (2010)
+! var_front_R04: Subroutine to compute presence of a front following Rodrigues et al.(2004)
 ! var_potevap_orPM: potential evapotranspiration following Penman-Monteith formulation implemented in ORCHIDEE
 ! var_psl_ecmwf: sea level pressure using ECMWF method following Mats Hamrud and Philippe Courtier [Pa]
@@ -1976 +1977 @@
   END SUBROUTINE Svar_tws_S11
 
+  SUBROUTINE var_front_R04(dx, dy, dt, tas, uas, vas, ddx, ddy, front)
+! Subroutine to compute presence of a front following 
+!   Rodrigues et al.(2004), Rev. Bras.  Geofis. 22, 135-151, doi: 10.1590/S0102-261X2004000200004
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: dx, dy, dt
+    REAL(r_k), DIMENSION(dx,dy), INTENT(in)              :: ddx, ddy
+    REAL(r_k), DIMENSION(dx,dy,dt), INTENT(in)           :: tas, uas, vas
+    INTEGER, DIMENSION(dx,dy,dt), INTENT(out)            :: front
+
+! Local
+    INTEGER                                              :: i, j, l
+    REAL, DIMENSION(dx,dy,dt)                            :: dt1uas, dt1vas, dt1tas, dc1wss, d2tas
+
+!!!!!!! Variables
+! tas: 2-m temperature [K]
+! dd[x/y]: real distance between grid points along x and y axes [m]
+! uas: 10m eastward wind speed [ms-1]
+! vas: 10m northward wind speed [ms-1]
+! front: presence of a front in the grid point [0: no, 1: yes]
+
+    fname = 'var_front_R04'
+
+    dt1uas = zeroRK
+    dt1vas = zeroRK
+    dt1tas = zeroRK
+    dc1wss = zeroRK
+    d2tas = zeroRK
+
+    ! 1-time-step derivatives
+    DO l=1,dt-1
+      dt1tas(:,:,l) = ABS(tas(:,:,l+1) - tas(:,:,l))
+      dt1vas(:,:,l) = ABS(vas(:,:,l+1) - vas(:,:,l))
+    END DO
+
+    ! First order curl
+    DO l=1,dt-1
+      CALL curl2D_1o(dx,dy,ddx,ddy,uas(:,:,l),vas(:,:,l),dc1wss(:,:,l))
+    END DO
+
+    ! 2-time-step centered derivatives
+    DO l=2,dt-1
+      d2tas(:,:,l) = ABS(tas(:,:,l+1) - tas(:,:,l-1))
+    END DO
+
+    front = 0
+    DO l=1,dt
+      DO i=1,dx
+        DO j=1,dy
+          IF ( (dc1wss(i,j,l) /= zeroRK) .AND. (d2tas(i,j,l) > 0.5) ) front(i,j,l) = 1
+        END DO
+      END DO
+    END DO
+ 
+    RETURN
+
+  END SUBROUTINE var_front_R04
+
 END MODULE module_ForDiagnosticsVars