/[lmdze]/trunk/libf/phylmd/moy_undefSTD.f
ViewVC logotype

Annotation of /trunk/libf/phylmd/moy_undefSTD.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 7 - (hide annotations)
Mon Mar 31 12:24:17 2008 UTC (16 years, 2 months ago) by guez
File size: 2258 byte(s) 
This revision is not in working order. Pending some moving of files.

Important changes. In the program "etat0_lim": ozone coefficients from
Mobidic are regridded in time instead of pressure ; consequences in
"etat0". In the program "gcm", ozone coefficients from Mobidic are
read once per day only for the current day and regridded in pressure ;
consequences in "o3_chem_m", "regr_pr_coefoz", "phytrac" and
"regr_pr_comb_coefoz_m".

NetCDF95 is a library and does not export NetCDF.

New variables "nag_gl_options", "nag_fcalls_options" and
"nag_cross_options" in "nag_tools.mk".

"check_coefoz.jnl" rewritten entirely for new version of
"coefoz_LMDZ.nc".

Target "obj_etat0_lim" moved from "GNUmakefile" to "nag_rules.mk".

Added some "intent" attributes in "calfis", "clmain", "clqh",
"cltrac", "cltracrn", "cvltr", "ini_undefSTD", "moy_undefSTD",
"nflxtr", "phystokenc", "phytrac", "readsulfate", "readsulfate_preind"
and "undefSTD".

In "dynetat0", "dynredem0" and "gcm", "phis" has rank 2 instead of
1. "phis" has assumed shape in "dynredem0".

Added module containing "dynredem0". Changed some calls with NetCDF
Fortran 77 interface to calls with NetCDF95 interface.

Replaced calls to "ssum" by calls to "sum" in "inigeom".

In "make.sh", new option "-c" to change compiler.

In "aaam_bud", argument "rjour" deleted.

In "physiq": renamed some variables; deleted variable "xjour".

In "phytrac": renamed some variables; new argument "lmt_pas".
1 guez 3 SUBROUTINE moy_undefSTD(nlevSTD,itap,
2     $ ecrit_day,ecrit_mth,ecrit_hf2mth,
3     $ tnondef,tsumSTD)
4     use dimens_m
5     use dimphy
6     IMPLICIT none
7     c moyenne d'1 var a 1 niveau de pression
8     c====================================================================
9     c
10     c I. Musat : 09.2004
11     c
12     c Moyenne - a des frequences differentes - des valeurs bien definies
13     c (.NE.1.E+20) des variables interpolees a un niveau de
14     c pression.
15     c 1) les variables de type "day" (nout=1) ou "mth" (nout=2) sont sommees
16     c tous les pas de temps de la physique
17     c
18     c 2) les variables de type "NMC" (nout=3) sont calculees a partir
19     c des valeurs instantannees toutes les 6 heures
20     c
21     c
22     c NB: mettre "inst(X)" dans le write_histXXX.h !
23     c====================================================================
24     integer jjmp1
25     parameter (jjmp1=jjm+1-1/jjm)
26     c
27     c variables Input
28 guez 7 INTEGER nlevSTD, klevSTD
29     integer, intent(in):: itap
30 guez 3 PARAMETER(klevSTD=17)
31     INTEGER ecrit_day, ecrit_mth, ecrit_hf2mth
32     c
33     c variables locales
34     INTEGER i, k, nout
35     PARAMETER(nout=3) !nout=1 day/nout=2 mth/nout=3 NMC
36     c
37     c variables Output
38     REAL tnondef(klon,klevSTD,nout)
39     REAL tsumSTD(klon,klevSTD,nout)
40     c
41     c calcul 1 fois par jour
42     c
43     IF(MOD(itap,ecrit_day).EQ.0) THEN
44     DO k=1, nlevSTD
45     DO i=1, klon
46     IF(tnondef(i,k,1).NE.FLOAT(ecrit_day)) THEN
47     tsumSTD(i,k,1)=tsumSTD(i,k,1)/
48     $ (FLOAT(ecrit_day)-tnondef(i,k,1))
49     ELSE
50     tsumSTD(i,k,1)=1.E+20
51     ENDIF !tnondef
52     ENDDO !i
53     ENDDO !k
54     ENDIF !MOD(itap,ecrit_day).EQ.0
55     c
56     c calcul 1 fois par mois
57     c
58     IF(MOD(itap,ecrit_mth).EQ.0) THEN
59     DO k=1, nlevSTD
60     DO i=1, klon
61     IF(tnondef(i,k,2).NE.FLOAT(ecrit_mth)) THEN
62     tsumSTD(i,k,2)=tsumSTD(i,k,2)/
63     $ (FLOAT(ecrit_mth)-tnondef(i,k,2))
64     ELSE
65     tsumSTD(i,k,2)=1.E+20
66     ENDIF !tnondef
67     c
68     IF(tnondef(i,k,3).NE.FLOAT(ecrit_hf2mth)) THEN
69     tsumSTD(i,k,3)=tsumSTD(i,k,3)/
70     $ (FLOAT(ecrit_hf2mth)-tnondef(i,k,3))
71     ELSE
72     tsumSTD(i,k,3)=1.E+20
73     ENDIF !tnondef
74     c
75     ENDDO !i
76     ENDDO !k
77     ENDIF !MOD
78     c
79     RETURN
80     END
  ViewVC Help
Powered by ViewVC 1.1.21