Changeset 51 in lmdz_wrf


Ignore:
Timestamp:
Jul 23, 2014, 5:06:27 PM (10 years ago)
Author:
lfita
Message:

Changing largest values

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/LMDZ_WRFmeas/WRFV3/lmdz/diagphy_mod.F90

    r50 r51  
    9292      fname = 'diagphy'
    9393      errmsg = 'ERROR -- error -- ERROR -- error'
    94       largest = 10.e3
     94      largest = 10.e4
    9595
    9696      pas=pas+1
     
    164164        PRINT *,'    fs_bound = stops-stopl - (ssols+ssoll)+ssens+sfront'
    165165        PRINT *,'    airetot= ',airetot
    166         IF (airetot .NE. airetot .OR. ABS(airetot) > largest) THEN
     166        IF (airetot .NE. airetot .OR. ABS(airetot) > largest*10.e6) THEN
    167167          varname = 'airephy'
    168           CALL check_var(fname, varname, airephy, klon, largest*10.e9, .FALSE.)
     168          CALL check_var(fname, varname, airephy, klon, largest*10.e6, .FALSE.)
    169169        END IF
    170170        PRINT *,'    stops= ',stops
     
    216216      END IF
    217217
    218       IF (fq_bound .NE. fq_bound .OR. ABS(fq_bound) > 100000.) THEN
     218      IF (fq_bound .NE. fq_bound .OR. ABS(fq_bound) > largest) THEN
    219219        PRINT *,TRIM(errmsg)
    220220        PRINT *,'  ' // TRIM(fname) // ': Wrong fq_bound= ',fs_bound,' !!!'
     
    223223        PRINT *,'    evap_tot= ',evap_tot
    224224        PRINT *,'    airetot= ',airetot
    225         IF (airetot .NE. airetot .OR. ABS(airetot) > largest) THEN
     225        IF (airetot .NE. airetot .OR. ABS(airetot) > largest*10.e5) THEN
    226226          varname = 'airephy'
    227           CALL check_var(fname, varname, airephy, klon, largest, .FALSE.)
     227          CALL check_var(fname, varname, airephy, klon, largest*10.e5, .FALSE.)
    228228        END IF
    229229        IF (evap_tot .NE. evap_tot .OR. ABS(evap_tot) > largest) THEN
     
    449449      fname = 'diagetpq'
    450450      errmsg = 'ERROR -- error -- ERROR -- error'
    451       largest = 10.e3
     451      largest = 10.e4
    452452!c======================================================================
    453453!C
     
    603603          PRINT *,'      h_vcol_tot = h_dair_tot+h_qw_tot+h_ql_tot+h_qs_tot'
    604604          PRINT *,'    airetot= ',airetot
    605           IF (airetot .NE. airetot .OR. ABS(airetot) > largest) THEN
     605          IF (airetot .NE. airetot .OR. ABS(airetot) > largest*10.e5) THEN
    606606            varname = 'airephy'
    607             CALL check_var(fname, varname, airephy, klon, largest, .FALSE.)
     607            CALL check_var(fname, varname, airephy, klon, largest*10.e5, .FALSE.)
    608608          END IF
    609609          PRINT *,'    h_dair_tot= ',h_dair_tot
    610           IF (h_dair_tot .NE. h_dair_tot .OR. ABS(h_dair_tot) > largest) THEN
     610          IF (h_dair_tot .NE. h_dair_tot .OR. ABS(h_dair_tot) > largest*10.e6) THEN
    611611            varname = 'zh_dair_col'
    612             CALL check_var(fname, varname, zh_dair_col, klon, largest, .FALSE.)
     612            CALL check_var(fname, varname, zh_dair_col, klon, largest*10.e6, .FALSE.)
    613613          END IF
    614614          PRINT *,'    h_qw_tot= ',h_qw_tot
    615           IF (h_qw_tot .NE. h_qw_tot .OR. ABS(h_qw_tot) > largest) THEN
     615          IF (h_qw_tot .NE. h_qw_tot .OR. ABS(h_qw_tot) > largest*10.e4) THEN
    616616            varname = 'zh_qw_col'
    617             CALL check_var(fname, varname, zh_qw_col, klon, largest, .FALSE.)
     617            CALL check_var(fname, varname, zh_qw_col, klon, largest*10.e4, .FALSE.)
    618618          END IF
    619619          PRINT *,'    h_ql_tot= ',h_ql_tot
    620           IF (h_ql_tot .NE. h_ql_tot .OR. ABS(h_ql_tot) > largest) THEN
     620          IF (h_ql_tot .NE. h_ql_tot .OR. ABS(h_ql_tot) > largest*10.e2) THEN
    621621            varname = 'zh_ql_col'
    622             CALL check_var(fname, varname, zh_ql_col, klon, largest, .FALSE.)
     622            CALL check_var(fname, varname, zh_ql_col, klon, largest*10.e2, .FALSE.)
    623623          END IF
    624624          PRINT *,'    h_qs_tot= ',h_qs_tot
Note: See TracChangeset for help on using the changeset viewer.