Changeset 2743 for trunk/MESOSCALE

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:

• LES/modif_mars/Registry.EM.newphys (modified) (1 diff)
• LES/modif_mars/module_initialize_les.F (modified) (4 diffs)
• WRFV2/phys/dynphy_wrf_titan_lmd_new/update_inputs_physiq_mod.F (modified) (3 diffs)
• WRFV2/phys/dynphy_wrf_titan_lmd_new/variables_mod.F (modified) (1 diff)
• WRFV2/phys/module_lmd_driver.F.new (modified) (1 diff)

trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM.newphys

@@ -108 +108 @@
 state  real  GSNOW  ikj  misc  1  -  -  "SNOW"     "SNOW kg/kg.s-1"           ""
 state  real  ALBEQ  ij misc  1  -  -  "ALBEQ"    "ALBEDO EQUIVALENT"   ""
-state  real  FLUXTOP_DN ij misc  1  -  - "FLUXTOP_DN"   "Incoming SW radiation at the TOA" "W.m-2"
-state  real  FLUXABS_SW ij misc  1  -  - "FLUXABS_SW"  "Absorbed SW flux" "W.m-2"
-state  real  FLUXTOP_LW ij misc  1  -  - "FLUXTOP_LW"  "Outgoing LW flux TOA" "W.m-2"
-state  real  FLUXSURF_SW ij misc  1  -  - "FLUXSURF_SW" "SW Surface flux" "W.m-2"
-state  real  FLUXSURF_LW ij misc  1  -  -  "FLUXSURF_LW" "LW Surface flux" "W.m-2"
+state  real  FLUXTOP_DN ij misc  1  -  h "FLUXTOP_DN"   "Incoming SW radiation at the TOA" "W.m-2"
+state  real  FLUXABS_SW ij misc  1  -  h "FLUXABS_SW"  "Absorbed SW flux" "W.m-2"
+state  real  FLUXTOP_LW ij misc  1  -  h "FLUXTOP_LW"  "Outgoing LW flux TOA" "W.m-2"
+state  real  FLUXSURF_SW ij misc  1  -  h "FLUXSURF_SW" "SW Surface flux" "W.m-2"
+state  real  FLUXSURF_LW ij misc  1  -  h  "FLUXSURF_LW" "LW Surface flux" "W.m-2"
 state  real  FLXGRD ij misc  1  -  h "FLXGRD" "heat flux from ground"  "W.m-2"
 state  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

@@ -176 +176 @@
    END SELECT
 
-   IF (planet == "mars") THEN
+   IF (planet == "mars" .or. planet == "titan") THEN
      stretch_grid = .false.
      !! FOR LES, set stretch to false
@@ -259 +259 @@
      !!MARS
      !!MARS
-     IF (planet .ne. "mars") THEN
+     IF (planet .ne. "mars" .and. planet .ne. "titan") THEN
        open(unit=12,file='levels',form='formatted',status='old')
        rewind(12)
@@ -276 +276 @@
 
    !! SPECIFIC FOR LES PBL MARS
-   IF (planet == "mars") THEN
+   IF (planet == "mars" .or. planet == "titan") THEN
      !!!MARS
      grid%znw(1)=1.000
@@ -286 +286 @@
      ENDDO
    ENDIF
+
+   open(unit=12,file='levels',form='formatted',status='old')
+   rewind(12)
+   DO k=1, kde
+     write(12,*) grid%znw(k)
+     write(6,*) 'for update_inputs_physiq_mod (e.g. Titan, generic)',k,grid%znw(k)
+   ENDDO
+   close(12)
+
 
 

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