Ignore:
Timestamp:
Jan 30, 2017, 11:24:47 AM (8 years ago)
Author:
Ehouarn Millour
Message:

Further convergence with planetary models:

  • move vertical_layers_mod.F90 to phy_common and call init_vertical_layers in inigeomphy_mod
  • add aps and bps (mid-layer coordinates) to vertical_layers_mod.F90
  • same adaptations for the 1D case

EM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/dyn1d/lmdz1d.F90

    r2716 r2786  
    4040   USE physiq_mod, ONLY: physiq
    4141   USE comvert_mod, ONLY: presnivs, ap, bp, dpres,nivsig, nivsigs, pa, &
    42                           preff
     42                          preff, aps, bps, pseudoalt, scaleheight
     43   USE vertical_layers_mod, ONLY: init_vertical_layers
    4344   USE temps_mod, ONLY: annee_ref, calend, day_end, day_ini, day_ref, &
    4445                        itau_dyn, itau_phy, start_time
     
    634635        call disvert0(pa,preff,ap,bp,dpres,presnivs,nivsigs,nivsig)
    635636        print *,'On utilise disvert0'
     637        aps(1:llm)=0.5*(ap(1:llm)+ap(2:llm+1))
     638        bps(1:llm)=0.5*(bp(1:llm)+bp(2:llm+1))
     639        scaleheight=8.
     640        pseudoalt(1:llm)=-scaleheight*log(presnivs(1:llm)/preff)
    636641      ELSE
    637642        call disvert()
     
    641646      ENDIF
    642647      ! initialize ap,bp, etc. in vertical_layers_mod
     648      CALL init_vertical_layers(llm,preff,scaleheight, &
     649                                ap,bp,aps,bps,presnivs,pseudoalt)
     650
    643651      sig_s=presnivs/preff
    644652      plev =ap+bp*psurf
Note: See TracChangeset for help on using the changeset viewer.