Changeset 4244 for dynamico_lmdz

Timestamp:

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

Author:

dubos

Message:

simple_physics : output LW fluxes

Location:

dynamico_lmdz/simple_physics

Files:

• config/DYNAMICO/TEST/field_def_physics.xml (modified) (1 diff)
• config/DYNAMICO/TEST/file_def_physics.xml (modified) (1 diff)
• phyparam/physics/radiative_lw.F90 (modified) (5 diffs)
• phyparam/physics/radiative_mod.F90 (modified) (1 diff)

dynamico_lmdz/simple_physics/config/DYNAMICO/TEST/field_def_physics.xml

@@ -39 +39 @@
         <field id="phyparam_swflux_down" name="swflux_down" long_name="Downward SW radiative flux" unit="W/m2" />
         <field id="phyparam_swflux_up"   name="swflux_up"   long_name="Upward SW radiative flux"   unit="W/m2" />
+        <field id="phyparam_lwflux_down" name="lwflux_down" long_name="Downward LW radiative flux" unit="W/m2" />
+        <field id="phyparam_lwflux_up"   name="lwflux_up"   long_name="Upward LW radiative flux"   unit="W/m2" />
       </field_group>
       

dynamico_lmdz/simple_physics/config/DYNAMICO/TEST/file_def_physics.xml

@@ -20 +20 @@
       <field field_ref="phyparam_swflux_down" />
       <field field_ref="phyparam_swflux_up" />
+      <field field_ref="phyparam_lwflux_down" />
+      <field field_ref="phyparam_lwflux_up" />
 
     </field_group>

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
     !-----------------------------------------------------------------------
 

dynamico_lmdz/simple_physics/phyparam/physics/radiative_mod.F90

@@ -95 +95 @@
          &             pplev,ps_rad,tsurf,pt, &
          &             zfluxlw,zdtlw,         &
-         &             lverbose)
+         &             lverbose, lwrite)
 
     !    2.4 surface fluxes