Changeset 2743 for trunk/LMDZ.TITAN/libf

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/LMDZ.TITAN/libf/phytitan

Files:

• call_profilgases.F90 (modified) (3 diffs)
• vdifc.F (modified) (2 diffs)

trunk/LMDZ.TITAN/libf/phytitan/call_profilgases.F90

@@ -43 +43 @@
 
 !$OMP MASTER
-     
+#ifndef MESOSCALE     
      ! Load CH4 vertical profile from file 'profile.def'
      open(90,file='profile.def',status='old',form='formatted',iostat=ierr)
@@ -60 +60 @@
            endif          
         enddo        
-        
+#else        
+        gfrac(igas_CH4,:) = 1.443E-02
+        print*,"MESOSCALE. setting constant CH4",gfrac(igas_CH4,1)
+#endif
         ! Then set H2 (fixed) and N2 (what remains)
 
         gfrac(igas_H2,:) = 1.0E-3
         gfrac(igas_N2,:) = 1.0 - ( gfrac(igas_CH4,:) + gfrac(igas_H2,:) )
-        
+
+
+#ifndef MESOSCALE        
      else
         write(*,*) 'Cannot find required file "profile.def"'
@@ -72 +77 @@
 
      close(90)
+#endif
 !$OMP END MASTER
 !$OMP BARRIER

trunk/LMDZ.TITAN/libf/phytitan/vdifc.F

@@ -12 +12 @@
       use comcstfi_mod, only: g, r, cpp, rcp
       use callkeys_mod, only: tracer,nosurf
+      use turb_mod, only: ustar
 
       implicit none
@@ -203 +204 @@
          zcdv(ig)=zcdv_true(ig)*sqrt(zu2)
          zcdh(ig)=zcdh_true(ig)*sqrt(zu2)
+         ustar(ig)=sqrt(zcdv_true(ig))*sqrt(zu2)
       ENDDO