641 |
ENDDO |
ENDDO |
642 |
ENDIF |
ENDIF |
643 |
|
|
644 |
|
#ifdef ALLOW_GEOTHERMAL_FLUX |
645 |
|
IF ( usingZCoords ) THEN |
646 |
|
DO j=1,sNy |
647 |
|
DO i=1,sNx |
648 |
|
IF ( k.EQ.kLowC(i,j,bi,bj) ) THEN |
649 |
|
gT_arr(i,j)=gT_arr(i,j) |
650 |
|
& + geothermalFlux(i,j,bi,bj) |
651 |
|
& *recip_Cp*mass2rUnit |
652 |
|
& *recip_drF(k)*_recip_hFacC(i,j,k,bi,bj) |
653 |
|
ENDIF |
654 |
|
ENDDO |
655 |
|
ENDDO |
656 |
|
ENDIF |
657 |
|
#endif /* ALLOW_GEOTHERMAL_FLUX */ |
658 |
|
|
659 |
#ifdef SHORTWAVE_HEATING |
#ifdef SHORTWAVE_HEATING |
660 |
C Penetrating SW radiation |
C Penetrating SW radiation |
661 |
c IF ( usePenetratingSW ) THEN |
c IF ( usePenetratingSW ) THEN |