75 |
#ifdef SHORTWAVE_HEATING |
#ifdef SHORTWAVE_HEATING |
76 |
& -Qsw(i,j,bi,bj)+ |
& -Qsw(i,j,bi,bj)+ |
77 |
#endif |
#endif |
78 |
& (surfaceTendencyT(i,j,bi,bj)+surfaceTendencyTice(I,J,bi,bj)) |
& (surfaceForcingT(i,j,bi,bj)+surfaceForcingTice(I,J,bi,bj)) |
79 |
& *HeatCapacity_Cp*recip_horiVertRatio*rhoConst |
& *HeatCapacity_Cp*recip_horiVertRatio*rhoConst |
|
& *drF(k)*hFacC(i,j,k,bi,bj) |
|
80 |
ENDDO |
ENDDO |
81 |
ENDDO |
ENDDO |
82 |
#ifdef NONLIN_FRSURF |
#ifdef NONLIN_FRSURF |
98 |
DO j=1,sNy |
DO j=1,sNy |
99 |
DO i=1,sNx |
DO i=1,sNx |
100 |
tmpFld(i,j)= |
tmpFld(i,j)= |
101 |
& surfaceTendencyS(i,j,bi,bj)* |
& surfaceForcingS(i,j,bi,bj)* |
102 |
& recip_horiVertRatio*rhoConst*drF(k)*hFacC(i,j,k,bi,bj) |
& recip_horiVertRatio*rhoConst |
103 |
ENDDO |
ENDDO |
104 |
ENDDO |
ENDDO |
105 |
#ifdef NONLIN_FRSURF |
#ifdef NONLIN_FRSURF |