Changeset 1012


Ignore:
Timestamp:
Jul 24, 2013, 6:04:38 PM (11 years ago)
Author:
aslmd
Message:

LMDZ.COMMON. No need for ok_strato to activate sponge layer at model top. Default iflag_top_bound is 0 and no longer 1. Modified doc accordingly.

Location:
trunk
Files:
8 edited

Legend:

Unmodified
Added
Removed

  • trunk/DOC/documentation/top_bound.tex

    r1010 r1012  
    6565\section{Pratical aspects in the code}
    6666
    67 The sponge layer is applied at the upper boundary when the \textsf{ok\_strato}
    68 flag is set to {\em True} in \textsf{gcm.def}
    69 (this parameter also controls the application of a second step in the
    70 determination of vertical variation of coefficients for
    71 the horizontal dissipation, see \textsf{inidissip.F} and
    72 \textsf{disspi\_horiz.pdf} document).
     67%The sponge layer is applied at the upper boundary when the \textsf{ok\_strato}
     68%flag is set to {\em True} in \textsf{gcm.def}
     69%(this parameter also controls the application of a second step in the
     70%determination of vertical variation of coefficients for
     71%the horizontal dissipation, see \textsf{inidissip.F} and
     72%\textsf{disspi\_horiz.pdf} document).
    7373
    7474The tendencies for the upper boundary sponge layer are computed separately in

  • trunk/LMDZ.COMMON/libf/dyn3d/conf_gcm.F

    r1010 r1012  
    347347!                   iflag_top_bound=1 for sponge over 4 topmost layers
    348348!                   iflag_top_bound=2 for sponge from top to ~1% of top layer pressure
    349        iflag_top_bound=1
     349       iflag_top_bound=0
    350350       CALL getin('iflag_top_bound',iflag_top_bound)
    351351

  • trunk/LMDZ.COMMON/libf/dyn3d/leapfrog.F

    r1010 r1012  
    504504c      Couche superieure :
    505505c      -------------------
    506          IF (ok_strato) THEN
     506         IF (iflag_top_bound > 0) THEN
    507507           CALL top_bound(vcov,ucov,teta,masse,dtphys,dutop)
    508508           dutop(:,:)=dutop(:,:)/dtphys   ! convert to a tendency in (m/s)/s

  • trunk/LMDZ.COMMON/libf/dyn3d/top_bound.F

    r1010 r1012  
    7878      INTEGER j,l
    7979     
    80       if (iflag_top_bound.eq.0) return
    81 
    8280      if (first) then
    8381         if (iflag_top_bound.eq.1) then

  • trunk/LMDZ.COMMON/libf/dyn3dpar/conf_gcm.F

    r1010 r1012  
    374374!                   iflag_top_bound=1 for sponge over 4 topmost layers
    375375!                   iflag_top_bound=2 for sponge from top to ~1% of top layer pressure
    376        iflag_top_bound=1
     376       iflag_top_bound=0
    377377       CALL getin('iflag_top_bound',iflag_top_bound)
    378378

  • trunk/LMDZ.COMMON/libf/dyn3dpar/leapfrog_p.F

    r1010 r1012  
    10011001c      Couche superieure :
    10021002c      -------------------
    1003          IF (ok_strato) THEN
     1003         IF (iflag_top_bound > 0) THEN
    10041004           CALL top_bound_p(vcov,ucov,teta,masse,dtphys,
    10051005     $                       dutop)

  • trunk/LMDZ.COMMON/libf/dyn3dpar/top_bound_p.F

    r1010 r1012  
    7676      INTEGER j,l,jjb,jje
    7777
    78 
    79       if (iflag_top_bound == 0) return
    8078
    8179      if (first) then
Note: See TracChangeset for help on using the changeset viewer.