source: LMDZ4/branches/LMDZ4_AR5/libf/phylmd/write_histmthNMC.h @ 5759

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

Ajouts CFMIP2/CMIP5

  • 6eme fichier de sortie "stations" histstn.nc qui necessite 2 fichiers: PARAM/npCFMIP_param.data contenant le nombre de points (120 pour simulations AMIP, 73 pour aqua) PARAM/pointlocations.txt contenat 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

IM
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.8 KB
RevLine 
[524]1!
2! $Header$
3!
[1374]4      IF (ok_histNMC(1)) THEN
[524]5c
[644]6       ndex3d = 0
7       itau_w = itau_phy + itap
8ccc
9c  Champs interpolles sur des niveaux de pression du NMC
[524]10c
[1352]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
[644]15ccc
[1534]16      CALL histwrite_phy(nid_mthnmc,lNMC,"tnondef",itau_w,
17     $tnondef(:,:,1))
[524]18c
[1534]19      CALL histwrite_phy(nid_mthnmc,lNMC,"ta",itau_w,
20     $twriteSTD(:,:,1))
[524]21c
[1534]22      CALL histwrite_phy(nid_mthnmc,lNMC,"zg",itau_w,
23     $phiwriteSTD(:,:,1))
[524]24c
[1534]25       CALL histwrite_phy(nid_mthnmc,lNMC,"hus",itau_w,
26     $qwriteSTD(:,:,1))
[524]27c
[1534]28       CALL histwrite_phy(nid_mthnmc,lNMC,"hur",itau_w,
29     $rhwriteSTD(:,:,1))
[644]30c
[1534]31       CALL histwrite_phy(nid_mthnmc,lNMC,"ua",itau_w,
32     $uwriteSTD(:,:,1))
[644]33c
[1534]34       CALL histwrite_phy(nid_mthnmc,lNMC,"va",itau_w,
35     $vwriteSTD(:,:,1))
[644]36c
[1534]37       CALL histwrite_phy(nid_mthnmc,lNMC,"wap",itau_w,
38     $wwriteSTD(:,:,1))
[644]39c
40       DO k=1, nlevSTD
41        DO i=1, klon
[1352]42         IF(tnondef(i,k,1).NE.missing_val) THEN
[1398]43          zx_tmp_fiNC(i,k) = (100.*tnondef(i,k,1))/freq_moyNMC(1)
[644]44         ELSE
[1352]45          zx_tmp_fiNC(i,k) = missing_val
[524]46         ENDIF
[644]47        ENDDO
48       ENDDO !k=1, nlevSTD
[524]49c
[1534]50      CALL histwrite_phy(nid_mthnmc,lNMC,"psbg",itau_w,
51     $zx_tmp_fiNC)
[524]52c
[1534]53      CALL histwrite_phy(nid_mthnmc,lNMC,"uv",itau_w,
54     $uvsumSTD(:,:,1))
[524]55c
[1534]56      CALL histwrite_phy(nid_mthnmc,lNMC,"vq",itau_w,
57     $vqsumSTD(:,:,1))
[524]58c
[1534]59      CALL histwrite_phy(nid_mthnmc,lNMC,"vT",itau_w,
60     $vTsumSTD(:,:,1))
[524]61c
[1534]62      CALL histwrite_phy(nid_mthnmc,lNMC,"wq",itau_w,
63     $wqsumSTD(:,:,1))
[524]64c
[1534]65      CALL histwrite_phy(nid_mthnmc,lNMC,"vphi",itau_w,
66     $vphisumSTD(:,:,1))
[524]67c
[1534]68       CALL histwrite_phy(nid_mthnmc,lNMC,"wT",itau_w,
69     $wTsumSTD(:,:,1))
[524]70c
[1534]71       CALL histwrite_phy(nid_mthnmc,lNMC,"uxu",itau_w,
72     $u2sumSTD(:,:,1))
[644]73c
[1534]74       CALL histwrite_phy(nid_mthnmc,lNMC,"vxv",itau_w,
75     $v2sumSTD(:,:,1))
[644]76c
[1534]77       CALL histwrite_phy(nid_mthnmc,lNMC,"TxT",itau_w,
78     $T2sumSTD(:,:,1))
[644]79c
[1398]80       DO k=1, nlevSTD
81        DO i=1, klon
82         IF(O3sumSTD(i,k,1).NE.missing_val) THEN
83          zx_tmp_fiNC(i,k) = O3sumSTD(i,k,1) * 1.e+9
84         ELSE
85          zx_tmp_fiNC(i,k) = missing_val
86         ENDIF
87        ENDDO
88       ENDDO !k=1, nlevSTD
[1534]89       CALL histwrite_phy(nid_mthnmc,lNMC,"tro3",itau_w,
[1398]90     $ zx_tmp_fiNC)
91c
92       if (read_climoz == 2) THEN
93       DO k=1, nlevSTD
94        DO i=1, klon
95         IF(O3daysumSTD(i,k,1).NE.missing_val) THEN
96          zx_tmp_fiNC(i,k) = O3daysumSTD(i,k,1) * 1.e+9
97         ELSE
98          zx_tmp_fiNC(i,k) = missing_val
99         ENDIF
100        ENDDO
101       ENDDO !k=1, nlevSTD
102c
[1534]103        CALL histwrite_phy(nid_mthnmc,lNMC,"tro3_daylight",
104     $itau_w, zx_tmp_fiNC)
[1398]105       endif
106c
[524]107      if (ok_sync) then
[766]108c$OMP MASTER
[1352]109        call histsync(nid_mthnmc)
[766]110c$OMP END MASTER
[524]111      endif
[1374]112c
113      ENDIF !(ok_histNMC(1)) THEN
Note: See TracBrowser for help on using the repository browser.