Changeset 2989 for LMDZ5/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Sep 18, 2017, 5:32:12 PM (7 years ago)

Author:

musat

Message:

The 4co2 output fields are now managed by the XIOS server via
file*lmdz*xml files and the ok_4xCO2atm logical flag.
ok_4xCO2atm is initialized to TRUE if

• one of the 4co2 variables (rsut4cO2, rlut4co2, etc) are asked for

by XIOS sefrver or

• if one of the perturbed GES concentration (RCO2_per, RCH4_per, etc)

is different from the actual (RCO2, RCH4, etc) GES concentration.

File:

• LMDZ5/trunk/libf/phylmd/physiq_mod.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -1610 +1610 @@
          IF (xios_field_is_active("dryod550_"//name_aero_tau(naero))) dryaod_diag=.TRUE. 
        ENDDO
+       !
+       !--setting up ok_4xCO2atm to TRUE in XIOS case 
+       IF (xios_field_is_active("rsut4co2").OR.xios_field_is_active("rlut4co2").OR. & 
+           xios_field_is_active("rsutcs4co2").OR.xios_field_is_active("rlutcs4co2").OR. &
+           xios_field_is_active("rsu4co2").OR.xios_field_is_active("rsucs4co2").OR. &
+           xios_field_is_active("rsd4co2").OR.xios_field_is_active("rsdcs4co2").OR. &
+           xios_field_is_active("rlu4co2").OR.xios_field_is_active("rlucs4co2").OR. &
+           xios_field_is_active("rld4co2").OR.xios_field_is_active("rldcs4co2")) &
+           ok_4xCO2atm=.TRUE. 
 #endif 
 
@@ -3725 +3734 @@
           IF (debut) swaero_diag = .FALSE.
           IF (debut) dryaod_diag = .FALSE.
+          !--IM 15/09/2017 here we return ok_4xCO2atm to FALSE
+          !--as for swaero_diag, see above
+          IF (debut) ok_4xCO2atm = .FALSE.
 #endif
           !
@@ -3731 +3743 @@
           !IM Par defaut on a les taux perturbes egaux aux taux actuels
           !
+          IF (RCO2_per.NE.RCO2_act.OR. &
+              RCH4_per.NE.RCH4_act.OR. &
+              RN2O_per.NE.RN2O_act.OR. &
+              RCFC11_per.NE.RCFC11_act.OR. &
+              RCFC12_per.NE.RCFC12_act) ok_4xCO2atm =.TRUE. 
+   ! 
           IF (ok_4xCO2atm) THEN
-             IF (RCO2_per.NE.RCO2_act.OR.RCH4_per.NE.RCH4_act.OR.     &
-                 RN2O_per.NE.RN2O_act.OR.RCFC11_per.NE.RCFC11_act.OR. &
-                 RCFC12_per.NE.RCFC12_act) THEN
                 !
                 RCO2 = RCO2_per
@@ -3782 +3797 @@
                      ZLWFT0_i, ZFLDN0, ZFLUP0, &
                      ZSWFT0_i, ZFSDN0, ZFSUP0)
-             endif
           endif
           !