Reorder SNICAR SSP test data arrays to be consistent with source code

eclare108213 · Aug 26, 2022 · e569c4c · e569c4c
1 parent 4e2d90b
commit e569c4c
Showing 1 changed file with 9 additions and 6 deletions.
diff --git a/columnphysics/icepack_shortwave.F90 b/columnphysics/icepack_shortwave.F90
@@ -194,53 +194,56 @@ subroutine icepack_init_radiation
             allocate(ssp_sasymmdr(nspint_5bd,nmbrad_snicar)) ! snow asymmetry factor, direct
             allocate(ssp_sasymmdf(nspint_5bd,nmbrad_snicar)) ! snow asymmetry factor, diffuse
 
+            ! tcraig, these data statements are not consistent with the array index order
+            ! reshape from nmbrad,nspint to nspint,nmbrad with order = 2,1
+
             ssp_snwextdr = reshape((/ &
                46.27374983, 24.70286257, 6.54918455, 2.6035624,  1.196168,   &
                46.56827715, 24.81790668, 6.56181227, 2.60604155, 1.19682614, &
                46.76114033, 24.90468677, 6.56950323, 2.60780948, 1.19737512, &
                46.9753992,  24.9495155,  6.57687695, 2.60937383, 1.19774008, &
                47.48598349, 25.14100194, 6.59708024, 2.61372576, 1.19897351 /), &
-               (/nspint_5bd,nmbrad_snicar/))
+               (/nmbrad_snicar,nspint_5bd/), order=(/ 2, 1 /))
 
             ssp_snwextdf = reshape((/ &
                46.26936158, 24.70165487, 6.54903637, 2.60353051, 1.19615825, &
                46.56628244, 24.81707286, 6.56164279, 2.60601584, 1.1968169,  &
                46.75501968, 24.90175807, 6.5693214,  2.60775795, 1.19736237, &
                46.95368476, 24.94497414, 6.57612007, 2.60924059, 1.19770981, &
                47.29620774, 25.0713585,  6.5891698,  2.61198929, 1.19850197 /), &
-               (/nspint_5bd,nmbrad_snicar/))
+               (/nmbrad_snicar,nspint_5bd/), order=(/ 2, 1 /))
 
             ssp_snwalbdr = reshape((/ &
                0.99999581, 0.99999231, 0.99997215, 0.99993093, 0.99985157, &
                0.99987892, 0.99978212, 0.99923115, 0.99817946, 0.99628374, &
                0.99909682, 0.99835523, 0.99418008, 0.98592085, 0.97062632, &
                0.99308867, 0.98972448, 0.97261528, 0.93947995, 0.88207117, &
                0.93322249, 0.89645776, 0.75765643, 0.6272883,  0.55435033 /), &
-               (/nspint_5bd,nmbrad_snicar/))
+               (/nmbrad_snicar,nspint_5bd/), order=(/ 2, 1 /))
 
             ssp_snwalbdf = reshape((/ &
                0.99999586, 0.9999924,  0.99997249, 0.99993178, 0.9998534,  &
                0.99988441, 0.99979202, 0.99926745, 0.99826865, 0.99647135, &
                0.99910997, 0.99837683, 0.99426171, 0.98610431, 0.97097643, &
                0.99476731, 0.9916114,  0.97441209, 0.94127464, 0.88440735, &
                0.9419837,  0.90613207, 0.77042376, 0.63887161, 0.5566098  /), &
-               (/nspint_5bd,nmbrad_snicar/))
+               (/nmbrad_snicar,nspint_5bd/), order=(/ 2, 1 /))
 
             ssp_sasymmdr = reshape((/ &
                0.88503036, 0.88778473, 0.89049023, 0.89112501, 0.89136157, &
                0.88495225, 0.88905367, 0.89315385, 0.89433673, 0.8950027,  &
                0.88440201, 0.88928256, 0.89503166, 0.89762648, 0.90045378, &
                0.88590371, 0.89350221, 0.90525156, 0.91314567, 0.92157748, &
                0.9033537,  0.91778261, 0.94615574, 0.96323447, 0.97167644 /), &
-               (/nspint_5bd,nmbrad_snicar/))
+               (/nmbrad_snicar,nspint_5bd/), order=(/ 2, 1 /))
 
             ssp_sasymmdf = reshape((/ &
                0.88500571, 0.88773868, 0.89042496, 0.8910553,  0.89128949, &
                0.88495225, 0.88905367, 0.89315385, 0.89433673, 0.8950027 , &
                0.88441433, 0.88929133, 0.89500611, 0.89756091, 0.90035504, &
                0.88495554, 0.89201556, 0.90204619, 0.90914885, 0.91769988, &
                0.89620237, 0.90998944, 0.94126152, 0.96209938, 0.9726631  /), &
-               (/nspint_5bd,nmbrad_snicar/))
+               (/nmbrad_snicar,nspint_5bd/), order=(/ 2, 1 /))
 
          else
             call icepack_warnings_setabort(.true.,__FILE__,__LINE__)