source: LMDZ4/trunk/libf/phylmd/write_histhfNMC.h @ 1351

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

Add new output files for variables on standard pressure levels at daily and high
frequency (6hr)
IM

  • Property svn:executable set to *
File size: 1.8 KB
Line 
1!
2! $Header$
3!
4       ndex3d = 0
5       itau_w = itau_phy + itap
6ccc
7c  Champs interpolles sur des niveaux de pression du NMC
8c
9c     PARAMETER(nout=3) 
10c nout=1 : in=pdtphys,    out=mth
11c nout=2 : in=pdtphys,    out=day
12c nout=3 : in=pdtphys,    out=hf
13ccc
14       CALL histwrite_phy(nid_hfnmc,"tnondef",itau_w,tnondef(:,:,3))
15c
16       CALL histwrite_phy(nid_hfnmc,"ta",itau_w,twriteSTD3)
17c
18       CALL histwrite_phy(nid_hfnmc,"zg",itau_w,phiwriteSTD3)
19c
20       CALL histwrite_phy(nid_hfnmc,"hus",itau_w,qwriteSTD3)
21c
22       CALL histwrite_phy(nid_hfnmc,"hur",itau_w,rhwriteSTD3)
23c
24       CALL histwrite_phy(nid_hfnmc,"ua",itau_w,uwriteSTD3)
25c
26       CALL histwrite_phy(nid_hfnmc,"va",itau_w,vwriteSTD3)
27c
28       CALL histwrite_phy(nid_hfnmc,"wap",itau_w,wwriteSTD3)
29c
30       DO k=1, nlevSTD
31        DO i=1, klon
32         IF(tnondef(i,k,3).NE.missing_val) THEN
33          zx_tmp_fiNC(i,k) = (100.*tnondef(i,k,3))/ecrit_hf2mth
34         ELSE
35          zx_tmp_fiNC(i,k) = missing_val
36         ENDIF
37        ENDDO
38       ENDDO !k=1, nlevSTD
39c
40       CALL histwrite_phy(nid_hfnmc,"psbg",itau_w,zx_tmp_fiNC)
41c
42       CALL histwrite_phy(nid_hfnmc,"uv",itau_w,uvsumSTD(:,:,3))
43c
44       CALL histwrite_phy(nid_hfnmc,"vq",itau_w,vqsumSTD(:,:,3))
45c
46       CALL histwrite_phy(nid_hfnmc,"vT",itau_w,vTsumSTD(:,:,3))
47c
48       CALL histwrite_phy(nid_hfnmc,"wq",itau_w,wqsumSTD(:,:,3))
49c
50       CALL histwrite_phy(nid_hfnmc,"vphi",itau_w,vphisumSTD(:,:,3))
51c
52       CALL histwrite_phy(nid_hfnmc,"wT",itau_w,wTsumSTD(:,:,3))
53c
54       CALL histwrite_phy(nid_hfnmc,"uxu",itau_w,u2sumSTD(:,:,3))
55c
56       CALL histwrite_phy(nid_hfnmc,"vxv",itau_w,v2sumSTD(:,:,3))
57c
58       CALL histwrite_phy(nid_hfnmc,"TxT",itau_w,T2sumSTD(:,:,3))
59c
60c     ENDIF !type_run
61c
62      if (ok_sync) then
63c$OMP MASTER
64        call histsync(nid_hfnmc)
65c$OMP END MASTER
66      endif
67
68c     ENDIF
Note: See TracBrowser for help on using the repository browser.