Changeset 2987 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

Jul 3, 2023, 5:22:09 PM (17 months ago)

Author:

emillour

Message:

Fixed the hardcoded value for the altitude of topmost level

In physiq_mod, the variable zzlev(:,nlayer+1) was hardcoded to 1e7. This is not enough when simulating the atmosphere of inflated hot Jupiter. Thus, to compute this value, we now use the thickness of the level below the top one to compute the altitude of topmost level. This variable is used only in the turbulent diffusion (we use 1/thickness) and the sedimentation (use of thickness).

File:

• trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90 (modified) (2 diffs)

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

@@ -845 +845 @@
 
       zzlev(1:ngrid,1)=0.
-      zzlev(1:ngrid,nlayer+1)=1.e7 ! Dummy top of last layer above 10000 km...
 
       do l=2,nlayer
@@ -854 +853 @@
          enddo
       enddo     
+
+      !Altitude of top interface (nlayer+1), using the thicknesss of the level below the top one. LT22
+      
+      zzlev(1:ngrid,nlayer+1) = 2*zzlev(1:ngrid,nlayer)-zzlev(1:ngrid,nlayer-1)
 
       ! Compute potential temperature