![(please configure the [header_logo] section in trac.ini)](/LMDZ_WRF/chrome/site/your_project_logo.png){kind=logo}
Changeset 178 in lmdz_wrf
- Timestamp:
- Oct 9, 2014, 5:25:55 PM (10 years ago)
- Location:
- trunk/WRFV3
- Files:
-
- 4 edited
-
lmdz/physiq.F90 (modified) (4 diffs)
-
lmdz/thermcell_main.F90 (modified) (8 diffs)
-
lmdz/thermcell_plume.F90 (modified) (2 diffs)
-
run/namelist.input.backup (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/WRFV3/lmdz/physiq.F90
r176 r178 2476 2476 !c================================================================================================= 2477 2477 2478 PRINT *,' Lluis end calc itap: ', itap,' it_wape_prescr: ',& 2479 it_wape_prescr,' k ale ale_wake alp alp_bl ale ale_bl _____' 2480 DO i=1, klev 2481 PRINT *, i, ale(i), ale_wake(i), alp(i),alp_bl(i), ale(i), & 2482 ale_bl(i) 2483 END DO 2478 2484 2479 2485 !c sb, oct02: … … 2791 2797 !ccc fin nrlmd le 10/04/2012 2792 2798 2799 PRINT *,' Lluis before the itap: ', itap,' it_wape_prescr: ',& 2800 it_wape_prescr,' k ale ale_wake alp alp_bl ale ale_bl _____' 2801 DO i=1, klev 2802 PRINT *, i, ale(i), ale_wake(i), alp(i),alp_bl(i), ale(i), & 2803 ale_bl(i) 2804 END DO 2805 2793 2806 if (iflag_thermals>=1) then 2794 2807 call calltherm(pdtphys & … … 2811 2824 !ccc fin nrlmd le 10/04/2012 & 2812 2825 & ,zqla,ztva ) 2826 2827 PRINT *,' Lluis after the itap: ', itap,' it_wape_prescr: ',& 2828 it_wape_prescr,' k ale ale_wake alp alp_bl ale ale_bl _____' 2829 DO i=1, klev 2830 PRINT *, i, ale(i), ale_wake(i), alp(i),alp_bl(i), ale(i), & 2831 ale_bl(i) 2832 END DO 2813 2833 2814 2834 !ccc nrlmd le 10/04/2012 … … 4285 4305 4286 4306 ! Lluis 4287 PRINT *,' Lluis: ',klev,' UBOUNDS: ',UBOUND(t_seri), UBOUND(u_seri), & 4288 UBOUND(d_q_con), UBOUND(d_t_con) 4289 PRINT *,' Lluis llp ',llp,' itap: ',itap,' zlev t_seri u_seri d_q_con d_t_con_____' 4290 DO i=1,klev 4291 PRINT *,i,t_seri(llp,i), u_seri(llp,i), d_q_con(llp,i), d_t_con(llp,i) 4292 END DO 4307 PRINT *,' Lluis: ',klev,' UBOUNDS: ',UBOUND(t_seri), & 4308 UBOUND(u_seri), UBOUND(d_q_con), UBOUND(d_t_con) 4309 PRINT *,' Lluis llp ',llp,' itap: ',itap,' zlev t_seri ' // & 4310 'u_seri d_q_con d_t_con_____' 4311 DO i=1,klev 4312 PRINT *,i,t_seri(llp,i), u_seri(llp,i), d_q_con(llp,i), & 4313 d_t_con(llp,i) 4314 END DO 4293 4315 4294 4316 -
trunk/WRFV3/lmdz/thermcell_main.F90
r142 r178 253 253 ! L. Fita, LMD July 2014. Initializing variables. 254 254 ! Some not initializated according to values: iflag_trig_bl, iflag_clos_bl 255 zdthladj = 0. 256 pbl_tke_max0 = 0. 257 fraca0 = 0. 258 w_conv = 0. 259 w0 = 0. 260 therm_tke_max0 = 0. 261 env_tke_max0 = 0. 262 alp_bl_det = 0. 263 alp_bl_fluct_m = 0. 264 alp_bl_fluct_tke = 0. 265 alp_bl_conv = 0. 266 alp_bl_stat = 0. 267 interp = 0. 268 klcl = 0 255 PRINT *,' Lluis debut thermcell_main i Alp_bl Ale_bl __' 256 DO l=1, ngrid 257 PRINT *, l, Alp_bl(l), Ale_bl(l) 258 END DO 259 260 !! zdthladj = 0. 261 !! pbl_tke_max0 = 0. 262 !! fraca0 = 0. 263 !! w_conv = 0. 264 !! w0 = 0. 265 !! therm_tke_max0 = 0. 266 !! env_tke_max0 = 0. 267 !! alp_bl_det = 0. 268 !! alp_bl_fluct_m = 0. 269 !! alp_bl_fluct_tke = 0. 270 !! alp_bl_conv = 0. 271 !! alp_bl_stat = 0. 272 !! interp = 0. 273 !! klcl = 0 269 274 270 275 !----------------------------------------------------------------------- … … 487 492 ! print*,'THERM thermcell_main iflag_thermals_ed=',iflag_thermals_ed 488 493 if (iflag_thermals_ed<=9) then 494 PRINT *,' Lluis going to thermcell_plume!!!!' 489 495 ! print*,'THERM NOUVELLE/NOUVELLE Arnaud' 490 496 CALL thermcell_plume(itap,ngrid,nlay,ptimestep,ztv,zthl,po,zl,rhobarz,& … … 495 501 496 502 elseif (iflag_thermals_ed>9) then 503 PRINT *,' Lluis going to thermcellV1_plume!!!!' 497 504 ! print*,'THERM RIO et al 2010, version d Arnaud' 498 505 CALL thermcellV1_plume(itap,ngrid,nlay,ptimestep,ztv,zthl,po,zl,rhobarz,& … … 503 510 504 511 endif 512 PRINT *,' Lluis in thermcell_main iflag_thermals_ed: ', & 513 iflag_thermals_ed, ' ig=528 ',zw2(528,:) 505 514 506 515 if (prt_level.ge.1) print*,'apres thermcell_plume ',lev_out … … 601 610 & detr,zqla,lev_out,lunout1,igout) 602 611 !IM 060508 & detr,zqla,zmax,lev_out,lunout,igout) 612 PRINT *,' Lluis after thermcell_flux2 ig=528 zw2: ',zw2(528,:) 603 613 604 614 if (prt_level.ge.1) print*,'thermcell_main apres thermcell_flux' … … 745 755 if (prt_level.ge.10)write(lunout,*) & 746 756 & 'WARNING thermcell_main wth2=0. si zw2 > 1.e-10' 757 PRINT *,' Lluis in thermcell_main ig=528 l wth3 zf2 fraca zw2 ___' 747 758 do l=1,nlay 748 759 do ig=1,ngrid … … 758 769 wth3(ig,l)=zf2*(1-2.*fraca(ig,l))/(1-fraca(ig,l)) & 759 770 & *zw2(ig,l)*zw2(ig,l)*zw2(ig,l) 771 IF (ig == 528) PRINT *, l, wth3(ig,l), zf2, & 772 fraca(ig,l), zw2(ig,l) 773 ! PRINT *, ig, l, wth3(ig,l), zf2, & 774 ! fraca(ig,l), zw2(ig,l) 760 775 q2(ig,l)=zf2*(zqta(ig,l)*1000.-po(ig,l)*1000.)**2 761 776 !test: on calcul q2/po=ratqsc … … 1001 1016 enddo 1002 1017 enddo 1018 PRINT *,' Lluis in thermcell_main i Alp_bl Srhobarz Swth3 ' //& 1019 'Ale_bl Szw2 ______', alp_bl_k 1020 DO l=1, ngrid 1021 PRINT *, l, Alp_bl(l), SUM(rhobarz(l,:)), SUM(wth3(l,:)), & 1022 Ale_bl(l), SUM(zw2(l,:)) 1023 END DO 1003 1024 1004 1025 !test:calcul de la ponderation des couches pour KE -
trunk/WRFV3/lmdz/thermcell_plume.F90
r142 r178 191 191 endif 192 192 enddo 193 PRINT *,' Lluis in themcell_plume ig actvie l zw2 ztv zlev pphi ___' 194 DO ig=1,ngrid 195 PRINT *,ig, active(ig), 1, zw2(ig,1), ztv(ig,1), zlev(ig,1), pphi(ig,1), & 196 2, zw2(ig,2), ztv(ig,2), zlev(ig,2), pphi(ig,2) 197 END DO 193 198 ! 194 199 !============================================================================== … … 645 650 endif 646 651 enddo 652 PRINT *,' Lluis in themcellV1_plume ig actvie l zw2 ztv zlev pphi ___' 653 DO ig=1,ngrid 654 PRINT *,ig, active(ig), 1, zw2(ig,1), ztv(ig,1), zlev(ig,1), pphi(ig,1), & 655 2, zw2(ig,2), ztv(ig,2), zlev(ig,2), pphi(ig,2) 656 END DO 647 657 ! 648 658 -
trunk/WRFV3/run/namelist.input.backup
r175 r178 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.
