Home
last modified time | relevance | path

Searched refs:kproma (Results 1 – 24 of 24) sorted by relevance

/dports/math/libpgmath/flang-d07daf3/test/ncar_kernels/PSRAD_lrtm_Bangalore/src/
H A Dmo_rrtm_coeffs.f90115 forfac(1:kproma,1:klev) = forfac(1:kproma,1:klev) * colh2o(1:kproma,1:klev)
116 selffac(1:kproma,1:klev) = selffac(1:kproma,1:klev) * colh2o(1:kproma,1:klev)
126 rat_h2oo3(1:kproma, jk) = chi_mls(1,jp(1:kproma, jk))/chi_mls(3,jp(1:kproma, jk))
144 scaleminor(1:kproma,1:klev) = play(1:kproma,1:klev)/tlay(1:kproma,1:klev)
230 fp(1:kproma,jk) = 5._wp *(preflog(jp(1:kproma,jk)) - plog(1:kproma,jk))
253 water(1:kproma,1:klev) = wkl(1:kproma,1,1:klev)/coldry(1:kproma,1:klev)
264 fac10(1:kproma,1:klev) = compfp(1:kproma,1:klev) * ft(1:kproma,1:klev)
265 fac00(1:kproma,1:klev) = compfp(1:kproma,1:klev) * (1._wp - ft(1:kproma,1:klev))
266 fac11(1:kproma,1:klev) = fp(1:kproma,1:klev) * ft1(1:kproma,1:klev)
267 fac01(1:kproma,1:klev) = fp(1:kproma,1:klev) * (1._wp - ft1(1:kproma,1:klev))
[all …]
H A Dmo_lrtm_solver.f9089 odepth(1:kproma,jk) = max(0._wp, secdiff(1:kproma) * tau(1:kproma,jk))
100 tfn(1:kproma) = tautrans(odepth(:,jk), kproma)
101 dplnkup(1:kproma,jk) = levplnk(1:kproma,jk) - layplnk(1:kproma,jk)
102 dplnkdn(1:kproma,jk) = levplnk(1:kproma,jk-1) - layplnk(1:kproma,jk)
103 …bbup(1:kproma,jk) = weights(1:kproma,jk) * (layplnk(1:kproma,jk) + dplnkup(1:kproma,jk) * tfn(1:kp…
104 …bbdn(1:kproma,jk) = weights(1:kproma,jk) * (layplnk(1:kproma,jk) + dplnkdn(1:kproma,jk) * tfn(1:kp…
113 trans(1:kproma,jk) = transmit(odepth(:,jk), kproma)
115 …rad_dn(1:kproma,jk-1) = rad_dn(1:kproma,jk) + (bbdn(1:kproma,jk) - rad_dn(1:kproma,jk)) * trans(1:…
120 …rad_up(1:kproma, 0) = weights(1:kproma, 1) * surfemis(1:kproma) * surfplanck(1:kproma) …
121 surfemis(1:kproma)) * rad_dn(1:kproma, 0)
[all …]
H A Dmo_lrtm_driver.f90249 DO jl = 1, kproma
303 wbrodl(1:kproma,1:klev) = coldry(1:kproma,1:klev) - sum(wkl(1:kproma,2:,1:klev), dim=2)
386 …tautot(1:kproma,:,:) = taut(1:kproma,:,:) + smp_tau(1:kproma,:,:) ! All-sky optical depth. Mask fo…
402 planklay(1:kproma,1:klev,ig) = planckfunction(tlay(1:kproma,1:klev ),ig)
403 planklev(1:kproma,0:klev,ig) = planckfunction(tlev(1:kproma,1:klev+1),ig)
409 …pwvcm(1:kproma) = ((amw * sum(wkl(1:kproma,1,1:klev), dim=2)) / (amd * sum(…
415 secdiff(1:kproma) = find_secdiff(ibs(1:kproma, ig), pwvcm(1:kproma))
429 …uflx(1:kproma,0:klev) = uflx (1:kproma,0:klev) + zgpfu(1:kproma,0:kle…
430 …dflx(1:kproma,0:klev) = dflx (1:kproma,0:klev) + zgpfd(1:kproma,0:kle…
440 … uflxc(1:kproma,0:klev) = uflxc(1:kproma,0:klev) + zgpcu(1:kproma,0:klev) * gpt_scaling
[all …]
H A Dmo_cld_sampling.f9027 SUBROUTINE sample_cld_state(kproma, kbdim, klev, ksamps, rnseeds, i_overlap, cld_frac, cldy) argument
31 INTEGER, intent(in) :: kproma !< numbers of columns, levels, samples local
41 cldy(1:kproma,1:klev,1) = cld_frac(1:kproma,1:klev) > 0._wp
48 CALL get_random(kproma, kbdim, rnseeds, cldy(:,jk,1), rank(:,jk,js))
54 …rank(1:kproma,jk,js) = merge(rank(1:kproma,jk+1,js), …
55 … rank(1:kproma,jk,js) * (1._wp - cld_frac(1:kproma,jk+1)), &
56 rank(1:kproma,jk+1,js) > 1._wp - cld_frac(1:kproma,jk+1))
67 CALL get_random(kproma, kbdim, rnseeds, rank(:, 1, js))
68 rank(1:kproma,2:klev,js) = spread(rank(1:kproma,1,js), dim=2, ncopies=(klev-1))
77 CALL get_random(kproma, kbdim, rnseeds, cldy(:,jk,1), rank(:,jk,js))
[all …]
H A Dmo_spec_sampling.f90106 FUNCTION get_gpoint_set(kproma, kbdim, strategy, seeds)
107 INTEGER, intent(in) :: kproma local
111 INTEGER, dimension(kproma, strategy%num_gpts_ts) :: get_gpoint_set
125 get_gpoint_set(1:kproma,:) = spread(strategy%teams(:, 1), dim = 1, ncopies = kproma)
131 CALL get_random(kproma, kbdim, seeds, rn)
132 team(1:kproma) = min(int(rn(1:kproma) * num_teams) + 1, num_teams)
133 DO jl = 1, kproma
142 CALL get_random(kproma, kbdim, seeds, rn)
143 team(1:kproma) = min(int(rn(1:kproma) * num_teams) + 1, num_teams)
144 get_gpoint_set(1:kproma, it) = strategy%teams(1, team(1:kproma))
H A Dmo_random_numbers.f9062 SUBROUTINE kissvec_all(kproma, kbdim, seed, harvest) argument
64 INTEGER, intent(in ) :: kproma local
69 CALL kissvec(kproma, kbdim, seed, mask, harvest)
73 SUBROUTINE kissvec(kproma, kbdim, seed, mask, harvest) argument
75 INTEGER, intent(in ) :: kproma local
79 INTEGER(KIND=i8) :: kiss(kproma)
81 DO jk = 1, kproma
H A Dkernel_driver.f9050 INTEGER :: kproma variable
92 READ(UNIT=kgen_unit) kproma
96 call psrad_interface(kbdim, klev, nb_sw, kproma, ktrac, tk_sfc, kgen_unit)
H A Dmo_psrad_interface.f90109 SUBROUTINE psrad_interface(kbdim, klev, nb_sw, kproma, ktrac, tk_sfc, kgen_unit) argument
140 INTEGER, intent(in) :: kproma local
327 …CALL lrtm(kproma, kbdim, klev, pm_fl_vr, pm_sfc, tk_fl_vr, tk_hl_vr, tk_sfc, wkl_vr, wx_vr, col_dr…
337 …CALL lrtm(kproma, kbdim, klev, pm_fl_vr, pm_sfc, tk_fl_vr, tk_hl_vr, tk_sfc, wkl_vr, wx_vr, col_dr…
/dports/math/libpgmath/flang-d07daf3/test/ncar_kernels/PSRAD_lrtm_codereview/src/
H A Dmo_rrtm_coeffs.f90115 forfac(1:kproma,1:klev) = forfac(1:kproma,1:klev) * colh2o(1:kproma,1:klev)
116 selffac(1:kproma,1:klev) = selffac(1:kproma,1:klev) * colh2o(1:kproma,1:klev)
126 rat_h2oo3(1:kproma, jk) = chi_mls(1,jp(1:kproma, jk))/chi_mls(3,jp(1:kproma, jk))
144 scaleminor(1:kproma,1:klev) = play(1:kproma,1:klev)/tlay(1:kproma,1:klev)
230 fp(1:kproma,jk) = 5._wp *(preflog(jp(1:kproma,jk)) - plog(1:kproma,jk))
253 water(1:kproma,1:klev) = wkl(1:kproma,1,1:klev)/coldry(1:kproma,1:klev)
264 fac10(1:kproma,1:klev) = compfp(1:kproma,1:klev) * ft(1:kproma,1:klev)
265 fac00(1:kproma,1:klev) = compfp(1:kproma,1:klev) * (1._wp - ft(1:kproma,1:klev))
266 fac11(1:kproma,1:klev) = fp(1:kproma,1:klev) * ft1(1:kproma,1:klev)
267 fac01(1:kproma,1:klev) = fp(1:kproma,1:klev) * (1._wp - ft1(1:kproma,1:klev))
[all …]
H A Dmo_lrtm_solver.f9089 odepth(1:kproma,jk) = max(0._wp, secdiff(1:kproma) * tau(1:kproma,jk))
100 tfn(1:kproma) = tautrans(odepth(:,jk), kproma)
101 dplnkup(1:kproma,jk) = levplnk(1:kproma,jk) - layplnk(1:kproma,jk)
102 dplnkdn(1:kproma,jk) = levplnk(1:kproma,jk-1) - layplnk(1:kproma,jk)
103 …bbup(1:kproma,jk) = weights(1:kproma,jk) * (layplnk(1:kproma,jk) + dplnkup(1:kproma,jk) * tfn(1:kp…
104 …bbdn(1:kproma,jk) = weights(1:kproma,jk) * (layplnk(1:kproma,jk) + dplnkdn(1:kproma,jk) * tfn(1:kp…
113 trans(1:kproma,jk) = transmit(odepth(:,jk), kproma)
115 …rad_dn(1:kproma,jk-1) = rad_dn(1:kproma,jk) + (bbdn(1:kproma,jk) - rad_dn(1:kproma,jk)) * trans(1:…
120 …rad_up(1:kproma, 0) = weights(1:kproma, 1) * surfemis(1:kproma) * surfplanck(1:kproma) …
121 surfemis(1:kproma)) * rad_dn(1:kproma, 0)
[all …]
H A Dmo_cld_sampling.f9027 SUBROUTINE sample_cld_state(kproma, kbdim, klev, ksamps, rnseeds, i_overlap, cld_frac, cldy) argument
31 INTEGER, intent(in) :: kproma !< numbers of columns, levels, samples local
41 cldy(1:kproma,1:klev,1) = cld_frac(1:kproma,1:klev) > 0._wp
48 CALL get_random(kproma, kbdim, rnseeds, cldy(:,jk,1), rank(:,jk,js))
54 …rank(1:kproma,jk,js) = merge(rank(1:kproma,jk+1,js), …
55 … rank(1:kproma,jk,js) * (1._wp - cld_frac(1:kproma,jk+1)), &
56 rank(1:kproma,jk+1,js) > 1._wp - cld_frac(1:kproma,jk+1))
67 CALL get_random(kproma, kbdim, rnseeds, rank(:, 1, js))
68 rank(1:kproma,2:klev,js) = spread(rank(1:kproma,1,js), dim=2, ncopies=(klev-1))
77 CALL get_random(kproma, kbdim, rnseeds, cldy(:,jk,1), rank(:,jk,js))
[all …]
H A Dmo_lrtm_driver.f90251 DO jl = 1, kproma
305 wbrodl(1:kproma,1:klev) = coldry(1:kproma,1:klev) - sum(wkl(1:kproma,2:,1:klev), dim=2)
406 …tautot(1:kproma,:,:) = taut(1:kproma,:,:) + smp_tau(1:kproma,:,:) ! All-sky optical depth. Mask fo…
422 planklay(1:kproma,1:klev,ig) = planckfunction(tlay(1:kproma,1:klev ),ig)
423 planklev(1:kproma,0:klev,ig) = planckfunction(tlev(1:kproma,1:klev+1),ig)
429 …pwvcm(1:kproma) = ((amw * sum(wkl(1:kproma,1,1:klev), dim=2)) / (amd * sum(…
435 secdiff(1:kproma) = find_secdiff(ibs(1:kproma, ig), pwvcm(1:kproma))
449 …uflx(1:kproma,0:klev) = uflx (1:kproma,0:klev) + zgpfu(1:kproma,0:kle…
450 …dflx(1:kproma,0:klev) = dflx (1:kproma,0:klev) + zgpfd(1:kproma,0:kle…
460 … uflxc(1:kproma,0:klev) = uflxc(1:kproma,0:klev) + zgpcu(1:kproma,0:klev) * gpt_scaling
[all …]
H A Dmo_spec_sampling.f90106 FUNCTION get_gpoint_set(kproma, kbdim, strategy, seeds)
107 INTEGER, intent(in) :: kproma local
111 INTEGER, dimension(kproma, strategy%num_gpts_ts) :: get_gpoint_set
125 get_gpoint_set(1:kproma,:) = spread(strategy%teams(:, 1), dim = 1, ncopies = kproma)
131 CALL get_random(kproma, kbdim, seeds, rn)
132 team(1:kproma) = min(int(rn(1:kproma) * num_teams) + 1, num_teams)
133 DO jl = 1, kproma
142 CALL get_random(kproma, kbdim, seeds, rn)
143 team(1:kproma) = min(int(rn(1:kproma) * num_teams) + 1, num_teams)
144 get_gpoint_set(1:kproma, it) = strategy%teams(1, team(1:kproma))
H A Dmo_random_numbers.f9062 SUBROUTINE kissvec_all(kproma, kbdim, seed, harvest) argument
64 INTEGER, intent(in ) :: kproma local
69 CALL kissvec(kproma, kbdim, seed, mask, harvest)
73 SUBROUTINE kissvec(kproma, kbdim, seed, mask, harvest) argument
75 INTEGER, intent(in ) :: kproma local
79 INTEGER(KIND=i8) :: kiss(kproma)
81 DO jk = 1, kproma
H A Dkernel_driver.f9050 INTEGER :: kproma variable
92 READ(UNIT=kgen_unit) kproma
96 call psrad_interface(kbdim, klev, nb_sw, kproma, ktrac, tk_sfc, kgen_unit)
H A Dmo_psrad_interface.f90109 SUBROUTINE psrad_interface(kbdim, klev, nb_sw, kproma, ktrac, tk_sfc, kgen_unit) argument
140 INTEGER, intent(in) :: kproma local
327 …CALL lrtm(kproma, kbdim, klev, pm_fl_vr, pm_sfc, tk_fl_vr, tk_hl_vr, tk_sfc, wkl_vr, wx_vr, col_dr…
337 …CALL lrtm(kproma, kbdim, klev, pm_fl_vr, pm_sfc, tk_fl_vr, tk_hl_vr, tk_sfc, wkl_vr, wx_vr, col_dr…
/dports/math/libpgmath/flang-d07daf3/test/ncar_kernels/PSRAD_lrtm/src/
H A Dmo_rrtm_coeffs.f90115 forfac(1:kproma,1:klev) = forfac(1:kproma,1:klev) * colh2o(1:kproma,1:klev)
116 selffac(1:kproma,1:klev) = selffac(1:kproma,1:klev) * colh2o(1:kproma,1:klev)
126 rat_h2oo3(1:kproma, jk) = chi_mls(1,jp(1:kproma, jk))/chi_mls(3,jp(1:kproma, jk))
144 scaleminor(1:kproma,1:klev) = play(1:kproma,1:klev)/tlay(1:kproma,1:klev)
230 fp(1:kproma,jk) = 5._wp *(preflog(jp(1:kproma,jk)) - plog(1:kproma,jk))
253 water(1:kproma,1:klev) = wkl(1:kproma,1,1:klev)/coldry(1:kproma,1:klev)
264 fac10(1:kproma,1:klev) = compfp(1:kproma,1:klev) * ft(1:kproma,1:klev)
265 fac00(1:kproma,1:klev) = compfp(1:kproma,1:klev) * (1._wp - ft(1:kproma,1:klev))
266 fac11(1:kproma,1:klev) = fp(1:kproma,1:klev) * ft1(1:kproma,1:klev)
267 fac01(1:kproma,1:klev) = fp(1:kproma,1:klev) * (1._wp - ft1(1:kproma,1:klev))
[all …]
H A Dmo_lrtm_solver.f9089 odepth(1:kproma,jk) = max(0._wp, secdiff(1:kproma) * tau(1:kproma,jk))
100 tfn(1:kproma) = tautrans(odepth(:,jk), kproma)
101 dplnkup(1:kproma,jk) = levplnk(1:kproma,jk) - layplnk(1:kproma,jk)
102 dplnkdn(1:kproma,jk) = levplnk(1:kproma,jk-1) - layplnk(1:kproma,jk)
103 …bbup(1:kproma,jk) = weights(1:kproma,jk) * (layplnk(1:kproma,jk) + dplnkup(1:kproma,jk) * tfn(1:kp…
104 …bbdn(1:kproma,jk) = weights(1:kproma,jk) * (layplnk(1:kproma,jk) + dplnkdn(1:kproma,jk) * tfn(1:kp…
113 trans(1:kproma,jk) = transmit(odepth(:,jk), kproma)
115 …rad_dn(1:kproma,jk-1) = rad_dn(1:kproma,jk) + (bbdn(1:kproma,jk) - rad_dn(1:kproma,jk)) * trans(1:…
120 …rad_up(1:kproma, 0) = weights(1:kproma, 1) * surfemis(1:kproma) * surfplanck(1:kproma) …
121 surfemis(1:kproma)) * rad_dn(1:kproma, 0)
[all …]
H A Dmo_lrtm_driver.f90228 DO jl = 1, kproma
282 wbrodl(1:kproma,1:klev) = coldry(1:kproma,1:klev) - sum(wkl(1:kproma,2:,1:klev), dim=2)
314 …tautot(1:kproma,:,:) = taut(1:kproma,:,:) + smp_tau(1:kproma,:,:) ! All-sky optical depth. Mask fo…
330 planklay(1:kproma,1:klev,ig) = planckfunction(tlay(1:kproma,1:klev ),ig)
331 planklev(1:kproma,0:klev,ig) = planckfunction(tlev(1:kproma,1:klev+1),ig)
337 …pwvcm(1:kproma) = ((amw * sum(wkl(1:kproma,1,1:klev), dim=2)) / (amd * sum(…
343 secdiff(1:kproma) = find_secdiff(ibs(1:kproma, ig), pwvcm(1:kproma))
357 …uflx(1:kproma,0:klev) = uflx (1:kproma,0:klev) + zgpfu(1:kproma,0:kle…
358 …dflx(1:kproma,0:klev) = dflx (1:kproma,0:klev) + zgpfd(1:kproma,0:kle…
368 … uflxc(1:kproma,0:klev) = uflxc(1:kproma,0:klev) + zgpcu(1:kproma,0:klev) * gpt_scaling
[all …]
H A Dmo_cld_sampling.f9027 SUBROUTINE sample_cld_state(kproma, kbdim, klev, ksamps, rnseeds, i_overlap, cld_frac, cldy) argument
31 INTEGER, intent(in) :: kproma !< numbers of columns, levels, samples local
41 cldy(1:kproma,1:klev,1) = cld_frac(1:kproma,1:klev) > 0._wp
48 CALL get_random(kproma, kbdim, rnseeds, cldy(:,jk,1), rank(:,jk,js))
54 …rank(1:kproma,jk,js) = merge(rank(1:kproma,jk+1,js), …
55 … rank(1:kproma,jk,js) * (1._wp - cld_frac(1:kproma,jk+1)), &
56 rank(1:kproma,jk+1,js) > 1._wp - cld_frac(1:kproma,jk+1))
67 CALL get_random(kproma, kbdim, rnseeds, rank(:, 1, js))
68 rank(1:kproma,2:klev,js) = spread(rank(1:kproma,1,js), dim=2, ncopies=(klev-1))
77 CALL get_random(kproma, kbdim, rnseeds, cldy(:,jk,1), rank(:,jk,js))
[all …]
H A Dmo_spec_sampling.f90106 FUNCTION get_gpoint_set(kproma, kbdim, strategy, seeds)
107 INTEGER, intent(in) :: kproma local
111 INTEGER, dimension(kproma, strategy%num_gpts_ts) :: get_gpoint_set
125 get_gpoint_set(1:kproma,:) = spread(strategy%teams(:, 1), dim = 1, ncopies = kproma)
131 CALL get_random(kproma, kbdim, seeds, rn)
132 team(1:kproma) = min(int(rn(1:kproma) * num_teams) + 1, num_teams)
133 DO jl = 1, kproma
142 CALL get_random(kproma, kbdim, seeds, rn)
143 team(1:kproma) = min(int(rn(1:kproma) * num_teams) + 1, num_teams)
144 get_gpoint_set(1:kproma, it) = strategy%teams(1, team(1:kproma))
H A Dmo_random_numbers.f9062 SUBROUTINE kissvec_all(kproma, kbdim, seed, harvest) argument
64 INTEGER, intent(in ) :: kproma local
69 CALL kissvec(kproma, kbdim, seed, mask, harvest)
73 SUBROUTINE kissvec(kproma, kbdim, seed, mask, harvest) argument
75 INTEGER, intent(in ) :: kproma local
79 INTEGER(KIND=i8) :: kiss(kproma)
81 DO jk = 1, kproma
H A Dkernel_driver.f9050 INTEGER :: kproma variable
92 READ(UNIT=kgen_unit) kproma
96 call psrad_interface(kbdim, klev, nb_sw, kproma, ktrac, tk_sfc, kgen_unit)
H A Dmo_psrad_interface.f90109 SUBROUTINE psrad_interface(kbdim, klev, nb_sw, kproma, ktrac, tk_sfc, kgen_unit) argument
140 INTEGER, intent(in) :: kproma local
327 …CALL lrtm(kproma, kbdim, klev, pm_fl_vr, pm_sfc, tk_fl_vr, tk_hl_vr, tk_sfc, wkl_vr, wx_vr, col_dr…
337 …CALL lrtm(kproma, kbdim, klev, pm_fl_vr, pm_sfc, tk_fl_vr, tk_hl_vr, tk_sfc, wkl_vr, wx_vr, col_dr…