/[MITgcm]/MITgcm/pkg/exf/exf_bulk_largeyeager04.F
ViewVC logotype

Diff of /MITgcm/pkg/exf/exf_bulk_largeyeager04.F

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

revision 1.2 by jmc, Mon May 7 19:35:20 2007 UTC revision 1.3 by jmc, Mon May 14 14:03:50 2007 UTC
# Line 178  C    take U_N = del.U ; stability from d Line 178  C    take U_N = del.U ; stability from d
178  C     in this case ustress and vstress are defined a u and v points  C     in this case ustress and vstress are defined a u and v points
179  C     respectively, and we need to average to mass points to avoid  C     respectively, and we need to average to mass points to avoid
180  C     windStress = 0 near coasts or other boundaries  C     windStress = 0 near coasts or other boundaries
181               windStress = sqrt(0.5 _d 0*              windStress = sqrt( 0.5 _d 0*
182       &                   (ustress(i,  j,bi,bj)*ustress(i  ,j,bi,bj)       &                   (ustress(i,  j,bi,bj)*ustress(i  ,j,bi,bj)
183       &                   +ustress(i+1,j,bi,bj)*ustress(i+1,j,bi,bj)       &                   +ustress(i+1,j,bi,bj)*ustress(i+1,j,bi,bj)
184       &                   +vstress(i,j,  bi,bj)*vstress(i,j  ,bi,bj)       &                   +vstress(i,j,  bi,bj)*vstress(i,j  ,bi,bj)
185       &                   +vstress(i,j+1,bi,bj)*vstress(i,j+1,bi,bj)       &                   +vstress(i,j+1,bi,bj)*vstress(i,j+1,bi,bj)
186       &                   )    )       &                   )   )
187               ustar = sqrt(windStess/atmrho)              ustar = sqrt(windStress/atmrho)
188    c           tau   = windStress/ustar
189                tau   = sqrt(windStress*atmrho)
190  #endif /* ALLOW_ATM_WIND */  #endif /* ALLOW_ATM_WIND */
191    
192  C--  initial guess for exchange other coefficients:  C--  initial guess for exchange other coefficients:
# Line 247  C-   Shift wind speed using old coeffici Line 249  C-   Shift wind speed using old coeffici
249                usn = ws/(exf_one + rdn/karman*(zwln-psimh) )                usn = ws/(exf_one + rdn/karman*(zwln-psimh) )
250                usm = MAX(usn, umin)                usm = MAX(usn, umin)
251    
252  C-   Update the 10m, neutral stability transfer coefficients  C-   Update the 10m, neutral stability transfer coefficients (momentum)
253                tmpbulk = cdrag_1/usm + cdrag_2 + cdrag_3*usm                tmpbulk = cdrag_1/usm + cdrag_2 + cdrag_3*usm
254                rdn = SQRT(tmpbulk)                rdn = SQRT(tmpbulk)
255                rd = rdn/(exf_one + rdn*(zwln-psimh)/karman)                rd = rdn/(exf_one + rdn*(zwln-psimh)/karman)
256                ustar = rd*wsm                ustar = rd*wsm
257    C-   Coeff:
258                  tau   = atmrho*rd*ws
259  #endif  #endif
260    
261  C-   Update the 10m, neutral stability transfer coefficients  C-   Update the 10m, neutral stability transfer coefficients (sens&evap)
262                rhn = (exf_one-stable)*cstanton_1 + stable*cstanton_2                rhn = (exf_one-stable)*cstanton_1 + stable*cstanton_2
263                ren = cDalton                ren = cDalton
264    
# Line 278  CADJ STORE cw(i,j,bi ,bj)     = comlev1_ Line 282  CADJ STORE cw(i,j,bi ,bj)     = comlev1_
282  CADJ STORE sw(i,j,bi,bj)     = comlev1_exf_1, key = ikey_1  CADJ STORE sw(i,j,bi,bj)     = comlev1_exf_1, key = ikey_1
283  #endif  #endif
284    
 C-   Coeff:  
             tau           = atmrho*rd*ws  
   
285  C-   Turbulent Fluxes  C-   Turbulent Fluxes
286              hs(i,j,bi,bj) = atmcp*tau*tstar              hs(i,j,bi,bj) = atmcp*tau*tstar
287              hl(i,j,bi,bj) = flamb*tau*qstar              hl(i,j,bi,bj) = flamb*tau*qstar
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22