MODULE module_diagnostics_driver CONTAINS SUBROUTINE diagnostics_driver ( grid, config_flags, & moist, chem, tracer, scalar, & th_phy, pi_phy, p_phy, & p8w, t8w, dz8w, & curr_secs, curr_secs2, & diag_flag, & ids, ide, jds, jde, kds, kde, & ims, ime, jms, jme, kms, kme, & ips, ipe, jps, jpe, kps, kpe, & imsx, imex, jmsx, jmex, kmsx, kmex, & ipsx, ipex, jpsx, jpex, kpsx, kpex, & imsy, imey, jmsy, jmey, kmsy, kmey, & ipsy, ipey, jpsy, jpey, kpsy, kpey ) USE module_state_description, ONLY: num_moist, num_chem, num_tracer, num_scalar, & SKIP_PRESS_DIAGS, SKIP_Z_DIAGS, & DO_TRAD_FIELDS, & DO_SOLAR_OUTPUT, & PARAM_FIRST_SCALAR, & P_QG, P_QH, P_QV, P_QC, P_QI, P_QS, & P_QNG, P_QH, P_QNH, P_QR, P_QNR, P_QVOLG, & F_QV, F_QC, F_QI, F_QS, & KESSLERSCHEME, LINSCHEME, SBU_YLINSCHEME, WSM3SCHEME, WSM5SCHEME, & WSM6SCHEME, ETAMPNEW, THOMPSON, THOMPSONAERO, THOMPSONGH, & MORR_TWO_MOMENT, GSFCGCESCHEME, WDM5SCHEME, WDM6SCHEME, & MILBRANDT2MOM , CAMMGMPSCHEME, FAST_KHAIN_LYNN_SHPUND, FULL_KHAIN_LYNN, & NSSL_2MOM, MORR_TM_AERO USE module_driver_constants, ONLY: max_plevs, max_zlevs USE module_model_constants, ONLY: g, R_v, R_d, Cp, T0, RCP USE module_domain, ONLY : domain, domain_clock_get, domain_get_current_time USE module_configure, ONLY : grid_config_rec_type, & model_config_rec USE module_streams USE module_utility, ONLY : WRFU_Time USE module_lightning_driver, ONLY : lightning_driver USE module_diag_misc, ONLY : diagnostic_output_calc, noahmp_output_calc USE module_diag_nwp, ONLY : diagnostic_output_nwp USE module_diag_cl, ONLY : clwrf_output_calc USE module_diag_pld, ONLY : pld USE module_diag_zld, ONLY : zld USE module_diag_afwa, ONLY : afwa_diagnostics_driver USE module_diag_rasm, ONLY : mean_output_calc, diurnalcycle_output_calc USE module_diag_hailcast, ONLY : hailcast_diagnostic_driver USE module_diag_solar, ONLY : solar_diag USE module_trad_fields, ONLY : trad_fields IMPLICIT NONE TYPE ( domain ), INTENT(INOUT) :: grid TYPE ( grid_config_rec_type ), INTENT(IN) :: config_flags REAL , DIMENSION(ims:ime,kms:kme,jms:jme,num_moist ) , INTENT(IN) :: moist REAL , DIMENSION(ims:ime,kms:kme,jms:jme,num_chem ) , INTENT(IN) :: chem REAL , DIMENSION(ims:ime,kms:kme,jms:jme,num_tracer) , INTENT(IN) :: tracer REAL , DIMENSION(ims:ime,kms:kme,jms:jme,num_scalar) , INTENT(IN) :: scalar REAL , DIMENSION(ims:ime,kms:kme,jms:jme) , INTENT(IN) :: th_phy , & p_phy , & pi_phy , & dz8w , & p8w , & t8w REAL(8) :: curr_secs REAL :: curr_secs2 LOGICAL :: diag_flag INTEGER , INTENT(IN) :: ids, ide, jds, jde, kds, kde, & ims, ime, jms, jme, kms, kme, & ips, ipe, jps, jpe, kps, kpe INTEGER , INTENT(IN) :: imsx,imex,jmsx,jmex,kmsx,kmex, & ipsx,ipex,jpsx,jpex,kpsx,kpex, & imsy,imey,jmsy,jmey,kmsy,kmey, & ipsy,ipey,jpsy,jpey,kpsy,kpey CHARACTER (LEN=1000) :: diag_message INTEGER :: ij INTEGER :: k_start, k_end TYPE(WRFU_Time) :: currentTime CALL wrf_debug ( 100 , '--> TOP OF DIAGNOSTICS PACKAGE' ) k_start = kps k_end = kpe-1 CALL update_phys_fields ( grid, config_flags, moist, & ids, ide, jds, jde, kds, kde, & ims, ime, jms, jme, kms, kme, & ips, ipe, jps, jpe, kps, kpe ) LIGHTNING: IF ( config_flags%lightning_option /= 0 ) THEN CALL wrf_debug ( 100 , '--> CALL DIAGNOSTICS PACKAGE: LIGHTNING_DRIVER' ) CALL lightning_driver ( & curr_secs, grid%dt, grid%dx, grid%dy, & grid%xlat, grid%xlong, grid%xland, grid%ht, & grid%t_phy, p_phy, grid%rho, & grid%u_phy, grid%v_phy, grid%w_2, & th_phy, pi_phy,dz8w, & grid%z, moist, & grid%ktop_deep, grid%refl_10cm, & domain_get_current_time( grid ), & config_flags%lightning_option, & config_flags%lightning_dt, & config_flags%lightning_start_seconds, & grid%ltngacttime, & config_flags%flashrate_factor, & config_flags%iccg_method, & config_flags%iccg_prescribed_num, & config_flags%iccg_prescribed_den, & grid%iccg_in_num, grid%iccg_in_den, & config_flags%cellcount_method, & config_flags%cldtop_adjustment, & ids, ide, jds, jde, kds, kde, & ims, ime, jms, jme, kms, kme, & ips, ipe, jps, jpe, kps, kpe, & grid%ic_flashcount, grid%ic_flashrate, & grid%cg_flashcount, grid%cg_flashrate, & grid%lpi & ) END IF LIGHTNING HAILCAST: IF ( config_flags%hailcast_opt /= 0 ) THEN IF ( ( config_flags%history_interval == 0 ) ) THEN WRITE (diag_message , * ) & "HAILCAST Error : No 'history_interval' defined in namelist" CALL wrf_error_fatal3("",241,& diag_message ) END IF !$OMP PARALLEL DO & !$OMP PRIVATE ( ij ) DO ij = 1 , grid%num_tiles CALL wrf_debug ( 100 , & '--> CALL DIAGNOSTICS PACKAGE: HAILCAST_DIAGNOSTIC_DRIVER' ) CALL hailcast_diagnostic_driver ( grid , config_flags & ,moist, grid%rho & ,ids, ide, jds, jde, kds, kde & ,ims, ime, jms, jme, kms, kme & ,ips, ipe, jps, jpe, kps, kpe & ,ITS=grid%i_start(ij),ITE=grid%i_end(ij) & ,JTS=grid%j_start(ij),JTE=grid%j_end(ij) & ,K_START=k_start,K_END=k_end & ,dt=grid%dt, itimestep=grid%itimestep & ,haildt=grid%haildt, curr_secs=curr_secs2 & ,haildtacttime=grid%haildtacttime ) END DO !$OMP END PARALLEL DO END IF HAILCAST TRADITIONAL_FIELDS: IF ( config_flags%diag_nwp2 == do_trad_fields ) THEN !$OMP PARALLEL DO & !$OMP PRIVATE ( ij ) DO ij = 1 , grid%num_tiles CALL wrf_debug ( 100 , '--> CALL DIAGNOSTICS PACKAGE: TRAD_FIELDS' ) CALL trad_fields ( & U=grid%u_2 & ,V=grid%v_2 & ,W=grid%w_2 & ,t=grid%th_phy_m_t0 & ,qv=moist(:,:,:,P_QV) & ,zp=grid%ph_2 & ,zb=grid%phb & ,pp=grid%p & ,pb=grid%pb & ,p=grid%p_hyd & ,pw=grid%p_hyd_w & ,msfux=grid%msfux & ,msfuy=grid%msfuy & ,msfvx=grid%msfvx & ,msfvy=grid%msfvy & ,msftx=grid%msftx & ,msfty=grid%msfty & ,f=grid%f & ,e=grid%e & ,sina=grid%sina & ,cosa=grid%cosa & ,rho=grid%rho & ,dz8w=dz8w & ,qc=moist(:,:,:,P_QC) & ,rainc=grid%rainc & ,rainnc=grid%rainnc & ,snownc=grid%snownc & ,graupelnc=grid%graupelnc & ,hailnc=grid%hailnc & ,ht=grid%ht & ,use_theta_m=config_flags%use_theta_m & ,sealevelp=grid%sealevelp & ,temperature=grid%temperature & ,pressure=grid%pressure & ,geoheight=grid%geoheight & ,umet=grid%umet & ,vmet=grid%vmet & ,speed=grid%speed & ,dir=grid%dir & ,psfc=grid%psfc & ,rain=grid%rain & ,liqrain=grid%liqrain & ,tpw=grid%tpw & ,potential_t=grid%potential_t & ,rh=grid%rh & ,IDS=ids,IDE=ide, JDS=jds,JDE=jde, KDS=kds,KDE=kde & ,IMS=ims,IME=ime, JMS=jms,JME=jme, KMS=kms,KME=kme & ,IPS=ips,IPE=ipe, JPS=jps,JPE=jpe, KPS=kps,KPE=kpe & ,ITS=grid%i_start(ij),ITE=grid%i_end(ij) & ,JTS=grid%j_start(ij),JTE=grid%j_end(ij) & ,kts=kps,kte=kpe ) END DO !$OMP END PARALLEL DO END IF TRADITIONAL_FIELDS CALL wrf_debug ( 100 , '--> CALL DIAGNOSTICS PACKAGE: ACCUMULATED AND BUCKET DIAGNOSTICS' ) CALL diagnostic_output_calc( & DPSDT=grid%dpsdt ,DMUDT=grid%dmudt & ,P8W=p8w ,PK1M=grid%pk1m & ,MU_2=grid%mu_2 ,MU_2M=grid%mu_2m & ,RAINCV=grid%raincv ,RAINNCV=grid%rainncv & ,RAINC=grid%rainc ,RAINNC=grid%rainnc & ,I_RAINC=grid%i_rainc ,I_RAINNC=grid%i_rainnc & ,HFX=grid%hfx ,SFCEVP=grid%sfcevp ,LH=grid%lh & ,T2=grid%t2 & ,DT=grid%dt & ,XTIME=grid%xtime & ,ACSWUPT=grid%acswupt ,ACSWUPTC=grid%acswuptc & ,ACSWDNT=grid%acswdnt ,ACSWDNTC=grid%acswdntc & ,ACSWUPB=grid%acswupb ,ACSWUPBC=grid%acswupbc & ,ACSWDNB=grid%acswdnb ,ACSWDNBC=grid%acswdnbc & ,ACLWUPT=grid%aclwupt ,ACLWUPTC=grid%aclwuptc & ,ACLWDNT=grid%aclwdnt ,ACLWDNTC=grid%aclwdntc & ,ACLWUPB=grid%aclwupb ,ACLWUPBC=grid%aclwupbc & ,ACLWDNB=grid%aclwdnb ,ACLWDNBC=grid%aclwdnbc & ,I_ACSWUPT=grid%i_acswupt ,I_ACSWUPTC=grid%i_acswuptc & ,I_ACSWDNT=grid%i_acswdnt ,I_ACSWDNTC=grid%i_acswdntc & ,I_ACSWUPB=grid%i_acswupb ,I_ACSWUPBC=grid%i_acswupbc & ,I_ACSWDNB=grid%i_acswdnb ,I_ACSWDNBC=grid%i_acswdnbc & ,I_ACLWUPT=grid%i_aclwupt ,I_ACLWUPTC=grid%i_aclwuptc & ,I_ACLWDNT=grid%i_aclwdnt ,I_ACLWDNTC=grid%i_aclwdntc & ,I_ACLWUPB=grid%i_aclwupb ,I_ACLWUPBC=grid%i_aclwupbc & ,I_ACLWDNB=grid%i_aclwdnb ,I_ACLWDNBC=grid%i_aclwdnbc & ,athmp=grid%athmpten ,aqvmp=grid%aqvmpten & ,athcu=grid%athcuten ,aqvcu=grid%aqvcuten & ,aucu=grid%aucuten ,avcu=grid%avcuten & ,athsh=grid%athshten ,aqvsh=grid%aqvshten & ,aush=grid%aushten ,avsh=grid%avshten & ,athbl=grid%athblten ,aqvbl=grid%aqvblten & ,aubl=grid%aublten ,avbl=grid%avblten & ,athralw=grid%athratenlw ,athrasw=grid%athratensw & ,h_diabatic=grid%h_diabatic ,qv_diabatic=grid%qv_diabatic & ,rthcuten=grid%rthcuten ,rqvcuten=grid%rqvcuten & ,rucuten=grid%rucuten ,rvcuten=grid%rvcuten & ,rthshten=grid%rthshten ,rqvshten=grid%rqvshten & ,rushten=grid%rushten ,rvshten=grid%rvshten & ,rthblten=grid%rthblten ,rqvblten=grid%rqvblten & ,rublten=grid%rublten ,rvblten=grid%rvblten & ,rthratenlw=grid%rthratenlw ,rthratensw=grid%rthratensw & ,DIAG_PRINT=config_flags%diag_print & ,BUCKET_MM=config_flags%bucket_mm & ,BUCKET_J =config_flags%bucket_J & ,SNOWNCV=grid%snowncv, SNOW_ACC_NC=grid%snow_acc_nc & ,PREC_ACC_C=grid%prec_acc_c & ,PREC_ACC_NC=grid%prec_acc_nc & ,PREC_ACC_DT=config_flags%prec_acc_dt & ,CURR_SECS2=curr_secs2 & ,HISTORY_INTERVAL=grid%history_interval & ,ITIMESTEP=grid%itimestep & ,cu_used=grid%cu_used & ,shcu_used=grid%shcu_used & ,acc_phy_tend=grid%acc_phy_tend & ,IDS=ids,IDE=ide, JDS=jds,JDE=jde, KDS=kds,KDE=kde & ,IMS=ims,IME=ime, JMS=jms,JME=jme, KMS=kms,KME=kme & ,IPS=ips,IPE=ipe, JPS=jps,JPE=jpe, KPS=kps,KPE=kpe & ,I_START=grid%i_start,I_END=min(grid%i_end, ide-1) & ,J_START=grid%j_start,J_END=min(grid%j_end, jde-1) & ,KTS=k_start, KTE=min(k_end,kde-1) & ,NUM_TILES=grid%num_tiles & ) NWPDIAGS: IF ( config_flags%nwp_diagnostics == 1 ) THEN CALL wrf_debug ( 100 , '--> CALL DIAGNOSTICS PACKAGE: NWP DIAGNOSTICS' ) mp_select: SELECT CASE(config_flags%mp_physics) CASE (LINSCHEME, WSM6SCHEME, WDM6SCHEME, GSFCGCESCHEME) CALL diagnostic_output_nwp( & config_flags=config_flags, & U=grid%u_2 ,V=grid%v_2 & ,TEMP=grid%t_phy ,P8W=p8w & ,DT=grid%dt ,SBW=config_flags%spec_bdy_width & ,XTIME=grid%xtime & ,MPHYSICS_OPT=config_flags%mp_physics & ,GSFCGCE_HAIL=config_flags%gsfcgce_hail & ,GSFCGCE_2ICE=config_flags%gsfcgce_2ice & ,MPUSE_HAIL=config_flags%hail_opt & ,NSSL_ALPHAH=config_flags%nssl_alphah & ,NSSL_ALPHAHL=config_flags%nssl_alphahl & ,NSSL_CNOH=config_flags%nssl_cnoh & ,NSSL_CNOHL=config_flags%nssl_cnohl & ,NSSL_RHO_QH=config_flags%nssl_rho_qh & ,NSSL_RHO_QHL=config_flags%nssl_rho_qhl & ,CURR_SECS2=curr_secs2 & ,NWP_DIAGNOSTICS=config_flags%nwp_diagnostics & ,DIAGFLAG=diag_flag & ,HISTORY_INTERVAL=grid%history_interval & ,ITIMESTEP=grid%itimestep & ,U10=grid%u10,V10=grid%v10,W=grid%w_2 & ,WSPD10MAX=grid%wspd10max & ,UP_HELI_MAX=grid%up_heli_max & ,W_UP_MAX=grid%w_up_max,W_DN_MAX=grid%w_dn_max & ,ZNW=grid%znw,W_COLMEAN=grid%w_colmean & ,NUMCOLPTS=grid%numcolpts,W_MEAN=grid%w_mean & ,GRPL_MAX=grid%grpl_max,GRPL_COLINT=grid%grpl_colint & ,REFD_MAX=grid%refd_max & ,refl_10cm=grid%refl_10cm & ,HAIL_MAXK1=grid%hail_maxk1,HAIL_MAX2D=grid%hail_max2d & ,QG_CURR=moist(ims,kms,jms,P_QG) & ,RHO=grid%rho,PH=grid%ph_2,PHB=grid%phb,G=g & ,IDS=ids,IDE=ide, JDS=jds,JDE=jde, KDS=kds,KDE=kde & ,IMS=ims,IME=ime, JMS=jms,JME=jme, KMS=kms,KME=kme & ,IPS=ips,IPE=ipe, JPS=jps,JPE=jpe, KPS=kps,KPE=kpe & ,I_START=grid%i_start,I_END=min(grid%i_end, ide-1) & ,J_START=grid%j_start,J_END=min(grid%j_end, jde-1) & ,KTS=k_start, KTE=min(k_end,kde-1) & ,NUM_TILES=grid%num_tiles & ,MAX_TIME_STEP=grid%max_time_step & ,ADAPTIVE_TS=config_flags%use_adaptive_time_step & ) CASE (THOMPSON, THOMPSONAERO) CALL diagnostic_output_nwp( & config_flags=config_flags, & U=grid%u_2 ,V=grid%v_2 & ,TEMP=grid%t_phy ,P8W=p8w & ,DT=grid%dt ,SBW=config_flags%spec_bdy_width & ,XTIME=grid%xtime & ,MPHYSICS_OPT=config_flags%mp_physics & ,GSFCGCE_HAIL=config_flags%gsfcgce_hail & ,GSFCGCE_2ICE=config_flags%gsfcgce_2ice & ,MPUSE_HAIL=config_flags%hail_opt & ,NSSL_ALPHAH=config_flags%nssl_alphah & ,NSSL_ALPHAHL=config_flags%nssl_alphahl & ,NSSL_CNOH=config_flags%nssl_cnoh & ,NSSL_CNOHL=config_flags%nssl_cnohl & ,NSSL_RHO_QH=config_flags%nssl_rho_qh & ,NSSL_RHO_QHL=config_flags%nssl_rho_qhl & ,CURR_SECS2=curr_secs2 & ,NWP_DIAGNOSTICS=config_flags%nwp_diagnostics & ,DIAGFLAG=diag_flag & ,HISTORY_INTERVAL=grid%history_interval & ,ITIMESTEP=grid%itimestep & ,U10=grid%u10,V10=grid%v10,W=grid%w_2 & ,WSPD10MAX=grid%wspd10max & ,UP_HELI_MAX=grid%up_heli_max & ,W_UP_MAX=grid%w_up_max,W_DN_MAX=grid%w_dn_max & ,ZNW=grid%znw,W_COLMEAN=grid%w_colmean & ,NUMCOLPTS=grid%numcolpts,W_MEAN=grid%w_mean & ,GRPL_MAX=grid%grpl_max,GRPL_COLINT=grid%grpl_colint & ,REFD_MAX=grid%refd_max & ,refl_10cm=grid%refl_10cm & ,HAIL_MAXK1=grid%hail_maxk1,HAIL_MAX2D=grid%hail_max2d & ,QG_CURR=moist(ims,kms,jms,P_QG) & ,QR_CURR=moist(ims,kms,jms,P_QR) & ,NR_CURR=scalar(ims,kms,jms,P_QNR) & ,RHO=grid%rho,PH=grid%ph_2,PHB=grid%phb,G=g & ,IDS=ids,IDE=ide, JDS=jds,JDE=jde, KDS=kds,KDE=kde & ,IMS=ims,IME=ime, JMS=jms,JME=jme, KMS=kms,KME=kme & ,IPS=ips,IPE=ipe, JPS=jps,JPE=jpe, KPS=kps,KPE=kpe & ,I_START=grid%i_start,I_END=min(grid%i_end, ide-1) & ,J_START=grid%j_start,J_END=min(grid%j_end, jde-1) & ,KTS=k_start, KTE=min(k_end,kde-1) & ,NUM_TILES=grid%num_tiles & ,MAX_TIME_STEP=grid%max_time_step & ,ADAPTIVE_TS=config_flags%use_adaptive_time_step & ) CASE (THOMPSONGH) CALL diagnostic_output_nwp( & config_flags=config_flags, & U=grid%u_2 ,V=grid%v_2 & ,TEMP=grid%t_phy ,P8W=p8w & ,DT=grid%dt ,SBW=config_flags%spec_bdy_width & ,XTIME=grid%xtime & ,MPHYSICS_OPT=config_flags%mp_physics & ,GSFCGCE_HAIL=config_flags%gsfcgce_hail & ,GSFCGCE_2ICE=config_flags%gsfcgce_2ice & ,MPUSE_HAIL=config_flags%hail_opt & ,NSSL_ALPHAH=config_flags%nssl_alphah & ,NSSL_ALPHAHL=config_flags%nssl_alphahl & ,NSSL_CNOH=config_flags%nssl_cnoh & ,NSSL_CNOHL=config_flags%nssl_cnohl & ,NSSL_RHO_QH=config_flags%nssl_rho_qh & ,NSSL_RHO_QHL=config_flags%nssl_rho_qhl & ,CURR_SECS2=curr_secs2 & ,NWP_DIAGNOSTICS=config_flags%nwp_diagnostics & ,DIAGFLAG=diag_flag & ,HISTORY_INTERVAL=grid%history_interval & ,ITIMESTEP=grid%itimestep & ,U10=grid%u10,V10=grid%v10,W=grid%w_2 & ,WSPD10MAX=grid%wspd10max & ,UP_HELI_MAX=grid%up_heli_max & ,W_UP_MAX=grid%w_up_max,W_DN_MAX=grid%w_dn_max & ,ZNW=grid%znw,W_COLMEAN=grid%w_colmean & ,NUMCOLPTS=grid%numcolpts,W_MEAN=grid%w_mean & ,GRPL_MAX=grid%grpl_max,GRPL_COLINT=grid%grpl_colint & ,REFD_MAX=grid%refd_max & ,refl_10cm=grid%refl_10cm & ,HAIL_MAXK1=grid%hail_maxk1,HAIL_MAX2D=grid%hail_max2d & ,QG_CURR=moist(ims,kms,jms,P_QG) & ,QR_CURR=moist(ims,kms,jms,P_QR) & ,NR_CURR=scalar(ims,kms,jms,P_QNR) & ,NG_CURR=scalar(ims,kms,jms,P_QNG) & ,QVOLG_CURR=scalar(ims,kms,jms,P_QVOLG) & ,RHO=grid%rho,PH=grid%ph_2,PHB=grid%phb,G=g & ,IDS=ids,IDE=ide, JDS=jds,JDE=jde, KDS=kds,KDE=kde & ,IMS=ims,IME=ime, JMS=jms,JME=jme, KMS=kms,KME=kme & ,IPS=ips,IPE=ipe, JPS=jps,JPE=jpe, KPS=kps,KPE=kpe & ,I_START=grid%i_start,I_END=min(grid%i_end, ide-1) & ,J_START=grid%j_start,J_END=min(grid%j_end, jde-1) & ,KTS=k_start, KTE=min(k_end,kde-1) & ,NUM_TILES=grid%num_tiles & ,MAX_TIME_STEP=grid%max_time_step & ,ADAPTIVE_TS=config_flags%use_adaptive_time_step & ) CASE (MORR_TWO_MOMENT, MORR_TM_AERO) CALL diagnostic_output_nwp( & config_flags=config_flags, & U=grid%u_2 ,V=grid%v_2 & ,TEMP=grid%t_phy ,P8W=p8w & ,DT=grid%dt ,SBW=config_flags%spec_bdy_width & ,XTIME=grid%xtime & ,MPHYSICS_OPT=config_flags%mp_physics & ,GSFCGCE_HAIL=config_flags%gsfcgce_hail & ,GSFCGCE_2ICE=config_flags%gsfcgce_2ice & ,MPUSE_HAIL=config_flags%morr_rimed_ice & ,NSSL_ALPHAH=config_flags%nssl_alphah & ,NSSL_ALPHAHL=config_flags%nssl_alphahl & ,NSSL_CNOH=config_flags%nssl_cnoh & ,NSSL_CNOHL=config_flags%nssl_cnohl & ,NSSL_RHO_QH=config_flags%nssl_rho_qh & ,NSSL_RHO_QHL=config_flags%nssl_rho_qhl & ,CURR_SECS2=curr_secs2 & ,NWP_DIAGNOSTICS=config_flags%nwp_diagnostics & ,DIAGFLAG=diag_flag & ,HISTORY_INTERVAL=grid%history_interval & ,ITIMESTEP=grid%itimestep & ,U10=grid%u10,V10=grid%v10,W=grid%w_2 & ,WSPD10MAX=grid%wspd10max & ,UP_HELI_MAX=grid%up_heli_max & ,W_UP_MAX=grid%w_up_max,W_DN_MAX=grid%w_dn_max & ,ZNW=grid%znw,W_COLMEAN=grid%w_colmean & ,NUMCOLPTS=grid%numcolpts,W_MEAN=grid%w_mean & ,GRPL_MAX=grid%grpl_max,GRPL_COLINT=grid%grpl_colint & ,REFD_MAX=grid%refd_max & ,refl_10cm=grid%refl_10cm & ,HAIL_MAXK1=grid%hail_maxk1,HAIL_MAX2D=grid%hail_max2d & ,QG_CURR=moist(ims,kms,jms,P_QG) & ,NG_CURR=scalar(ims,kms,jms,P_QNG) & ,RHO=grid%rho,PH=grid%ph_2,PHB=grid%phb,G=g & ,IDS=ids,IDE=ide, JDS=jds,JDE=jde, KDS=kds,KDE=kde & ,IMS=ims,IME=ime, JMS=jms,JME=jme, KMS=kms,KME=kme & ,IPS=ips,IPE=ipe, JPS=jps,JPE=jpe, KPS=kps,KPE=kpe & ,I_START=grid%i_start,I_END=min(grid%i_end, ide-1) & ,J_START=grid%j_start,J_END=min(grid%j_end, jde-1) & ,KTS=k_start, KTE=min(k_end,kde-1) & ,NUM_TILES=grid%num_tiles & ,MAX_TIME_STEP=grid%max_time_step & ,ADAPTIVE_TS=config_flags%use_adaptive_time_step & ) CASE (MILBRANDT2MOM, NSSL_2MOM) CALL diagnostic_output_nwp( & config_flags=config_flags, & U=grid%u_2 ,V=grid%v_2 & ,TEMP=grid%t_phy ,P8W=p8w & ,DT=grid%dt ,SBW=config_flags%spec_bdy_width & ,XTIME=grid%xtime & ,MPHYSICS_OPT=config_flags%mp_physics & ,GSFCGCE_HAIL=config_flags%gsfcgce_hail & ,GSFCGCE_2ICE=config_flags%gsfcgce_2ice & ,MPUSE_HAIL=config_flags%hail_opt & ,NSSL_ALPHAH=config_flags%nssl_alphah & ,NSSL_ALPHAHL=config_flags%nssl_alphahl & ,NSSL_CNOH=config_flags%nssl_cnoh & ,NSSL_CNOHL=config_flags%nssl_cnohl & ,NSSL_RHO_QH=config_flags%nssl_rho_qh & ,NSSL_RHO_QHL=config_flags%nssl_rho_qhl & ,CURR_SECS2=curr_secs2 & ,NWP_DIAGNOSTICS=config_flags%nwp_diagnostics & ,DIAGFLAG=diag_flag & ,HISTORY_INTERVAL=grid%history_interval & ,ITIMESTEP=grid%itimestep & ,U10=grid%u10,V10=grid%v10,W=grid%w_2 & ,WSPD10MAX=grid%wspd10max & ,UP_HELI_MAX=grid%up_heli_max & ,W_UP_MAX=grid%w_up_max,W_DN_MAX=grid%w_dn_max & ,ZNW=grid%znw,W_COLMEAN=grid%w_colmean & ,NUMCOLPTS=grid%numcolpts,W_MEAN=grid%w_mean & ,GRPL_MAX=grid%grpl_max,GRPL_COLINT=grid%grpl_colint & ,REFD_MAX=grid%refd_max & ,refl_10cm=grid%refl_10cm & ,HAIL_MAXK1=grid%hail_maxk1,HAIL_MAX2D=grid%hail_max2d & ,QG_CURR=moist(ims,kms,jms,P_QG) & ,NG_CURR=scalar(ims,kms,jms,P_QNG) & ,QH_CURR=moist(ims,kms,jms,P_QH) & ,NH_CURR=scalar(ims,kms,jms,P_QNH) & ,RHO=grid%rho,PH=grid%ph_2,PHB=grid%phb,G=g & ,IDS=ids,IDE=ide, JDS=jds,JDE=jde, KDS=kds,KDE=kde & ,IMS=ims,IME=ime, JMS=jms,JME=jme, KMS=kms,KME=kme & ,IPS=ips,IPE=ipe, JPS=jps,JPE=jpe, KPS=kps,KPE=kpe & ,I_START=grid%i_start,I_END=min(grid%i_end, ide-1) & ,J_START=grid%j_start,J_END=min(grid%j_end, jde-1) & ,KTS=k_start, KTE=min(k_end,kde-1) & ,NUM_TILES=grid%num_tiles & ,MAX_TIME_STEP=grid%max_time_step & ,ADAPTIVE_TS=config_flags%use_adaptive_time_step & ) CASE DEFAULT CALL diagnostic_output_nwp( & config_flags=config_flags, & U=grid%u_2 ,V=grid%v_2 & ,TEMP=grid%t_phy ,P8W=p8w & ,DT=grid%dt ,SBW=config_flags%spec_bdy_width & ,XTIME=grid%xtime & ,MPHYSICS_OPT=config_flags%mp_physics & ,GSFCGCE_HAIL=config_flags%gsfcgce_hail & ,GSFCGCE_2ICE=config_flags%gsfcgce_2ice & ,MPUSE_HAIL=config_flags%hail_opt & ,NSSL_ALPHAH=config_flags%nssl_alphah & ,NSSL_ALPHAHL=config_flags%nssl_alphahl & ,NSSL_CNOH=config_flags%nssl_cnoh & ,NSSL_CNOHL=config_flags%nssl_cnohl & ,NSSL_RHO_QH=config_flags%nssl_rho_qh & ,NSSL_RHO_QHL=config_flags%nssl_rho_qhl & ,CURR_SECS2=curr_secs2 & ,NWP_DIAGNOSTICS=config_flags%nwp_diagnostics & ,DIAGFLAG=diag_flag & ,HISTORY_INTERVAL=grid%history_interval & ,ITIMESTEP=grid%itimestep & ,U10=grid%u10,V10=grid%v10,W=grid%w_2 & ,WSPD10MAX=grid%wspd10max & ,UP_HELI_MAX=grid%up_heli_max & ,W_UP_MAX=grid%w_up_max,W_DN_MAX=grid%w_dn_max & ,ZNW=grid%znw,W_COLMEAN=grid%w_colmean & ,NUMCOLPTS=grid%numcolpts,W_MEAN=grid%w_mean & ,GRPL_MAX=grid%grpl_max,GRPL_COLINT=grid%grpl_colint & ,REFD_MAX=grid%refd_max & ,refl_10cm=grid%refl_10cm & ,HAIL_MAXK1=grid%hail_maxk1,HAIL_MAX2D=grid%hail_max2d & ,QG_CURR=moist(ims,kms,jms,P_QG) & ,RHO=grid%rho,PH=grid%ph_2,PHB=grid%phb,G=g & ,IDS=ids,IDE=ide, JDS=jds,JDE=jde, KDS=kds,KDE=kde & ,IMS=ims,IME=ime, JMS=jms,JME=jme, KMS=kms,KME=kme & ,IPS=ips,IPE=ipe, JPS=jps,JPE=jpe, KPS=kps,KPE=kpe & ,I_START=grid%i_start,I_END=min(grid%i_end, ide-1) & ,J_START=grid%j_start,J_END=min(grid%j_end, jde-1) & ,KTS=k_start, KTE=min(k_end,kde-1) & ,NUM_TILES=grid%num_tiles & ,MAX_TIME_STEP=grid%max_time_step & ,ADAPTIVE_TS=config_flags%use_adaptive_time_step & ) END SELECT mp_select END IF NWPDIAGS CLIMATE_DIAGS : IF ( config_flags%output_diagnostics == 1 ) THEN IF ( ( config_flags%auxhist3_interval == 0 ) ) THEN WRITE (diag_message , * ) & "CLWRF: ERROR -- error -- ERROR -- error : NO 'auxhist3_interval' has been defined in 'namelist.input'" CALL wrf_error_fatal3("",749,& diag_message ) END IF CALL wrf_debug ( 100 , '--> CALL DIAGNOSTICS PACKAGE: CLIMATE DIAGNOSTICS' ) CALL clwrf_output_calc( & is_restart=config_flags%restart & ,clwrfH=config_flags%auxhist3_interval & ,T2=grid%t2, Q2=grid%q2, U10=grid%u10, V10=grid%v10 & ,SKINTEMP=grid%tsk & ,T2CLMIN=grid%t2min, T2CLMAX=grid%t2max & ,TT2CLMIN=grid%tt2min, TT2CLMAX=grid%tt2max & ,T2CLMEAN=grid%t2mean, T2CLSTD=grid%t2std & ,Q2CLMIN=grid%q2min, Q2CLMAX=grid%q2max & ,TQ2CLMIN=grid%tq2min, TQ2CLMAX=grid%tq2max & ,Q2CLMEAN=grid%q2mean, Q2CLSTD=grid%q2std & ,U10CLMAX=grid%u10max, V10CLMAX=grid%v10max & ,SPDUV10CLMAX=grid%spduv10max & ,TSPDUV10CLMAX=grid%tspduv10max & ,U10CLMEAN=grid%u10mean, V10CLMEAN=grid%v10mean & ,SPDUV10CLMEAN=grid%spduv10mean & ,U10CLSTD=grid%u10std, V10CLSTD=grid%v10std & ,SPDUV10CLSTD=grid%spduv10std & ,RAINCCLMAX=grid%raincvmax & ,RAINNCCLMAX=grid%rainncvmax & ,TRAINCCLMAX=grid%traincvmax & ,TRAINNCCLMAX=grid%trainncvmax & ,RAINCCLMEAN=grid%raincvmean & ,RAINNCCLMEAN=grid%rainncvmean & ,RAINCCLSTD=grid%raincvstd & ,RAINNCCLSTD=grid%rainncvstd & ,SKINTEMPCLMIN=grid%skintempmin & ,SKINTEMPCLMAX=grid%skintempmax & ,TSKINTEMPCLMIN=grid%tskintempmin & ,TSKINTEMPCLMAX=grid%tskintempmax & ,SKINTEMPCLMEAN=grid%skintempmean & ,SKINTEMPCLSTD=grid%skintempstd & ,RAINCV=grid%raincv ,RAINNCV=grid%rainncv & ,DT=grid%dt & ,XTIME=grid%xtime,CURR_SECS2=curr_secs2 & ,NSTEPS=grid%nsteps & ,IDS=ids,IDE=ide, JDS=jds,JDE=jde, KDS=kds,KDE=kde & ,IMS=ims,IME=ime, JMS=jms,JME=jme, KMS=kms,KME=kme & ,IPS=ips,IPE=ipe, JPS=jps,JPE=jpe, KPS=kps,KPE=kpe & ,I_START=grid%i_start,I_END=min(grid%i_end, ide-1) & ,J_START=grid%j_start,J_END=min(grid%j_end, jde-1) & ,KTS=k_start, KTE=k_end & ,NUM_TILES=grid%num_tiles & ) END IF CLIMATE_DIAGS PL_DIAGNOSTICS : IF ( config_flags%p_lev_diags .NE. SKIP_PRESS_DIAGS ) THEN TIME_TO_DO_PL_DIAGS : IF ( ( ( MOD(NINT(curr_secs2+grid%dt),NINT(config_flags%p_lev_interval)) .EQ. 0 ) ) .OR. & ( config_flags%use_adaptive_time_step ) ) THEN !$OMP PARALLEL DO & !$OMP PRIVATE ( ij ) DO ij = 1 , grid%num_tiles CALL wrf_debug ( 100 , '--> CALL DIAGNOSTICS PACKAGE: PRESSURE LEVEL DIAGNOSTICS' ) CALL pld ( & U=grid%u_2 & ,V=grid%v_2 & ,W=grid%w_2 & ,t=grid%th_phy_m_t0 & ,qv=moist(:,:,:,P_QV) & ,zp=grid%ph_2 & ,zb=grid%phb & ,pp=grid%p & ,pb=grid%pb & ,p=grid%p_hyd & ,pw=grid%p_hyd_w & ,msfux=grid%msfux & ,msfuy=grid%msfuy & ,msfvx=grid%msfvx & ,msfvy=grid%msfvy & ,msftx=grid%msftx & ,msfty=grid%msfty & ,f=grid%f & ,e=grid%e & ,use_tot_or_hyd_p=config_flags%use_tot_or_hyd_p & ,extrap_below_grnd=config_flags%extrap_below_grnd & ,missing=config_flags%p_lev_missing & ,num_press_levels=config_flags%num_press_levels & ,max_press_levels=max_plevs & ,press_levels=model_config_rec%press_levels & ,p_pl = grid%p_pl & ,u_pl = grid%u_pl & ,v_pl = grid%v_pl & ,t_pl = grid%t_pl & ,rh_pl = grid%rh_pl & ,ght_pl= grid%ght_pl & ,s_pl = grid%s_pl & ,td_pl = grid%td_pl & ,q_pl = grid%q_pl & ,IDS=ids,IDE=ide, JDS=jds,JDE=jde, KDS=kds,KDE=kde & ,IMS=ims,IME=ime, JMS=jms,JME=jme, KMS=kms,KME=kme & ,ITS=grid%i_start(ij),ITE=grid%i_end(ij) & ,JTS=grid%j_start(ij),JTE=grid%j_end(ij) & ,KTS=k_start,KTE=k_end+1 ) END DO !$OMP END PARALLEL DO END IF TIME_TO_DO_PL_DIAGS END IF PL_DIAGNOSTICS ZL_DIAGNOSTICS : IF ( config_flags%z_lev_diags .NE. SKIP_Z_DIAGS ) THEN TIME_TO_DO_ZL_DIAGS : IF ( ( ( MOD(NINT(curr_secs2+grid%dt),NINT(config_flags%z_lev_interval)) .EQ. 0 ) ) .OR. & ( config_flags%use_adaptive_time_step ) ) THEN !$OMP PARALLEL DO & !$OMP PRIVATE ( ij ) DO ij = 1 , grid%num_tiles CALL wrf_debug ( 100 , '--> CALL DIAGNOSTICS PACKAGE: HEIGHT LEVEL AND AGL DIAGNOSTICS' ) CALL zld ( & U=grid%u_2 & ,V=grid%v_2 & ,W=grid%w_2 & ,t=grid%th_phy_m_t0 & ,qv=moist(:,:,:,P_QV) & ,zp=grid%ph_2 & ,zb=grid%phb & ,pp=grid%p & ,pb=grid%pb & ,p=grid%p_hyd & ,pw=grid%p_hyd_w & ,msfux=grid%msfux & ,msfuy=grid%msfuy & ,msfvx=grid%msfvx & ,msfvy=grid%msfvy & ,msftx=grid%msftx & ,msfty=grid%msfty & ,f=grid%f & ,e=grid%e & ,ht=grid%ht & ,use_tot_or_hyd_p=config_flags%use_tot_or_hyd_p & ,extrap_below_grnd=config_flags%extrap_below_grnd & ,missing=config_flags%z_lev_missing & ,num_z_levels=config_flags%num_z_levels & ,max_z_levels=max_zlevs & ,z_levels=model_config_rec%z_levels & ,z_zl = grid%z_zl & ,u_zl = grid%u_zl & ,v_zl = grid%v_zl & ,t_zl = grid%t_zl & ,rh_zl = grid%rh_zl & ,ght_zl= grid%ght_zl & ,s_zl = grid%s_zl & ,td_zl = grid%td_zl & ,q_zl = grid%q_zl & ,p_zl = grid%p_zl & ,IDS=ids,IDE=ide, JDS=jds,JDE=jde, KDS=kds,KDE=kde & ,IMS=ims,IME=ime, JMS=jms,JME=jme, KMS=kms,KME=kme & ,ITS=grid%i_start(ij),ITE=grid%i_end(ij) & ,JTS=grid%j_start(ij),JTE=grid%j_end(ij) & ,KTS=k_start,KTE=k_end+1 ) END DO !$OMP END PARALLEL DO END IF TIME_TO_DO_ZL_DIAGS END IF ZL_DIAGNOSTICS AFWA_DIAGS : IF ( config_flags%afwa_diag_opt == 1 ) THEN IF ( ( config_flags%history_interval == 0 ) ) THEN WRITE (diag_message , * ) & "AFWA Diagnostics Error : No 'history_interval' defined in namelist" CALL wrf_error_fatal3("",956,& diag_message ) END IF !$OMP PARALLEL DO & !$OMP PRIVATE ( ij ) DO ij = 1 , grid%num_tiles CALL wrf_debug ( 100 , '--> CALL DIAGNOSTICS PACKAGE: AFWA DIAGNOSTICS' ) CALL afwa_diagnostics_driver ( grid , config_flags & ,moist & ,scalar & ,chem & ,th_phy , pi_phy , p_phy & ,grid%u_phy , grid%v_phy & ,dz8w , p8w , t8w , grid%rho & ,ids, ide, jds, jde, kds, kde & ,ims, ime, jms, jme, kms, kme & ,ips, ipe, jps, jpe, kps, kpe & ,ITS=grid%i_start(ij),ITE=grid%i_end(ij) & ,JTS=grid%j_start(ij),JTE=grid%j_end(ij) & ,K_START=k_start,K_END=k_end ) END DO !$OMP END PARALLEL DO ENDIF AFWA_DIAGS RASM_DIAGS_MEAN : IF ( config_flags%mean_diag == 1 ) THEN CALL wrf_debug ( 100 , '--> CALL DIAGNOSTICS PACKAGE: RASM DIAGNOSTICS - MEAN' ) CALL domain_clock_get ( grid, current_time=currentTime) CALL mean_output_calc( & is_restart=config_flags%restart & ,CURRENTTIME=currentTime & ,stats_interval=config_flags%mean_interval & ,output_freq=config_flags%mean_freq & ,run_days=config_flags%run_days & ,DT=grid%dt, XTIME=grid%xtime & ,PSFC=grid%psfc, PSFC_MEAN=grid%psfc_mean & ,TSK=grid%tsk, TSK_MEAN=grid%tsk_mean & ,PMSL_MEAN=grid%pmsl_mean & ,T2=grid%t2, T2_MEAN=grid%t2_mean & ,T=grid%th_phy_m_t0, P=grid%p, PB=grid%pb & ,MOIST=grid%moist(:,:,:,P_QV), HT=grid%ht & ,TH2=grid%th2, TH2_MEAN=grid%th2_mean & ,Q2=grid%q2, Q2_MEAN=grid%q2_mean & ,U10=grid%u10, U10_MEAN=grid%u10_mean & ,V10=grid%v10, V10_MEAN=grid%v10_mean & ,HFX=grid%hfx, HFX_MEAN=grid%hfx_mean & ,LH=grid%lh, LH_MEAN=grid%lh_mean & ,SWDNB=grid%swdnb, SWDNB_MEAN=grid%swdnb_mean & ,GLW=grid%glw , GLW_MEAN=grid%glw_mean & ,LWUPB=grid%lwupb, LWUPB_MEAN=grid%lwupb_mean & ,SWUPB=grid%swupb, SWUPB_MEAN=grid%swupb_mean & ,SWUPT=grid%swupt, SWUPT_MEAN=grid%swupt_mean & ,SWDNT=grid%swdnt, SWDNT_MEAN=grid%swdnt_mean & ,LWUPT=grid%lwupt, LWUPT_MEAN=grid%lwupt_mean & ,LWDNT=grid%lwdnt, LWDNT_MEAN=grid%lwdnt_mean & ,AVGOUTALARM=grid%alarms(AUXHIST5_ALARM) & ,AVGOUTDATESTR=grid%OUTDATE_MEAN & ,NSTEPS=grid%NSTEPS_MEAN & ,IDS=ids,IDE=ide, JDS=jds,JDE=jde, KDS=kds,KDE=kde & ,IMS=ims,IME=ime, JMS=jms,JME=jme, KMS=kms,KME=kme & ,IPS=ips,IPE=ipe, JPS=jps,JPE=jpe, KPS=kps,KPE=kpe & ,I_START=grid%i_start,I_END=min(grid%i_end, ide-1) & ,J_START=grid%j_start,J_END=min(grid%j_end, jde-1) & ,NUM_TILES=grid%num_tiles & ) END IF RASM_DIAGS_MEAN RASM_DIAGS_DIURNAL : IF ( config_flags%diurnal_diag == 1 ) THEN CALL wrf_debug ( 100 , '--> CALL DIAGNOSTICS PACKAGE: RASM DIAGNOSTICS - DIURNAL' ) CALL domain_clock_get ( grid, current_time=currentTime) CALL diurnalcycle_output_calc( & is_restart=config_flags%restart & ,CURRENTTIME=currentTime & ,DT=grid%dt, XTIME=grid%xtime & ,PSFC=grid%psfc, PSFC_DTMP=grid%psfc_dtmp & ,TSK=grid%tsk, TSK_DTMP=grid%tsk_dtmp & ,T2=grid%t2, T2_DTMP=grid%t2_dtmp & ,T=grid%th_phy_m_t0, P=grid%p, PB=grid%pb & ,MOIST=grid%moist(:,:,:,P_QV) & ,TH2=grid%th2, TH2_DTMP=grid%th2_dtmp & ,Q2=grid%q2, Q2_DTMP=grid%q2_dtmp & ,U10=grid%u10, U10_DTMP=grid%u10_dtmp & ,V10=grid%v10, V10_DTMP=grid%v10_dtmp & ,HFX=grid%hfx, HFX_DTMP=grid%hfx_dtmp & ,LH=grid%lh, LH_DTMP=grid%lh_dtmp & ,SWDNB=grid%swdnb, SWDNB_DTMP=grid%swdnb_dtmp & ,GLW=grid%glw, GLW_DTMP=grid%glw_dtmp & ,LWUPB=grid%lwupb, LWUPB_DTMP=grid%lwupb_dtmp & ,SWUPB=grid%swupb, SWUPB_DTMP=grid%swupb_dtmp & ,SWUPT=grid%swupt, SWUPT_DTMP=grid%swupt_dtmp & ,SWDNT=grid%swdnt, SWDNT_DTMP=grid%swdnt_dtmp & ,LWUPT=grid%lwupt, LWUPT_DTMP=grid%lwupt_dtmp & ,LWDNT=grid%lwdnt, LWDNT_DTMP=grid%lwdnt_dtmp & ,AVGOUTALARM=grid%alarms(AUXHIST6_ALARM) & ,DIURNOUTDATESTR=grid%OUTDATE_DIURN & ,AVG_NSTEPS=grid%NSTEPSMEAN_DIURN & ,DIURNAL_NSTEPS=grid%NSTEPS_DIURN & ,PSFC_DIURN=grid%PSFC_DIURN & ,TSK_DIURN=grid%TSK_DIURN, T2_DIURN=grid%T2_DIURN & ,TH2_DIURN=grid%TH2_DIURN, Q2_DIURN=grid%Q2_DIURN & ,U10_DIURN=grid%U10_DIURN, V10_DIURN=grid%V10_DIURN & ,HFX_DIURN=grid%HFX_DIURN, LH_DIURN=grid%LH_DIURN & ,SWDNB_DIURN=grid%SWDNB_DIURN, GLW_DIURN=grid%GLW_DIURN & ,LWUPB_DIURN=grid%LWUPB_DIURN, SWUPB_DIURN=grid%SWUPB_DIURN & ,SWUPT_DIURN=grid%SWUPT_DIURN, SWDNT_DIURN=grid%SWDNT_DIURN & ,LWUPT_DIURN=grid%LWUPT_DIURN, LWDNT_DIURN=grid%LWDNT_DIURN & ,IDS=ids, IDE=ide, JDS=jds, JDE=jde, KDS=kds, KDE=kde & ,IMS=ims, IME=ime, JMS=jms, JME=jme, KMS=kms, KME=kme & ,IPS=ips, IPE=ipe, JPS=jps, JPE=jpe, KPS=kps, KPE=kpe & ,I_START=grid%i_start, I_END=min(grid%i_end, ide-1) & ,J_START=grid%j_start, J_END=min(grid%j_end, jde-1) & ,NUM_TILES=grid%num_tiles & ) END IF RASM_DIAGS_DIURNAL NOAHMP_DIAGS : IF ( config_flags%sf_surface_physics == 4 .and. config_flags%noahmp_output == 3 ) THEN CALL wrf_debug ( 100 , '--> CALL DIAGNOSTICS PACKAGE: NOAHMP_DIAGS' ) CALL domain_clock_get ( grid, current_time=currentTime) CALL noahmp_output_calc( & ims=ims, ime=ime, jms=jms, jme=jme ,& i_start=grid%i_start, i_end=min(grid%i_end, ide-1) ,& j_start=grid%j_start, j_end=min(grid%j_end, jde-1) ,& num_tiles=grid%num_tiles ,& dt=grid%dt ,& noahmp_acc_dt=config_flags%noahmp_acc_dt ,& opt_run=config_flags%opt_run ,& stepwtd=grid%stepwtd ,& itimestep=grid%itimestep ,& curr_secs2=curr_secs2 ,& acsagb=grid%acsagb , acsagv=grid%acsagv , sag=grid%sagxy ,& acirg=grid%acirg , irg=grid%irgxy ,& acshg=grid%acshg , shg=grid%shgxy ,& acevg=grid%acevg , evg=grid%evgxy ,& acghv=grid%acghv , ghv=grid%ghvxy ,& acpahg=grid%acpahg , pahg=grid%pahgxy ,& acirb=grid%acirb , irb=grid%irbxy ,& acshb=grid%acshb , shb=grid%shbxy ,& acevb=grid%acevb , evb=grid%evbxy ,& acghb=grid%acghb , ghb=grid%ghbxy ,& acpahb=grid%acpahb , pahb=grid%pahbxy ,& acsav=grid%acsav , sav=grid%savxy ,& acirc=grid%acirc , ircmp=grid%ircxy ,& acshc=grid%acshc , shc=grid%shcxy ,& acevc=grid%acevc , evc=grid%evcxy ,& actr=grid%actr , tr=grid%trxy ,& acpahv=grid%acpahv , pahv=grid%pahvxy ,& acswdnlsm=grid%acswdnlsm , swdown=grid%swnorm ,& acswuplsm=grid%acswuplsm , albedo=grid%albedo ,& aclwdnlsm=grid%aclwdnlsm , emiss=grid%emiss , glw=grid%glw ,& aclwuplsm=grid%aclwuplsm , trad=grid%tradxy ,& acshflsm=grid%acshflsm ,aclhflsm=grid%aclhflsm ,& acghflsm=grid%acghflsm ,acpahlsm=grid%acpahlsm ,& acints=grid%acints , qints=grid%qintsxy ,& acintr=grid%acintr , qintr=grid%qintrxy ,& acdrips=grid%acdrips , qdrips=grid%qdripsxy ,& acdripr=grid%acdripr , qdripr=grid%qdriprxy ,& acthros=grid%acthros , qthros=grid%qthrosxy ,& acthror=grid%acthror , qthror=grid%qthrorxy ,& acsnsub=grid%acsnsub , qsnsub=grid%qsnsubxy ,& acsnfro=grid%acsnfro , qsnfro=grid%qsnfroxy ,& acsubc=grid%acsubc , qsubc=grid%qsubcxy ,& acfroc=grid%acfroc , qfroc=grid%qfrocxy ,& acevac=grid%acevac , qevac=grid%qevacxy ,& acdewc=grid%acdewc , qdewc=grid%qdewcxy ,& acfrzc=grid%acfrzc , qfrzc=grid%qfrzcxy ,& acmeltc=grid%acmeltc , qmeltc=grid%qmeltcxy ,& acsnbot=grid%acsnbot , qsnbot=grid%qsnbotxy ,& acsnmelt=grid%acsnmelt , qmelt=grid%qmeltxy ,& acponding=grid%acponding , ponding=grid%pondingxy ,& acrainlsm=grid%acrainlsm , rainlsm=grid%rainlsm ,& acsnowlsm=grid%acsnowlsm , snowlsm=grid%snowlsm ,& acrainsnow=grid%acrainsnow, isnow=grid%isnowxy ,& acrunsb=grid%acrunsb , runsb=grid%runsbxy ,& acrunsf=grid%acrunsf , runsf=grid%runsfxy ,& acecan=grid%acecan , ecan=grid%ecanxy ,& acetran=grid%acetran , etran=grid%etranxy ,& acedir=grid%acedir , edir=grid%edirxy ,& acetlsm=grid%acetlsm ,& acqlat=grid%acqlat , qlat=grid%qlatxy ,& acqrf=grid%acqrf , qrf=grid%qrfxy ,& acrech=grid%acrech ,deeprech=grid%deeprechxy ,& acqspring=grid%acqspring , qspring=grid%qspringxy ,& fpice=grid%fpicexy , fveg=grid%fvegxy ,& aceflxb=grid%aceflxb , eflxb=grid%eflxbxy ,& accanhs=grid%accanhs , canhs=grid%canhsxy ) END IF NOAHMP_DIAGS SOLAR_FIELDS: IF (config_flags%solar_diagnostics == do_solar_output) THEN CALL wrf_debug (100 , '--> CALL DIAGNOSTICS PACKAGE: SOLAR_DIAG') CALL solar_diag (dt=grid%dt, & rho=grid%rho, dz8w=dz8w, ph=grid%ph_2, phb=grid%phb, & cldfrac3d=grid%cldfra, coszen=grid%coszen, swdnb=grid%swdnb, & swdnt=grid%swdnt, param_first_scalar=param_first_scalar, & p_qc=p_qc, p_qi=p_qi, p_qs=p_qs, qv=moist(ims,kms,jms,P_QV), & qc=moist(ims,kms,jms,P_QC), qi=moist(ims,kms,jms,P_QI), & qs=moist(ims,kms,jms,P_QS), qc_tot=grid%qc_tot, & qi_tot=grid%qi_tot, & has_reqc=grid%has_reqc, has_reqi=grid%has_reqi, & has_reqs=grid%has_reqs, f_qv=f_qv, f_qc=f_qc, f_qi=f_qi, & f_qs=f_qs, re_cloud=grid%re_cloud, re_ice=grid%re_ice, & re_snow=grid%re_snow, cldfrac2d=grid%cldfrac2d, & wvp2d=grid%wvp, lwp2d=grid%lwp, iwp2d=grid%iwp, & swp2d=grid%swp, wp2d_sum=grid%wp_sum, lwp2d_tot=grid%lwp_tot, & iwp2d_tot=grid%iwp_tot, wp2d_tot_sum=grid%wp_tot_sum, & re_cloud_path=grid%re_qc, re_ice_path=grid%re_qi, & re_snow_path=grid%re_qs, re_cloud_path_tot=grid%re_qc_tot, & re_ice_path_tot=grid%re_qi_tot, tau_qc=grid%tau_qc, & tau_qi=grid%tau_qi, tau_qs=grid%tau_qs, & tau_qc_tot=grid%tau_qc_tot, tau_qi_tot=grid%tau_qi_tot, & cbase=grid%cbaseht, ctop=grid%ctopht, & cbase_tot=grid%cbaseht_tot, ctop_tot=grid%ctopht_tot, & clrnidx=grid%clrnidx, sza=grid%sza, ghi_accum=grid%ghi_accum, & ids=ids, ide=ide, jds=jds, jde=jde, kds=kds, kde=kde, & ims=ims, ime=ime, jms=jms, jme=jme, kms=kms, kme=kme, & ips=ips, ipe=ipe, jps=jps, jpe=jpe, kps=kps, kpe=kpe, & i_start=grid%i_start, i_end=min(grid%i_end, ide-1), & j_start=grid%j_start, j_end=min(grid%j_end, jde-1), & kts=k_start, kte=min(k_end,kde-1), num_tiles=grid%num_tiles ) END IF SOLAR_FIELDS END SUBROUTINE diagnostics_driver SUBROUTINE update_phys_fields ( grid, config_flags, moist, & ids, ide, jds, jde, kds, kde, & ims, ime, jms, jme, kms, kme, & ips, ipe, jps, jpe, kps, kpe ) USE module_domain, ONLY : domain USE module_configure, ONLY : grid_config_rec_type USE module_state_description, ONLY: num_moist, P_Qv USE module_model_constants, ONLY: g, R_v, R_d, Cp, T0, RCP TYPE ( domain ), INTENT(INOUT) :: grid TYPE ( grid_config_rec_type ), INTENT(IN) :: config_flags REAL , DIMENSION(ims:ime,kms:kme,jms:jme,num_moist ) , INTENT(IN) :: moist INTEGER , INTENT(IN) :: ids, ide, jds, jde, kds, kde, & ims, ime, jms, jme, kms, kme, & ips, ipe, jps, jpe, kps, kpe INTEGER :: i, j, k IF ( config_flags%use_theta_m .EQ. 1 ) THEN DO j = jps, MIN(jde-1, jpe) DO k = kps, kde-1 DO i = ips, MIN(ide-1, ipe) grid%th_phy_m_t0(i,k,j) = (grid%t_2(i,k,j)+T0) / & (1.+(R_v/R_d)*moist(i,k,j,P_Qv)) - T0 END DO END DO END DO ELSE DO j = jps, MIN(jde-1, jpe) DO k = kps, kde-1 DO i = ips, MIN(ide-1, ipe) grid%th_phy_m_t0(i,k,j) = grid%t_2(i,k,j) END DO END DO END DO END IF END SUBROUTINE update_phys_fields END MODULE module_diagnostics_driver