Changeset 888 for trunk/LMDZ.TITAN/libf/phytitan

Timestamp:

Feb 18, 2013, 2:56:41 PM (12 years ago)

Author:

slebonnois

Message:

SL: small modifications to the tools, to Venus default .def files and to outputs (including forgotten modifications linked to the 1D); + bug corrections in phytitan

Location:

trunk/LMDZ.TITAN/libf/phytitan

Files:

• cooling.F (modified) (2 diffs)
• ini_histday.h (modified) (1 diff)
• ini_histins.h (modified) (1 diff)
• ini_histmth.h (modified) (1 diff)
• optci_1pt_3.F (modified) (1 diff)
• phytrac.F (modified) (4 diffs)

trunk/LMDZ.TITAN/libf/phytitan/cooling.F

@@ -137 +137 @@
       save RHOP,UBARI,WNOI,DWNI
 
+      REAL effg    ! effg est une fonction(z en m)
+
 c-----------------------------------------------------------------------
 
@@ -321 +323 @@
 C TURN THE Q'S INTO TIMESCALES.....
           DO 3300 ig=1,NG
-             eff_g = RG*(RA/(RA+Z(ig,J)))**2 ! 10% DIFF AT 1 MBAR
-             COLDEN = RHOP*(PRESS(ig,J+1)-PRESS(ig,J))/eff_g
+          COLDEN = RHOP*(PRESS(ig,J+1)-PRESS(ig,J))/effg(Z(ig,J))
 c            Q0(J) = (COLDEN * CSUBP )/Q0(J)
              Q0(ig,J) = Q0(ig,J) / (COLDEN*CSUBP) 

trunk/LMDZ.TITAN/libf/phytitan/ini_histday.h

@@ -11 +11 @@
          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlond,zx_lon)
          DO i = 1, iim
-            zx_lon(i,1) = rlond(i+1)
-            zx_lon(i,jjmp1) = rlond(i+1)
+            zx_lon(i,1) = rlond(i+jjmp1-jjm)
+            zx_lon(i,jjmp1) = rlond(i+jjmp1-jjm)
          ENDDO
          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlatd,zx_lat)

trunk/LMDZ.TITAN/libf/phytitan/ini_histins.h

@@ -10 +10 @@
          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlond,zx_lon)
          DO i = 1, iim
-            zx_lon(i,1) = rlond(i+1)
-            zx_lon(i,jjmp1) = rlond(i+1)
+            zx_lon(i,1) = rlond(i+jjmp1-jjm)
+            zx_lon(i,jjmp1) = rlond(i+jjmp1-jjm)
          ENDDO
          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlatd,zx_lat)

trunk/LMDZ.TITAN/libf/phytitan/ini_histmth.h

@@ -11 +11 @@
          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlond,zx_lon)
          DO i = 1, iim
-            zx_lon(i,1) = rlond(i+1)
-            zx_lon(i,jjmp1) = rlond(i+1)
+            zx_lon(i,1) = rlond(i+jjmp1-jjm)
+            zx_lon(i,jjmp1) = rlond(i+jjmp1-jjm)
          ENDDO
          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlatd,zx_lat)

trunk/LMDZ.TITAN/libf/phytitan/optci_1pt_3.F

@@ -16 +16 @@
      &     pi, sigma, prod,reali,fhvis
 
-      integer k, j,inq
+      integer k, j,inq,kgas
 
       real tbar, pbar, bmu, coef1, effg, taeros, taeroscat, cbar, 
      &     qext, qsct, qabs, qbar, xmono, xrule, deltaz, tnuext, 
      &     tnuscat, cnbar, qextc, qsctc, qabsc, qbarc, taugas, pnn, 
-     &     pcc, pcn, phn, kgas, u, ig, tau2, tlimit
+     &     pcc, pcn, phn, u, ig, tau2, tlimit
 
       PARAMETER(NLAYER=llm,NLEVEL=NLAYER+1)

trunk/LMDZ.TITAN/libf/phytitan/phytrac.F

@@ -712 +712 @@
 c    On a quand meme le droit de produire des traceurs dans la cellule.
 c    On considere donc que la valeur de sortie 3D correspond a la valeur de sortie 2D.
-c    Cela permet aussi entre autre d'eviter les NaN pour les traceurs des nuages !
+c    Cela permet aussi entre autre d eviter les NaN pour les traceurs des nuages !
 c    (au dessus de la tropo pas de nuages donc qaer(nrad+1:ntype*nrad) = 0 !!!)
                  IF (zqaer0(j,l,iq).lt.1e-100) THEN
@@ -788 +788 @@
 c----------------------
 c La microphysique avec nuages doit se faire obligatoirement en 3D.  (FAUX ACTUELLEMENT)
-c Rien n'empeche de faire la chimie en 2D. Cependant pour prendre en compte la 
+c Rien n empeche de faire la chimie en 2D. Cependant pour prendre en compte la 
 c condensation due a la microfi (en 3D) on recalcule la tendance finale pour
 c les especes concernees (CH4, C2H6 pour le moment). 
@@ -852 +852 @@
 
 c--------------------------------------------------
-c CALCUL DU FLUX DE CHALEUR LATENTE D'EVAPORATION 
+c CALCUL DU FLUX DE CHALEUR LATENTE D EVAPORATION 
 c DU METHANE 
 c--------------------------------------------------
@@ -941 +941 @@
 c      OCCCLD
 c      Calcul le nombre d'occurence d'un nuage 
-c      d'opacité comprise en kmin et kmax
+c      d opacité comprise en kmin et kmax
 c          k        kmin            kmax
 c          1   0.0000000      0.10000000