Changeset 4244 for dynamico_lmdz/simple_physics/phyparam/physics/radiative_lw.F90

Timestamp:

Jun 13, 2020, 1:59:26 PM (4 years ago)

Author:

dubos

Message:

simple_physics : output LW fluxes

File:

• dynamico_lmdz/simple_physics/phyparam/physics/radiative_lw.F90 (modified) (5 diffs)

dynamico_lmdz/simple_physics/phyparam/physics/radiative_lw.F90

@@ -18 +18 @@
        pp,ps_rad,ptsurf,pt,              &
        pfluxir,pdtlw,                    &
-       lwrite)
+       lverbose, lwrite)
     USE phys_const, ONLY : cpp, g
+    USE writefield_mod, ONLY : writefield
     !=======================================================================
     !
@@ -56 +57 @@
     REAL,    INTENT(IN)  ::  ptsurf(ngrid),pt(ngrid,nlayer),pp(ngrid,nlayer+1)
     REAL,    INTENT(OUT) ::  pdtlw(ngrid,nlayer),pfluxir(ngrid)
-    LOGICAL, INTENT(IN)  :: lwrite
+    LOGICAL, INTENT(IN)  :: lwrite, lverbose
 
     !   variables locales:
@@ -86 +87 @@
           ENDDO
        ENDDO
-       IF(lwrite) THEN
+       IF(lverbose) THEN
           DO ilev=1,nlayer
              WRITELOG(*,*) ' up(',ilev,')  =  ',zup(ngrid/2+1,ilev)
@@ -205 +206 @@
     !   ------------------------
 
-    IF (lwrite) THEN
+    IF (lverbose) THEN
        WRITELOG(*,*) 'Diagnostique rayonnement thermique'
        WRITELOG(*,*) 'temperature     ', &
@@ -219 +220 @@
     ENDIF
 
+    IF(lwrite) THEN
+       CALL writefield('lwflux_up', 'Upward LW flux', 'W/m2', zfluxup)
+       CALL writefield('lwflux_down', 'Downward LW flux', 'W/m2', zfluxdn)
+    END IF
     !-----------------------------------------------------------------------