source: LMDZ5/branches/testing/libf/phylmd/undefSTD.F90 @ 2314

Last change on this file since 2314 was 2314, checked in by musat, 9 years ago

Corrections missing_value pour compilation en mode debug
JG/IM

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 3.3 KB
Line 
1
2! $Id: undefSTD.F90 2314 2015-06-25 13:40:53Z musat $
3
4SUBROUTINE undefstd(itap, read_climoz)
5  USE netcdf
6  USE dimphy
7#ifdef CPP_IOIPSL
8  USE phys_state_var_mod
9#endif
10#ifdef CPP_XIOS
11  USE wxios, ONLY: missing_val
12#endif
13
14  IMPLICIT NONE
15  include "clesphys.h"
16#ifndef CPP_XIOS
17  REAL :: missing_val
18#endif
19
20  ! ====================================================================
21
22  ! I. Musat : 09.2004
23
24  ! Calcul * du nombre de pas de temps (FLOAT(ecrit_XXX)-tnondef))
25  ! ou la variable tlevSTD est bien definie (.NE.missing_val),
26  ! et
27  ! * de la somme de tlevSTD => tsumSTD
28
29  ! nout=1 !var. journaliere "day" moyenne sur tous les pas de temps
30  ! ! de la physique
31  ! nout=2 !var. mensuelle "mth" moyennee sur tous les pas de temps
32  ! ! de la physique
33  ! nout=3 !var. mensuelle "NMC" moyennee toutes les ecrit_hf
34
35
36  ! NB: mettre "inst(X)" dans le write_hist*NMC.h !
37  ! ====================================================================
38
39  ! ym#include "dimensions.h"
40  ! ym      integer jjmp1
41  ! ym      parameter (jjmp1=jjm+1-1/jjm)
42  ! ym#include "dimphy.h"
43  ! variables Input
44
45  ! INTEGER nlevSTD, klevSTD, itap
46  ! PARAMETER(klevSTD=17)
47  INTEGER itap
48  ! REAL dtime
49
50  ! variables locales
51  ! INTEGER i, k, nout, n
52  ! PARAMETER(nout=3) !nout=1 : day; =2 : mth; =3 : NMC
53  INTEGER i, k, n
54  ! REAL freq_calNMC(nout)
55  INTEGER read_climoz
56
57  ! variables Output
58  ! REAL tlevSTD(klon,klevSTD), tsumSTD(klon,klevSTD,nout)
59  ! LOGICAL oknondef(klon,klevSTD,nout)
60  ! REAL tnondef(klon,klevSTD,nout)
61
62! REAL missing_val
63
64! missing_val = nf90_fill_real
65#ifndef CPP_XIOS
66      missing_val=missing_val_nf90
67#endif
68
69  DO n = 1, nout
70
71
72    ! calcul variables tous les freq_calNMC(n)/dtime pas de temps
73    ! de la physique
74
75    IF (mod(itap,nint(freq_calnmc(n)/dtime))==0) THEN
76      DO k = 1, nlevstd
77        DO i = 1, klon
78          IF (tlevstd(i,k)==missing_val) THEN
79            ! IF(oknondef(i,k,n)) THEN
80            tnondef(i, k, n) = tnondef(i, k, n) + 1.
81            ! ENDIF !oknondef(i,k)
82
83          ELSE IF (tlevstd(i,k)/=missing_val) THEN
84            tsumstd(i, k, n) = tsumstd(i, k, n) + tlevstd(i, k)
85            usumstd(i, k, n) = usumstd(i, k, n) + ulevstd(i, k)
86            vsumstd(i, k, n) = vsumstd(i, k, n) + vlevstd(i, k)
87            wsumstd(i, k, n) = wsumstd(i, k, n) + wlevstd(i, k)
88            phisumstd(i, k, n) = phisumstd(i, k, n) + philevstd(i, k)
89            qsumstd(i, k, n) = qsumstd(i, k, n) + qlevstd(i, k)
90            rhsumstd(i, k, n) = rhsumstd(i, k, n) + rhlevstd(i, k)
91            uvsumstd(i, k, n) = uvsumstd(i, k, n) + uvstd(i, k)
92            vqsumstd(i, k, n) = vqsumstd(i, k, n) + vqstd(i, k)
93            vtsumstd(i, k, n) = vtsumstd(i, k, n) + vtstd(i, k)
94            wqsumstd(i, k, n) = wqsumstd(i, k, n) + wqstd(i, k)
95            vphisumstd(i, k, n) = vphisumstd(i, k, n) + vphistd(i, k)
96            wtsumstd(i, k, n) = wtsumstd(i, k, n) + wtstd(i, k)
97            u2sumstd(i, k, n) = u2sumstd(i, k, n) + u2std(i, k)
98            v2sumstd(i, k, n) = v2sumstd(i, k, n) + v2std(i, k)
99            t2sumstd(i, k, n) = t2sumstd(i, k, n) + t2std(i, k)
100            o3sumstd(i, k, n) = o3sumstd(i, k, n) + o3std(i, k)
101            IF (read_climoz==2) o3daysumstd(i, k, n) = o3daysumstd(i, k, n) + &
102              o3daystd(i, k)
103
104          END IF
105        END DO !i
106      END DO !k
107
108    END IF !MOD(itap,NINT(freq_calNMC(n)/dtime)).EQ.0
109
110  END DO !n
111
112  RETURN
113END SUBROUTINE undefstd
Note: See TracBrowser for help on using the repository browser.