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

Last change on this file since 1377 was 1374, checked in by musat, 15 years ago

Add logical flags ok_histNMC(3) to manage output of the
standard pressure files hist*NMC.nc
IM

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