Changeset 1539 for LMDZ5/trunk/libf/phylmd/phystokenc.F90

Timestamp:

Jun 9, 2011, 12:13:33 AM (13 years ago)

Author:

musat

Message:

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

File:

• LMDZ5/trunk/libf/phylmd/phystokenc.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/phystokenc.F90

@@ -125 +125 @@
   LOGICAL ok_sync
   CHARACTER(len=12) :: nvar
+  logical, parameter :: lstokenc=.FALSE.
 !
 !======================================================================
@@ -164 +165 @@
      
      ! Write field phis and aire only once
-     CALL histwrite_phy(physid,"phis",itap,pphis)
-     CALL histwrite_phy(physid,"aire",itap,paire)
-     CALL histwrite_phy(physid,"longitudes",itap,rlon)
-     CALL histwrite_phy(physid,"latitudes",itap,rlat)
+     CALL histwrite_phy(physid,lstokenc,"phis",itap,pphis)
+     CALL histwrite_phy(physid,lstokenc,"aire",itap,paire)
+     CALL histwrite_phy(physid,lstokenc,"longitudes",itap,rlon)
+     CALL histwrite_phy(physid,lstokenc,"latitudes",itap,rlat)
 
   END IF
@@ -288 +289 @@
 
      ! write fields
-     CALL histwrite_phy(physid,"t",itap,t)
-     CALL histwrite_phy(physid,"mfu",itap,mfu)
-     CALL histwrite_phy(physid,"mfd",itap,mfd)
-     CALL histwrite_phy(physid,"en_u",itap,en_u)
-     CALL histwrite_phy(physid,"de_u",itap,de_u)
-     CALL histwrite_phy(physid,"en_d",itap,en_d)
-     CALL histwrite_phy(physid,"de_d",itap,de_d)
-     CALL histwrite_phy(physid,"coefh",itap,coefh)      
-     CALL histwrite_phy(physid,"fm_th",itap,fm_therm)
-     CALL histwrite_phy(physid,"en_th",itap,entr_therm)
-     CALL histwrite_phy(physid,"frac_impa",itap,frac_impa)
-     CALL histwrite_phy(physid,"frac_nucl",itap,frac_nucl)
-     CALL histwrite_phy(physid,"pyu1",itap,pyu1)
-     CALL histwrite_phy(physid,"pyv1",itap,pyv1)
-     CALL histwrite_phy(physid,"ftsol1",itap,pftsol(:,1))
-     CALL histwrite_phy(physid,"ftsol2",itap,pftsol(:,2))
-     CALL histwrite_phy(physid,"ftsol3",itap,pftsol(:,3))
-     CALL histwrite_phy(physid,"ftsol4",itap,pftsol(:,4))
-     CALL histwrite_phy(physid,"psrf1",itap,ppsrf(:,1))
-     CALL histwrite_phy(physid,"psrf2",itap,ppsrf(:,2))
-     CALL histwrite_phy(physid,"psrf3",itap,ppsrf(:,3))
-     CALL histwrite_phy(physid,"psrf4",itap,ppsrf(:,4))
-     CALL histwrite_phy(physid,"sh",itap,sh)
-     CALL histwrite_phy(physid,"da",itap,da)
-     CALL histwrite_phy(physid,"mp",itap,mp)
-     CALL histwrite_phy(physid,"upwd",itap,upwd)
-     CALL histwrite_phy(physid,"dnwd",itap,dnwd)
+     CALL histwrite_phy(physid,lstokenc,"t",itap,t)
+     CALL histwrite_phy(physid,lstokenc,"mfu",itap,mfu)
+     CALL histwrite_phy(physid,lstokenc,"mfd",itap,mfd)
+     CALL histwrite_phy(physid,lstokenc,"en_u",itap,en_u)
+     CALL histwrite_phy(physid,lstokenc,"de_u",itap,de_u)
+     CALL histwrite_phy(physid,lstokenc,"en_d",itap,en_d)
+     CALL histwrite_phy(physid,lstokenc,"de_d",itap,de_d)
+     CALL histwrite_phy(physid,lstokenc,"coefh",itap,coefh)     
+     CALL histwrite_phy(physid,lstokenc,"fm_th",itap,fm_therm)
+     CALL histwrite_phy(physid,lstokenc,"en_th",itap,entr_therm)
+     CALL histwrite_phy(physid,lstokenc,"frac_impa",itap,frac_impa)
+     CALL histwrite_phy(physid,lstokenc,"frac_nucl",itap,frac_nucl)
+     CALL histwrite_phy(physid,lstokenc,"pyu1",itap,pyu1)
+     CALL histwrite_phy(physid,lstokenc,"pyv1",itap,pyv1)
+     CALL histwrite_phy(physid,lstokenc,"ftsol1",itap,pftsol(:,1))
+     CALL histwrite_phy(physid,lstokenc,"ftsol2",itap,pftsol(:,2))
+     CALL histwrite_phy(physid,lstokenc,"ftsol3",itap,pftsol(:,3))
+     CALL histwrite_phy(physid,lstokenc,"ftsol4",itap,pftsol(:,4))
+     CALL histwrite_phy(physid,lstokenc,"psrf1",itap,ppsrf(:,1))
+     CALL histwrite_phy(physid,lstokenc,"psrf2",itap,ppsrf(:,2))
+     CALL histwrite_phy(physid,lstokenc,"psrf3",itap,ppsrf(:,3))
+     CALL histwrite_phy(physid,lstokenc,"psrf4",itap,ppsrf(:,4))
+     CALL histwrite_phy(physid,lstokenc,"sh",itap,sh)
+     CALL histwrite_phy(physid,lstokenc,"da",itap,da)
+     CALL histwrite_phy(physid,lstokenc,"mp",itap,mp)
+     CALL histwrite_phy(physid,lstokenc,"upwd",itap,upwd)
+     CALL histwrite_phy(physid,lstokenc,"dnwd",itap,dnwd)
 
 
@@ -328 +329 @@
         nvar='phi_lev'//trim(nvar)
         
-        CALL histwrite_phy(physid,nvar,itap,phi(:,:,k))
+        CALL histwrite_phy(physid,lstokenc,nvar,itap,phi(:,:,k))
      END DO