Ignore:
Timestamp:
Mar 16, 2021, 11:28:16 AM (4 years ago)
Author:
aboissinot
Message:

Commit the last changes in the thermal plume model (more english comments).

Location:
trunk/LMDZ.GENERIC/libf/phystd
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.GENERIC/libf/phystd/thermcell_flux.F90

    r2177 r2480  
    177177            print *, '- - - - - - - - - - - - - - - - - - - - - - - - - - - - -'
    178178            DO k=nlay,1,-1
    179                print *, 'k, fm ', k+1, fm(igout,k+1)
     179               print *, 'fm,w ', fm(igout,k+1), zw2(igout,k+1)
    180180               print *, 'entr,detr', entr(igout,k), detr(igout,k)
    181181            ENDDO
    182             print *, 'k, fm ', 1, fm(igout,1)
     182            print *, 'fm,w ', fm(igout,1), zw2(igout,1)
    183183            print *, '---------------------------------------------------------'
    184184            CALL abort
    … …  
    222222                  fm(ig,l+1) = fm(ig,l) + entr(ig,l)
    223223                  entr(ig,l+1) = entr(ig,l+1) + ddd
     224               ELSE
    224225!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    225226! AB: Simulation abort (try to reduce the physical time step).
    226227!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    227                ELSE
    228                   entr(ig,l) = entr(ig,l) + eee
    229                   igout = ig
    230                   lout = l
    231                   labort_physic = .true.
     228!                  entr(ig,l) = entr(ig,l) + eee
     229!                  igout = ig
     230!                  lout = l
     231!                  labort_physic = .true.
    232232!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    233233! AB: We can renormalize the plume mass fluxes. I think it does not work.
    234234!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    235 !               ELSE
    236235!                  fact = max(fact, eee0 / emax)
    237 !                  entr(ig,l) = eee0
    238 !                  ncorecfact = ncorecfact + 1
     236                  fact = eee0 / emax
     237                  entr(ig,l) = eee0
     238                  ncorecfact = ncorecfact + 1
    239239!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    240240! AB: The renormalization can be just applied in the local plume.
    241241!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    242 !                  DO k=lmin(ig),lmax(ig)
    243 !                     entr(ig,k) = entr(ig,k) * emax / eee0
    244 !                     detr(ig,k) = detr(ig,k) * emax / eee0
    245 !                     fm(ig,k) = fm(ig,k) * emax / eee0
    246 !                  ENDDO
     242                  DO k=lmin(ig),lmax(ig)
     243                     entr(ig,k) = entr(ig,k) * emax / eee0
     244                     detr(ig,k) = detr(ig,k) * emax / eee0
     245                     fm(ig,k) = fm(ig,k) * emax / eee0
     246                  ENDDO
    247247               ENDIF
    248248            ENDIF
    … …  
    264264            print *, '- - - - - - - - - - - - - - - - - - - - - - - - - - - - -'
    265265            DO k=nlay,1,-1
    266                print *, 'fm ', fm(igout,k+1)
     266               print *, 'fm,w ', fm(igout,k+1), zw2(igout,k+1)
    267267               print *, 'entr,detr', entr(igout,k), detr(igout,k)
    268268            ENDDO
    269             print *, 'fm ', fm(igout,1)
     269            print *, 'fm,w ', fm(igout,1), zw2(igout,1)
    270270            print *, '---------------------------------------------------------'
    271271            CALL abort

  • trunk/LMDZ.GENERIC/libf/phystd/thermcell_height.F90

    r2177 r2480  
    22!
    33!
    4 SUBROUTINE thermcell_height(ngrid,nlay,lmin,lmix,lmax,zlev,                   &
    5                             zmin,zmix,zmax,zw2,wmax,f_star)
     4SUBROUTINE thermcell_height(ngrid,nlay,lmin,lmix,lmax,                        &
     5                            zlev,zmin,zmix,zmax,zw2,wmax,f_star)
    66     
    77     

  • trunk/LMDZ.GENERIC/libf/phystd/thermcell_main.F90

    r2232 r2480  
    383383     
    384384!-------------------------------------------------------------------------------
    385 ! Calcul du transport vertical de la temperature potentielle
     385! Potential temperature
    386386!-------------------------------------------------------------------------------
    387387     
    … …  
    396396     
    397397!-------------------------------------------------------------------------------
    398 ! Calcul du transport vertical des traceurs
     398! Tracers
    399399!-------------------------------------------------------------------------------
    400400     
    … …  
    405405     
    406406!-------------------------------------------------------------------------------
    407 ! Calcul du transport vertical du moment horizontal
     407! Horizontal momentum
    408408!-------------------------------------------------------------------------------
    409409     

  • trunk/LMDZ.GENERIC/libf/phystd/thermcell_mod.F90

    r2144 r2480  
    1313      REAL,SAVE :: r_aspect_thermals             !  1.0      Aspect ratio of the thermals (width / height)
    1414      REAL,SAVE :: tau_thermals                  !  0.       Relaxation time
    15       REAL,SAVE :: betalpha                      !  0.9      - between 0 (e=d) and 1 (rho*fraca=cst)
     15      REAL,SAVE :: betalpha                      !  0.9      - between 0 (e=d=0) and 1 (rho*fraca=cst)
    1616      REAL,SAVE :: afact                         !  2./3.    - buoyancy acceleration efficiency, between 0 and 1
    1717      REAL,SAVE :: fact_epsilon                  !  2.e-3    - turbulence and friction deceleration
Note: See TracChangeset for help on using the changeset viewer.