Changeset 2645

Timestamp:

Oct 4, 2016, 5:49:13 PM (8 years ago)

Author:

musat

Message:

Add histwrite for scalar outputs in LMDZ.
The scalars (GES, rsun*, solaire and Earths' orbital params)
outputed before in paramLMDZ_phy.nc are now in
histmth.nc file.

Location:

LMDZ5/trunk/libf/phylmd

Files:

• iophy.F90 (modified) (2 diffs)
• phys_output_write_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/iophy.F90

@@ -17 +17 @@
 #ifdef CPP_XIOS
   INTERFACE histwrite_phy
-    MODULE PROCEDURE histwrite2d_phy,histwrite3d_phy,histwrite2d_phy_old,histwrite3d_phy_old,histwrite2d_xios,histwrite3d_xios
+    MODULE PROCEDURE histwrite2d_phy,histwrite3d_phy,histwrite2d_phy_old,histwrite3d_phy_old,histwrite2d_xios,histwrite3d_xios,histwrite0d_xios
   END INTERFACE
 #else
@@ -1331 +1331 @@
   IF (prt_level >= 10) write(lunout,*)'End histrwrite3d_xios ',field_name
   END SUBROUTINE histwrite3d_xios
+
+  SUBROUTINE histwrite0d_xios(field_name, field)
+  USE xios, only: xios_send_scalar 
+  IMPLICIT NONE
+
+    CHARACTER(LEN=*), INTENT(IN) :: field_name
+    REAL, INTENT(IN) :: field ! --> scalar
+
+!$OMP MASTER
+   CALL xios_send_scalar(field_name, field)
+!$OMP END MASTER
+
+  END SUBROUTINE histwrite0d_xios
 #endif
 end module iophy

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -312 +312 @@
     USE phys_cal_mod, only : mth_len
 
+#ifdef CPP_RRTM
+    USE YOESW, ONLY : RSUN
+#endif
 
     IMPLICIT NONE
@@ -355 +358 @@
 #endif
     REAL, PARAMETER :: un_jour=86400.
+    INTEGER ISW
+    CHARACTER*1 ch1
 
     ! On calcul le nouveau tau:
@@ -403 +408 @@
        CALL histwrite_phy(o_contfracOR, pctsrf(:,is_ter))
        CALL histwrite_phy(o_aireTER, paire_ter)
+!
+#ifdef CPP_XIOS
+       CALL histwrite_phy("R_ecc",R_ecc)
+       CALL histwrite_phy("R_peri",R_peri)
+       CALL histwrite_phy("R_incl",R_incl)
+       CALL histwrite_phy("solaire",solaire)
+!
+#ifdef CPP_RRTM
+      IF (iflag_rrtm.EQ.1) THEN
+        DO ISW=1, NSW
+          WRITE(ch1,'(i1)') ISW
+!         zx_tmp_0d=RSUN(ISW)
+!         CALL histwrite_phy("rsun"//ch1,zx_tmp_0d)
+          CALL histwrite_phy("rsun"//ch1,RSUN(ISW))
+        ENDDO
+      ENDIF
+#endif
+!
+       CALL histwrite_phy("co2_ppm",co2_ppm)
+       CALL histwrite_phy("CH4_ppb",CH4_ppb)
+       CALL histwrite_phy("N2O_ppb",N2O_ppb)
+       CALL histwrite_phy("CFC11_ppt",CFC11_ppt)
+       CALL histwrite_phy("CFC12_ppt",CFC12_ppt)
+!
+#endif
+
 !!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! Simulateur AIRS