Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_adding_ica_sw.F90
- Timestamp:
- Aug 2, 2024, 2:12:03 PM (7 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad/radiation/radiation_adding_ica_sw.F90
r4773 r5158 84 84 ! half-level by working down through the atmosphere 85 85 flux_dn_direct(:,1) = incoming_toa 86 dojlev = 1,nlev86 DO jlev = 1,nlev 87 87 flux_dn_direct(:,jlev+1) = flux_dn_direct(:,jlev)*trans_dir_dir(:,jlev) 88 88 end do … … 100 100 ! Added for DWD (2020) 101 101 !NEC$ outerloop_unroll(8) 102 dojlev = nlev,1,-1102 DO jlev = nlev,1,-1 103 103 ! Next loop over columns. We could do this by indexing the 104 104 ! entire inner dimension as follows, e.g. for the first line: … … 107 107 ! routine by a factor of 2! Rather, we do it with an explicit 108 108 ! loop. 109 dojcol = 1,ncol109 DO jcol = 1,ncol 110 110 ! Lacis and Hansen (1974) Eq 33, Shonk & Hogan (2008) Eq 10: 111 111 inv_denominator(jcol,jlev) = 1.0_jprb / (1.0_jprb-albedo(jcol,jlev+1)*reflectance(jcol,jlev)) … … 132 132 ! Added for DWD (2020) 133 133 !NEC$ outerloop_unroll(8) 134 dojlev = 1,nlev135 dojcol = 1,ncol134 DO jlev = 1,nlev 135 DO jcol = 1,ncol 136 136 ! Shonk & Hogan (2008) Eq 14 (after simplification): 137 137 flux_dn_diffuse(jcol,jlev+1) &
Note: See TracChangeset
for help on using the changeset viewer.