Changeset 1072 for trunk/LMDZ.TITAN/libf

Timestamp:

Oct 15, 2013, 1:10:02 PM (11 years ago)

Author:

slebonnois

Message:

SL: toujours des corrections oubliees pour 1056

Location:

trunk/LMDZ.TITAN/libf/phytitan

Files:

• gasses.F (modified) (1 diff)
• physiq.F (modified) (2 diffs)
• phytrac.F (modified) (6 diffs)

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

@@ -42 +42 @@
   9   FORMAT(///' ALT   CH4      C2H2        C2H6: MASS MIXING RATIOS')
       DO 103 J=1,NLAYER
-      WRITE (6,10)Z(J),GAS1(J),C2H2(J),C2H6(J)
+      WRITE (6,*)Z(J),GAS1(J),C2H2(J),C2H6(J)
+c     WRITE (6,10)Z(J),GAS1(J),C2H2(J),C2H6(J)
 103   CONTINUE
  10   FORMAT(1X,F6.2,1P3E9.1)

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

@@ -1126 +1126 @@
          tr_seri(:,:,1:nmicro) = tr_seri(:,:,1:nmicro)
      .                        + d_tr_mph(:,:,1:nmicro)*dtime
-c        call WriteField_phy('physiq_d_tr_mph01',
-c    .                          d_tr_mph(:,:,1),klev)
-c        call WriteField_phy('physiq_d_tr_mph10',
-c    .                          d_tr_mph(:,:,10),klev)
-        endif
+
 c       PAS ELEGANT mais je n'ai pas trouve d'autres solutions :
 c       Il semblerait qu'il y ait un probleme lorsque les tendances de traceurs
@@ -1144 +1140 @@
           ENDDO
         ENDDO
+
+        endif
 
 c condensation:

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

@@ -44 +44 @@
      &      evapch4,occcld_m,occcld,satch4,satc2h6,satc2h2,rmcloud
       USE moyzon_mod
+      USE write_field_phy
       IMPLICIT none
 #include "dimensions.h"
@@ -184 +185 @@
         DO iq=1,nmicro
          qaer(:,:,iq) = zqfibar(:,:,iq)*delpbar(:,:)/RG
+        DO l=1,klev
+          DO i = 1, klon
+              if (qaer(i,l,iq).lt.0.) then
+        print*,"NEGS ICI ICI !!!!",qaer(i,l,iq),i,l,iq
+                qaer(i,l,iq)=0.
+c          stop
+              endif
+              if (delpbar(i,l).lt.0.) then
+        print*,"NEGS DELP ICI !!!!",i,l,iq,delpbar(i,l)
+           stop
+              endif
+          ENDDO
+        ENDDO
          qaer0(:,:,iq)= tr_seri(:,:,iq)*delp(:,:)/RG
          qaer0bar(:,:,iq) = qaer(:,:,iq)
@@ -455 +469 @@
 c   ---------------------------------
 
+c     call WriteField_phy('phytrac_qaer01',qaer(:,:,1),klev)
+
+c===============================
 c TRACEURS MICROPHYSIQUES
 c                                        
+c===============================
 c ---> pas de microphysique
        IF (microfi.eq.0) THEN
@@ -463 +481 @@
          ENDDO
        ENDIF
+c===============================
 c ---> microphysique 2D
-      IF (microfi.eq.1) THEN
-         DO iq=1,nmicro
-           DO l=1,llm
-             DO i=1,klon
+c     IF (microfi.eq.1) THEN
+c        DO iq=1,nmicro
+c          DO l=1,llm
+c            DO i=1,klon
 c  on repasse le champ de traceurs en 3D (pas les tendances)
 c qaer est ce qui entre dans muphy, donc la moyenne zonale
@@ -481 +500 @@
 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 (qaer0bar(i,l,iq).gt.1e-100) THEN
-                   qaer(i,l,iq) = qaer0(i,l,iq) *
-     &             qaer(i,l,iq)/qaer0bar(i,l,iq)
-               ENDIF
+c              IF (qaer0bar(i,l,iq).gt.1e-100) THEN
+c                  qaer(i,l,iq) = qaer0(i,l,iq) *
+c    &             qaer(i,l,iq)/qaer0bar(i,l,iq)
+c              ENDIF
 c        La tendance correspond a (qaer-qaer0)/ptimestep
-               d_tr_mph(i,l,iq) = (qaer(i,l,iq)-qaer0(i,l,iq))/
+c              d_tr_mph(i,l,iq) = (qaer(i,l,iq)-qaer0(i,l,iq))/
+c    &                            ptimestep
+c            ENDDO
+c          ENDDO
+c        ENDDO
+c ---> microphysique 3D
+c      ELSEIF(microfi.gt.1) THEN 
+c        DO iq=1,nmicro
+c          d_tr_mph(:,:,iq)=(qaer(:,:,iq)-qaer0(:,:,iq))/ptimestep
+c        ENDDO
+c      ENDIF   ! microfi
+
+c      DO iq=1,nmicro
+c  Traceurs microphysiques: passage en intensif: n/m^2 --> n/kg
+c        d_tr_mph(:,:,iq) = d_tr_mph(:,:,iq)*RG/delp(:,:)
+c      ENDDO
+
+c===============================
+c TOUT CE QUI EST AU-DESSUS NE MARCHE PAS: PLEIN DE NEGS...
+c CA MARCHE EN 3D, MAIS PAS EN MOY ZON...
+c===============================
+
+c ---> microphysique 2D
+      IF (microfi.eq.1) THEN
+         DO iq=1,nmicro
+           DO l=1,llm
+             DO i=1,klon
+c ici, qaer correspond a la moy zonale modifiee par la microphys.
+c  Traceurs microphysiques: passage en intensif: n/m^2 --> n/kg
+c en mettant ceci:
+               d_tr_mph(i,l,iq) = (qaer(i,l,iq)*RG/delpbar(i,l)
+     &                            -qaer0(i,l,iq)*RG/delp(i,l))/
      &                            ptimestep
+c on remplace le champ 3D initial par la valeur modifiee de sa moyenne zonale
+c => on remet un champ uniforme en zonal...
              ENDDO
            ENDDO
@@ -495 +547 @@
          DO iq=1,nmicro
            d_tr_mph(:,:,iq)=(qaer(:,:,iq)-qaer0(:,:,iq))/ptimestep
+c  Traceurs microphysiques: passage en intensif: n/m^2 --> n/kg
+           d_tr_mph(:,:,iq) = d_tr_mph(:,:,iq)*RG/delp(:,:)
          ENDDO
        ENDIF   ! microfi
 
-       DO iq=1,nmicro
-c  Traceurs microphysiques: passage en intensif: n/m^2 --> n/kg
-         d_tr_mph(:,:,iq) = d_tr_mph(:,:,iq)*RG/delp(:,:)
-       ENDDO
+c===============================
 
 c AUTRES TRACEURS