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

Last change on this file since 1085 was 776, checked in by Laurent Fairhead, 18 years ago

Suite du merge entre la version et la HEAD: quelques modifications
de Yann sur le
LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.2 KB
Line 
1!
2! $Header$
3!
4      IF (ok_mensuel) THEN
5c
6       ndex3d = 0
7       itau_w = itau_phy + itap
8ccc
9c  Champs interpolles sur des niveaux de pression du NMC
10c
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
15cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,tsumSTD(:,:,2),
16cym     $      zx_tmp_NC)
17       CALL histwrite_phy(nid_nmc,"temp",itau_w,tsumSTD(:,:,2))
18c
19cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,phisumSTD(:,:,2),
20cym     $     zx_tmp_NC)
21       CALL histwrite_phy(nid_nmc,"phi",itau_w,phisumSTD(:,:,2))
22c
23cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,qsumSTD(:,:,2),
24cym     $     zx_tmp_NC)
25       CALL histwrite_phy(nid_nmc,"q",itau_w,qsumSTD(:,:,2))
26c
27cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,rhsumSTD(:,:,2),
28cym     $     zx_tmp_NC)
29       CALL histwrite_phy(nid_nmc,"rh",itau_w,rhsumSTD(:,:,2))
30c
31cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,usumSTD(:,:,2),
32cym     $     zx_tmp_NC)
33       CALL histwrite_phy(nid_nmc,"u",itau_w,usumSTD(:,:,2))
34c
35cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vsumSTD(:,:,2),
36cym     $     zx_tmp_NC)
37       CALL histwrite_phy(nid_nmc,"v",itau_w,vsumSTD(:,:,2))
38ccc
39      ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
40ccc
41cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,tsumSTD(:,:,3),
42cym     $     zx_tmp_NC)
43       CALL histwrite_phy(nid_nmc,"temp",itau_w,tsumSTD(:,:,3))
44c
45cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,phisumSTD(:,:,3),
46cym     $     zx_tmp_NC)
47       CALL histwrite_phy(nid_nmc,"phi",itau_w,phisumSTD(:,:,3))
48c
49cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,qsumSTD(:,:,3),
50cym     $     zx_tmp_NC)
51       CALL histwrite_phy(nid_nmc,"q",itau_w,qsumSTD(:,:,3))
52c
53cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,rhsumSTD(:,:,3),
54cym     $     zx_tmp_NC)
55       CALL histwrite_phy(nid_nmc,"rh",itau_w,rhsumSTD(:,:,3))
56c
57cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,usumSTD(:,:,3),
58cym     $     zx_tmp_NC)
59       CALL histwrite_phy(nid_nmc,"u",itau_w,usumSTD(:,:,3))
60c
61cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vsumSTD(:,:,3),
62cym     $     zx_tmp_NC)
63       CALL histwrite_phy(nid_nmc,"v",itau_w,vsumSTD(:,:,3))
64c
65cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,wsumSTD(:,:,3),
66cym     $     zx_tmp_NC)
67       CALL histwrite_phy(nid_nmc,"w",itau_w,wsumSTD(:,:,3))
68c
69       DO k=1, nlevSTD
70        DO i=1, klon
71         IF(tnondef(i,k,3).NE.1.E+20) THEN
72          zx_tmp_fiNC(i,k) = (100.*tnondef(i,k,3))/ecrit_hf2mth
73         ELSE
74          zx_tmp_fiNC(i,k) = 1.E+20
75         ENDIF
76        ENDDO
77       ENDDO !k=1, nlevSTD
78c
79cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,zx_tmp_fiNC,zx_tmp_NC)
80       CALL histwrite_phy(nid_nmc,"psbg",itau_w,zx_tmp_fiNC)
81c
82cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,uvsumSTD(:,:,3),
83cym     $     zx_tmp_NC)
84       CALL histwrite_phy(nid_nmc,"uv",itau_w,uvsumSTD(:,:,3))
85c
86cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vqsumSTD(:,:,3),
87cym     $     zx_tmp_NC)
88       CALL histwrite_phy(nid_nmc,"vq",itau_w,vqsumSTD(:,:,3))
89c
90cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vTsumSTD(:,:,3),
91cym     $     zx_tmp_NC)
92       CALL histwrite_phy(nid_nmc,"vT",itau_w,vTsumSTD(:,:,3))
93c
94cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1, wqsumSTD(:,:,3),
95cym     $     zx_tmp_NC)
96       CALL histwrite_phy(nid_nmc,"wq",itau_w,wqsumSTD(:,:,3))
97c
98cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vphisumSTD(:,:,3),
99cym     $     zx_tmp_NC)
100       CALL histwrite_phy(nid_nmc,"vphi",itau_w,vphisumSTD(:,:,3))
101c
102cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,wTsumSTD(:,:,3),
103cym     $     zx_tmp_NC)
104       CALL histwrite_phy(nid_nmc,"wT",itau_w,wTsumSTD(:,:,3))
105c
106cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,u2sumSTD(:,:,3),
107cym     $     zx_tmp_NC)
108       CALL histwrite_phy(nid_nmc,"uxu",itau_w,u2sumSTD(:,:,3))
109c
110cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,v2sumSTD(:,:,3),
111cym     $     zx_tmp_NC)
112       CALL histwrite_phy(nid_nmc,"vxv",itau_w,v2sumSTD(:,:,3))
113c
114cym       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,T2sumSTD(:,:,3),
115cym     $     zx_tmp_NC)
116       CALL histwrite_phy(nid_nmc,"TxT",itau_w,T2sumSTD(:,:,3))
117c
118      ENDIF !type_run
119c
120      if (ok_sync) then
121c$OMP MASTER
122        call histsync(nid_nmc)
123c$OMP END MASTER
124      endif
125
126      ENDIF
Note: See TracBrowser for help on using the repository browser.