Changeset 33 in lmdz_wrf

Timestamp:

Jul 23, 2014, 12:44:01 PM (10 years ago)

Author:

lfita

Message:

Adding subroutine 'check_var' to check consistency of a given variable

File:

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

branches/LMDZ_WRFmeas/WRFV3/lmdz/diagphy_mod.F90

@@ -81 +81 @@
 ! L. Fita, LMD July 2014
       CHARACTER(LEN=50)                                  :: errmsg, fname
+      LOGICAL                                            :: found
 
       fname = 'diagphy'
@@ -185 +186 @@
       end SUBROUTINE diagphy
 
+!L. Fita, LMD 2004. Check variable
+SUBROUTINE check_var(funcn, varn, var, sizev, bigvalue, stoprun)
+!  Subroutine to check the consistency of a variable
+
+  IMPLICIT NONE
+
+  INTEGER, INTENT(IN)                                    :: sizev
+  CHARACTER(LEN=50), INTENT(IN)                          :: funcn, varn
+  REAL, DIMENSION(sizev), INTENT(IN)                     :: var
+  REAL, INTENT(IN)                                       :: bigvalue
+  LOGICAL, INTENT(IN)                                    :: stoprun
+
+! Local
+  INTEGER                                                :: i, wrongi
+  CHARACTER(LEN=50)                                      :: errmsg
+  LOGICAL                                                :: found
+  REAL, DIMENSION(sizev)                                 :: wrongvalues
+  INTEGER, DIMENSION(sizev)                              :: wronggridpt
+
+!!!!!!! Variables
+! funcn: at which functino of part of the program variable is checked
+! varn: name of the variable
+! var: variable to check
+! sizev: sixe of the variable
+! bigvalue: biggest attenaible value for the variable
+! stoprun: Should the run stop if it founds a problem?
+
+  errmsg = 'ERROR -- error -- ERROR -- error'
+
+  found = .FALSE.
+  wrongi = 0
+  DO i=1,sizev
+    IF (var(i) /= var(i) .OR. var(i) > 100000. ) THEN
+      IF (wrongi == 0) found = .TRUE.
+      wrongi = wrongi + 1
+      wrongvalues(wrongi) = var(i)
+      wronggridpt(wrongi) = i
+    END IF
+  END DO
+
+  IF (found) THEN
+    PRINT *,TRIM(errmsg)
+    PRINT *,"  at '" // TRIM(funcn) // "' variable '" //TRIM(varn)//                 &
+      "' is wrong at i value___'
+    DO i=1,wrongi
+      PRINT *,wronggridpt(i), wrongvalues(i)
+    END DO
+    IF (stoprun) THEN
+      STOP
+    END IF
+  END IF
+
+  RETURN
+
+END SUBROUTINE check_var
+
 !C======================================================================
       SUBROUTINE diagetpq(airephy,tit,iprt,idiag,idiag2,dtime                        &