| 33 |
DO j=1, sNy |
DO j=1, sNy |
| 34 |
DO i=1, sNx |
DO i=1, sNx |
| 35 |
|
|
|
CJRS ask J-M about next lines |
|
| 36 |
icFlxAtm(i,j,1,1) = -qneti_2D(i,j) |
icFlxAtm(i,j,1,1) = -qneti_2D(i,j) |
| 37 |
icFlxSw(i,j,1,1) = - solarinc_2D(i,j) !at this point, incident SW |
icFlxSw(i,j,1,1) = - solarinc_2D(i,j) !at this point, incident SW |
| 38 |
dFdT(i,j,1,1) = - dFdT_ice_2D(i,j) |
dFdT(i,j,1,1) = - dFdT_ice_2D(i,j) |
| 39 |
|
|
| 40 |
C Note converted to ice model sign convention below |
C Note converion to ice model sign convention below |
| 41 |
sum_prcIce(i,j) = sum_prcIce(i,j) - precipi_2D(i,j)*rhofw |
sum_prcIce(i,j) = sum_prcIce(i,j) - precipi_2D(i,j)*rhofw |
| 42 |
sum_evapIce(i,j) = sum_evapIce(i,j) - evapi_2D(i,j)*rhofw |
sum_evapIce(i,j) = sum_evapIce(i,j) - evapi_2D(i,j)*rhofw |
| 43 |
IF (Tair_2D(i,j) .LE. 0. _d 0) |
IF (Tair_2D(i,j) .LE. 0. _d 0) |
| 44 |
& sum_snowPrc(i,j) = sum_snowPrc(i,j) - precipi_2D(i,j)*rhofw |
& sum_snowPrc(i,j) = sum_snowPrc(i,j) - precipi_2D(i,j)*rhofw |
| 45 |
|
|
| 46 |
C IF (iceMask(i,j,1,1) .GT. 0. _d 0) THEN |
C IF (iceMask(i,j,1,1) .GT. 0. _d 0) THEN |
| 47 |
C PRINT *,'Pass seaice fluxes at, T :',i,j,Tair_2D(i,j) |
C PRINT *,'Pass seaice fluxes at, T :',i,j,Tair_2D(i,j) |
| 48 |
C PRINT *,precipi_2D(i,j),sum_snowPrc(i,j) |
C PRINT *,precipi_2D(i,j),sum_snowPrc(i,j) |