Ignore:
Timestamp:
Aug 28, 2018, 4:22:24 PM (6 years ago)
Author:
jleconte
Message:

28/08/2018 == JL

We now shift the radiative model top from p=0 to the middle of the last physical layer. This is done by changing pmid and plevrad in callcorrk and some corrections need to be done in gfluxv.
This seems to get rid of the aratic temperature behavior in the last two layers of the model (especially on the night side on synchronous planets).
Additional speedup corrections have been made in gfluxi that change nothing to the result.
Finally, if aerosols are present in the last layer we must account for them. Provides better upper boundary condition in the IR. They must however be put to zero in the sw (see optcv and changes in last commit.)
This has been done for water ice in aeropacity, but same correction should probably be done for other aerosol types.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/README

    r1987 r1988  
    13781378Tauaero and tauray are set to 0 (a small value for rayleigh because the code crashes otherwise) in the 4 first semilayers in optcv, but not optci.
    13791379This solves random variations of the sw heating at the model top.
     1380
     1381== 28/08/2018 == JL
     1382We now shift the radiative model top from p=0 to the middle of the last physical layer. This is done by changing pmid and plevrad in callcorrk and some corrections need to be done in gfluxv.
     1383This seems to get rid of the aratic temperature behavior in the last two layers of the model (especially on the night side on synchronous planets).
     1384Additional speedup corrections have been made in gfluxi that change nothing to the result.
     1385Finally, if aerosols are present in the last layer we must account for them. Provides better upper boundary condition in the IR. They must however be put to zero in the sw (see optcv and changes in last commit.)
     1386This has been done for water ice in aeropacity, but same correction should probably be done for other aerosol types.
     1387
Note: See TracChangeset for help on using the changeset viewer.