Changeset 1682 for trunk/LMDZ.GENERIC/libf/phystd
- Timestamp:
- Mar 31, 2017, 11:31:38 AM (8 years ago)
- Location:
- trunk/LMDZ.GENERIC/libf/phystd
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified trunk/LMDZ.GENERIC/libf/phystd/inifis_mod.F90 ¶
r1677 r1682 9 9 prad,pg,pr,pcpp) 10 10 11 use init_print_control_mod, only: init_print_control 11 12 use radinc_h, only: ini_radinc_h, naerkind 12 13 use radcommon_h, only: ini_radcommon_h … … 76 77 REAL SSUM 77 78 79 ! Initialize flags lunout, prt_level, debug (in print_control_mod) 80 CALL init_print_control 81 78 82 ! initialize constants in comcstfi_mod 79 83 rad=prad -
TabularUnified trunk/LMDZ.GENERIC/libf/phystd/ocean_slab_mod.F90 ¶
r1397 r1682 94 94 IF (error /= 0) THEN 95 95 abort_message='Pb allocation tmp_tslab' 96 CALL abort_ gcm(modname,abort_message,1)96 CALL abort_physic(modname,abort_message,1) 97 97 ENDIF 98 98 tmp_tslab(:,:) = tslab_rst(:,:) … … 100 100 IF (error /= 0) THEN 101 101 abort_message='Pb allocation tmp_tslab_loc' 102 CALL abort_ gcm(modname,abort_message,1)102 CALL abort_physic(modname,abort_message,1) 103 103 ENDIF 104 104 tmp_tslab_loc(:,:) = tslab_rst(:,:) … … 107 107 IF (error /= 0) THEN 108 108 abort_message='Pb allocation tmp_seaice' 109 CALL abort_ gcm(modname,abort_message,1)109 CALL abort_physic(modname,abort_message,1) 110 110 ENDIF 111 111 tmp_seaice(:) = seaice_rst(:) … … 114 114 IF (error /= 0) THEN 115 115 abort_message='Pb allocation tmp_pctsrf_slab' 116 CALL abort_ gcm(modname,abort_message,1)116 CALL abort_physic(modname,abort_message,1) 117 117 ENDIF 118 118 tmp_pctsrf_slab(:) = pctsrf_rst(:) … … 122 122 IF (error /= 0) THEN 123 123 abort_message='Pb allocation tmp_radsol' 124 CALL abort_ gcm(modname,abort_message,1)124 CALL abort_physic(modname,abort_message,1) 125 125 ENDIF 126 126 … … 128 128 IF (error /= 0) THEN 129 129 abort_message='Pb allocation tmp_flux_o' 130 CALL abort_ gcm(modname,abort_message,1)130 CALL abort_physic(modname,abort_message,1) 131 131 ENDIF 132 132 … … 134 134 IF (error /= 0) THEN 135 135 abort_message='Pb allocation tmp_flux_g' 136 CALL abort_ gcm(modname,abort_message,1)136 CALL abort_physic(modname,abort_message,1) 137 137 ENDIF 138 138 … … 141 141 IF (error /= 0) THEN 142 142 abort_message='Pb allocation slab_bils' 143 CALL abort_ gcm(modname,abort_message,1)143 CALL abort_physic(modname,abort_message,1) 144 144 ENDIF 145 145 slab_bils(:) = 0.0 … … 148 148 IF (error /= 0) THEN 149 149 abort_message='Pb allocation dt_hdiff' 150 CALL abort_ gcm(modname,abort_message,1)150 CALL abort_physic(modname,abort_message,1) 151 151 ENDIF 152 152 dt_hdiff = 0.0 … … 155 155 IF (error /= 0) THEN 156 156 abort_message='Pb allocation dt_hdiff' 157 CALL abort_ gcm(modname,abort_message,1)157 CALL abort_physic(modname,abort_message,1) 158 158 ENDIF 159 159 dt_ekman = 0.0 … … 163 163 IF (error /= 0) THEN 164 164 abort_message='Pb allocation lmt_bils' 165 CALL abort_ gcm(modname,abort_message,1)165 CALL abort_physic(modname,abort_message,1) 166 166 ENDIF 167 167 lmt_bils(:) = 0.0 … … 170 170 IF (error /= 0) THEN 171 171 abort_message='Pb allocation slabh' 172 CALL abort_ gcm(modname,abort_message,1)172 CALL abort_physic(modname,abort_message,1) 173 173 ENDIF 174 174 slabh(1)=50. -
TabularUnified trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90 ¶
r1669 r1682 48 48 use callkeys_mod 49 49 use vertical_layers_mod, only: presnivs, pseudoalt 50 use mod_phys_lmdz_omp_data, ONLY: is_omp_master 50 51 #ifdef CPP_XIOS 51 52 use xios_output_mod, only: initialize_xios_output, & 52 53 update_xios_timestep, & 53 54 send_xios_field 55 use wxios, only: wxios_context_init, xios_context_finalize 54 56 #endif 55 57 implicit none … … 529 531 endif 530 532 533 #ifdef CPP_XIOS 534 ! Initialize XIOS context 535 write(*,*) "physiq: call wxios_context_init" 536 CALL wxios_context_init 537 #endif 531 538 532 539 ! Read 'startfi.nc' file. … … 735 742 presnivs,pseudoalt) 736 743 #endif 744 write(*,*) "physiq: end of firstcall" 737 745 endif ! end of 'firstcall' 738 746 … … 1893 1901 end if 1894 1902 1895 1896 endif ! end of 'lastcall' 1903 endif ! end of 'lastcall' 1897 1904 1898 1905 … … 2181 2188 CALL send_xios_field("u",zu) 2182 2189 CALL send_xios_field("v",zv) 2183 2190 2191 if (lastcall.and.is_omp_master) then 2192 write(*,*) "physiq: call xios_context_finalize" 2193 call xios_context_finalize 2194 endif 2184 2195 #endif 2185 2196 -
TabularUnified trunk/LMDZ.GENERIC/libf/phystd/xios_output_mod.F90 ¶
r1626 r1682 26 26 mpi_size, mpi_rank, klon_mpi, & 27 27 is_sequential, is_south_pole_dyn 28 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, klon_glo 28 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, klon_glo, grid_type, unstructured 29 29 USE print_control_mod, ONLY: lunout, prt_level 30 30 USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat 31 31 USE regular_lonlat_mod, ONLY: lon_reg, lat_reg 32 32 USE nrtype, ONLY: pi 33 USE wxios 33 #ifdef CPP_XIOS 34 USE xios 35 #endif 36 USE wxios, ONLY: wxios_domain_param, wxios_domain_param_unstructured, wxios_closedef 34 37 IMPLICIT NONE 35 38 … … 60 63 ! 2. Declare horizontal domain 61 64 ! Set values for the mask: 62 IF (mpi_rank == 0) THEN 63 data_ibegin = 0 64 ELSE 65 data_ibegin = ii_begin - 1 66 END IF 67 68 IF (mpi_rank == mpi_size-1) THEN 69 data_iend = nbp_lon 70 ELSE 71 data_iend = ii_end + 1 72 END IF 73 74 if (prt_level>=10) then 75 write(lunout,*) "initialize_xios_output: mpirank=",mpi_rank," iibegin=",ii_begin , " ii_end=",ii_end," jjbegin=",jj_begin," jj_nb=",jj_nb," jj_end=",jj_end 76 write(lunout,*) "initialize_xios_output: mpirank=",mpi_rank," nbp_lon=",nbp_lon," nbp_lat=",nbp_lat 77 write(lunout,*) "initialize_xios_output: mpirank=",mpi_rank," data_ibegin=",data_ibegin," data_iend=",data_iend 78 write(lunout,*) "initialize_xios_output: mpirank=",mpi_rank," data_ibegin=",data_ibegin," data_iend=",data_iend 79 write(lunout,*) "initialize_xios_output: mpirank=",mpi_rank," is_south_pole=",is_south_pole_dyn 80 endif 81 65 ! IF (mpi_rank == 0) THEN 66 ! data_ibegin = 0 67 ! ELSE 68 ! data_ibegin = ii_begin - 1 69 ! END IF 70 71 ! IF (mpi_rank == mpi_size-1) THEN 72 ! data_iend = nbp_lon 73 ! ELSE 74 ! data_iend = ii_end + 1 75 ! END IF 76 77 ! if (prt_level>=10) then 78 ! write(lunout,*) "initialize_xios_output: mpirank=",mpi_rank," iibegin=",ii_begin , " ii_end=",ii_end," jjbegin=",jj_begin," jj_nb=",jj_nb," jj_end=",jj_end 79 ! write(lunout,*) "initialize_xios_output: mpirank=",mpi_rank," nbp_lon=",nbp_lon," nbp_lat=",nbp_lat 80 ! write(lunout,*) "initialize_xios_output: mpirank=",mpi_rank," data_ibegin=",data_ibegin," data_iend=",data_iend 81 ! write(lunout,*) "initialize_xios_output: mpirank=",mpi_rank," data_ibegin=",data_ibegin," data_iend=",data_iend 82 ! write(lunout,*) "initialize_xios_output: mpirank=",mpi_rank," is_south_pole=",is_south_pole_dyn 83 ! endif 84 85 !$OMP END MASTER 86 !$OMP BARRIER 82 87 ! Initialize the XIOS domain coreesponding to this process: 83 88 if (prt_level>=10) write(lunout,*) "initialize_xios_output: call wxios_domain_param" 84 CALL wxios_domain_param("dom_glo", is_sequential, nbp_lon, jj_nb, nbp_lon, nbp_lat, & 85 1, nbp_lon, ii_begin, ii_end, jj_begin, jj_end, & 86 klon_mpi+2*(nbp_lon-1), data_ibegin, data_iend, & 87 lat_reg*(180./pi), lon_reg*(180./pi), & 88 is_south_pole_dyn,mpi_rank) 89 89 ! CALL wxios_domain_param("dom_glo", is_sequential, nbp_lon, jj_nb, nbp_lon, nbp_lat, & 90 ! 1, nbp_lon, ii_begin, ii_end, jj_begin, jj_end, & 91 ! klon_mpi+2*(nbp_lon-1), data_ibegin, data_iend, & 92 ! lat_reg*(180./pi), lon_reg*(180./pi), & 93 ! is_south_pole_dyn,mpi_rank) 94 95 IF (grid_type==unstructured) THEN 96 CALL wxios_domain_param_unstructured("dom_glo") 97 ELSE 98 CALL wxios_domain_param("dom_glo") 99 ENDIF 100 101 !$OMP MASTER 90 102 ! 3. Declare calendar and time step 91 103 if (prt_level>=10) write(lunout,*) "initialize_xios_output: build calendar"
Note: See TracChangeset
for help on using the changeset viewer.