source: LMDZ4/branches/LMDZ4_AR5/libf/phylmd/write_histdayNMC.h @ 1595

Last change on this file since 1595 was 1534, checked in by musat, 13 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:executable set to *
File size: 1.7 KB
Line 
1!
2! $Header$
3!
4      IF (ok_histNMC(2)) 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      IF(lev_histdayNMC.EQ.nlevSTD) THEN
17       CALL histwrite_phy(nid_daynmc,lNMC,"tnondef",itau_w,
18     $tnondef(:,:,2))
19       CALL histwrite_phy(nid_daynmc,lNMC,"ta",itau_w,
20     $twriteSTD(:,:,2))
21       CALL histwrite_phy(nid_daynmc,lNMC,"zg",itau_w,
22     $phiwriteSTD(:,:,2))
23       CALL histwrite_phy(nid_daynmc,lNMC,"hus",itau_w,
24     $qwriteSTD(:,:,2))
25       CALL histwrite_phy(nid_daynmc,lNMC,"hur",itau_w,
26     $rhwriteSTD(:,:,2))
27       CALL histwrite_phy(nid_daynmc,lNMC,"ua",itau_w,
28     $uwriteSTD(:,:,2))
29       CALL histwrite_phy(nid_daynmc,lNMC,"va",itau_w,
30     $vwriteSTD(:,:,2))
31       CALL histwrite_phy(nid_daynmc,lNMC,"wap",itau_w,
32     $wwriteSTD(:,:,2))
33      ELSE IF(lev_histdayNMC.EQ.nlevSTD8) THEN
34       CALL histwrite_phy(nid_daynmc,lNMC,"tnondef",itau_w,
35     $tnondefSTD8)
36       CALL histwrite_phy(nid_daynmc,lNMC,"ta",itau_w,
37     $twriteSTD8)
38       CALL histwrite_phy(nid_daynmc,lNMC,"zg",itau_w,
39     $phiwriteSTD8)
40       CALL histwrite_phy(nid_daynmc,lNMC,"hus",itau_w,
41     $qwriteSTD8)
42       CALL histwrite_phy(nid_daynmc,lNMC,"hur",itau_w,
43     $rhwriteSTD8)
44       CALL histwrite_phy(nid_daynmc,lNMC,"ua",itau_w,
45     $uwriteSTD8)
46       CALL histwrite_phy(nid_daynmc,lNMC,"va",itau_w,
47     $vwriteSTD8)
48       CALL histwrite_phy(nid_daynmc,lNMC,"wap",itau_w,
49     $wwriteSTD8)
50      ENDIF
51c
52      if (ok_sync) then
53c$OMP MASTER
54        call histsync(nid_daynmc)
55c$OMP END MASTER
56      endif
57c
58      ENDIF ! (ok_histNMC(2)) THEN
Note: See TracBrowser for help on using the repository browser.