16 |
C !DESCRIPTION: \bv |
C !DESCRIPTION: \bv |
17 |
C *==========================================================* |
C *==========================================================* |
18 |
C | S/R THSICE_MAP_EXF |
C | S/R THSICE_MAP_EXF |
19 |
C | Interface S/R : map Precip, Snow and shortwave fluxes |
C | Interface S/R : map Precip, Snow and shortwave fluxes |
20 |
C | from pkg EXF to thsice variables |
C | from pkg EXF to thsice variables |
21 |
C *==========================================================* |
C *==========================================================* |
22 |
C \ev |
C \ev |
40 |
C iceMsk :: sea-ice fraction: no ice=0, grid all ice 1 [] |
C iceMsk :: sea-ice fraction: no ice=0, grid all ice 1 [] |
41 |
C totPrc :: Total Precipitation (including run-off) [kg/m2/s] |
C totPrc :: Total Precipitation (including run-off) [kg/m2/s] |
42 |
C snowPrc :: Snow Precipitation [kg/m2/s] |
C snowPrc :: Snow Precipitation [kg/m2/s] |
43 |
C flxSW :: Net short-wave surface flux (+=down) [W/m2] |
C flxSW :: Downward short-wave surface flux (+=down) [W/m2] |
44 |
C iMin,iMax :: range of indices of computation domain |
C iMin,iMax :: range of indices of computation domain |
45 |
C jMin,jMax :: range of indices of computation domain |
C jMin,jMax :: range of indices of computation domain |
46 |
C bi,bj :: current tile indices |
C bi,bj :: current tile indices |
72 |
#else |
#else |
73 |
STOP 'ABNORMAL END: S/R THSICE_MAP_EXF: ALLOW_RUNOFF undef' |
STOP 'ABNORMAL END: S/R THSICE_MAP_EXF: ALLOW_RUNOFF undef' |
74 |
#endif |
#endif |
|
CML flxSW (i,j) = -Qsw(i,j,bi,bj) |
|
75 |
#ifdef ALLOW_DOWNWARD_RADIATION |
#ifdef ALLOW_DOWNWARD_RADIATION |
76 |
flxSW (i,j) = swdown(i,j,bi,bj) |
flxSW (i,j) = swdown(i,j,bi,bj) |
77 |
#else |
#else |
78 |
STOP 'ABNORMAL END: S/R THSICE_MAP_EXF: DOWNWARD_RADIATION undef' |
STOP 'ABNORMAL END: S/R THSICE_MAP_EXF: DOWNWARD_RADIATION undef' |
79 |
#endif |
#endif |
80 |
ENDDO |
ENDDO |
81 |
ENDDO |
ENDDO |
82 |
|
|
83 |
#ifdef ALLOW_ATM_TEMP |
#ifdef ALLOW_ATM_TEMP |
84 |
IF ( snowPrecipFile .NE. ' ' ) THEN |
IF ( snowPrecipFile .NE. ' ' ) THEN |
85 |
DO j = jMin, jMax |
DO j = jMin, jMax |
88 |
ENDDO |
ENDDO |
89 |
ENDDO |
ENDDO |
90 |
ELSE |
ELSE |
91 |
C If specific snow precipitiation is now available, use |
C If specific snow precipitiation is now available, use |
92 |
C precipitation when ever the air temperature is below 0 degC |
C precipitation when ever the air temperature is below 0 degC |
93 |
DO j = jMin, jMax |
DO j = jMin, jMax |
94 |
DO i = iMin, iMax |
DO i = iMin, iMax |