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

Diff of /MITgcm/pkg/aim_v23/phy_suflux_ocean.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 Jun 24 23:43:11 2004 UTC
# Line 9  C     !INTERFACE: Line 9  C     !INTERFACE:
9        SUBROUTINE SUFLUX_OCEAN(        SUBROUTINE SUFLUX_OCEAN(
10       I                   PSA, FMASK,       I                   PSA, FMASK,
11       I                   Tsurf, SSR, SLRD,       I                   Tsurf, SSR, SLRD,
12       I                   T0, Q0, CDENVV,       I                   T1, T0, Q0, DENVV,
13       O                   SHF, EVAP, SLRU,       O                   SHF, EVAP, SLRU,
14       I                   bi,bj,myThid)       I                   bi,bj,myThid)
15    
# Line 43  C    FMASK  :: fractional land-sea mask Line 43  C    FMASK  :: fractional land-sea mask
43  C    Tsurf  :: surface temperature        (2-dim)  C    Tsurf  :: surface temperature        (2-dim)
44  C    SSR    :: sfc sw radiation (net flux)     (2-dim)  C    SSR    :: sfc sw radiation (net flux)     (2-dim)
45  C    SLRD   :: sfc lw radiation (downward flux)(2-dim)  C    SLRD   :: sfc lw radiation (downward flux)(2-dim)
46    C    T1     :: near-surface air temperature (from Pot.temp)
47  C    T0     :: near-surface air temperature    (2-dim)  C    T0     :: near-surface air temperature    (2-dim)
48  C    Q0     :: near-surface sp. humidity [g/kg](2-dim)  C    Q0     :: near-surface sp. humidity [g/kg](2-dim)
49  C    CDENVV :: sensible heat flux coefficient  (2-dim)  C    DENVV  :: surface flux (sens,lat.) coeff. (=Rho*|V|) [kg/m2/s]
50  C--   Output:  C--   Output:
51  C    SHF    :: sensible heat flux              (2-dim)  C    SHF    :: sensible heat flux              (2-dim)
52  C    EVAP   :: evaporation [g/(m^2 s)]         (2-dim)  C    EVAP   :: evaporation [g/(m^2 s)]         (2-dim)
# Line 56  C    myThid :: Thread number for this in Line 57  C    myThid :: Thread number for this in
57  C--  C--
58        _RL  PSA(NGP)        _RL  PSA(NGP)
59        _RL  FMASK(NGP), Tsurf(NGP)        _RL  FMASK(NGP), Tsurf(NGP)
60        _RL  SSR(NGP), SLRD(NGP), T0(NGP), Q0(NGP), CDENVV(NGP)        _RL  SSR(NGP), SLRD(NGP)
61          _RL  T1(NGP), T0(NGP), Q0(NGP), DENVV(NGP)
62    
63        _RL  SHF(NGP), EVAP(NGP), SLRU(NGP)        _RL  SHF(NGP), EVAP(NGP), SLRU(NGP)
64    
# Line 66  CEOP Line 68  CEOP
68  #ifdef ALLOW_AIM  #ifdef ALLOW_AIM
69    
70  C-- Local variables:  C-- Local variables:
71    C    CDENVV :: surf. heat flux (sens.,lat.) coeff including stability effect
72          _RL CDENVV(NGP), RDTH, FSSEA
73        _RL QSAT0(NGP)        _RL QSAT0(NGP)
74        _RL QDUMMY(1), RDUMMY(1), TS4        _RL QDUMMY(1), RDUMMY(1), TS4
75        INTEGER J        INTEGER J
# Line 78  C     2.1 Wind stress Line 82  C     2.1 Wind stress
82    
83  C     2.2 Sensible heat flux (from clim. TS over land)  C     2.2 Sensible heat flux (from clim. TS over land)
84    
85    C     Stability correction
86    
87          RDTH = FSTAB/DTHETA
88    
89          DO J=1,NGP
90            FSSEA =1.+MIN(DTHETA,MAX(-DTHETA, Tsurf(J)-T1(J)))*RDTH
91            CDENVV(J)=CHS*DENVV(J)*FSSEA
92          ENDDO
93    
94        DO J=1,NGP        DO J=1,NGP
95          SHF(J) = CDENVV(J)*CP*(Tsurf(J)-T0(J))          SHF(J) = CDENVV(J)*CP*(Tsurf(J)-T0(J))
96        ENDDO        ENDDO
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22