Changeset 2328

Timestamp:

Jul 16, 2015, 2:09:53 PM (10 years ago)

Author:

Laurent Fairhead

Message:

Special option for Dice 1D case with a new value of iflag_radia=2.
It allows to force LWUP computation via ground temperature
MPL

Location:

LMDZ5/trunk/libf/phylmd

Files:

• dyn1d/compar1d.h (moved) (moved from LMDZ5/trunk/libf/phylmd/compar1d.h)
• dyn1d/flux_arp.h (moved) (moved from LMDZ5/trunk/libf/phylmd/flux_arp.h)
• dyn1d/perturb_radlwsw.F90 (added)
• physiq.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -578 +578 @@
   !  REAL zxsnow(klon)
   REAL zxsnow_dummy(klon)
+  REAL zsav_tsol(klon)
   !
   REAL dist, rmu0(klon), fract(klon)
@@ -3527 +3528 @@
      endif
 
+!
+!sonia :   If Iflag_radia >=2, pertubation of some variables input to radiation
+!(DICE)
+!
+      IF (iflag_radia .ge. 2) THEN
+        zsav_tsol (:) = zxtsol(:)
+        call perturb_radlwsw(zxtsol,iflag_radia)
+      ENDIF
+
      IF (aerosol_couple) THEN 
 #ifdef INCA
@@ -3660 +3670 @@
      ENDIF ! aerosol_couple
      itaprad = 0
+!
+!  If Iflag_radia >=2, reset pertubed variables
+!
+      IF (iflag_radia .ge. 2) THEN
+        zxtsol(:) = zsav_tsol (:)
+      ENDIF
   ENDIF ! MOD(itaprad,radpas)
   itaprad = itaprad + 1