343 |
C-- Return atmospheric fluxes in evpAtm & flxSW (same sign and units): |
C-- Return atmospheric fluxes in evpAtm & flxSW (same sign and units): |
344 |
evpAtm(i,j) = frwAtm |
evpAtm(i,j) = frwAtm |
345 |
flxSW (i,j) = flxAtm(i,j) |
flxSW (i,j) = flxAtm(i,j) |
346 |
|
|
347 |
|
#ifdef ATMOSPHERIC_LOADING |
348 |
|
C-- Compute Sea-Ice Loading (= mass of sea-ice + snow / area unit) |
349 |
|
sIceLoad(i,j,bi,bj) = ( snowHeight(i,j,bi,bj)*rhos |
350 |
|
& + iceHeight(i,j,bi,bj)*rhoi |
351 |
|
& )*iceMask(i,j,bi,bj) |
352 |
|
#endif |
353 |
|
|
354 |
ENDDO |
ENDDO |
355 |
ENDDO |
ENDDO |
356 |
|
|