Changeset 699 for LMDZ4/trunk/libf/phylmd/radlwsw.F

Timestamp:

May 9, 2006, 11:46:14 AM (18 years ago)

Author:

Laurent Fairhead

Message:

En fait le bug sur l'ozone n'en etait pas un! on revient a la situation precedente MPL/JLD
LF

File:

• LMDZ4/trunk/libf/phylmd/radlwsw.F (modified) (2 diffs)

LMDZ4/trunk/libf/phylmd/radlwsw.F

@@ -238 +238 @@
 c wo:    cm.atm (epaisseur en cm dans la situation standard)
 c POZON: kg/kg
-         IF (bug_ozone) then
-           POZON(i,k) = MAX(wo(iof+i,k),1.0e-12)*RG/46.6968
+         POZON(i,k) = MAX(wo(iof+i,k),1.0e-12)*RG/46.6968
      .               /(paprs(iof+i,k)-paprs(iof+i,k+1))
      .               *(paprs(iof+i,1)/101325.0)
-         ELSE
-c le calcul qui suit est maintenant fait dans ozonecm (MPL)
-           POZON(i,k) = wo(i,k)
-         ENDIF
          PCLDLD(i,k) = cldfra(iof+i,k)*cldemi(iof+i,k)
          PCLDLU(i,k) = cldfra(iof+i,k)*cldemi(iof+i,k)
@@ -555 +550 @@
       DO JL = 1, KDLON
          ZCLDSW0(JL,JK) = 0.0
-         IF (bug_ozone) then
-           ZOZ(JL,JK) = POZON(JL,JK)*46.6968/RG
+         ZOZ(JL,JK) = POZON(JL,JK)*46.6968/RG
      .               *PDP(JL,JK)*(101325.0/PPSOL(JL))
-         ELSE
-c        Correction MPL 100505
-           ZOZ(JL,JK) = POZON(JL,JK)*RMD/RMO3*46.6968/RG*PDP(JL,JK)
-         ENDIF           
       ENDDO
       ENDDO