Changeset 2743 for trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2

Timestamp:

Jul 5, 2022, 11:56:27 PM (3 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/WRFV2/phys

Files:

• dynphy_wrf_titan_lmd_new/update_inputs_physiq_mod.F (modified) (3 diffs)
• dynphy_wrf_titan_lmd_new/variables_mod.F (modified) (1 diff)
• module_lmd_driver.F.new (modified) (1 diff)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_titan_lmd_new/update_inputs_physiq_mod.F

@@ -187 +187 @@
    use comgeomfi_h, only: totarea, totarea_planet
    use planete_mod, only: ini_planete_mod
+   ! in LMD-WRF interface
+   USE variables_mod, only: ptopwrf
 
    INTEGER, INTENT(IN) :: ims,ime,jms,jme
@@ -198 +200 @@
    REAL, DIMENSION(nlayer+1) :: znw
    REAL*8, DIMENSION(nlayer+1) :: apdyn,bpdyn
-   REAL :: ptop
    REAL*8 :: PP0
    DO j = jps,jpe
@@ -280 +281 @@
    close(12)
    
-   ptop=0.5
-   apdyn=ptop*(1-znw)
+   apdyn=ptopwrf*(1-znw)
    bpdyn=znw
    PP0=p0

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_titan_lmd_new/variables_mod.F

@@ -7 +7 @@
 REAL*8 :: zdt_split !ptimestep !zdt_split ! time step over which the physics are evaluated
 REAL*8 :: phour_ini ! start time (fraction of day) of the run 0=<phour_ini<1
+REAL*8 :: ptopwrf ! pressure at top of model in WRF [should be equal to grid%p_top]
 REAL*8,DIMENSION(:,:),ALLOCATABLE :: zplev_omp !pplev !zplev_omp(klon,llm+1) ! interlayer pressure (Pa)
 REAL*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

@@ -264 +264 @@
 wappel_phys = RADT
 zdt_split = dt*wappel_phys            ! physical timestep (s)
+ptopwrf = ptop                        ! top pressure in WRF coordinates (in variables_mod)
 ngrid=(ipe-ips+1)*(jpe-jps+1)         ! size of the horizontal grid
 nlayer = kpe-kps+1                    ! number of vertical layers: nlayermx