Changeset 2735 for trunk/LMDZ.GENERIC
- Timestamp:
- Jun 21, 2022, 11:24:14 AM (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.GENERIC/libf/phystd/xios_output_mod.F90
r2632 r2735 11 11 12 12 INTERFACE send_xios_field 13 MODULE PROCEDURE histwrite0d_xios,histwrite2d_xios,histwrite3d_xios 13 MODULE PROCEDURE histwrite0d_xios,histwrite2d_xios,histwrite3d_xios!,histwrite1d_xios 14 14 END INTERFACE 15 15 … … 18 18 19 19 SUBROUTINE initialize_xios_output(day,timeofday,dtphys,daysec,& 20 yearday,presnivs,pseudoalt )20 yearday,presnivs,pseudoalt,wnoi,wnov) 21 21 ! USE mod_phys_lmdz_para, only: gather, bcast, & 22 22 ! jj_nb, jj_begin, jj_end, ii_begin, ii_end, & … … 44 44 REAL,INTENT(IN) :: presnivs(:) ! vertical grid approximate pressure (Pa) 45 45 REAL,INTENT(IN) :: pseudoalt(:) ! vertical grid approximate altitude (km) 46 46 REAL,INTENT(IN) :: wnoi(:) ! Array of wavenumbers at the spectral interval centers for the infrared. 47 real,intent(in) :: wnov (:) !Array of wavenumbers at the spectral interval centers for the visible. 47 48 48 49 INTEGER :: data_ibegin, data_iend … … 61 62 CALL xios_set_axis_attr("altitude", n_glo=size(pseudoalt), value=pseudoalt,& 62 63 unit="km",positive="up") 63 64 if (prt_level >=10) write(lunout,*) "initialize_xios_output: call xios_set_axis_attr for IR_Wavenumber" 65 write(lunout,*) "writing IR_Wavenumber now in initialize_xios_output" 66 call xios_set_axis_attr("IR_Wavenumber",n_glo=size(wnoi),value=wnoi, & 67 unit="cm^-1",positive="up") 68 if (prt_level >=10) write(lunout,*) "initialize_xios_output: call xios_set_axis_attr for VI_Wavenumber" 69 write(lunout,*) "writing VI_Wavenumber now in initialize_xios_output" 70 call xios_set_axis_attr("VI_Wavenumber",n_glo=size(wnov),value=wnov, & 71 unit="cm^-1",positive="up") 64 72 ! 2. Declare horizontal domain 65 73 ! Set values for the mask: … … 213 221 USE print_control_mod, ONLY: prt_level, lunout 214 222 USE mod_grid_phy_lmdz, ONLY: nbp_lon 223 USE radinc_h ,only: L_NSPECTI,L_NSPECTV 215 224 IMPLICIT NONE 216 225 … … 222 231 223 232 IF (prt_level >= 10) WRITE(lunout,*)'Begin histrwrite2d_xios ',trim(field_name) 224 233 if ((size(field) .eq. L_NSPECTI) .or. (size(field) .eq. L_NSPECTV)) then 234 !$OMP MASTER 235 ! only for spectral stuff: IR_Bandwidth and VI_Bandwidth 236 call xios_send_field(field_name,field) 237 !$OMP END MASTER 238 return 239 endif 225 240 IF (SIZE(field)/=klon) CALL abort_physic('iophy::histwrite2d_xios','Field first DIMENSION not equal to klon',1) 226 241
Note: See TracChangeset
for help on using the changeset viewer.