source: LMDZ4/trunk/libf/phylmd/undefSTD.F @ 1398

Last change on this file since 1398 was 1398, checked in by musat, 14 years ago

Last corrections for CMIP5:

  • Add O3 at standard level files histmthNMC.nc
  • Add positive attribute "down" for vertical axes for all output files
  • Replace "inst" by "ave" for hist*NMC.nc files to have time_counter and bounds for time axis (Marie-Alice's hint)
  • Correct units for vertical axes : mb instead of hPa
  • Add mass flux at the bottom of clouds
  • Comment non initialized variables (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) for the output files
  • Geopotential field phy850, phi700, phi500, etc are modified to "geopotential height and are called z850, z700, z500, etc
  • Meaning of specific humidity outputs - ovapinit and ovap - were interchanged
  • Fields albs, albslw become alb1, alb2 in output files
  • Correct title for rugs_* fields
  • Correct units for pbase and ptop are Pa (not mb)
  • Correct ndayrain field

FH/JLD/JYG/MAF/IM

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 3.0 KB
Line 
1!
2! $Header$
3!
4      SUBROUTINE undefSTD(itap,freq_calNMC)
5      USE netcdf
6      USE dimphy
7      USE phys_state_var_mod ! Variables sauvegardees de la physique
8      IMPLICIT none
9c
10c====================================================================
11c
12c I. Musat : 09.2004
13c
14c Calcul * du nombre de pas de temps (FLOAT(ecrit_XXX)-tnondef))
15c          ou la variable tlevSTD est bien definie (.NE.missing_val),
16c et
17c        * de la somme de tlevSTD => tsumSTD
18c
19c nout=1 !var. journaliere "day" moyenne sur tous les pas de temps
20c        ! de la physique
21c nout=2 !var. mensuelle "mth" moyennee sur tous les pas de temps
22c        ! de la physique
23c nout=3 !var. mensuelle "NMC" moyennee toutes les ecrit_hf
24c
25c
26c NB: mettre "inst(X)" dans le write_hist*NMC.h !
27c====================================================================
28c
29cym#include "dimensions.h"
30cym      integer jjmp1
31cym      parameter (jjmp1=jjm+1-1/jjm)
32cym#include "dimphy.h"
33c variables Input
34c
35c     INTEGER nlevSTD, klevSTD, itap
36c     PARAMETER(klevSTD=17)
37      INTEGER itap
38c     REAL dtime
39c
40c variables locales
41c     INTEGER i, k, nout, n
42c     PARAMETER(nout=3) !nout=1 : day; =2 : mth; =3 : NMC
43      INTEGER i, k, n
44      REAL freq_calNMC(nout)
45c
46c variables Output
47c     REAL tlevSTD(klon,klevSTD), tsumSTD(klon,klevSTD,nout)
48c     LOGICAL oknondef(klon,klevSTD,nout)
49c     REAL tnondef(klon,klevSTD,nout)
50c
51      REAL missing_val
52c
53      missing_val=nf90_fill_real
54c
55      DO n=1, nout
56c
57c
58c calcul variables tous les freq_calNMC(n)/dtime pas de temps
59c de la physique
60c
61       IF(MOD(itap,NINT(freq_calNMC(n)/dtime)).EQ.0) THEN
62        DO k=1, nlevSTD
63         DO i=1, klon
64          IF(tlevSTD(i,k).EQ.missing_val) THEN
65c          IF(oknondef(i,k,n)) THEN         
66            tnondef(i,k,n)=tnondef(i,k,n)+1.
67c          ENDIF !oknondef(i,k)
68c
69          ELSE IF(tlevSTD(i,k).NE.missing_val) THEN
70           tsumSTD(i,k,n)=tsumSTD(i,k,n)+tlevSTD(i,k)
71           usumSTD(i,k,n)=usumSTD(i,k,n)+ulevSTD(i,k)
72           vsumSTD(i,k,n)=vsumSTD(i,k,n)+vlevSTD(i,k)
73           wsumSTD(i,k,n)=wsumSTD(i,k,n)+wlevSTD(i,k)
74           phisumSTD(i,k,n)=phisumSTD(i,k,n)+philevSTD(i,k)
75           qsumSTD(i,k,n)=qsumSTD(i,k,n)+qlevSTD(i,k)
76           rhsumSTD(i,k,n)=rhsumSTD(i,k,n)+rhlevSTD(i,k)
77           uvsumSTD(i,k,n)=uvsumSTD(i,k,n)+uvSTD(i,k)
78           vqsumSTD(i,k,n)=vqsumSTD(i,k,n)+vqSTD(i,k)
79           vTsumSTD(i,k,n)=vTsumSTD(i,k,n)+vTSTD(i,k)
80           wqsumSTD(i,k,n)=wqsumSTD(i,k,n)+wqSTD(i,k)
81           vphisumSTD(i,k,n)=vphisumSTD(i,k,n)+vphiSTD(i,k)
82           wTsumSTD(i,k,n)=wTsumSTD(i,k,n)+wTSTD(i,k)
83           u2sumSTD(i,k,n)=u2sumSTD(i,k,n)+u2STD(i,k)
84           v2sumSTD(i,k,n)=v2sumSTD(i,k,n)+v2STD(i,k)
85           T2sumSTD(i,k,n)=T2sumSTD(i,k,n)+T2STD(i,k)
86           O3sumSTD(i,k,n)=O3sumSTD(i,k,n)+O3STD(i,k)
87           O3daysumSTD(i,k,n)=O3daysumSTD(i,k,n)+O3daySTD(i,k)
88          ENDIF
89         ENDDO !i
90        ENDDO !k
91c
92       ENDIF !MOD(itap,NINT(freq_calNMC(n)/dtime)).EQ.0
93c
94      ENDDO !n
95c
96      RETURN
97      END 
Note: See TracBrowser for help on using the repository browser.