source: LMDZ4/trunk/libf/phylmd/write_histmthNMC.h @ 693

Last change on this file since 693 was 684, checked in by lmdzadmin, 19 years ago

Ajout nouveaux flags (ip_ebil_phy, ok_slab_sicOBS) ds physiq.def +
Lecture freq. ecriture en nombre de jours
IM
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.7 KB
RevLine 
[524]1!
2! $Header$
3!
4      IF (ok_mensuel) THEN
5c
[644]6       ndex3d = 0
7       itau_w = itau_phy + itap
8ccc
9c  Champs interpolles sur des niveaux de pression du NMC
[524]10c
[644]11c     PARAMETER(nout=3) !nout=1 : day; =2 : mth; =3 : NMC
12ccc
13      IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN
14ccc
15       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,tsumSTD(:,:,2),
16     $      zx_tmp_3d)
17       CALL histwrite(nid_nmc,"temp",itau_w,zx_tmp_3d,
18     $       iim*jjmp1*nlevSTD,ndex3d)
[524]19c
[644]20       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,phisumSTD(:,:,2),
21     $     zx_tmp_3d)
22       CALL histwrite(nid_nmc,"phi",itau_w,zx_tmp_3d,
23     $       iim*jjmp1*nlevSTD,ndex3d)
[524]24c
[644]25       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,qsumSTD(:,:,2),
26     $     zx_tmp_3d)
27       CALL histwrite(nid_nmc,"q",itau_w,zx_tmp_3d,
28     $       iim*jjmp1*nlevSTD,ndex3d)
[524]29c
[644]30       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,rhsumSTD(:,:,2),
31     $     zx_tmp_3d)
32       CALL histwrite(nid_nmc,"rh",itau_w,zx_tmp_3d,
33     $       iim*jjmp1*nlevSTD,ndex3d)
[524]34c
[644]35       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,usumSTD(:,:,2),
36     $     zx_tmp_3d)
37       CALL histwrite(nid_nmc,"u",itau_w,zx_tmp_3d,
38     $       iim*jjmp1*nlevSTD,ndex3d)
39c
40       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vsumSTD(:,:,2),
41     $     zx_tmp_3d)
42       CALL histwrite(nid_nmc,"v",itau_w,zx_tmp_3d,
43     $       iim*jjmp1*nlevSTD,ndex3d)
44ccc
45      ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
46ccc
47       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,tsumSTD(:,:,3),
48     $     zx_tmp_3d)
49       CALL histwrite(nid_nmc,"temp",itau_w,zx_tmp_3d,
50     $       iim*jjmp1*nlevSTD,ndex3d)
51c
52       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,phisumSTD(:,:,3),
53     $     zx_tmp_3d)
54       CALL histwrite(nid_nmc,"phi",itau_w,zx_tmp_3d,
55     $       iim*jjmp1*nlevSTD,ndex3d)
56c
57       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,qsumSTD(:,:,3),
58     $     zx_tmp_3d)
59       CALL histwrite(nid_nmc,"q",itau_w,zx_tmp_3d,
60     $       iim*jjmp1*nlevSTD,ndex3d)
61c
62       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,rhsumSTD(:,:,3),
63     $     zx_tmp_3d)
64       CALL histwrite(nid_nmc,"rh",itau_w,zx_tmp_3d,
65     $       iim*jjmp1*nlevSTD,ndex3d)
66c
67       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,usumSTD(:,:,3),
68     $     zx_tmp_3d)
69       CALL histwrite(nid_nmc,"u",itau_w,zx_tmp_3d,
70     $       iim*jjmp1*nlevSTD,ndex3d)
71c
72       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vsumSTD(:,:,3),
73     $     zx_tmp_3d)
74       CALL histwrite(nid_nmc,"v",itau_w,zx_tmp_3d,
75     $       iim*jjmp1*nlevSTD,ndex3d)
76c
77       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,wsumSTD(:,:,3),
78     $     zx_tmp_3d)
79       CALL histwrite(nid_nmc,"w",itau_w,zx_tmp_3d,
80     $       iim*jjmp1*nlevSTD,ndex3d)
81c
82       DO k=1, nlevSTD
83        DO i=1, klon
84         IF(tnondef(i,k,3).NE.1.E+20) THEN
[684]85          zx_tmp_fi3d(i,k) = (100.*tnondef(i,k,3))/ecrit_hf2mth
[644]86         ELSE
87          zx_tmp_fi3d(i,k) = 1.E+20
[524]88         ENDIF
[644]89        ENDDO
90       ENDDO !k=1, nlevSTD
[524]91c
[644]92       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
93       CALL histwrite(nid_nmc,"psbg",itau_w,zx_tmp_3d,
94     $       iim*jjmp1*nlevSTD,ndex3d)
[524]95c
[644]96       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,uvsumSTD(:,:,3),
97     $     zx_tmp_3d)
98       CALL histwrite(nid_nmc,"uv",itau_w,zx_tmp_3d,
99     $       iim*jjmp1*nlevSTD,ndex3d)
[524]100c
[644]101       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vqsumSTD(:,:,3),
102     $     zx_tmp_3d)
103       CALL histwrite(nid_nmc,"vq",itau_w,zx_tmp_3d,
104     $       iim*jjmp1*nlevSTD,ndex3d)
[524]105c
[644]106       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vTsumSTD(:,:,3),
107     $     zx_tmp_3d)
108       CALL histwrite(nid_nmc,"vT",itau_w,zx_tmp_3d,
109     $       iim*jjmp1*nlevSTD,ndex3d)
[524]110c
[644]111       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1, wqsumSTD(:,:,3),
112     $     zx_tmp_3d)
113       CALL histwrite(nid_nmc,"wq",itau_w,zx_tmp_3d,
114     $       iim*jjmp1*nlevSTD,ndex3d)
[524]115c
[644]116       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vphisumSTD(:,:,3),
117     $     zx_tmp_3d)
118       CALL histwrite(nid_nmc,"vphi",itau_w,zx_tmp_3d,
119     $       iim*jjmp1*nlevSTD,ndex3d)
[524]120c
[644]121       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,wTsumSTD(:,:,3),
122     $     zx_tmp_3d)
123       CALL histwrite(nid_nmc,"wT",itau_w,zx_tmp_3d,
124     $       iim*jjmp1*nlevSTD,ndex3d)
[524]125c
[644]126       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,u2sumSTD(:,:,3),
127     $     zx_tmp_3d)
128       CALL histwrite(nid_nmc,"uxu",itau_w,zx_tmp_3d,
129     $       iim*jjmp1*nlevSTD,ndex3d)
130c
131       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,v2sumSTD(:,:,3),
132     $     zx_tmp_3d)
133       CALL histwrite(nid_nmc,"vxv",itau_w,zx_tmp_3d,
134     $       iim*jjmp1*nlevSTD,ndex3d)
135c
136       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,T2sumSTD(:,:,3),
137     $     zx_tmp_3d)
138       CALL histwrite(nid_nmc,"TxT",itau_w,zx_tmp_3d,
139     $       iim*jjmp1*nlevSTD,ndex3d)
140c
141      ENDIF !type_run
142c
[524]143      if (ok_sync) then
144        call histsync(nid_nmc)
145      endif
146
147      ENDIF
Note: See TracBrowser for help on using the repository browser.