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

Contents of /trunk/libf/phylmd/undefSTD.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 7 - (show annotations)
Mon Mar 31 12:24:17 2008 UTC (16 years, 1 month ago) by guez
File size: 2364 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 SUBROUTINE undefSTD(nlevSTD,itap,tlevSTD,
2 $ ecrit_hf,
3 $ oknondef,tnondef,tsumSTD)
4 use dimens_m
5 use dimphy
6 IMPLICIT none
7 c
8 c====================================================================
9 c
10 c I. Musat : 09.2004
11 c
12 c Calcul * du nombre de pas de temps (FLOAT(ecrit_XXX)-tnondef))
13 c ou la variable tlevSTD est bien definie (.NE.1.E+20),
14 c et
15 c * de la somme de tlevSTD => tsumSTD
16 c
17 c nout=1 !var. journaliere "day" moyenne sur tous les pas de temps
18 c ! de la physique
19 c nout=2 !var. mensuelle "mth" moyennee sur tous les pas de temps
20 c ! de la physique
21 c nout=3 !var. mensuelle "NMC" moyennee toutes les 6heures
22 c
23 c
24 c NB: mettre "inst(X)" dans le write_histXXX.h !
25 c====================================================================
26 c
27 integer jjmp1
28 parameter (jjmp1=jjm+1-1/jjm)
29 c variables Input
30 INTEGER nlevSTD, klevSTD
31 integer, intent(in):: itap
32 PARAMETER(klevSTD=17)
33 INTEGER, intent(in):: ecrit_hf
34 c
35 c variables locales
36 INTEGER i, k, nout
37 PARAMETER(nout=3) !nout=1 : day; =2 : mth; =3 : NMC
38 c
39 c variables Output
40 REAL tlevSTD(klon,klevSTD), tsumSTD(klon,klevSTD,nout)
41 LOGICAL oknondef(klon,klevSTD,nout)
42 REAL tnondef(klon,klevSTD,nout)
43 c
44 c calcul variables tous les pas de temps de la physique
45 c
46 DO k=1, nlevSTD
47 DO i=1, klon
48 IF(tlevSTD(i,k).EQ.1E+20) THEN
49 IF(oknondef(i,k,1)) THEN
50 tnondef(i,k,1)=tnondef(i,k,1)+1.
51 ENDIF !oknondef(i,k)
52 c
53 IF(oknondef(i,k,2)) THEN
54 tnondef(i,k,2)=tnondef(i,k,2)+1.
55 ENDIF !oknondef(i,k)
56 c
57 ELSE IF(tlevSTD(i,k).NE.1E+20) THEN
58 tsumSTD(i,k,1)=tsumSTD(i,k,1)+tlevSTD(i,k)
59 tsumSTD(i,k,2)=tsumSTD(i,k,2)+tlevSTD(i,k)
60 ENDIF
61 ENDDO !i
62 ENDDO !k
63 c
64 c calcul variables toutes les 6h
65 c
66 IF(MOD(itap,ecrit_hf).EQ.0) THEN
67 c
68 DO k=1, nlevSTD
69 DO i=1, klon
70 IF(tlevSTD(i,k).EQ.1E+20) THEN
71 IF(oknondef(i,k,3)) THEN
72 tnondef(i,k,3)=tnondef(i,k,3)+1.
73 ENDIF !oknondef(i,k)
74 c
75 ELSE IF(tlevSTD(i,k).NE.1E+20) THEN
76 tsumSTD(i,k,3)=tsumSTD(i,k,3)+tlevSTD(i,k)
77 ENDIF
78 ENDDO !i
79 ENDDO !k
80
81 ENDIF !MOD(itap,ecrit_hf).EQ.0
82 c
83 RETURN
84 END
  ViewVC Help
Powered by ViewVC 1.1.21