source: LMDZ5/branches/LMDZ5_AR5/libf/phylmd/write_histhfNMC.h

Last change on this file was 1539, checked in by musat, 13 years ago

Ajouts CFMIP2/CMIP5

  • 6eme fichier de sortie "stations" histstn.nc qui necessite 2 fichiers (voir DefLists?): npCFMIP_param.data(_*) contenant le nombre de points (120 pour simulations AMIP, 73 pour aqua) pointlocations.txt(_*) contenant le numero, les coordonnees (lon,lat) et le nom de chaque station
  • flag LOGICAL dans tous les appels histwrite_phy pour pouvoir sortir le fichier histstn.nc

NB: 1) les flags de type phys_ que l'on met dans le physiq.def_L* pour ajouter plus de sorties

necessitent dorenavant 6 valeurs, la 6eme correspondant au fichier histstn.nc

2) par defaut le fichier histstn.nc ne sort pas; pour le sortir ajouter les lignes suivantes

dans physiq.def_L*

### Type de fichier : global (n) ou stations (y)
phys_out_filestations = n n n n n y

  • introduction de 2 jeux de flags pour les taux des GES; taux actuels avec suffixes _act, taux futurs avec "_per" avec 2 appels au rayonnement si taux "_per" different des taux "_act" (utiles pour diags. CFMIP 4CO2)
  • flags "betaCRF" pour calculs CRF pour experiences sensibilite proprietes optiques eau liquide nuageuse avec initialisations par defaut; sinon besoin de fichier beta_crf.data

Ajout flag LOGICAL lCOSP necessaire pour sortir un fichier stations
IM

  • Property svn:executable set to *
File size: 2.1 KB
Line 
1!
2! $Header$
3!
4      IF (ok_histNMC(3)) 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) 
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,lNMC,"tnondef",itau_w,
17     $tnondef(:,:,3))
18c
19       CALL histwrite_phy(nid_hfnmc,lNMC,"ta",itau_w,
20     $twriteSTD3)
21c
22       CALL histwrite_phy(nid_hfnmc,lNMC,"zg",itau_w,
23     $phiwriteSTD3)
24c
25       CALL histwrite_phy(nid_hfnmc,lNMC,"hus",itau_w,
26     $qwriteSTD3)
27c
28       CALL histwrite_phy(nid_hfnmc,lNMC,"hur",itau_w,
29     $rhwriteSTD3)
30c
31       CALL histwrite_phy(nid_hfnmc,lNMC,"ua",itau_w,
32     $uwriteSTD3)
33c
34       CALL histwrite_phy(nid_hfnmc,lNMC,"va",itau_w,
35     $vwriteSTD3)
36c
37       CALL histwrite_phy(nid_hfnmc,lNMC,"wap",itau_w,
38     $wwriteSTD3)
39c
40       IF (1.EQ.0) THEN
41c
42       DO k=1, nlevSTD
43        DO i=1, klon
44         IF(tnondef(i,k,3).NE.missing_val) THEN
45          zx_tmp_fiNC(i,k) = (100.*tnondef(i,k,3))/freq_moyNMC(3)
46         ELSE
47          zx_tmp_fiNC(i,k) = missing_val
48         ENDIF
49        ENDDO
50       ENDDO !k=1, nlevSTD
51c
52       CALL histwrite_phy(nid_hfnmc,lNMC,"psbg",itau_w,
53     $zx_tmp_fiNC)
54c
55       CALL histwrite_phy(nid_hfnmc,lNMC,"uv",itau_w,
56     $uvsumSTD(:,:,3))
57c
58       CALL histwrite_phy(nid_hfnmc,lNMC,"vq",itau_w,
59     $vqsumSTD(:,:,3))
60c
61       CALL histwrite_phy(nid_hfnmc,lNMC,"vT",itau_w,
62     $vTsumSTD(:,:,3))
63c
64       CALL histwrite_phy(nid_hfnmc,lNMC,"wq",itau_w,
65     $wqsumSTD(:,:,3))
66c
67      CALL histwrite_phy(nid_hfnmc,lNMC,"vphi",itau_w,
68     $vphisumSTD(:,:,3))
69c
70       CALL histwrite_phy(nid_hfnmc,lNMC,"wT",itau_w,
71     $wTsumSTD(:,:,3))
72c
73       CALL histwrite_phy(nid_hfnmc,lNMC,"uxu",itau_w,
74     $u2sumSTD(:,:,3))
75c
76       CALL histwrite_phy(nid_hfnmc,lNMC,"vxv",itau_w,
77     $v2sumSTD(:,:,3))
78c
79       CALL histwrite_phy(nid_hfnmc,lNMC,"TxT",itau_w,
80     $T2sumSTD(:,:,3))
81c
82       ENDIF !(1.EQ.0) THEN
83c
84      if (ok_sync) then
85c$OMP MASTER
86        call histsync(nid_hfnmc)
87c$OMP END MASTER
88      endif
89c
90      ENDIF !      (ok_histNMC(3)) THEN
Note: See TracBrowser for help on using the repository browser.