![(please configure the [header_logo] section in trac.ini)](/LMDZ_WRF/chrome/site/your_project_logo.png){kind=logo}
Changeset 43 in lmdz_wrf
- Timestamp:
- Jul 23, 2014, 3:06:33 PM (10 years ago)
- Location:
- branches/LMDZ_WRFmeas/WRFV3
- Files:
-
- 3 edited
-
dyn_em/solve_em.F (modified) (1 diff)
-
lmdz/physiq.F90 (modified) (1 diff)
-
run/namelist.input.backup (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/LMDZ_WRFmeas/WRFV3/dyn_em/solve_em.F
r26 r43 859 859 im2 = ims + (ime - ims) / 2 860 860 jm2 = jms + (jme - jms) / 2 861 862 IF (grid%t_2(im2,iz,jm2) /= grid%t_2(im2,iz,jm2) .OR. ABS(grid%t_2(im2,iz,jm2)) > 10000. ) THEN 863 PRINT *,TRIM(errmsg) 864 WRITE(wrF_err_message,*)'solve_em: wrong T value=', & 865 grid%t_2(im2,iz,jm2),' at: ', im2,', ', iz,', ', jm2,' !!!' 861 DO iz = kms, kme 862 863 IF (grid%t_2(im2,iz,jm2) /= grid%t_2(im2,iz,jm2) .OR. ABS(grid%t_2(im2,iz,jm2)) > 10000. ) THEN 864 PRINT *,TRIM(errmsg) 865 WRITE(wrF_err_message,*)'solve_em: wrong T value=', & 866 grid%t_2(im2,iz,jm2),' at: ', im2,', ', iz,', ', jm2,' !!!' 866 867 #ifdef DM_PARALLEL 867 CALL wrf_error_fatal(TRIM(wrf_err_message))868 CALL wrf_error_fatal(TRIM(wrf_err_message)) 868 869 #else 869 PRINT *,TRIM(wrf_err_message) 870 STOP 871 #endif 872 END IF 870 PRINT *,TRIM(wrf_err_message) 871 STOP 872 #endif 873 END IF 874 END DO 873 875 ! Checking for NaNs (should not be necessary but....) 874 876 IF (grid%psfc(im2,jm2) /= grid%psfc(im2,jm2) .OR. ABS(grid%psfc(im2,jm2)) > 1000000. ) THEN -
branches/LMDZ_WRFmeas/WRFV3/lmdz/physiq.F90
r6 r43 2168 2168 zdtime=dtime*REAL(radpas) ! pas de temps du rayonnement (s) 2169 2169 CALL zenang(zlongi,jH_cur,zdtime,rlat,rlon,rmu0,fract) 2170 PRINT *,' Lluis: radpas: ',radpas,' fract: ',fract(llp), &2171 ' zdtime: ',zdtime,' jH_cur: ',jH_cur,' rmu0: ',rmu0(llp),&2172 ' dtime: ',dtime,' nbapp_rad: ',nbapp_rad2173 PRINT *,' Lluis @rlat@: ',rlat2174 PRINT *,' Lluis @rlon@: ',rlon2175 PRINT *,' Lluis @rmu0@: ',rmu02176 PRINT *,' Lluis @fract@: ',fract2177 2170 ELSE 2178 2171 CALL angle(zlongi, rlat, fract, rmu0) 2179 2172 ENDIF 2180 2173 ENDIF 2181 PRINT *,' Lluis zenang_an rlat: ',rlat(llp), ' rlon:',rlon(llp), &2182 ' rmu0: ',rmu0(llp),' fract: ',fract(llp),' jH_cur: ',jH_cur, &2183 ' cycle_diurne: ',cycle_diurne2184 2174 2185 2175 if (mydebug) then -
branches/LMDZ_WRFmeas/WRFV3/run/namelist.input.backup
r28 r43 1 1 &time_control 2 2 run_days = 0, 3 run_hours = 0,3 run_hours = 12, 4 4 run_minutes = 0, 5 5 run_seconds = 0, 6 start_year = 1995,7 start_month = 01, 8 start_day = 13,9 start_hour = 00,10 start_minute = 00, 11 start_second = 00, 12 end_year = 1995,13 end_month = 01, 14 end_day = 18,15 end_hour = 00,16 end_minute = 00, 17 end_second = 00, 18 history_interval = 30,19 in terval_seconds = 2160020 input_from_file = .true.,21 frames_per_outfile = 240,6 start_year = 2000, 2000, 2000, 7 start_month = 01, 01, 01, 8 start_day = 24, 24, 24, 9 start_hour = 12, 12, 12, 10 start_minute = 00, 00, 00, 11 start_second = 00, 00, 00, 12 end_year = 2000, 2000, 2000, 13 end_month = 01, 01, 01, 14 end_day = 25, 25, 25, 15 end_hour = 12, 12, 12, 16 end_minute = 00, 00, 00, 17 end_second = 00, 00, 00, 18 interval_seconds = 21600 19 input_from_file = .true.,.true.,.true., 20 history_interval = 180, 60, 60, 21 frames_per_outfile = 1000, 1000, 1000, 22 22 restart = .false., 23 restart_interval _h = 120,23 restart_interval = 5000, 24 24 io_form_history = 2 25 25 io_form_restart = 2 26 26 io_form_input = 2 27 27 io_form_boundary = 2 28 auxinput4_inname = 'wrflowinp_d<domain>'29 auxinput4_interval = 360,30 io_form_auxinput4 = 231 override_restart_timers = .true.32 28 debug_level = 0 33 29 / 34 30 35 31 &domains 36 time_step = 180 37 time_step_fract_num = 0 038 time_step_fract_den = 1 32 time_step = 180, 33 time_step_fract_num = 0, 34 time_step_fract_den = 1, 39 35 max_dom = 1, 40 e_we = 80,41 e_sn = 90,42 e_vert = 39,43 p_top_requested = 1000,44 num_metgrid_levels = 38,36 e_we = 74, 112, 94, 37 e_sn = 61, 97, 91, 38 e_vert = 28, 28, 28, 39 p_top_requested = 5000, 40 num_metgrid_levels = 27, 45 41 num_metgrid_soil_levels = 4, 46 dx = 20000.,47 dy = 20000.,48 grid_id = 1, 49 parent_id = 0, 50 i_parent_start = 1, 51 j_parent_start = 1, 52 parent_grid_ratio = 1, 53 parent_time_step_ratio = 1, 42 dx = 30000, 10000, 3333.33, 43 dy = 30000, 10000, 3333.33, 44 grid_id = 1, 2, 3, 45 parent_id = 0, 1, 2, 46 i_parent_start = 1, 31, 30, 47 j_parent_start = 1, 17, 30, 48 parent_grid_ratio = 1, 3, 3, 49 parent_time_step_ratio = 1, 3, 3, 54 50 feedback = 1, 55 51 smooth_option = 0 56 bottom_top_plus1 = 39,57 max_lidar_locs = 20,58 52 / 59 53 60 54 &physics 61 mp_physics = 0,62 ra_lw_physics = 0,63 ra_sw_physics = 0,64 radt = 120,65 sf_sfclay_physics = 0,66 sf_surface_physics = 0,67 bl_pbl_physics = 0,68 bldt = 0, 69 cu_physics = 0,70 cudt = 0,55 mp_physics = 3, 3, 3, 56 ra_lw_physics = 1, 1, 1, 57 ra_sw_physics = 1, 1, 1, 58 radt = 30, 30, 30, 59 sf_sfclay_physics = 1, 1, 1, 60 sf_surface_physics = 2, 2, 2, 61 bl_pbl_physics = 1, 1, 1, 62 bldt = 0, 0, 0, 63 cu_physics = 1, 1, 0, 64 cudt = 5, 5, 5, 71 65 isfflx = 1, 72 66 ifsnow = 0, 73 67 icloud = 1, 74 68 surface_input_source = 1, 75 num_soil_layers = 5, 76 mp_zero_out = 0, 77 sst_update = 1 69 num_soil_layers = 4, 70 sf_urban_physics = 0, 0, 0, 78 71 / 79 72 … … 82 75 83 76 &dynamics 84 diff_opt = 0, 85 km_opt = 0, 86 damp_opt = 3, 77 w_damping = 0, 78 diff_opt = 1, 79 km_opt = 4, 80 diff_6th_opt = 0, 0, 0, 81 diff_6th_factor = 0.12, 0.12, 0.12, 87 82 base_temp = 290. 88 zdamp = 5000., 89 dampcoef = 0.2, 90 khdif = 0, 91 kvdif = 0, 92 non_hydrostatic = .true., 93 moist_adv_opt = 0, 94 scalar_adv_opt = 0, 95 chem_adv_opt = 0, 96 tke_adv_opt = 0, 97 fft_filter_lat = 75., 98 w_damping = 1, 83 damp_opt = 0, 84 zdamp = 5000., 5000., 5000., 85 dampcoef = 0.2, 0.2, 0.2 86 khdif = 0, 0, 0, 87 kvdif = 0, 0, 0, 88 non_hydrostatic = .true., .true., .true., 89 moist_adv_opt = 1, 1, 1, 90 scalar_adv_opt = 1, 1, 1, 99 91 / 100 92 101 93 &bdy_control 102 spec_bdy_width= 5,103 spec_zone= 1,104 relax_zone= 4,105 specified = .true.,106 nested = .false.,94 spec_bdy_width = 5, 95 spec_zone = 1, 96 relax_zone = 4, 97 specified = .true., .false.,.false., 98 nested = .false., .true., .true., 107 99 / 108 100
Note: See TracChangeset
for help on using the changeset viewer.
