Changeset 1202

Timestamp:

Jul 7, 2009, 4:48:23 PM (15 years ago)

Author:

Laurent Fairhead

Message:

Probleme dans le test pour savoir si les latitudes sont inversees ou non
FH/LF

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol.F90 (modified) (2 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/readaerosol.F90

@@ -243 +243 @@
        CALL check_err( nf90_get_var(ncid, varid, lon_src(:)) )
        
-       IF (minval(ABS(lon_src - io_lon)) > 0.001) THEN
+       IF (maxval(ABS(lon_src - io_lon)) > 0.001) THEN
           WRITE(lunout,*) 'Problem in longitudes read from file : ',TRIM(fname)
           WRITE(lunout,*) 'longitudes in file ', TRIM(fname),' : ', lon_src
@@ -260 +260 @@
        END DO
 
-       IF (minval(ABS(lat_src - io_lat)) < 0.001) THEN
+       IF (maxval(ABS(lat_src - io_lat)) < 0.001) THEN
           ! Latitudes are the same
           invert_lat=.FALSE.
-       ELSE IF (minval(ABS(lat_src_inv - io_lat)) < 0.001) THEN
+       ELSE IF (maxval(ABS(lat_src_inv - io_lat)) < 0.001) THEN
           ! Inverted source latitudes correspond to model latitudes
           WRITE(lunout,*) 'latitudes will be inverted for file : ',TRIM(fname)