@@ -3556,17 +3556,27 @@ end subroutine fv_getDivergence
35563556subroutine fv_getUpdraftHelicity (uh25 )
35573557 use constants_mod, only: fms_grav= >grav
35583558 use fv_diagnostics_mod, only: get_vorticity, updraft_helicity
3559- real (FVPRC), intent (OUT ) :: uh25(FV_Atm(1 )% bd% isc:FV_Atm(1 )% bd% iec,FV_Atm(1 )% bd% jsc:FV_Atm(1 )% bd% jec)
3559+ ! made this REAL4
3560+ real (REAL4), intent (OUT ) :: uh25(FV_Atm(1 )% bd% isc:FV_Atm(1 )% bd% iec,FV_Atm(1 )% bd% jsc:FV_Atm(1 )% bd% jec)
3561+
3562+ ! made an intermediate output of FVPRC
3563+ real (FVPRC) :: uh25_tmp(FV_Atm(1 )% bd% isc:FV_Atm(1 )% bd% iec,FV_Atm(1 )% bd% jsc:FV_Atm(1 )% bd% jec)
35603564 integer :: sphum= 1
35613565 real (FVPRC) :: vort(FV_Atm(1 )% bd% isc:FV_Atm(1 )% bd% iec,FV_Atm(1 )% bd% jsc:FV_Atm(1 )% bd% jec,FV_Atm(1 )% npz)
3566+ ! introduced these two variables for the literals
3567+ real (FVPRC) :: z_bot, z_top
3568+
3569+ z_bot = 2.e3
3570+ z_top = 5.e3
35623571 call get_vorticity(FV_Atm(1 )% bd% isc, FV_Atm(1 )% bd% iec, FV_Atm(1 )% bd% jsc, FV_Atm(1 )% bd% jec, &
35633572 FV_Atm(1 )% bd% isd, FV_Atm(1 )% bd% ied, FV_Atm(1 )% bd% jsd, FV_Atm(1 )% bd% jed, &
35643573 FV_Atm(1 )% npz, FV_Atm(1 )% u, FV_Atm(1 )% v, vort, &
35653574 FV_Atm(1 )% gridstruct% dx, FV_Atm(1 )% gridstruct% dy, FV_Atm(1 )% gridstruct% rarea)
35663575 call updraft_helicity(FV_Atm(1 )% bd% isc, FV_Atm(1 )% bd% iec, FV_Atm(1 )% bd% jsc, FV_Atm(1 )% bd% jec, FV_Atm(1 )% ng, FV_Atm(1 )% npz, &
3567- zvir, sphum, uh25 , &
3576+ zvir, sphum, uh25_tmp , &
35683577 FV_Atm(1 )% w, vort, FV_Atm(1 )% delz, FV_Atm(1 )% q, &
3569- FV_Atm(1 )% flagstruct% hydrostatic, FV_Atm(1 )% pt, FV_Atm(1 )% peln, FV_Atm(1 )% phis, fms_grav, 2.e3 , 5.e3 )
3578+ FV_Atm(1 )% flagstruct% hydrostatic, FV_Atm(1 )% pt, FV_Atm(1 )% peln, FV_Atm(1 )% phis, fms_grav, z_bot, z_top)
3579+ uh25 = uh25_tmp
35703580end subroutine fv_getUpdraftHelicity
35713581
35723582subroutine fv_getEPV (pt , vort , ua , va , epv )
0 commit comments