Changeset 1355 in lmdz_wrf

Timestamp:

Nov 25, 2016, 9:57:52 PM (8 years ago)

Author:

lfita

Message:

Fixing some issues of the `reprojection' modules
Generalizing the use of `module_generic'
Definition of errormsg' and warnmsg' in 'module_generic.F90'

Location:

trunk/tools

Files:

• interpolate.F90 (modified) (1 diff)
• module_ForInterpolate.F90 (modified) (18 diffs)
• module_generic.F90 (modified) (1 diff)

trunk/tools/interpolate.F90

@@ -225 +225 @@
   REAL(r_k), ALLOCATABLE, DIMENSION(:,:)                 :: fractionlon, fractionlat
   INTEGER                                                :: fracdx, fracdy
-  CHARACTER(LEN=50)                                      :: fname, errormsg
+  CHARACTER(LEN=50)                                      :: fname
 
 !!!!!!! Variables

trunk/tools/module_ForInterpolate.F90

@@ -10 +10 @@
 MODULE module_ForInterpolate
 
+  USE module_generic
+
   CONTAINS
 
@@ -15 +17 @@
   Nperx, Npery, ilonlat, mindiffLl)
 ! Function to search a given value from a coarser version of the data
-
-  USE module_generic
 
   IMPLICIT NONE
@@ -206 +206 @@
   iv, lonv, latv, per, Nperx, Npery, mindiff, ilonlat, mindiffLl)
 ! Function to search a given value from a coarser version of the data
-
-  USE module_generic
 
   IMPLICIT NONE
@@ -437 +435 @@
 ! Function to search a given value from a coarser version of the data
 
-  USE module_generic
-
   IMPLICIT NONE
 
@@ -489 +485 @@
 !   approche from percentages of the whole domain
 
-  USE module_generic
-
   IMPLICIT NONE
 
@@ -592 +586 @@
 ! Subroutine which finds the closest grid point within a projection throughout a first guest 
 !   and then whole domain approche from percentages of the whole domain 
-
-  USE module_generic
 
   IMPLICIT NONE
@@ -724 +716 @@
 ! Subroutine which provides the indices for a given interpolation of a projection
 
-  USE module_generic
-
   IMPLICIT NONE
 
@@ -738 +728 @@
   INTEGER                                                :: i,j,iv, ix, iy
   REAL(r_k)                                              :: mindiffLl, dist
-  REAL(r_k)                                              :: inclon, inclat, maxdiffprojlonlat
+  REAL(r_k)                                              :: inclon, inclat, maxdiffprojlonlat,        &
+    maxdiffinlonlat
   REAL(r_k), DIMENSION(idimx,idimy)                      :: difflonlat
+  REAL(r_k), DIMENSION(idimx,idimy)                      :: idifflon, idifflat
   REAL(r_k), DIMENSION(pdimx,pdimy)                      :: difflon, difflat
   REAL(r_k), DIMENSION(2)                                :: extremelon, extremelat, ipos
@@ -763 +755 @@
   extremelat = (/ MINVAL(projlat), MAXVAL(projlat) /)
  
+  ! Maximum distance between grid points in input projection
+  idifflon = 0.
+  idifflat = 0.
+  idifflon(1:idimx-1,:) = inlonv(2:idimx,:)-inlonv(1:idimx-1,:)
+  idifflat(:,1:idimy-1) = inlatv(:,2:idimy)-inlatv(:,1:idimy-1)
+  maxdiffinlonlat = MAXVAL(SQRT(idifflon**2. + idifflat**2.))
   ! Maximum distance between grid points in target projection
   difflon = 0.
@@ -770 +768 @@
   maxdiffprojlonlat = MAXVAL(SQRT(difflon**2. + difflat**2.))
 
+  IF (maxdiffinlonlat > maxdiffprojlonlat) THEN
+    PRINT *,TRIM(warnmsg)
+    PRINT *,'  ' //TRIM(fname)// '; input resolution: ', maxdiffinlonlat, ' is coarser than target:', &
+      maxdiffprojlonlat, ' !!'
+  END IF
+
   ! Using case outside loop to be more efficient
   SELECT CASE(TRIM(intkind))
@@ -781 +785 @@
           difflonlat = SQRT((projlon(i,j)-inlonv)**2. + (projlat(i,j)-inlatv)**2.)
           mindiffLl = MINVAL(difflonlat)
-          IF (mindiffLl > maxdiffprojlonlat) THEN
+          IF ( (mindiffLl > maxdiffprojlonlat) .AND. (mindiffLl > maxdiffinlonlat)) THEN
             outLlw(3,:,i,j) = 0.
             outLlw(3,:,i,j) = -1.
@@ -852 +856 @@
           ! We do not want that values larger that the maximum distance between target grid points
 !          PRINT *,i,j,':',mindiffLl,'maxdiffLl:',maxdiffprojlonlat
-          IF (mindiffLl .gt. maxdiffprojlonlat) THEN
+          IF ((mindiffLl .gt. maxdiffprojlonlat) .AND. (mindiffLl > maxdiffinlonlat)) THEN
 !            PRINT *,'  ' // TRIM(fname) // ': reprojected minimum distance to nearest grid point:',   &
 !              mindiffLl, ' larger than the maximum distance between grid points in target projection!!'
@@ -870 +874 @@
   idimy, pdimx, pdimy)
 ! Subroutine to interpolate a 2D variable
-
-  USE module_generic
 
   IMPLICIT NONE
@@ -949 +951 @@
   idimy, pdimx, pdimy, d3)
 ! Subroutine to interpolate a 3D variable
-
-  USE module_generic
 
   IMPLICIT NONE
@@ -1032 +1032 @@
   idimy, pdimx, pdimy, d3, d4)
 ! Subroutine to interpolate a 4D variable
-
-  USE module_generic
 
   IMPLICIT NONE
@@ -1119 +1117 @@
   idimy, pdimx, pdimy, d3, d4, d5)
 ! Subroutine to interpolate a 5D variable
-
-  USE module_generic
 
   IMPLICIT NONE
@@ -1211 +1207 @@
 ! Subroutine which finds the closest grid point within a projection
 
-  USE module_generic
-
   IMPLICIT NONE
 
@@ -1309 +1303 @@
   Ninpts)
 ! Subroutine which finds the closest grid point within a projection with 1D longitudes and latitudes
-
-  USE module_generic
 
   IMPLICIT NONE

trunk/tools/module_generic.F90

@@ -8 +8 @@
   ! Fill value at 64 bits
   REAL(r_k)                                              :: fillval64 = 1.e20
+  CHARACTER(len=50), PARAMETER                           :: errormsg = 'ERROR -- error -- ERROR -- error'
+  CHARACTER(len=50), PARAMETER                           :: warnmsg = 'WARNING -- warning -- WARNING -- warning'
 
   CONTAINS