Changeset 2735 for trunk/LMDZ.GENERIC/libf/phystd/xios_output_mod.F90

Timestamp:

Jun 21, 2022, 11:24:14 AM (2 years ago)

Author:

aslmd

Message:

Created the spectral axes for XspecIR and VI

File:

• trunk/LMDZ.GENERIC/libf/phystd/xios_output_mod.F90 (modified) (6 diffs)

trunk/LMDZ.GENERIC/libf/phystd/xios_output_mod.F90

@@ -11 +11 @@
 
  INTERFACE send_xios_field
-    MODULE PROCEDURE histwrite0d_xios,histwrite2d_xios,histwrite3d_xios
+    MODULE PROCEDURE histwrite0d_xios,histwrite2d_xios,histwrite3d_xios!,histwrite1d_xios
  END INTERFACE
  
@@ -18 +18 @@
 
   SUBROUTINE initialize_xios_output(day,timeofday,dtphys,daysec,&
-                                    yearday,presnivs,pseudoalt)
+                                    yearday,presnivs,pseudoalt,wnoi,wnov)
 !  USE mod_phys_lmdz_para, only: gather, bcast, &
 !                                jj_nb, jj_begin, jj_end, ii_begin, ii_end, &
@@ -44 +44 @@
   REAL,INTENT(IN) :: presnivs(:) ! vertical grid approximate pressure (Pa)
   REAL,INTENT(IN) :: pseudoalt(:) ! vertical grid approximate altitude (km)
-  
+  REAL,INTENT(IN) :: wnoi(:) ! Array of wavenumbers at the spectral interval centers for the infrared.
+  real,intent(in) :: wnov (:) !Array of wavenumbers at the spectral interval centers for the visible.
   
   INTEGER :: data_ibegin, data_iend
@@ -61 +62 @@
     CALL xios_set_axis_attr("altitude", n_glo=size(pseudoalt), value=pseudoalt,&
                             unit="km",positive="up")
-    
+    if (prt_level >=10) write(lunout,*) "initialize_xios_output: call xios_set_axis_attr for IR_Wavenumber"
+    write(lunout,*) "writing IR_Wavenumber now in initialize_xios_output"
+    call xios_set_axis_attr("IR_Wavenumber",n_glo=size(wnoi),value=wnoi, &
+                            unit="cm^-1",positive="up")
+    if (prt_level >=10) write(lunout,*) "initialize_xios_output: call xios_set_axis_attr for VI_Wavenumber"
+    write(lunout,*) "writing VI_Wavenumber now in initialize_xios_output"
+    call xios_set_axis_attr("VI_Wavenumber",n_glo=size(wnov),value=wnov, &
+                            unit="cm^-1",positive="up")
     ! 2. Declare horizontal domain
     ! Set values for the mask:
@@ -213 +221 @@
   USE print_control_mod, ONLY: prt_level, lunout
   USE mod_grid_phy_lmdz, ONLY: nbp_lon
+  USE radinc_h ,only: L_NSPECTI,L_NSPECTV
   IMPLICIT NONE
 
@@ -222 +231 @@
 
     IF (prt_level >= 10) WRITE(lunout,*)'Begin histrwrite2d_xios ',trim(field_name)
-
+    if ((size(field) .eq. L_NSPECTI) .or. (size(field) .eq. L_NSPECTV)) then 
+!$OMP MASTER 
+      ! only for spectral stuff: IR_Bandwidth and VI_Bandwidth
+      call xios_send_field(field_name,field)
+!$OMP END MASTER
+      return
+    endif
     IF (SIZE(field)/=klon) CALL abort_physic('iophy::histwrite2d_xios','Field first DIMENSION not equal to klon',1)