91 |
C |
C |
92 |
tsurfLoc (I,J) = MIN(273.16 _d 0+MAX_TICE,TSURF(I,J,bi,bj)) |
tsurfLoc (I,J) = MIN(273.16 _d 0+MAX_TICE,TSURF(I,J,bi,bj)) |
93 |
C Is this necessary? |
C Is this necessary? |
94 |
|
#ifdef ALLOW_ATM_TEMP |
95 |
atempLoc (I,J) = MAX(273.16 _d 0+MIN_ATEMP,ATEMP(I,J,bi,bj)) |
atempLoc (I,J) = MAX(273.16 _d 0+MIN_ATEMP,ATEMP(I,J,bi,bj)) |
96 |
|
#else /* ALLOW_ATM_TEMP */ |
97 |
|
STOP 'ABNORMAL END: S/R SEAICE_BUDGET_OCE: ATM_TEMP undef' |
98 |
|
#endif /* ALLOW_ATM_TEMP */ |
99 |
|
#ifdef ALLOW_DOWNWARD_RADIATION |
100 |
lwdownLoc(I,J) = MAX(MIN_LWDOWN,LWDOWN(I,J,bi,bj)) |
lwdownLoc(I,J) = MAX(MIN_LWDOWN,LWDOWN(I,J,bi,bj)) |
101 |
|
#else |
102 |
|
STOP |
103 |
|
& 'ABNORMAL END: S/R SEAICE_BUDGET_OCE: DOWNWARD_RADIATION undef' |
104 |
|
#endif |
105 |
ENDDO |
ENDDO |
106 |
ENDDO |
ENDDO |
107 |
#endif /* SEAICE_EXTERNAL_FLUXES */ |
#endif /* SEAICE_EXTERNAL_FLUXES */ |
115 |
SWHeatFlux (I,J) = Qsw(I,J,bi,bj) |
SWHeatFlux (I,J) = Qsw(I,J,bi,bj) |
116 |
#else /* SEAICE_EXTERNAL_FLUXES undefined */ |
#else /* SEAICE_EXTERNAL_FLUXES undefined */ |
117 |
ALB(I,J)=SEAICE_waterAlbedo |
ALB(I,J)=SEAICE_waterAlbedo |
118 |
|
#if defined(ALLOW_DOWNWARD_RADIATION) && defined(ALLOW_ATM_TEMP) |
119 |
A1(I,J)=(ONE-ALB(I,J))*SWDOWN(I,J,bi,bj) |
A1(I,J)=(ONE-ALB(I,J))*SWDOWN(I,J,bi,bj) |
120 |
& +lwdownLoc(I,J)*0.97 _d 0 |
& +lwdownLoc(I,J)*0.97 _d 0 |
121 |
& +D1*UG(I,J)*atempLoc(I,J)+D1W*UG(I,J)*AQH(I,J,bi,bj) |
& +D1*UG(I,J)*atempLoc(I,J)+D1W*UG(I,J)*AQH(I,J,bi,bj) |
126 |
& -D3*(tsurfLoc(I,J)**4) |
& -D3*(tsurfLoc(I,J)**4) |
127 |
netHeatFlux(I,J)=-A1(I,J)-A2(I,J) |
netHeatFlux(I,J)=-A1(I,J)-A2(I,J) |
128 |
SWHeatFlux (I,J)=-(ONE-ALB(I,J))*SWDOWN(I,J,bi,bj) |
SWHeatFlux (I,J)=-(ONE-ALB(I,J))*SWDOWN(I,J,bi,bj) |
129 |
|
#else |
130 |
|
STOP |
131 |
|
& 'ABNORMAL END: S/R SEAICE_BUDGET_OCE: DOWNWARD_RADIATION undef' |
132 |
|
#endif |
133 |
#endif /* SEAICE_EXTERNAL_FLUXES */ |
#endif /* SEAICE_EXTERNAL_FLUXES */ |
134 |
ENDDO |
ENDDO |
135 |
ENDDO |
ENDDO |