Changeset 1534 for LMDZ4/branches/LMDZ4_AR5/libf/phylmd/phystokenc.F

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/phystokenc.F (modified) (4 diffs)

LMDZ4/branches/LMDZ4_AR5/libf/phylmd/phystokenc.F

@@ -123 +123 @@
       logical,save :: first=.true.
 c$OMP THREADPRIVATE(first)
+      logical, parameter :: lstokenc=.FALSE.
 c
 c   Couche limite:
@@ -169 +170 @@
       ndex2d = 0
       ndex3d = 0
-      i=itap 
+      i=itap
 cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d)
-      CALL histwrite_phy(physid,"phis",i,pphis)
+      CALL histwrite_phy(physid,lstokenc,"phis",i,pphis)
 c
       i=itap
 cym      CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d)
-      CALL histwrite_phy(physid,"aire",i,paire)
+      CALL histwrite_phy(physid,lstokenc,"aire",i,paire)
 
       iadvtr=iadvtr+1
@@ -282 +283 @@
 ccccc
 cym         CALL gr_fi_ecrit(klev,klon,iim,jjm+1, t, zx_tmp_3d)
-         CALL histwrite_phy(physid,"t",itap,t)
+         CALL histwrite_phy(physid,lstokenc,"t",itap,t)
 
 cym         CALL gr_fi_ecrit(klev,klon,iim,jjm+1, mfu, zx_tmp_3d)
-      CALL histwrite_phy(physid,"mfu",itap,mfu)
+      CALL histwrite_phy(physid,lstokenc,"mfu",itap,mfu)
 cym     CALL gr_fi_ecrit(klev,klon,iim,jjm+1, mfd, zx_tmp_3d)
-      CALL histwrite_phy(physid,"mfd",itap,mfd)
+      CALL histwrite_phy(physid,lstokenc,"mfd",itap,mfd)
 cym        CALL gr_fi_ecrit(klev,klon,iim,jjm+1, en_u, zx_tmp_3d)
-      CALL histwrite_phy(physid,"en_u",itap,en_u)
+      CALL histwrite_phy(physid,lstokenc,"en_u",itap,en_u)
 cym        CALL gr_fi_ecrit(klev,klon,iim,jjm+1, de_u, zx_tmp_3d)
-      CALL histwrite_phy(physid,"de_u",itap,de_u)
+      CALL histwrite_phy(physid,lstokenc,"de_u",itap,de_u)
 cym        CALL gr_fi_ecrit(klev,klon,iim,jjm+1, en_d, zx_tmp_3d)
-      CALL histwrite_phy(physid,"en_d",itap,en_d)
+      CALL histwrite_phy(physid,lstokenc,"en_d",itap,en_d)
 cym        CALL gr_fi_ecrit(klev,klon,iim,jjm+1, de_d, zx_tmp_3d)       
-      CALL histwrite_phy(physid,"de_d",itap,de_d)
+      CALL histwrite_phy(physid,lstokenc,"de_d",itap,de_d)
 cym        CALL gr_fi_ecrit(klev,klon,iim,jjm+1, coefh, zx_tmp_3d)         
-      CALL histwrite_phy(physid,"coefh",itap,coefh)     
+      CALL histwrite_phy(physid,lstokenc,"coefh",itap,coefh)    
 
 c ajou...
@@ -307 +308 @@
 
 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, fm_therm1, zx_tmp_3d)
-      CALL histwrite_phy(physid,"fm_th",itap,fm_therm1)
+      CALL histwrite_phy(physid,lstokenc,"fm_th",itap,fm_therm1)
 c
 cym      CALL gr_fi_ecrit(klev,klon,iim,jjm+1, entr_therm, zx_tmp_3d)
-      CALL histwrite_phy(physid,"en_th",itap,entr_therm)
+      CALL histwrite_phy(physid,lstokenc,"en_th",itap,entr_therm)
 cccc
 cym       CALL gr_fi_ecrit(klev,klon,iim,jjm+1,frac_impa,zx_tmp_3d)
-        CALL histwrite_phy(physid,"frac_impa",itap,frac_impa)
+        CALL histwrite_phy(physid,lstokenc,"frac_impa",itap,
+     $frac_impa)
 
 cym        CALL gr_fi_ecrit(klev,klon,iim,jjm+1,frac_nucl,zx_tmp_3d)
-        CALL histwrite_phy(physid,"frac_nucl",itap,frac_nucl)
+        CALL histwrite_phy(physid,lstokenc,"frac_nucl",itap,
+     $frac_nucl)
  
 cym        CALL gr_fi_ecrit(1, klon,iim,jjm+1, pyu1,zx_tmp_2d)
-      CALL histwrite_phy(physid,"pyu1",itap,pyu1)
+      CALL histwrite_phy(physid,lstokenc,"pyu1",itap,pyu1)
         
 cym     CALL gr_fi_ecrit(1, klon,iim,jjm+1, pyv1,zx_tmp_2d)
-      CALL histwrite_phy(physid,"pyv1",itap,pyv1)
+      CALL histwrite_phy(physid,lstokenc,"pyv1",itap,pyv1)
         
 cym     CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol1, zx_tmp_2d)
-      CALL histwrite_phy(physid,"ftsol1",itap,pftsol1)
+      CALL histwrite_phy(physid,lstokenc,"ftsol1",itap,pftsol1)
 cym         CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol2, zx_tmp_2d)
-      CALL histwrite_phy(physid,"ftsol2",itap,pftsol2)
+      CALL histwrite_phy(physid,lstokenc,"ftsol2",itap,pftsol2)
 cym          CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol3, zx_tmp_2d)
-      CALL histwrite_phy(physid,"ftsol3",itap,pftsol3)
+      CALL histwrite_phy(physid,lstokenc,"ftsol3",itap,pftsol3)
 cym         CALL gr_fi_ecrit(1,klon,iim,jjm+1, pftsol4, zx_tmp_2d)
-      CALL histwrite_phy(physid,"ftsol4",itap,pftsol4)
+      CALL histwrite_phy(physid,lstokenc,"ftsol4",itap,pftsol4)
 
 cym        CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf1, zx_tmp_2d)
-      CALL histwrite_phy(physid,"psrf1",itap,ppsrf1)
+      CALL histwrite_phy(physid,lstokenc,"psrf1",itap,ppsrf1)
 cym        CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf2, zx_tmp_2d)
-      CALL histwrite_phy(physid,"psrf2",itap,ppsrf2)
+      CALL histwrite_phy(physid,lstokenc,"psrf2",itap,ppsrf2)
 cym        CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf3, zx_tmp_2d)
-      CALL histwrite_phy(physid,"psrf3",itap,ppsrf3)
+      CALL histwrite_phy(physid,lstokenc,"psrf3",itap,ppsrf3)
 cym        CALL gr_fi_ecrit(1,klon,iim,jjm+1, ppsrf4, zx_tmp_2d)
-      CALL histwrite_phy(physid,"psrf4",itap,ppsrf4)
+      CALL histwrite_phy(physid,lstokenc,"psrf4",itap,ppsrf4)
 
 c$OMP MASTER