55 |
INTEGER myThid |
INTEGER myThid |
56 |
CEOP |
CEOP |
57 |
|
|
|
#ifdef ALLOW_THSICE |
|
58 |
#ifdef ALLOW_EXF |
#ifdef ALLOW_EXF |
59 |
|
|
60 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
64 |
|
|
65 |
DO j = jMin, jMax |
DO j = jMin, jMax |
66 |
DO i = iMin, iMax |
DO i = iMin, iMax |
67 |
|
#ifdef ALLOW_ATM_TEMP |
68 |
totPrc(i,j) = precip(i,j,bi,bj)*rhoConstFresh |
totPrc(i,j) = precip(i,j,bi,bj)*rhoConstFresh |
69 |
|
#endif |
70 |
#if defined (ALLOW_RUNOFF) || defined (ALLOW_SEAICE) |
#if defined (ALLOW_RUNOFF) || defined (ALLOW_SEAICE) |
71 |
totPrc(i,j) = totPrc(i,j) + runoff(i,j,bi,bj)*rhoConstFresh |
totPrc(i,j) = totPrc(i,j) + runoff(i,j,bi,bj)*rhoConstFresh |
72 |
#endif |
#endif |
73 |
CML flxSW (i,j) = -Qsw(i,j,bi,bj) |
CML flxSW (i,j) = -Qsw(i,j,bi,bj) |
74 |
|
#ifdef ALLOW_DOWNWARD_RADIATION |
75 |
flxSW (i,j) = swdown(i,j,bi,bj) |
flxSW (i,j) = swdown(i,j,bi,bj) |
76 |
|
#else |
77 |
|
STOP 'ABNORMAL END: S/R THSICE_MAP_EXF: DOWNWARD_RADIATION undef' |
78 |
|
#endif |
79 |
ENDDO |
ENDDO |
80 |
ENDDO |
ENDDO |
81 |
|
|
82 |
|
#ifdef ALLOW_ATM_TEMP |
83 |
IF ( snowPrecipFile .NE. ' ' ) THEN |
IF ( snowPrecipFile .NE. ' ' ) THEN |
84 |
DO j = jMin, jMax |
DO j = jMin, jMax |
85 |
DO i = iMin, iMax |
DO i = iMin, iMax |
99 |
ENDDO |
ENDDO |
100 |
ENDDO |
ENDDO |
101 |
ENDIF |
ENDIF |
102 |
|
#else /* ALLOW_ATM_TEMP */ |
103 |
|
STOP 'ABNORMAL END: S/R THSICE_MAP_EXF: ATM_TEMP undef' |
104 |
|
#endif /* ALLOW_ATM_TEMP */ |
105 |
|
|
106 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
107 |
|
|
108 |
#endif /* ALLOW_EXF */ |
#endif /* ALLOW_EXF */ |
|
#endif /* ALLOW_THSICE */ |
|
109 |
|
|
110 |
RETURN |
RETURN |
111 |
END |
END |