/[MITgcm]/MITgcm/pkg/aim_v23/phy_suflux_sice.F
ViewVC logotype

Diff of /MITgcm/pkg/aim_v23/phy_suflux_sice.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by jmc, Thu Mar 11 14:33:19 2004 UTC revision 1.2 by jmc, Thu Apr 8 00:14:09 2004 UTC
# Line 64  C    Evp0   :: evaporation computed over Line 64  C    Evp0   :: evaporation computed over
64  C    dEvp   :: evaporation derivative relative to surf. temp  C    dEvp   :: evaporation derivative relative to surf. temp
65  C    Slr0   :: upward long wave radiation over freezing surf.  C    Slr0   :: upward long wave radiation over freezing surf.
66  C    dSlr   :: upward long wave rad. derivative relative to surf. temp  C    dSlr   :: upward long wave rad. derivative relative to surf. temp
67  C    sFlx   :: net surface flux (+=down) function of surf. temp Ts:  C    sFlx   :: net heat flux (+=down) except SW, function of surf. temp Ts:
68  C              0: Flux(Ts=0.oC) ; 1: Flux(Ts^n) ; 2: d.Flux/d.Ts(Ts^n)  C              0: Flux(Ts=0.oC) ; 1: Flux(Ts^n) ; 2: d.Flux/d.Ts(Ts^n)
69  C    TSFC   :: surface temperature (clim.)     (2-dim)  C    TSFC   :: surface temperature (clim.)     (2-dim)
70  C    TSKIN  :: skin surface temperature        (2-dim)  C    TSKIN  :: skin surface temperature        (2-dim)
# Line 143  C     2.4 Emission of lw radiation from Line 143  C     2.4 Emission of lw radiation from
143  C--   Compute net surface heat flux and its derivative ./. surf. temp.  C--   Compute net surface heat flux and its derivative ./. surf. temp.
144        DO J=1,NGP        DO J=1,NGP
145          sFlx(J,0)= sFlx(J,0)          sFlx(J,0)= sFlx(J,0)
146       &           - ALHC*Evp0(J) - EMISloc*Slr0(J) + SLRD(J) + SSR(J)       &           - ALHC*Evp0(J) - EMISloc*Slr0(J) + SLRD(J)
147          sFlx(J,1)= sFlx(J,1)          sFlx(J,1)= sFlx(J,1)
148       &           - ALHC*EVAP(J) - EMISloc*SLRU(J) + SLRD(J) + SSR(J)       &           - ALHC*EVAP(J) - EMISloc*SLRU(J) + SLRD(J)
149          sFlx(J,2)= sFlx(J,2)          sFlx(J,2)= sFlx(J,2)
150       &           - ALHC*dEvp(J) - EMISloc*dSlr(J)       &           - ALHC*dEvp(J) - EMISloc*dSlr(J)
151        ENDDO        ENDDO
152        IF ( aim_energPrecip ) THEN        IF ( aim_energPrecip ) THEN
153  C-     Evap of snow: substract Latent Heat of freezing from heatFlux  C-     Evap of snow/ice: substract Latent Heat of freezing from heatFlux
154         DO J=1,NGP         DO J=1,NGP
         IF ( EnPrec(J) .LT. 0. ) THEN  
155           sFlx(J,0) = sFlx(J,0) - ALHF*Evp0(J)           sFlx(J,0) = sFlx(J,0) - ALHF*Evp0(J)
156           sFlx(J,1) = sFlx(J,1) - ALHF*EVAP(J)           sFlx(J,1) = sFlx(J,1) - ALHF*EVAP(J)
157           sFlx(J,2) = sFlx(J,2) - ALHF*dEvp(J)           sFlx(J,2) = sFlx(J,2) - ALHF*dEvp(J)
         ENDIF  
158         ENDDO         ENDDO
159        ENDIF        ENDIF
160    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22