Ignore:
Timestamp:
Jul 5, 2022, 11:56:27 PM (2 years ago)
Author:
aslmd
Message:

made the vertical discretization for Titan LES similar to that of Mars (simple linear law for znw and refinement close to the surface). no more need to prepare a file levels, just put ztop and number of levels in namelist. a file levels is saved at the end of the znw computation so that it can be read in update_inputs_physiq_mod in the case of Titan.

added MESOSCALE precompiling flags in call_profilgases to set constant methane abundance. therefore no need for profile.def file and the number of vertical levels is set by namelist -- no need to set the same number in both levels and profile.def

outputing ustar from vdifc

ending hardcoding of ptop in update_inputs_physiq_mod for Titan and generic LES/mesoscale. streamlining grod%ptop through module_lmd_driver as ptopwrf within variables_mod module.

adding flux outputs to check surface and atmosphere energy budgets

Location:
trunk/MESOSCALE/LMD_MM_MARS/SRC
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM.newphys

    r2208 r2743  
    108108state  real  GSNOW  ikj  misc  1  -  -  "SNOW"     "SNOW kg/kg.s-1"           ""
    109109state  real  ALBEQ  ij misc  1  -  -  "ALBEQ"    "ALBEDO EQUIVALENT"   ""
    110 state  real  FLUXTOP_DN ij misc  1  -  - "FLUXTOP_DN"   "Incoming SW radiation at the TOA" "W.m-2"
    111 state  real  FLUXABS_SW ij misc  1  -  - "FLUXABS_SW"  "Absorbed SW flux" "W.m-2"
    112 state  real  FLUXTOP_LW ij misc  1  -  - "FLUXTOP_LW"  "Outgoing LW flux TOA" "W.m-2"
    113 state  real  FLUXSURF_SW ij misc  1  -  - "FLUXSURF_SW" "SW Surface flux" "W.m-2"
    114 state  real  FLUXSURF_LW ij misc  1  -  -  "FLUXSURF_LW" "LW Surface flux" "W.m-2"
     110state  real  FLUXTOP_DN ij misc  1  -  h "FLUXTOP_DN"   "Incoming SW radiation at the TOA" "W.m-2"
     111state  real  FLUXABS_SW ij misc  1  -  h "FLUXABS_SW"  "Absorbed SW flux" "W.m-2"
     112state  real  FLUXTOP_LW ij misc  1  -  h "FLUXTOP_LW"  "Outgoing LW flux TOA" "W.m-2"
     113state  real  FLUXSURF_SW ij misc  1  -  h "FLUXSURF_SW" "SW Surface flux" "W.m-2"
     114state  real  FLUXSURF_LW ij misc  1  -  h  "FLUXSURF_LW" "LW Surface flux" "W.m-2"
    115115state  real  FLXGRD ij misc  1  -  h "FLXGRD" "heat flux from ground"  "W.m-2"
    116116state  real  LSCEZ ikj misc  1  -  - "LSCEZ" "heat from largescale" "W.m-2"
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F

    r2021 r2743  
    176176   END SELECT
    177177
    178    IF (planet == "mars") THEN
     178   IF (planet == "mars" .or. planet == "titan") THEN
    179179     stretch_grid = .false.
    180180     !! FOR LES, set stretch to false
     
    259259     !!MARS
    260260     !!MARS
    261      IF (planet .ne. "mars") THEN
     261     IF (planet .ne. "mars" .and. planet .ne. "titan") THEN
    262262       open(unit=12,file='levels',form='formatted',status='old')
    263263       rewind(12)
     
    276276
    277277   !! SPECIFIC FOR LES PBL MARS
    278    IF (planet == "mars") THEN
     278   IF (planet == "mars" .or. planet == "titan") THEN
    279279     !!!MARS
    280280     grid%znw(1)=1.000
     
    286286     ENDDO
    287287   ENDIF
     288
     289   open(unit=12,file='levels',form='formatted',status='old')
     290   rewind(12)
     291   DO k=1, kde
     292     write(12,*) grid%znw(k)
     293     write(6,*) 'for update_inputs_physiq_mod (e.g. Titan, generic)',k,grid%znw(k)
     294   ENDDO
     295   close(12)
     296
    288297
    289298
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_titan_lmd_new/update_inputs_physiq_mod.F

    r2742 r2743  
    187187   use comgeomfi_h, only: totarea, totarea_planet
    188188   use planete_mod, only: ini_planete_mod
     189   ! in LMD-WRF interface
     190   USE variables_mod, only: ptopwrf
    189191
    190192   INTEGER, INTENT(IN) :: ims,ime,jms,jme
     
    198200   REAL, DIMENSION(nlayer+1) :: znw
    199201   REAL*8, DIMENSION(nlayer+1) :: apdyn,bpdyn
    200    REAL :: ptop
    201202   REAL*8 :: PP0
    202203   DO j = jps,jpe
     
    280281   close(12)
    281282   
    282    ptop=0.5
    283    apdyn=ptop*(1-znw)
     283   apdyn=ptopwrf*(1-znw)
    284284   bpdyn=znw
    285285   PP0=p0
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_titan_lmd_new/variables_mod.F

    r2295 r2743  
    77REAL*8 :: zdt_split !ptimestep !zdt_split ! time step over which the physics are evaluated
    88REAL*8 :: phour_ini ! start time (fraction of day) of the run 0=<phour_ini<1
     9REAL*8 :: ptopwrf ! pressure at top of model in WRF [should be equal to grid%p_top]
    910REAL*8,DIMENSION(:,:),ALLOCATABLE :: zplev_omp !pplev !zplev_omp(klon,llm+1) ! interlayer pressure (Pa)
    1011REAL*8,DIMENSION(:,:),ALLOCATABLE :: zplay_omp !pplay !zplay_omp(klon,llm) ! mid-layer pressure (Pa)
  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F.new

    r2208 r2743  
    264264wappel_phys = RADT
    265265zdt_split = dt*wappel_phys            ! physical timestep (s)
     266ptopwrf = ptop                        ! top pressure in WRF coordinates (in variables_mod)
    266267ngrid=(ipe-ips+1)*(jpe-jps+1)         ! size of the horizontal grid
    267268nlayer = kpe-kps+1                    ! number of vertical layers: nlayermx
Note: See TracChangeset for help on using the changeset viewer.