Changeset 1534 for LMDZ4/branches/LMDZ4_AR5/libf/phylmd/write_histdayNMC.h

Timestamp:

Jun 3, 2011, 7:28:17 PM (13 years ago)

Author:

musat

Message:

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

File:

• LMDZ4/branches/LMDZ4_AR5/libf/phylmd/write_histdayNMC.h (modified) (1 diff)

LMDZ4/branches/LMDZ4_AR5/libf/phylmd/write_histdayNMC.h

@@ -15 +15 @@
 ccc
       IF(lev_histdayNMC.EQ.nlevSTD) THEN
-       CALL histwrite_phy(nid_daynmc,"tnondef",itau_w,tnondef(:,:,2))
-       CALL histwrite_phy(nid_daynmc,"ta",itau_w,twriteSTD(:,:,2))
-       CALL histwrite_phy(nid_daynmc,"zg",itau_w,phiwriteSTD(:,:,2))
-       CALL histwrite_phy(nid_daynmc,"hus",itau_w,qwriteSTD(:,:,2))
-       CALL histwrite_phy(nid_daynmc,"hur",itau_w,rhwriteSTD(:,:,2))
-       CALL histwrite_phy(nid_daynmc,"ua",itau_w,uwriteSTD(:,:,2))
-       CALL histwrite_phy(nid_daynmc,"va",itau_w,vwriteSTD(:,:,2))
-       CALL histwrite_phy(nid_daynmc,"wap",itau_w,wwriteSTD(:,:,2))
+       CALL histwrite_phy(nid_daynmc,lNMC,"tnondef",itau_w,
+     $tnondef(:,:,2))
+       CALL histwrite_phy(nid_daynmc,lNMC,"ta",itau_w,
+     $twriteSTD(:,:,2))
+       CALL histwrite_phy(nid_daynmc,lNMC,"zg",itau_w,
+     $phiwriteSTD(:,:,2))
+       CALL histwrite_phy(nid_daynmc,lNMC,"hus",itau_w,
+     $qwriteSTD(:,:,2))
+       CALL histwrite_phy(nid_daynmc,lNMC,"hur",itau_w,
+     $rhwriteSTD(:,:,2))
+       CALL histwrite_phy(nid_daynmc,lNMC,"ua",itau_w,
+     $uwriteSTD(:,:,2))
+       CALL histwrite_phy(nid_daynmc,lNMC,"va",itau_w,
+     $vwriteSTD(:,:,2))
+       CALL histwrite_phy(nid_daynmc,lNMC,"wap",itau_w,
+     $wwriteSTD(:,:,2))
       ELSE IF(lev_histdayNMC.EQ.nlevSTD8) THEN
-       CALL histwrite_phy(nid_daynmc,"tnondef",itau_w,tnondefSTD8)
-       CALL histwrite_phy(nid_daynmc,"ta",itau_w,twriteSTD8)
-       CALL histwrite_phy(nid_daynmc,"zg",itau_w,phiwriteSTD8)
-       CALL histwrite_phy(nid_daynmc,"hus",itau_w,qwriteSTD8)
-       CALL histwrite_phy(nid_daynmc,"hur",itau_w,rhwriteSTD8)
-       CALL histwrite_phy(nid_daynmc,"ua",itau_w,uwriteSTD8)
-       CALL histwrite_phy(nid_daynmc,"va",itau_w,vwriteSTD8)
-       CALL histwrite_phy(nid_daynmc,"wap",itau_w,wwriteSTD8)
+       CALL histwrite_phy(nid_daynmc,lNMC,"tnondef",itau_w,
+     $tnondefSTD8)
+       CALL histwrite_phy(nid_daynmc,lNMC,"ta",itau_w,
+     $twriteSTD8)
+       CALL histwrite_phy(nid_daynmc,lNMC,"zg",itau_w,
+     $phiwriteSTD8)
+       CALL histwrite_phy(nid_daynmc,lNMC,"hus",itau_w,
+     $qwriteSTD8)
+       CALL histwrite_phy(nid_daynmc,lNMC,"hur",itau_w,
+     $rhwriteSTD8)
+       CALL histwrite_phy(nid_daynmc,lNMC,"ua",itau_w,
+     $uwriteSTD8)
+       CALL histwrite_phy(nid_daynmc,lNMC,"va",itau_w,
+     $vwriteSTD8)
+       CALL histwrite_phy(nid_daynmc,lNMC,"wap",itau_w,
+     $wwriteSTD8)
       ENDIF
 c