Changeset 31 in lmdz_wrf

Timestamp:

Jul 23, 2014, 12:20:49 PM (10 years ago)

Author:

lfita

Message:

Adding checking for NaN

File:

• trunk/WRFV3/lmdz/diagphy_mod.F90 (modified) (2 diffs)

trunk/WRFV3/lmdz/diagphy_mod.F90

@@ -79 +79 @@
 !C
       integer pas
+! L. Fita, LMD July 2014
+      CHARACTER(LEN=50)                                  :: errmsg, fname
+
+      fname = 'diagphy'
+      errmsg = 'ERROR -- error -- ERROR -- error'
+
       save pas
       data pas/0/
@@ -146 +152 @@
        &    ,rain_fall_tot+snow_fall_tot
 
+! L. Fita, LMD July 2014. Checking for consistency
+      IF (fs_bound .NE. fs_bound .OR. ABS(fs_bound) > 100000.) THEN
+        PRINT *,TRIM(errmsg)
+        PRINT *,'  ' + TRIM(fname) + ': Wrong fs_bound= ',fs_bound,' !!!'
+        PRINT *,'    fs_bound: Heat flux at atm. boundaries'
+        PRINT *,'    fs_bound = stops-stopl - (ssols+ssoll)+ssens+sfront'
+        PRINT *,'    stops-stopl= ',stops-stopl
+        PRINT *,'    ssols= ',ssols
+        PRINT *,'    ssoll= ',ssoll
+        PRINT *,'    ssens= ',ssens
+        PRINT *,'    sfront= ',sfront
+        STOP
+      END IF
+
+      IF (fq_bound .NE. fq_bound .OR. ABS(fq_bound) > 100000.) THEN
+        PRINT *,TRIM(errmsg)
+        PRINT *,'  ' + TRIM(fname) + ': Wrong fs_bound= ',fs_bound,' !!!'
+        PRINT *,'    fq_bound: Watter flux at atm. boundaries'
+        PRINT *,'    fq_bound = evap_tot - rain_fall_tot -snow_fall_tot'
+        PRINT *,'    evap_tot= ',evap_tot
+        PRINT *,'    rain_fall_tot= ',rain_fall_tot
+        PRINT *,'    snow_fall_tot= ',snow_fall_tot
+        STOP
+      END IF
+
       return