29 |
#include "FFIELDS.h" |
#include "FFIELDS.h" |
30 |
#include "THSICE_PARAMS.h" |
#include "THSICE_PARAMS.h" |
31 |
#include "THSICE_VARS.h" |
#include "THSICE_VARS.h" |
|
#ifdef ALLOW_BULK_FORCE |
|
|
#include "BULKF.h" |
|
|
#endif |
|
32 |
|
|
33 |
C !INPUT/OUTPUT PARAMETERS: |
C !INPUT/OUTPUT PARAMETERS: |
34 |
C === Routine arguments === |
C === Routine arguments === |
127 |
flxSW (i,j) = 0. |
flxSW (i,j) = 0. |
128 |
snowPrc(i,j,bi,bj) = 0. _d 0 |
snowPrc(i,j,bi,bj) = 0. _d 0 |
129 |
siceAlb(i,j,bi,bj) = 0. _d 0 |
siceAlb(i,j,bi,bj) = 0. _d 0 |
|
#ifdef ALLOW_BULK_FORCE |
|
|
prcAtm(i,j) = ( rain(i,j,bi,bj)+runoff(i,j,bi,bj) )*rhofw |
|
|
flxSW (i,j) = solar(i,j,bi,bj) |
|
|
IF ( iceMask(i,j,bi,bj).GT.0. _d 0 |
|
|
& .AND. Tair(i,j,bi,bj).LE.Tf0kel ) THEN |
|
|
snowPrc(i,j,bi,bj) = rain(i,j,bi,bj)*rhofw |
|
|
ENDIF |
|
|
#endif |
|
130 |
ENDDO |
ENDDO |
131 |
ENDDO |
ENDDO |
132 |
|
#ifdef ALLOW_BULK_FORCE |
133 |
|
IF ( useBulkforce ) THEN |
134 |
|
CALL THSICE_GET_PRECIP( |
135 |
|
I iceMask, |
136 |
|
O prcAtm, snowPrc(1-Olx,1-Oly,bi,bj), flxSW, |
137 |
|
I iMin,iMax,jMin,jMax, bi,bj, myThid ) |
138 |
|
ENDIF |
139 |
|
#endif |
140 |
|
|
141 |
CALL THSICE_STEP_FWD( |
CALL THSICE_STEP_FWD( |
142 |
I bi, bj, iMin, iMax, jMin, jMax, |
I bi, bj, iMin, iMax, jMin, jMax, |