Changeset 973 for trunk/LMDZ.VENUS/libf/phyvenus/phyetat0.F90

Timestamp:

May 27, 2013, 9:47:37 AM (12 years ago)

Author:

slebonnois

Message:

EM+SL: bug corrections in Venus physics

File:

• trunk/LMDZ.VENUS/libf/phyvenus/phyetat0.F90 (modified) (3 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/phyetat0.F90

@@ -152 +152 @@
       PRINT*,'Rayonnement solaire au sol solsw:', xmin, xmax
 
-! Lecture rayonnement IF au sol:
+! Lecture rayonnement IR au sol:
 
       CALL get_field("sollw",sollw,found)
@@ -166 +166 @@
          xmax = MAX(sollw(i),xmax)
       ENDDO
-      PRINT*,'Rayonnement IF au sol sollw:', xmin, xmax
+      PRINT*,'Rayonnement IR au sol sollw:', xmin, xmax
 
 ! Lecture derive des flux:
@@ -183 +183 @@
       ENDDO
       PRINT*,'Derive des flux fder:', xmin, xmax
+
+! Lecture derive flux IR:
+
+      CALL get_field("dlw",dlw,found)
+      IF (.not.found) THEN
+         PRINT*, 'phyetat0: Le champ <dlw> est absent'
+         PRINT*, 'mis a zero'
+         dlw = 0.
+      ENDIF
+      xmin = 1.0E+20
+      xmax = -1.0E+20
+      DO i = 1, klon
+         xmin = MIN(dlw(i),xmin)
+         xmax = MAX(dlw(i),xmax)
+      ENDDO
+      PRINT*,'Derive flux IR dlw:', xmin, xmax
+
+! Lecture rayonnement IR vers le bas au sol:
+
+      CALL get_field("sollwdown",sollwdown,found)
+      IF (.not.found) THEN
+         PRINT*, 'phyetat0: Le champ <sollwdown> est absent'
+         PRINT*, 'mis a zero'
+         sollwdown = 0.
+      ENDIF
+      xmin = 1.0E+20
+      xmax = -1.0E+20
+      DO i = 1, klon
+         xmin = MIN(sollwdown(i),xmin)
+         xmax = MAX(sollwdown(i),xmax)
+      ENDDO
+      PRINT*,'Flux IR vers le bas au sol sollwdown:', xmin, xmax
 
 ! Lecture du rayonnement net au sol: