source: LMDZ6/branches/cirrus/libf/dyn3d_common/iso_verif_dyn.F @ 5409

Last change on this file since 5409 was 5203, checked in by Laurent Fairhead, 3 months ago

Merge with trunk revision 5202 before reintegration to trunk

File size: 3.4 KB
Line 
1        function iso_verif_noNaN_nostop(x,err_msg)
2        implicit none
3        ! si x est NaN, on affiche message
4        ! d'erreur et return 1 si erreur
5
6        ! input:
7        real x
8        character*(*) err_msg ! message d''erreur à afficher
9
10        ! output
11        real borne
12        parameter (borne=1e19)
13        integer iso_verif_noNaN_nostop
14
15        if ((x.gt.-borne).and.(x.lt.borne)) then
16                iso_verif_noNAN_nostop=0
17        else
18            write(*,*) 'erreur detectee par iso_verif_nonNaN:'
19            write(*,*) err_msg
20            write(*,*) 'x=',x
21            iso_verif_noNaN_nostop=1
22        endif     
23
24        return
25        end
26
27        function iso_verif_egalite_nostop
28     :           (a,b,err_msg)
29        implicit none
30        ! compare a et b. Si pas egal, on affiche message
31        ! d'erreur et stoppe
32        ! pour egalite, on verifie erreur absolue et arreur relative
33
34        ! input:
35        real a, b
36        character*(*) err_msg ! message d''erreur à afficher
37
38        ! locals
39        real errmax ! erreur maximale en absolu.
40        real errmaxrel ! erreur maximale en relatif autorisée
41        parameter (errmax=1e-8)
42        parameter (errmaxrel=1e-3)
43
44        ! output
45        integer iso_verif_egalite_nostop
46
47        iso_verif_egalite_nostop=0
48
49        if (abs(a-b).gt.errmax) then
50          if (abs((a-b)/max(max(abs(b),abs(a)),1e-18))
51     :            .gt.errmaxrel) then
52            write(*,*) 'erreur detectee par iso_verif_egalite:'
53            write(*,*) err_msg
54            write(*,*) 'a=',a
55            write(*,*) 'b=',b
56            iso_verif_egalite_nostop=1
57          endif
58        endif     
59       
60        return
61        end       
62
63
64        function iso_verif_aberrant_nostop
65     :           (x,iso,q,err_msg)
66#ifdef CPP_IOIPSL
67        USE IOIPSL, ONLY: getin
68#else
69        USE ioipsl_getincom, ONLY: getin
70#endif
71        USE infotrac, ONLY: isoName, getKey
72        implicit none
73       
74        ! input:
75        real x,q
76        integer iso ! 2=HDO, 1=O18
77        character*(*) err_msg ! message d''erreur à afficher
78
79        ! locals
80        real qmin,deltaD
81        real deltaDmax,deltaDmin,tnat
82        parameter (qmin=1e-11)
83        parameter (deltaDmax=200.0,deltaDmin=-999.9)
84        LOGICAL, SAVE :: ltnat1
85        LOGICAL, SAVE :: lFirst=.TRUE.
86
87        ! output
88        integer iso_verif_aberrant_nostop
89
90        IF(lFirst) THEN
91           ltnat1 = .TRUE.; CALL getin('tnateq1', ltnat1)
92           lFirst = .FALSE.
93        END IF
94        iso_verif_aberrant_nostop=0
95
96        ! verifier que HDO est raisonable
97         if (q.gt.qmin) then
98             IF(ltnat1) THEN
99                tnat = 1.0
100             ELSE IF(getKey('tnat', tnat, isoName(iso))) THEN
101                  err_msg = 'Missing isotopic parameter "tnat"'
102                  iso_verif_aberrant_nostop=1
103                  RETURN
104             END IF
105             deltaD=(x/q/tnat-1)*1000
106             if ((deltaD.gt.deltaDmax).or.(deltaD.lt.deltaDmin)) then
107                  write(*,*) 'erreur detectee par iso_verif_aberrant:'
108                  write(*,*) err_msg
109                  write(*,*) 'q=',q
110                  write(*,*) 'deltaD=',deltaD
111                  write(*,*) 'iso=',iso
112                  iso_verif_aberrant_nostop=1
113             endif !if ((deltaD.gt.deltaDmax).or.(deltaD.lt.deltaDmin)) then
114          endif !if (q(i,k,iq).gt.qmin) then
115
116       
117        return
118        end       
119
Note: See TracBrowser for help on using the repository browser.