Changeset 1095 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_mod.F90

Timestamp:

Feb 13, 2009, 10:52:29 AM (16 years ago)

Author:

jghattas

Message:

Ajoute des variables de sorties, necessaires à l'evaluation du code radiatif, qui avaient disparu avec la reecriture.
Anne C/ Abderrahmane I

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_mod.F90 (modified) (4 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_mod.F90

@@ -318 +318 @@
   type(ctrl_out) :: o_rugs         = ctrl_out((/ 10, 10, 10, 1, 1 /),'rugs')
 
+  type(ctrl_out) :: o_topswad      = ctrl_out((/ 4, 10, 10, 10, 10 /),'topswad')
+  type(ctrl_out) :: o_topswai      = ctrl_out((/ 4, 10, 10, 10, 10 /),'topswai')
+  type(ctrl_out) :: o_solswad      = ctrl_out((/ 4, 10, 10, 10, 10 /),'solswad')
+  type(ctrl_out) :: o_solswai      = ctrl_out((/ 4, 10, 10, 10, 10 /),'solswai')
 !!!!!!!!!!!!!!!!!!!!!! 3D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 
   type(ctrl_out) :: o_lwcon        = ctrl_out((/ 2, 5, 10, 10, 1 /),'lwcon')
@@ -429 +433 @@
                               ctetaSTD,dtime, presnivs, ok_veget, &
                               ocean, iflag_pbl,ok_mensuel,ok_journe, &
-                              ok_hf,ok_instan,ok_LES)   
+                              ok_hf,ok_instan,ok_LES,ok_ade,ok_aie)   
 
   USE iophy 
@@ -447 +451 @@
   integer                               :: nbteta, nlevSTD, radpas
   logical                               :: ok_mensuel, ok_journe, ok_hf, ok_instan
-  logical                               :: ok_LES
+  logical                               :: ok_LES,ok_ade,ok_aie
   real                                  :: dtime
   integer                               :: idayref
@@ -677 +681 @@
      END DO
 
+ IF (ok_ade) THEN
+  CALL histdef2d(iff,o_topswad%flag,o_topswad%name, "ADE at TOA", "W/m2")
+  CALL histdef2d(iff,o_solswad%flag,o_solswad%name, "ADE at SRF", "W/m2")
+ ENDIF
+
+ IF (ok_aie) THEN
+  CALL histdef2d(iff,o_topswai%flag,o_topswai%name, "AIE at TOA", "W/m2")
+  CALL histdef2d(iff,o_solswai%flag,o_solswai%name, "AIE at SFR", "W/m2")
+ ENDIF
+
+
  CALL histdef2d(iff,o_albs%flag,o_albs%name, "Surface albedo", "-")
  CALL histdef2d(iff,o_albslw%flag,o_albslw%name, "Surface albedo LW", "-")