Changeset 2743 for trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars

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/LES/modif_mars

Files:

• Registry.EM.newphys (modified) (1 diff)
• module_initialize_les.F (modified) (4 diffs)

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)
+