Changeset 138 in lmdz_wrf

Timestamp:

Jul 29, 2014, 11:36:55 AM (10 years ago)

Author:

lfita

Message:

Removing checkings of NaNs?

File:

• branches/LMDZ_WRFmeas/WRFV3/lmdz/physiq.F90 (modified) (5 diffs)

branches/LMDZ_WRFmeas/WRFV3/lmdz/physiq.F90

@@ -1299 +1299 @@
       integer iostat
 
-! Lluis
-      INTEGER                                            :: llp
-      CHARACTER(LEN=50)                                  :: lvarname, lfname
-      REAL                                               :: largest
-
-      llp = 734
-      lfname = 'physiq'
-      largest = 10.e5
-
 !c======================================================================
 ! Gestion calendrier : mise a jour du module phys_cal_mod
 !
       CALL phys_cal_update(jD_cur,jH_cur)
-      PRINT *,'  Lluis physiq: jD_cur: ',jD_cur, ' jH_cur: ',jH_cur, &
-        ' days_elapsed: ',days_elapsed
 
 !c======================================================================
@@ -2769 +2758 @@
 !c      detr_therm(:,:)=0.
 !c
-      lfname='physiq before thermals'
-      lvarname = 't_seri'
-      CALL check_var3D(lfname, lvarname, t_seri, klon, klev, largest, .FALSE.)
-      lvarname = 'q_seri'
-      CALL check_var3D(lfname, lvarname, q_seri, klon, klev, largest, .FALSE.)
-      lvarname = 'd_t_ajs'
-      CALL check_var3D(lfname, lvarname, d_t_ajs, klon, klev, largest, .FALSE.)
-      lvarname = 'd_q_ajs'
-      CALL check_var3D(lfname, lvarname, d_q_ajs, klon, klev, largest, .FALSE.)
 
       IF(prt_level>9)WRITE(lunout,*)                                                 &
@@ -2828 +2808 @@
        &      ,zqla,ztva )
 
-      lfname='physiq just after thermals'
-      lvarname = 't_seri'
-      CALL check_var3D(lfname, lvarname, t_seri, klon, klev, largest, .FALSE.)
-      lvarname = 'q_seri'
-      CALL check_var3D(lfname, lvarname, q_seri, klon, klev, largest, .FALSE.)
-      lvarname = 'd_t_ajs'
-      CALL check_var3D(lfname, lvarname, d_t_ajs, klon, klev, largest, .FALSE.)
-      lvarname = 'd_q_ajs'
-      CALL check_var3D(lfname, lvarname, d_q_ajs, klon, klev, largest, .FALSE.)
-
 !ccc nrlmd le 10/04/2012
 !c-----------Stochastic triggering-----------
@@ -2976 +2946 @@
 
       endif
-      lfname='after thermals'
-      lvarname = 't_seri'
-      CALL check_var3D(lfname, lvarname, t_seri, klon, klev, largest, .FALSE.)
-      lvarname = 'q_seri'
-      CALL check_var3D(lfname, lvarname, q_seri, klon, klev, largest, .FALSE.)
-      lvarname = 'd_t_ajsb'
-      CALL check_var3D(lfname, lvarname, d_t_ajsb, klon, klev, largest, .FALSE.)
-      lvarname = 'd_q_ajsb'
-      CALL check_var3D(lfname, lvarname, d_q_ajsb, klon, klev, largest, .FALSE.)
 
 !c
@@ -4384 +4345 @@
       RETURN
       END SUBROUTINE gr_fi_ecrit
-
-SUBROUTINE check_var(funcn, varn, var, sizev, bigvalue, stoprun)
-!  Subroutine to check the consistency of a variable
-!    * NaN value: by definition is variable /= variable
-!    * bigvalue: threshold for the variable
-
-  IMPLICIT NONE
-
-#include "dimensions.h"
-
-  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, xpt, ypt
-  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: size 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. ABS(var(i)) > bigvalue ) 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 in Nvalues= ",wrongi,' at i (x, y) value___', bigvalue
-    DO i=1,wrongi
-       ypt = INT(wronggridpt(i)/wiim) + 1
-       xpt = wronggridpt(i) - (ypt-1)*wiim
-      PRINT *,wronggridpt(i), '(',xpt,', ',ypt,')', wrongvalues(i)
-    END DO
-    IF (stoprun) THEN
-      STOP
-    END IF
-  END IF
-
-  RETURN
-
-END SUBROUTINE check_var
-
-SUBROUTINE check_var3D(funcn, varn, var, sizev, zsize, bigvalue, stoprun)
-!  Subroutine to check the consistency of a 3D LMDSZ - variable (klon, klev) !
-!    * NaN value: by definition is variable /= variable
-!    * bigvalue: threshold for the variable
-
-  IMPLICIT NONE
-
-#include "dimensions.h"
-
-  INTEGER, INTENT(IN)                                    :: sizev, zsize
-  CHARACTER(LEN=50), INTENT(IN)                          :: funcn, varn
-  REAL, DIMENSION(sizev,zsize), INTENT(IN)               :: var
-  REAL, INTENT(IN)                                       :: bigvalue
-  LOGICAL, INTENT(IN)                                    :: stoprun
-
-! Local
-  INTEGER                                                :: i, k, wrongi, xpt, ypt
-  CHARACTER(LEN=50)                                      :: errmsg
-  LOGICAL                                                :: found
-  REAL, DIMENSION(sizev*zsize)                           :: wrongvalues
-  INTEGER, DIMENSION(sizev*zsize,2)                      :: wronggridpt
-
-!!!!!!! Variables
-! funcn: at which functino of part of the program variable is checked
-! varn: name of the variable
-! var: variable to check
-! sizev: size of the variable
-! zsize: vertical size 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
-    DO k=1,zsize
-      IF (var(i,k) /= var(i,k) .OR. ABS(var(i,k)) > bigvalue ) THEN
-        IF (wrongi == 0) found = .TRUE.
-        wrongi = wrongi + 1
-        wrongvalues(wrongi) = var(i,k)
-        wronggridpt(wrongi,1) = i
-        wronggridpt(wrongi,2) = k
-      END IF
-    END DO
-  END DO
-
-  IF (found) THEN
-    PRINT *,TRIM(errmsg)
-    PRINT *,"  at '" // TRIM(funcn) // "' variable '" //TRIM(varn)//                 &
-      "' is wrong in Nvalues= ",wrongi,' at i (x,y) k value___'
-    DO i=1,wrongi
-       ypt = INT(wronggridpt(i,1)/wiim) + 1
-       xpt = wronggridpt(i,1) - (ypt-1)*wiim
-      PRINT *,wronggridpt(i,1), '(',xpt,', ',ypt,')', wronggridpt(i,2), wrongvalues(i)
-    END DO
-    IF (stoprun) THEN
-      STOP
-    END IF
-  END IF
-
-  RETURN
-
-END SUBROUTINE check_var3D
-