| 1 | !WRF:MEDIATION:IO |
|---|
| 2 | ! ---principal wrf output routine (called from routines in module_io_domain ) |
|---|
| 3 | SUBROUTINE output_wrf ( fid , grid , config_flags, switch , ierr ) |
|---|
| 4 | USE module_io |
|---|
| 5 | USE module_wrf_error |
|---|
| 6 | USE module_io_wrf |
|---|
| 7 | USE module_domain |
|---|
| 8 | USE module_state_description |
|---|
| 9 | USE module_configure |
|---|
| 10 | ! USE module_date_time |
|---|
| 11 | USE module_model_constants |
|---|
| 12 | USE module_utility |
|---|
| 13 | IMPLICIT NONE |
|---|
| 14 | #include <wrf_io_flags.h> |
|---|
| 15 | #include <wrf_status_codes.h> |
|---|
| 16 | TYPE(domain) :: grid |
|---|
| 17 | TYPE(grid_config_rec_type), INTENT(INOUT) :: config_flags |
|---|
| 18 | INTEGER, INTENT(IN) :: fid, switch |
|---|
| 19 | INTEGER, INTENT(INOUT) :: ierr |
|---|
| 20 | |
|---|
| 21 | ! Local data |
|---|
| 22 | INTEGER ids , ide , jds , jde , kds , kde , & |
|---|
| 23 | ims , ime , jms , jme , kms , kme , & |
|---|
| 24 | ips , ipe , jps , jpe , kps , kpe |
|---|
| 25 | |
|---|
| 26 | INTEGER , DIMENSION(3) :: domain_start , domain_end |
|---|
| 27 | INTEGER , DIMENSION(3) :: memory_start , memory_end |
|---|
| 28 | INTEGER , DIMENSION(3) :: patch_start , patch_end |
|---|
| 29 | INTEGER i,j |
|---|
| 30 | INTEGER julyr, julday, idt, iswater , map_proj |
|---|
| 31 | INTEGER filestate |
|---|
| 32 | LOGICAL dryrun |
|---|
| 33 | REAL gmt, cen_lat, cen_lon, bdyfrq , truelat1 , truelat2 , moad_cen_lat , stand_lon |
|---|
| 34 | INTEGER km_opt, diff_opt, damp_opt, & |
|---|
| 35 | mp_physics, ra_lw_physics, ra_sw_physics, sf_sfclay_physics, & |
|---|
| 36 | sf_surface_physics, bl_pbl_physics, cu_physics |
|---|
| 37 | REAL khdif, kvdif, swrad_scat |
|---|
| 38 | INTEGER ucmcall, w_damping, smooth_option, feedback, surface_input_source, sst_update |
|---|
| 39 | #if (EM_CORE == 1) |
|---|
| 40 | INTEGER grid_id , parent_id , i_parent_start , j_parent_start , parent_grid_ratio |
|---|
| 41 | INTEGER diff_6th_opt |
|---|
| 42 | REAL diff_6th_factor |
|---|
| 43 | INTEGER grid_fdda, gfdda_interval_m, gfdda_end_h, if_ramping, & |
|---|
| 44 | obs_nudge_opt, obs_nudge_wind, obs_nudge_temp, obs_nudge_mois, obs_nudge_pstr, obs_idynin, obs_ionf |
|---|
| 45 | REAL fgdt, guv, gt, gq, dtramp_min, & |
|---|
| 46 | obs_coef_wind, obs_coef_temp, obs_coef_mois, obs_coef_pstr, obs_dtramp, fdda_end |
|---|
| 47 | LOGICAL pd_moist, pd_scalar, pd_tke |
|---|
| 48 | #endif |
|---|
| 49 | CHARACTER (len=19) simulation_start_date |
|---|
| 50 | CHARACTER (len=len_current_date) current_date_save |
|---|
| 51 | INTEGER simulation_start_year , & |
|---|
| 52 | simulation_start_month , & |
|---|
| 53 | simulation_start_day , & |
|---|
| 54 | simulation_start_hour , & |
|---|
| 55 | simulation_start_minute , & |
|---|
| 56 | simulation_start_second |
|---|
| 57 | INTEGER rc |
|---|
| 58 | INTEGER :: io_form |
|---|
| 59 | LOGICAL, EXTERNAL :: multi_files |
|---|
| 60 | INTEGER, EXTERNAL :: use_package |
|---|
| 61 | INTEGER p_hr, p_min, p_sec, p_ms |
|---|
| 62 | |
|---|
| 63 | CHARACTER*256 message |
|---|
| 64 | CHARACTER*80 fname |
|---|
| 65 | CHARACTER*80 char_junk |
|---|
| 66 | INTEGER ibuf(1) |
|---|
| 67 | REAL rbuf(1) |
|---|
| 68 | TYPE(WRFU_TimeInterval) :: bdy_increment |
|---|
| 69 | TYPE(WRFU_Time) :: next_time, currentTime, startTime |
|---|
| 70 | CHARACTER*40 :: next_datestr |
|---|
| 71 | INTEGER :: start_year , start_month , start_day , start_hour , start_minute , start_second |
|---|
| 72 | LOGICAL :: adjust |
|---|
| 73 | |
|---|
| 74 | WRITE(wrf_err_message,*)'output_wrf: begin, fid = ',fid |
|---|
| 75 | CALL wrf_debug( 300 , wrf_err_message ) |
|---|
| 76 | |
|---|
| 77 | CALL wrf_inquire_filename ( fid , fname , filestate , ierr ) |
|---|
| 78 | IF ( ierr /= 0 ) THEN |
|---|
| 79 | WRITE(wrf_err_message,*)'module_io_wrf: output_wrf: wrf_inquire_filename Status = ',ierr |
|---|
| 80 | CALL wrf_error_fatal( wrf_err_message ) |
|---|
| 81 | ENDIF |
|---|
| 82 | |
|---|
| 83 | WRITE(wrf_err_message,*)'output_wrf: fid,filestate = ',fid,filestate |
|---|
| 84 | CALL wrf_debug( 300 , wrf_err_message ) |
|---|
| 85 | |
|---|
| 86 | ! io_form is used to determine if multi-file I/O is enabled and to |
|---|
| 87 | ! control writing of format-specific time-independent metadata |
|---|
| 88 | IF ( switch .EQ. model_input_only ) THEN |
|---|
| 89 | CALL nl_get_io_form_input( 1, io_form ) |
|---|
| 90 | ELSE IF ( switch .EQ. aux_model_input1_only ) THEN |
|---|
| 91 | CALL nl_get_io_form_auxinput1( 1, io_form ) |
|---|
| 92 | ELSE IF ( switch .EQ. aux_model_input2_only ) THEN |
|---|
| 93 | CALL nl_get_io_form_auxinput2( 1, io_form ) |
|---|
| 94 | ELSE IF ( switch .EQ. aux_model_input3_only ) THEN |
|---|
| 95 | CALL nl_get_io_form_auxinput3( 1, io_form ) |
|---|
| 96 | ELSE IF ( switch .EQ. aux_model_input4_only ) THEN |
|---|
| 97 | CALL nl_get_io_form_auxinput4( 1, io_form ) |
|---|
| 98 | ELSE IF ( switch .EQ. aux_model_input5_only ) THEN |
|---|
| 99 | CALL nl_get_io_form_auxinput5( 1, io_form ) |
|---|
| 100 | ELSE IF ( switch .EQ. aux_model_input6_only ) THEN |
|---|
| 101 | CALL nl_get_io_form_auxinput6( 1, io_form ) |
|---|
| 102 | ELSE IF ( switch .EQ. aux_model_input7_only ) THEN |
|---|
| 103 | CALL nl_get_io_form_auxinput7( 1, io_form ) |
|---|
| 104 | ELSE IF ( switch .EQ. aux_model_input8_only ) THEN |
|---|
| 105 | CALL nl_get_io_form_auxinput8( 1, io_form ) |
|---|
| 106 | ELSE IF ( switch .EQ. aux_model_input9_only ) THEN |
|---|
| 107 | CALL nl_get_io_form_auxinput9( 1, io_form ) |
|---|
| 108 | ELSE IF ( switch .EQ. aux_model_input10_only ) THEN |
|---|
| 109 | CALL nl_get_io_form_gfdda( 1, io_form ) |
|---|
| 110 | ELSE IF ( switch .EQ. aux_model_input11_only ) THEN |
|---|
| 111 | CALL nl_get_io_form_auxinput11( 1, io_form ) |
|---|
| 112 | |
|---|
| 113 | ELSE IF ( switch .EQ. history_only ) THEN |
|---|
| 114 | CALL nl_get_io_form_history( 1, io_form ) |
|---|
| 115 | ELSE IF ( switch .EQ. aux_hist1_only ) THEN |
|---|
| 116 | CALL nl_get_io_form_auxhist1( 1, io_form ) |
|---|
| 117 | ELSE IF ( switch .EQ. aux_hist2_only ) THEN |
|---|
| 118 | CALL nl_get_io_form_auxhist2( 1, io_form ) |
|---|
| 119 | ELSE IF ( switch .EQ. aux_hist3_only ) THEN |
|---|
| 120 | CALL nl_get_io_form_auxhist3( 1, io_form ) |
|---|
| 121 | ELSE IF ( switch .EQ. aux_hist4_only ) THEN |
|---|
| 122 | CALL nl_get_io_form_auxhist4( 1, io_form ) |
|---|
| 123 | ELSE IF ( switch .EQ. aux_hist5_only ) THEN |
|---|
| 124 | CALL nl_get_io_form_auxhist5( 1, io_form ) |
|---|
| 125 | ELSE IF ( switch .EQ. aux_hist6_only ) THEN |
|---|
| 126 | CALL nl_get_io_form_auxhist6( 1, io_form ) |
|---|
| 127 | ELSE IF ( switch .EQ. aux_hist7_only ) THEN |
|---|
| 128 | CALL nl_get_io_form_auxhist7( 1, io_form ) |
|---|
| 129 | ELSE IF ( switch .EQ. aux_hist8_only ) THEN |
|---|
| 130 | CALL nl_get_io_form_auxhist8( 1, io_form ) |
|---|
| 131 | ELSE IF ( switch .EQ. aux_hist9_only ) THEN |
|---|
| 132 | CALL nl_get_io_form_auxhist9( 1, io_form ) |
|---|
| 133 | ELSE IF ( switch .EQ. aux_hist10_only ) THEN |
|---|
| 134 | CALL nl_get_io_form_auxhist10( 1, io_form ) |
|---|
| 135 | ELSE IF ( switch .EQ. aux_hist11_only ) THEN |
|---|
| 136 | CALL nl_get_io_form_auxhist11( 1, io_form ) |
|---|
| 137 | |
|---|
| 138 | ELSE IF ( switch .EQ. restart_only ) THEN |
|---|
| 139 | CALL nl_get_io_form_restart( 1, io_form ) |
|---|
| 140 | ELSE IF ( switch .EQ. boundary_only ) THEN |
|---|
| 141 | CALL nl_get_io_form_boundary( 1, io_form ) |
|---|
| 142 | ELSE ! default: use history |
|---|
| 143 | CALL nl_get_io_form_history( 1, io_form ) |
|---|
| 144 | ENDIF |
|---|
| 145 | |
|---|
| 146 | dryrun = ( filestate .EQ. WRF_FILE_OPENED_NOT_COMMITTED ) |
|---|
| 147 | |
|---|
| 148 | WRITE(wrf_err_message,*)'output_wrf: dryrun = ',dryrun |
|---|
| 149 | CALL wrf_debug( 300 , wrf_err_message ) |
|---|
| 150 | |
|---|
| 151 | CALL get_ijk_from_grid ( grid , & |
|---|
| 152 | ids, ide, jds, jde, kds, kde, & |
|---|
| 153 | ims, ime, jms, jme, kms, kme, & |
|---|
| 154 | ips, ipe, jps, jpe, kps, kpe ) |
|---|
| 155 | |
|---|
| 156 | call nl_get_diff_opt ( 1, diff_opt ) |
|---|
| 157 | call nl_get_km_opt ( 1, km_opt ) |
|---|
| 158 | call nl_get_damp_opt ( 1, damp_opt ) |
|---|
| 159 | call nl_get_khdif ( grid%id, khdif ) |
|---|
| 160 | call nl_get_kvdif ( grid%id, kvdif ) |
|---|
| 161 | call nl_get_mp_physics ( grid%id, mp_physics ) |
|---|
| 162 | call nl_get_ra_lw_physics ( grid%id, ra_lw_physics ) |
|---|
| 163 | call nl_get_ra_sw_physics ( grid%id, ra_sw_physics ) |
|---|
| 164 | call nl_get_sf_sfclay_physics ( grid%id, sf_sfclay_physics ) |
|---|
| 165 | call nl_get_sf_surface_physics ( grid%id, sf_surface_physics ) |
|---|
| 166 | call nl_get_bl_pbl_physics ( grid%id, bl_pbl_physics ) |
|---|
| 167 | call nl_get_cu_physics ( grid%id, cu_physics ) |
|---|
| 168 | |
|---|
| 169 | ! add nml variables in 2.2 |
|---|
| 170 | call nl_get_surface_input_source ( 1 , surface_input_source ) |
|---|
| 171 | call nl_get_sst_update ( 1 , sst_update ) |
|---|
| 172 | call nl_get_feedback ( 1 , feedback ) |
|---|
| 173 | call nl_get_smooth_option ( 1 , smooth_option ) |
|---|
| 174 | call nl_get_swrad_scat ( 1 , swrad_scat ) |
|---|
| 175 | call nl_get_ucmcall ( 1 , ucmcall ) |
|---|
| 176 | call nl_get_w_damping ( 1 , w_damping ) |
|---|
| 177 | |
|---|
| 178 | #if (EM_CORE == 1) |
|---|
| 179 | CALL nl_get_pd_moist ( grid%id , pd_moist ) |
|---|
| 180 | CALL nl_get_pd_scalar ( grid%id , pd_scalar ) |
|---|
| 181 | CALL nl_get_pd_tke ( grid%id , pd_tke ) |
|---|
| 182 | CALL nl_get_diff_6th_opt ( grid%id , diff_6th_opt ) |
|---|
| 183 | CALL nl_get_diff_6th_factor ( grid%id , diff_6th_factor ) |
|---|
| 184 | CALL nl_get_grid_fdda ( grid%id , grid_fdda ) |
|---|
| 185 | CALL nl_get_gfdda_end_h( grid%id , gfdda_end_h ) |
|---|
| 186 | CALL nl_get_gfdda_interval_m ( grid%id , gfdda_interval_m ) |
|---|
| 187 | |
|---|
| 188 | IF ( grid_fdda == 1 ) THEN |
|---|
| 189 | CALL nl_get_fgdt ( grid%id , fgdt ) |
|---|
| 190 | CALL nl_get_guv ( grid%id , guv ) |
|---|
| 191 | CALL nl_get_gt ( grid%id , gt ) |
|---|
| 192 | CALL nl_get_gq ( grid%id , gq ) |
|---|
| 193 | CALL nl_get_if_ramping ( 1 , if_ramping ) |
|---|
| 194 | CALL nl_get_dtramp_min ( 1 , dtramp_min ) |
|---|
| 195 | ENDIF |
|---|
| 196 | |
|---|
| 197 | CALL nl_get_obs_nudge_opt ( grid%id , obs_nudge_opt ) |
|---|
| 198 | IF ( obs_nudge_opt == 1 ) THEN |
|---|
| 199 | CALL nl_get_fdda_end ( grid%id , fdda_end ) |
|---|
| 200 | CALL nl_get_obs_nudge_wind ( grid%id , obs_nudge_wind ) |
|---|
| 201 | CALL nl_get_obs_coef_wind ( grid%id , obs_coef_wind ) |
|---|
| 202 | CALL nl_get_obs_nudge_temp ( grid%id , obs_nudge_temp ) |
|---|
| 203 | CALL nl_get_obs_coef_temp ( grid%id , obs_coef_temp ) |
|---|
| 204 | CALL nl_get_obs_nudge_mois ( grid%id , obs_nudge_mois ) |
|---|
| 205 | CALL nl_get_obs_coef_mois ( grid%id , obs_coef_mois ) |
|---|
| 206 | CALL nl_get_obs_nudge_pstr ( grid%id , obs_nudge_pstr ) |
|---|
| 207 | CALL nl_get_obs_coef_pstr ( grid%id , obs_coef_pstr ) |
|---|
| 208 | CALL nl_get_obs_ionf ( 1 , obs_ionf ) |
|---|
| 209 | CALL nl_get_obs_idynin ( 1 , obs_idynin ) |
|---|
| 210 | CALL nl_get_obs_dtramp ( 1 , obs_dtramp ) |
|---|
| 211 | ENDIF |
|---|
| 212 | #endif |
|---|
| 213 | |
|---|
| 214 | ! julday and gmt can be set in namelist_03 for ideal.exe run |
|---|
| 215 | CALL nl_get_gmt (grid%id, gmt) |
|---|
| 216 | CALL nl_get_julyr (grid%id, julyr) |
|---|
| 217 | CALL nl_get_julday (grid%id, julday) |
|---|
| 218 | CALL nl_get_mminlu ( 1, char_junk(1:4) ) |
|---|
| 219 | CALL nl_get_iswater (grid%id, iswater ) |
|---|
| 220 | CALL nl_get_cen_lat ( grid%id , cen_lat ) |
|---|
| 221 | CALL nl_get_cen_lon ( grid%id , cen_lon ) |
|---|
| 222 | CALL nl_get_truelat1 ( grid%id , truelat1 ) |
|---|
| 223 | CALL nl_get_truelat2 ( grid%id , truelat2 ) |
|---|
| 224 | CALL nl_get_moad_cen_lat ( grid%id , moad_cen_lat ) |
|---|
| 225 | CALL nl_get_stand_lon ( grid%id , stand_lon ) |
|---|
| 226 | CALL nl_get_map_proj ( grid%id , map_proj ) |
|---|
| 227 | |
|---|
| 228 | #if (EM_CORE == 1) |
|---|
| 229 | CALL nl_get_parent_id ( grid%id , parent_id ) |
|---|
| 230 | CALL nl_get_i_parent_start ( grid%id , i_parent_start ) |
|---|
| 231 | CALL nl_get_j_parent_start ( grid%id , j_parent_start ) |
|---|
| 232 | CALL nl_get_parent_grid_ratio ( grid%id , parent_grid_ratio ) |
|---|
| 233 | #endif |
|---|
| 234 | |
|---|
| 235 | CALL domain_clockprint(150, grid, & |
|---|
| 236 | 'DEBUG output_wrf(): before call to domain_clock_get,') |
|---|
| 237 | CALL domain_clock_get( grid, current_time=currentTime, & |
|---|
| 238 | start_time=startTime, & |
|---|
| 239 | current_timestr=current_date ) |
|---|
| 240 | WRITE ( wrf_err_message , * ) 'output_wrf: begin, current_date=',current_date |
|---|
| 241 | CALL wrf_debug ( 300 , wrf_err_message ) |
|---|
| 242 | |
|---|
| 243 | WRITE( message , * ) "OUTPUT FROM " , TRIM(program_name) |
|---|
| 244 | CALL wrf_put_dom_ti_char ( fid , 'TITLE' , TRIM(message) , ierr ) |
|---|
| 245 | ! added grib-specific metadata: Todd Hutchinson 8/21/2005 |
|---|
| 246 | IF ( ( use_package( io_form ) == IO_GRIB1 ) .OR. & |
|---|
| 247 | ( use_package( io_form ) == IO_GRIB2 ) ) THEN |
|---|
| 248 | CALL wrf_put_dom_ti_char ( fid, 'PROGRAM_NAME', TRIM(program_name) , ierr ) |
|---|
| 249 | ENDIF |
|---|
| 250 | CALL nl_get_start_year(grid%id,start_year) |
|---|
| 251 | CALL nl_get_start_month(grid%id,start_month) |
|---|
| 252 | CALL nl_get_start_day(grid%id,start_day) |
|---|
| 253 | CALL nl_get_start_hour(grid%id,start_hour) |
|---|
| 254 | CALL nl_get_start_minute(grid%id,start_minute) |
|---|
| 255 | CALL nl_get_start_second(grid%id,start_second) |
|---|
| 256 | #ifdef PLANET |
|---|
| 257 | WRITE ( start_date , FMT = '(I4.4,"-",I5.5,"_",I2.2,":",I2.2,":",I2.2)' ) & |
|---|
| 258 | start_year,start_day,start_hour,start_minute,start_second |
|---|
| 259 | #else |
|---|
| 260 | WRITE ( start_date , FMT = '(I4.4,"-",I2.2,"-",I2.2,"_",I2.2,":",I2.2,":",I2.2)' ) & |
|---|
| 261 | start_year,start_month,start_day,start_hour,start_minute,start_second |
|---|
| 262 | #endif |
|---|
| 263 | CALL wrf_put_dom_ti_char ( fid , 'START_DATE', TRIM(start_date) , ierr ) |
|---|
| 264 | IF ( switch .EQ. model_input_only) THEN |
|---|
| 265 | CALL wrf_put_dom_ti_char ( fid , 'SIMULATION_START_DATE', TRIM(start_date) , ierr ) |
|---|
| 266 | ELSE IF ( ( switch .EQ. restart_only ) .OR. ( switch .EQ. history_only ) ) THEN |
|---|
| 267 | CALL nl_get_simulation_start_year ( 1, simulation_start_year ) |
|---|
| 268 | CALL nl_get_simulation_start_month ( 1, simulation_start_month ) |
|---|
| 269 | CALL nl_get_simulation_start_day ( 1, simulation_start_day ) |
|---|
| 270 | CALL nl_get_simulation_start_hour ( 1, simulation_start_hour ) |
|---|
| 271 | CALL nl_get_simulation_start_minute ( 1, simulation_start_minute ) |
|---|
| 272 | CALL nl_get_simulation_start_second ( 1, simulation_start_second ) |
|---|
| 273 | #ifdef PLANET |
|---|
| 274 | WRITE ( simulation_start_date , FMT = '(I4.4,"-",I5.5,"_",I2.2,":",I2.2,":",I2.2)' ) & |
|---|
| 275 | simulation_start_year,simulation_start_day,& |
|---|
| 276 | simulation_start_hour,simulation_start_minute,simulation_start_second |
|---|
| 277 | #else |
|---|
| 278 | WRITE ( simulation_start_date , FMT = '(I4.4,"-",I2.2,"-",I2.2,"_",I2.2,":",I2.2,":",I2.2)' ) & |
|---|
| 279 | simulation_start_year,simulation_start_month,simulation_start_day,& |
|---|
| 280 | simulation_start_hour,simulation_start_minute,simulation_start_second |
|---|
| 281 | #endif |
|---|
| 282 | CALL wrf_put_dom_ti_char ( fid , 'SIMULATION_START_DATE', TRIM(simulation_start_date) , ierr ) |
|---|
| 283 | END IF |
|---|
| 284 | |
|---|
| 285 | ibuf(1) = config_flags%e_we - config_flags%s_we + 1 |
|---|
| 286 | CALL wrf_put_dom_ti_integer ( fid , 'WEST-EAST_GRID_DIMENSION' , ibuf , 1 , ierr ) |
|---|
| 287 | |
|---|
| 288 | ibuf(1) = config_flags%e_sn - config_flags%s_sn + 1 |
|---|
| 289 | CALL wrf_put_dom_ti_integer ( fid , 'SOUTH-NORTH_GRID_DIMENSION' , ibuf , 1 , ierr ) |
|---|
| 290 | |
|---|
| 291 | ibuf(1) = config_flags%e_vert - config_flags%s_vert + 1 |
|---|
| 292 | CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_GRID_DIMENSION' , ibuf , 1 , ierr ) |
|---|
| 293 | |
|---|
| 294 | #if (EM_CORE == 1) |
|---|
| 295 | IF (grid%map_proj == 6) THEN |
|---|
| 296 | ! Global ... dx determined automatically |
|---|
| 297 | ! Don't use value from namelist ... used derived value instead |
|---|
| 298 | CALL wrf_put_dom_ti_real ( fid , 'DX' , grid%dx , 1 , ierr ) |
|---|
| 299 | CALL wrf_put_dom_ti_real ( fid , 'DY' , grid%dy , 1 , ierr ) |
|---|
| 300 | ELSE |
|---|
| 301 | CALL wrf_put_dom_ti_real ( fid , 'DX' , config_flags%dx , 1 , ierr ) |
|---|
| 302 | CALL wrf_put_dom_ti_real ( fid , 'DY' , config_flags%dy , 1 , ierr ) |
|---|
| 303 | END IF |
|---|
| 304 | #endif |
|---|
| 305 | |
|---|
| 306 | ! added this metadatum for H. Chuan, NCEP, 030417, JM |
|---|
| 307 | #if (NMM_CORE == 1) |
|---|
| 308 | CALL wrf_put_dom_ti_char ( fid , 'GRIDTYPE', 'E' , ierr ) |
|---|
| 309 | #endif |
|---|
| 310 | #if (EM_CORE == 1) |
|---|
| 311 | CALL wrf_put_dom_ti_char ( fid , 'GRIDTYPE', 'C' , ierr ) |
|---|
| 312 | #endif |
|---|
| 313 | #if (COAMPS_CORE == 1 ) |
|---|
| 314 | CALL wrf_put_dom_ti_char ( fid , 'GRIDTYPE', 'B' , ierr ) |
|---|
| 315 | #endif |
|---|
| 316 | |
|---|
| 317 | ! added these fields for W. Skamarock, 020402, JM |
|---|
| 318 | ibuf(1) = diff_opt |
|---|
| 319 | CALL wrf_put_dom_ti_integer ( fid , 'DIFF_OPT' , ibuf , 1 , ierr ) |
|---|
| 320 | ibuf(1) = km_opt |
|---|
| 321 | CALL wrf_put_dom_ti_integer ( fid , 'KM_OPT' , ibuf , 1 , ierr ) |
|---|
| 322 | ibuf(1) = damp_opt |
|---|
| 323 | CALL wrf_put_dom_ti_integer ( fid , 'DAMP_OPT' , ibuf , 1 , ierr ) |
|---|
| 324 | rbuf(1) = khdif |
|---|
| 325 | CALL wrf_put_dom_ti_real ( fid , 'KHDIF' , rbuf , 1 , ierr ) |
|---|
| 326 | rbuf(1) = kvdif |
|---|
| 327 | CALL wrf_put_dom_ti_real ( fid , 'KVDIF' , rbuf , 1 , ierr ) |
|---|
| 328 | ibuf(1) = mp_physics |
|---|
| 329 | CALL wrf_put_dom_ti_integer ( fid , 'MP_PHYSICS' , ibuf , 1 , ierr ) |
|---|
| 330 | ibuf(1) = ra_lw_physics |
|---|
| 331 | CALL wrf_put_dom_ti_integer ( fid , 'RA_LW_PHYSICS' , ibuf , 1 , ierr ) |
|---|
| 332 | ibuf(1) = ra_sw_physics |
|---|
| 333 | CALL wrf_put_dom_ti_integer ( fid , 'RA_SW_PHYSICS' , ibuf , 1 , ierr ) |
|---|
| 334 | ibuf(1) = sf_sfclay_physics |
|---|
| 335 | CALL wrf_put_dom_ti_integer ( fid , 'SF_SFCLAY_PHYSICS' , ibuf , 1 , ierr ) |
|---|
| 336 | ibuf(1) = sf_surface_physics |
|---|
| 337 | CALL wrf_put_dom_ti_integer ( fid , 'SF_SURFACE_PHYSICS' , ibuf , 1 , ierr ) |
|---|
| 338 | ibuf(1) = bl_pbl_physics |
|---|
| 339 | CALL wrf_put_dom_ti_integer ( fid , 'BL_PBL_PHYSICS' , ibuf , 1 , ierr ) |
|---|
| 340 | ibuf(1) = cu_physics |
|---|
| 341 | CALL wrf_put_dom_ti_integer ( fid , 'CU_PHYSICS' , ibuf , 1 , ierr ) |
|---|
| 342 | |
|---|
| 343 | ! added netcdf-specific metadata: |
|---|
| 344 | IF ( ( use_package( io_form ) == IO_NETCDF ) .OR. & |
|---|
| 345 | ( use_package( io_form ) == IO_PHDF5 ) .OR. & |
|---|
| 346 | ( use_package( io_form ) == IO_PNETCDF ) ) THEN |
|---|
| 347 | CALL wrf_put_dom_ti_integer ( fid, 'SURFACE_INPUT_SOURCE', surface_input_source , 1 , ierr ) |
|---|
| 348 | CALL wrf_put_dom_ti_integer ( fid, 'SST_UPDATE', sst_update , 1 , ierr ) |
|---|
| 349 | #if (EM_CORE == 1) |
|---|
| 350 | CALL wrf_put_dom_ti_integer ( fid, 'GRID_FDDA', grid_fdda , 1 , ierr ) |
|---|
| 351 | CALL wrf_put_dom_ti_integer ( fid, 'GFDDA_INTERVAL_M', gfdda_interval_m , 1 , ierr ) |
|---|
| 352 | CALL wrf_put_dom_ti_integer ( fid, 'GFDDA_END_H', gfdda_end_h , 1 , ierr ) |
|---|
| 353 | #endif |
|---|
| 354 | |
|---|
| 355 | IF ( switch .EQ. history_only ) THEN |
|---|
| 356 | CALL wrf_put_dom_ti_integer ( fid, 'UCMCALL', ucmcall , 1 , ierr ) |
|---|
| 357 | CALL wrf_put_dom_ti_integer ( fid, 'FEEDBACK', feedback , 1 , ierr ) |
|---|
| 358 | CALL wrf_put_dom_ti_integer ( fid, 'SMOOTH_OPTION', smooth_option , 1 , ierr ) |
|---|
| 359 | CALL wrf_put_dom_ti_real ( fid, 'SWRAD_SCAT', swrad_scat , 1 , ierr ) |
|---|
| 360 | CALL wrf_put_dom_ti_integer ( fid, 'W_DAMPING', w_damping , 1 , ierr ) |
|---|
| 361 | |
|---|
| 362 | #if (EM_CORE == 1) |
|---|
| 363 | CALL wrf_put_dom_ti_logical ( fid, 'PD_MOIST', pd_moist , 1 , ierr ) |
|---|
| 364 | CALL wrf_put_dom_ti_logical ( fid, 'PD_SCALAR', pd_scalar , 1 , ierr ) |
|---|
| 365 | CALL wrf_put_dom_ti_logical ( fid, 'PD_TKE', pd_tke , 1 , ierr ) |
|---|
| 366 | CALL wrf_put_dom_ti_integer ( fid, 'DIFF_6TH_OPT', diff_6th_opt , 1 , ierr ) |
|---|
| 367 | CALL wrf_put_dom_ti_real ( fid, 'DIFF_6TH_FACTOR', diff_6th_factor , 1 , ierr ) |
|---|
| 368 | |
|---|
| 369 | IF ( grid_fdda == 1 ) THEN |
|---|
| 370 | CALL wrf_put_dom_ti_real ( fid, 'FGDT', fgdt , 1 , ierr ) |
|---|
| 371 | CALL wrf_put_dom_ti_real ( fid, 'GUV', guv , 1 , ierr ) |
|---|
| 372 | CALL wrf_put_dom_ti_real ( fid, 'GT', gt , 1 , ierr ) |
|---|
| 373 | CALL wrf_put_dom_ti_real ( fid, 'GQ', gq , 1 , ierr ) |
|---|
| 374 | CALL wrf_put_dom_ti_integer ( fid, 'IF_RAMPING', if_ramping , 1 , ierr ) |
|---|
| 375 | CALL wrf_put_dom_ti_real ( fid, 'DTRAMP_MIN', dtramp_min , 1 , ierr ) |
|---|
| 376 | ENDIF |
|---|
| 377 | |
|---|
| 378 | CALL wrf_put_dom_ti_integer ( fid, 'OBS_NUDGE_OPT', obs_nudge_opt , 1 , ierr ) |
|---|
| 379 | IF ( obs_nudge_opt == 1 ) THEN |
|---|
| 380 | CALL wrf_put_dom_ti_real ( fid, 'FDDA_END', fdda_end , 1 , ierr ) |
|---|
| 381 | CALL wrf_put_dom_ti_integer ( fid, 'OBS_NUDGE_WIND', obs_nudge_wind , 1 , ierr ) |
|---|
| 382 | CALL wrf_put_dom_ti_real ( fid, 'OBS_COEF_WIND', obs_coef_wind , 1 , ierr ) |
|---|
| 383 | CALL wrf_put_dom_ti_integer ( fid, 'OBS_NUDGE_TEMP', obs_nudge_temp , 1 , ierr ) |
|---|
| 384 | CALL wrf_put_dom_ti_real ( fid, 'OBS_COEF_TEMP', obs_coef_temp , 1 , ierr ) |
|---|
| 385 | CALL wrf_put_dom_ti_integer ( fid, 'OBS_NUDGE_MOIS', obs_nudge_mois , 1 , ierr ) |
|---|
| 386 | CALL wrf_put_dom_ti_real ( fid, 'OBS_COEF_MOIS', obs_coef_mois , 1 , ierr ) |
|---|
| 387 | CALL wrf_put_dom_ti_integer ( fid, 'OBS_NUDGE_PSTR', obs_nudge_pstr , 1 , ierr ) |
|---|
| 388 | CALL wrf_put_dom_ti_real ( fid, 'OBS_COEF_PSTR', obs_coef_pstr , 1 , ierr ) |
|---|
| 389 | CALL wrf_put_dom_ti_integer ( fid, 'OBS_IONF', obs_ionf , 1 , ierr ) |
|---|
| 390 | CALL wrf_put_dom_ti_integer ( fid, 'OBS_IDYNIN', obs_idynin , 1 , ierr ) |
|---|
| 391 | CALL wrf_put_dom_ti_real ( fid, 'OBS_DTRAMP', obs_dtramp , 1 , ierr ) |
|---|
| 392 | ENDIF |
|---|
| 393 | #endif |
|---|
| 394 | ENDIF ! history_only |
|---|
| 395 | ENDIF |
|---|
| 396 | |
|---|
| 397 | ! added these fields for use by reassembly programs , 010831, JM |
|---|
| 398 | ! modified these fields so "patch" == "domain" when multi-file output |
|---|
| 399 | ! formats are not used. 051018, TBH |
|---|
| 400 | |
|---|
| 401 | ibuf(1) = MAX(ips,ids) |
|---|
| 402 | IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = ids |
|---|
| 403 | CALL wrf_put_dom_ti_integer ( fid , 'WEST-EAST_PATCH_START_UNSTAG' , ibuf , 1 , ierr ) |
|---|
| 404 | ibuf(1) = MIN(ipe,ide-1) |
|---|
| 405 | IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = ide - 1 |
|---|
| 406 | CALL wrf_put_dom_ti_integer ( fid , 'WEST-EAST_PATCH_END_UNSTAG' , ibuf , 1 , ierr ) |
|---|
| 407 | ibuf(1) = MAX(ips,ids) |
|---|
| 408 | IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = ids |
|---|
| 409 | CALL wrf_put_dom_ti_integer ( fid , 'WEST-EAST_PATCH_START_STAG' , ibuf , 1 , ierr ) |
|---|
| 410 | ibuf(1) = MIN(ipe,ide) |
|---|
| 411 | IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = ide |
|---|
| 412 | CALL wrf_put_dom_ti_integer ( fid , 'WEST-EAST_PATCH_END_STAG' , ibuf , 1 , ierr ) |
|---|
| 413 | ibuf(1) = MAX(jps,jds) |
|---|
| 414 | IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = jds |
|---|
| 415 | CALL wrf_put_dom_ti_integer ( fid , 'SOUTH-NORTH_PATCH_START_UNSTAG' , ibuf , 1 , ierr ) |
|---|
| 416 | ibuf(1) = MIN(jpe,jde-1) |
|---|
| 417 | IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = jde - 1 |
|---|
| 418 | CALL wrf_put_dom_ti_integer ( fid , 'SOUTH-NORTH_PATCH_END_UNSTAG' , ibuf , 1 , ierr ) |
|---|
| 419 | ibuf(1) = MAX(jps,jds) |
|---|
| 420 | IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = jds |
|---|
| 421 | CALL wrf_put_dom_ti_integer ( fid , 'SOUTH-NORTH_PATCH_START_STAG' , ibuf , 1 , ierr ) |
|---|
| 422 | ibuf(1) = MIN(jpe,jde) |
|---|
| 423 | IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = jde |
|---|
| 424 | CALL wrf_put_dom_ti_integer ( fid , 'SOUTH-NORTH_PATCH_END_STAG' , ibuf , 1 , ierr ) |
|---|
| 425 | |
|---|
| 426 | ibuf(1) = MAX(kps,kds) |
|---|
| 427 | IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = kds |
|---|
| 428 | CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_PATCH_START_UNSTAG' , ibuf , 1 , ierr ) |
|---|
| 429 | ibuf(1) = MIN(kpe,kde-1) |
|---|
| 430 | IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = kde - 1 |
|---|
| 431 | CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_PATCH_END_UNSTAG' , ibuf , 1 , ierr ) |
|---|
| 432 | ibuf(1) = MAX(kps,kds) |
|---|
| 433 | IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = kds |
|---|
| 434 | CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_PATCH_START_STAG' , ibuf , 1 , ierr ) |
|---|
| 435 | ibuf(1) = MIN(kpe,kde) |
|---|
| 436 | IF ( .NOT. multi_files ( io_form ) ) ibuf(1) = kde |
|---|
| 437 | CALL wrf_put_dom_ti_integer ( fid , 'BOTTOM-TOP_PATCH_END_STAG' , ibuf , 1 , ierr ) |
|---|
| 438 | #if (EM_CORE == 1) |
|---|
| 439 | ibuf(1) = grid%id |
|---|
| 440 | CALL wrf_put_dom_ti_integer ( fid , 'GRID_ID' , ibuf , 1 , ierr ) |
|---|
| 441 | ibuf(1) = parent_id |
|---|
| 442 | CALL wrf_put_dom_ti_integer ( fid , 'PARENT_ID' , ibuf , 1 , ierr ) |
|---|
| 443 | ibuf(1) = i_parent_start |
|---|
| 444 | CALL wrf_put_dom_ti_integer ( fid , 'I_PARENT_START' , ibuf , 1 , ierr ) |
|---|
| 445 | ibuf(1) = j_parent_start |
|---|
| 446 | CALL wrf_put_dom_ti_integer ( fid , 'J_PARENT_START' , ibuf , 1 , ierr ) |
|---|
| 447 | ibuf(1) = parent_grid_ratio |
|---|
| 448 | CALL wrf_put_dom_ti_integer ( fid , 'PARENT_GRID_RATIO' , ibuf , 1 , ierr ) |
|---|
| 449 | #endif |
|---|
| 450 | |
|---|
| 451 | ! end add 010831 JM |
|---|
| 452 | |
|---|
| 453 | #if (EM_CORE != 1) |
|---|
| 454 | CALL wrf_put_dom_ti_real ( fid , 'DX' , config_flags%dx , 1 , ierr ) |
|---|
| 455 | CALL wrf_put_dom_ti_real ( fid , 'DY' , config_flags%dy , 1 , ierr ) |
|---|
| 456 | #endif |
|---|
| 457 | ! Updated by T. Hutchinson to use adaptive time step. |
|---|
| 458 | CALL wrf_put_dom_ti_real ( fid , 'DT' , grid%dt , 1 , ierr ) |
|---|
| 459 | ! CALL wrf_put_dom_ti_real ( fid , 'DT' , config%dt , 1 , ierr ) |
|---|
| 460 | CALL wrf_put_dom_ti_real ( fid , 'CEN_LAT' , config_flags%cen_lat , 1 , ierr ) |
|---|
| 461 | CALL wrf_put_dom_ti_real ( fid , 'CEN_LON' , config_flags%cen_lon , 1 , ierr ) |
|---|
| 462 | CALL wrf_put_dom_ti_real ( fid , 'TRUELAT1', config_flags%truelat1, 1 , ierr ) |
|---|
| 463 | CALL wrf_put_dom_ti_real ( fid , 'TRUELAT2', config_flags%truelat2, 1 , ierr ) |
|---|
| 464 | CALL wrf_put_dom_ti_real ( fid , 'MOAD_CEN_LAT', config_flags%moad_cen_lat, 1 , ierr ) |
|---|
| 465 | CALL wrf_put_dom_ti_real ( fid , 'STAND_LON', config_flags%stand_lon, 1 , ierr ) |
|---|
| 466 | IF ( switch .NE. boundary_only ) THEN |
|---|
| 467 | #ifdef PLANET |
|---|
| 468 | ! When writing to restart files, use the values of the instantaneous |
|---|
| 469 | ! time for determining the values of JULYR, JULDAY, and GMT. If the |
|---|
| 470 | ! original values in config_flags are used, this assumes that the |
|---|
| 471 | ! restart simulation will start with an itimestep NE 0. If we use |
|---|
| 472 | ! the instantaneous time, we can start a restart simulation with a |
|---|
| 473 | ! different value of delta-t for timestep and still get the clocks |
|---|
| 474 | ! calendars (and orbital information!) correct. |
|---|
| 475 | ! |
|---|
| 476 | ! Current time is still defined from above call to WRF_UTIL_ClockGet |
|---|
| 477 | CALL WRFU_TimeGet( currentTime, YY=julyr, dayOfYear=julday, H=p_hr, M=p_min, S=p_sec, MS=p_ms, rc=rc) |
|---|
| 478 | WRITE(wrf_err_message,*)'output_wrf: julyr,julday,hr,min,sec,ms = ',julyr,julday,p_hr,p_min,p_sec,p_ms |
|---|
| 479 | CALL wrf_debug( 100 , wrf_err_message ) |
|---|
| 480 | gmt = REAL(p_hr)+REAL(p_min)/60.+REAL(p_sec)/3600.+REAL(p_ms)/3600000. |
|---|
| 481 | CALL wrf_put_dom_ti_real ( fid , 'GMT' , gmt , 1 , ierr ) |
|---|
| 482 | CALL wrf_put_dom_ti_integer ( fid , 'JULYR' , julyr , 1 , ierr ) |
|---|
| 483 | CALL wrf_put_dom_ti_integer ( fid , 'JULDAY' , julday , 1 , ierr ) |
|---|
| 484 | #else |
|---|
| 485 | CALL wrf_put_dom_ti_real ( fid , 'GMT' , config_flags%gmt , 1 , ierr ) |
|---|
| 486 | CALL wrf_put_dom_ti_integer ( fid , 'JULYR' , config_flags%julyr , 1 , ierr ) |
|---|
| 487 | CALL wrf_put_dom_ti_integer ( fid , 'JULDAY' , config_flags%julday , 1 , ierr ) |
|---|
| 488 | #endif |
|---|
| 489 | ENDIF |
|---|
| 490 | #if (NMM_CORE == 1) |
|---|
| 491 | write(0,*) 'MMINLU would be: ', MMINLU(1:4) |
|---|
| 492 | MMINLU(1:4)='USGS' |
|---|
| 493 | write(0,*) 'MMINLU now: ', MMINLU(1:4) |
|---|
| 494 | #endif |
|---|
| 495 | CALL wrf_put_dom_ti_integer ( fid , 'MAP_PROJ' , config_flags%map_proj , 1 , ierr ) |
|---|
| 496 | CALL wrf_put_dom_ti_char ( fid , 'MMINLU', mminlu(1:4) , ierr ) |
|---|
| 497 | CALL wrf_put_dom_ti_integer ( fid , 'ISWATER' , iswater , 1 , ierr ) |
|---|
| 498 | ! CALL wrf_put_dom_ti_integer ( fid , 'ISWATER' , config_flags%iswater , 1 , ierr ) |
|---|
| 499 | CALL wrf_put_dom_ti_integer ( fid , 'ISICE' , config_flags%isice , 1 , ierr ) |
|---|
| 500 | CALL wrf_put_dom_ti_integer ( fid , 'ISURBAN' , config_flags%isurban , 1 , ierr ) |
|---|
| 501 | CALL wrf_put_dom_ti_integer ( fid , 'ISOILWATER' , config_flags%isoilwater , 1 , ierr ) |
|---|
| 502 | ! added these fields for restarting of moving nests, JM |
|---|
| 503 | CALL wrf_put_dom_ti_integer ( fid , 'I_PARENT_START' , config_flags%i_parent_start , 1 , ierr ) |
|---|
| 504 | CALL wrf_put_dom_ti_integer ( fid , 'J_PARENT_START' , config_flags%j_parent_start , 1 , ierr ) |
|---|
| 505 | |
|---|
| 506 | |
|---|
| 507 | IF ( switch .EQ. boundary_only ) THEN |
|---|
| 508 | CALL WRFU_TimeIntervalSet( bdy_increment, S=NINT(config_flags%bdyfrq),rc=rc) |
|---|
| 509 | next_time = currentTime + bdy_increment |
|---|
| 510 | CALL wrf_timetoa ( next_time, next_datestr ) |
|---|
| 511 | CALL wrf_put_dom_td_char ( fid , 'THISBDYTIME' , current_date(1:19), current_date(1:19), ierr ) |
|---|
| 512 | CALL wrf_put_dom_td_char ( fid , 'NEXTBDYTIME' , current_date(1:19), next_datestr(1:19), ierr ) |
|---|
| 513 | ENDIF |
|---|
| 514 | |
|---|
| 515 | ! added grib2-specific metadata: Todd Hutchinson 8/21/2005 |
|---|
| 516 | IF ( use_package( io_form ) == IO_GRIB2 ) THEN |
|---|
| 517 | CALL wrf_put_dom_ti_integer ( fid , 'BACKGROUND_PROC_ID' , config_flags%background_proc_id , 1 , ierr ) |
|---|
| 518 | CALL wrf_put_dom_ti_integer ( fid , 'FORECAST_PROC_ID' , config_flags%forecast_proc_id , 1 , ierr ) |
|---|
| 519 | CALL wrf_put_dom_ti_integer ( fid , 'PRODUCTION_STATUS' , config_flags%production_status , 1 , ierr ) |
|---|
| 520 | CALL wrf_put_dom_ti_integer ( fid , 'COMPRESSION' , config_flags%compression , 1 , ierr ) |
|---|
| 521 | ENDIF |
|---|
| 522 | |
|---|
| 523 | CALL nl_get_adjust_output_times( grid%id, adjust ) |
|---|
| 524 | current_date_save = current_date |
|---|
| 525 | #if 1 |
|---|
| 526 | IF ( switch .EQ. model_input_only ) THEN |
|---|
| 527 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_inputout.inc' ) |
|---|
| 528 | CALL wrf_inputout( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 529 | ELSE IF ( switch .EQ. aux_model_input1_only ) THEN |
|---|
| 530 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput1out.inc' ) |
|---|
| 531 | CALL wrf_auxinput1out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 532 | ELSE IF ( switch .EQ. aux_model_input2_only ) THEN |
|---|
| 533 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput2out.inc' ) |
|---|
| 534 | CALL wrf_auxinput2out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 535 | ELSE IF ( switch .EQ. aux_model_input3_only ) THEN |
|---|
| 536 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput3out.inc' ) |
|---|
| 537 | CALL wrf_auxinput3out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 538 | ELSE IF ( switch .EQ. aux_model_input4_only ) THEN |
|---|
| 539 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput4out.inc' ) |
|---|
| 540 | CALL wrf_auxinput4out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 541 | ELSE IF ( switch .EQ. aux_model_input5_only ) THEN |
|---|
| 542 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput5out.inc' ) |
|---|
| 543 | CALL wrf_auxinput5out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 544 | ELSE IF ( switch .EQ. aux_model_input6_only ) THEN |
|---|
| 545 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput6out.inc' ) |
|---|
| 546 | CALL wrf_auxinput6out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 547 | ELSE IF ( switch .EQ. aux_model_input7_only ) THEN |
|---|
| 548 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput7out.inc' ) |
|---|
| 549 | CALL wrf_auxinput7out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 550 | ELSE IF ( switch .EQ. aux_model_input8_only ) THEN |
|---|
| 551 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput8out.inc' ) |
|---|
| 552 | CALL wrf_auxinput8out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 553 | ELSE IF ( switch .EQ. aux_model_input9_only ) THEN |
|---|
| 554 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput9out.inc' ) |
|---|
| 555 | CALL wrf_auxinput9out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 556 | ELSE IF ( switch .EQ. aux_model_input10_only ) THEN |
|---|
| 557 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput10out.inc' ) |
|---|
| 558 | CALL wrf_auxinput10out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 559 | ELSE IF ( switch .EQ. aux_model_input11_only ) THEN |
|---|
| 560 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxinput11out.inc' ) |
|---|
| 561 | CALL wrf_auxinput11out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 562 | ELSE IF ( switch .EQ. history_only ) THEN |
|---|
| 563 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_histout.inc' ) |
|---|
| 564 | IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( HISTORY_ALARM ), currentTime, startTime, current_date ) |
|---|
| 565 | CALL wrf_histout( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 566 | ELSE IF ( switch .EQ. aux_hist1_only ) THEN |
|---|
| 567 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist1out' ) |
|---|
| 568 | IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST1_ALARM ), currentTime, startTime, current_date ) |
|---|
| 569 | CALL wrf_auxhist1out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 570 | ELSE IF ( switch .EQ. aux_hist2_only ) THEN |
|---|
| 571 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist2out.inc' ) |
|---|
| 572 | IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST2_ALARM ), currentTime, startTime, current_date ) |
|---|
| 573 | CALL wrf_auxhist2out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 574 | ELSE IF ( switch .EQ. aux_hist3_only ) THEN |
|---|
| 575 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist3out.inc' ) |
|---|
| 576 | IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST3_ALARM ), currentTime, startTime, current_date ) |
|---|
| 577 | CALL wrf_auxhist3out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 578 | ELSE IF ( switch .EQ. aux_hist4_only ) THEN |
|---|
| 579 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist4out.inc' ) |
|---|
| 580 | IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST4_ALARM ), currentTime, startTime, current_date ) |
|---|
| 581 | CALL wrf_auxhist4out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 582 | ELSE IF ( switch .EQ. aux_hist5_only ) THEN |
|---|
| 583 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist5out.inc' ) |
|---|
| 584 | IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST5_ALARM ), currentTime, startTime, current_date ) |
|---|
| 585 | CALL wrf_auxhist5out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 586 | ELSE IF ( switch .EQ. aux_hist6_only ) THEN |
|---|
| 587 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist6out.inc' ) |
|---|
| 588 | IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST6_ALARM ), currentTime, startTime, current_date ) |
|---|
| 589 | CALL wrf_auxhist6out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 590 | ELSE IF ( switch .EQ. aux_hist7_only ) THEN |
|---|
| 591 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist7out.inc' ) |
|---|
| 592 | IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST7_ALARM ), currentTime, startTime, current_date ) |
|---|
| 593 | CALL wrf_auxhist7out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 594 | ELSE IF ( switch .EQ. aux_hist8_only ) THEN |
|---|
| 595 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist8out.inc' ) |
|---|
| 596 | IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST8_ALARM ), currentTime, startTime, current_date ) |
|---|
| 597 | CALL wrf_auxhist8out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 598 | ELSE IF ( switch .EQ. aux_hist9_only ) THEN |
|---|
| 599 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist9out.inc' ) |
|---|
| 600 | IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST9_ALARM ), currentTime, startTime, current_date ) |
|---|
| 601 | CALL wrf_auxhist9out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 602 | ELSE IF ( switch .EQ. aux_hist10_only ) THEN |
|---|
| 603 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist10out.inc' ) |
|---|
| 604 | IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST10_ALARM ), currentTime, startTime, current_date ) |
|---|
| 605 | CALL wrf_auxhist10out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 606 | ELSE IF ( switch .EQ. aux_hist11_only ) THEN |
|---|
| 607 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_auxhist11out.inc' ) |
|---|
| 608 | IF ( adjust ) CALL adjust_io_timestr( grid%io_intervals( AUXHIST11_ALARM ), currentTime, startTime, current_date ) |
|---|
| 609 | CALL wrf_auxhist11out( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 610 | ELSE IF ( switch .EQ. restart_only ) THEN |
|---|
| 611 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_restartout.inc' ) |
|---|
| 612 | CALL wrf_restartout( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 613 | ELSE IF ( switch .EQ. boundary_only ) THEN |
|---|
| 614 | CALL wrf_debug ( 300 , 'output_wrf: calling code in wrf_bdyout.inc' ) |
|---|
| 615 | CALL wrf_bdyout( fid , grid , config_flags, switch, dryrun, ierr ) |
|---|
| 616 | ENDIF |
|---|
| 617 | #else |
|---|
| 618 | CALL wrf_message ( "ALL I/O DISABLED IN share/module_io_wrf.F") |
|---|
| 619 | #endif |
|---|
| 620 | current_date = current_date_save |
|---|
| 621 | |
|---|
| 622 | IF ( .NOT. dryrun ) THEN |
|---|
| 623 | CALL wrf_debug ( 300 , 'output_wrf: calling wrf_iosync ' ) |
|---|
| 624 | CALL wrf_iosync ( fid , ierr ) |
|---|
| 625 | CALL wrf_debug ( 300 , 'output_wrf: back from wrf_iosync ' ) |
|---|
| 626 | ENDIF |
|---|
| 627 | |
|---|
| 628 | WRITE(wrf_err_message,*)'output_wrf: end, fid = ',fid |
|---|
| 629 | CALL wrf_debug( 300 , wrf_err_message ) |
|---|
| 630 | |
|---|
| 631 | RETURN |
|---|
| 632 | END SUBROUTINE output_wrf |
|---|
