Changeset 1505

Timestamp:

Apr 7, 2011, 4:15:05 PM (13 years ago)

Author:

Ehouarn Millour

Message:

Added possibility to have the base of the atmosphere heated (to mimic a constant uniform heat flux from below), set in .def file (parameter ihf = ... , in W/m2); only active if in "Newtonian mode" (iflag_phys=2) and if planet_type=="giant".
EM

Location:

LMDZ5/trunk/libf

Files:

• dyn3d/comconst.h (modified) (2 diffs)
• dyn3d/conf_planete.F90 (modified) (1 diff)
• dyn3d/iniacademic.F90 (modified) (1 diff)
• dyn3d/leapfrog.F (modified) (1 diff)
• dyn3dpar/comconst.h (modified) (2 diffs)
• dyn3dpar/conf_planete.F90 (modified) (1 diff)
• dyn3dpar/iniacademic.F90 (modified) (1 diff)
• dyn3dpar/leapfrog_p.F (modified) (1 diff)

LMDZ5/trunk/libf/dyn3d/comconst.h

@@ -11 +11 @@
      &                   ,dissip_factz,dissip_deltaz,dissip_zref        &
      &                   ,tau_top_bound,                                &
-     & daylen,year_day,molmass
+     & daylen,year_day,molmass, ihf
 
 
@@ -35 +35 @@
       REAL molmass ! (g/mol) molar mass of the atmosphere
 
+      REAL ihf ! (W/m2) Intrinsic heat flux (for giant planets)
 
 !-----------------------------------------------------------------------

LMDZ5/trunk/libf/dyn3d/conf_planete.F90

@@ -64 +64 @@
 CALL getin('omeg',omeg)
 
+! Intrinsic heat flux (default: none) (only used if planet_type="giant")
+ihf = 0.
+call getin('ihf',ihf)
+
 END SUBROUTINE conf_planete

LMDZ5/trunk/libf/dyn3d/iniacademic.F90

@@ -178 +178 @@
            tetajl(j,l)=teta0-delt_y*ddsin*ddsin+eps*ddsin &
                 -delt_z*(1.-ddsin*ddsin)*log(zsig)
+           if (planet_type=="giant") then
+             tetajl(j,l)=teta0+(delt_y*                   &
+                ((sin(rlatu(j)*3.14159*eps+0.0001))**2)   &
+                / ((rlatu(j)*3.14159*eps+0.0001)**2))     &
+                -delt_z*log(zsig)
+           endif
            ! Profil stratospherique isotherme (+vortex)
            w_pv=(1.-tanh((rlatu(j)-phi_pv)/dphi_pv))/2.

LMDZ5/trunk/libf/dyn3d/leapfrog.F

@@ -437 +437 @@
           ENDDO
         ENDDO ! of DO l=1,llm 
-          call friction(ucov,vcov,dtvr)
+        
+        if (planet_type.eq."giant") then
+          ! add an intrinsic heat flux at the base of the atmosphere
+          teta(:,1)=teta(:,1)+dtvr*aire(:)*ihf/cpp/masse(:,1)
+        endif
+
+        call friction(ucov,vcov,dtvr)
         
         ! Sponge layer (if any)

LMDZ5/trunk/libf/dyn3dpar/comconst.h

@@ -11 +11 @@
      &                   ,dissip_factz,dissip_deltaz,dissip_zref        &
      &                   ,tau_top_bound,                                &
-     & daylen,year_day,molmass
+     & daylen,year_day,molmass, ihf
 
 
@@ -35 +35 @@
       REAL molmass ! (g/mol) molar mass of the atmosphere
 
+      REAL ihf ! (W/m2) Intrinsic heat flux (for giant planets)
 
 !-----------------------------------------------------------------------

LMDZ5/trunk/libf/dyn3dpar/conf_planete.F90

@@ -64 +64 @@
 CALL getin('omeg',omeg)
 
+! Intrinsic heat flux (default: none) (only used if planet_type="giant")
+ihf = 0.
+call getin('ihf',ihf)
+
 END SUBROUTINE conf_planete

LMDZ5/trunk/libf/dyn3dpar/iniacademic.F90

@@ -178 +178 @@
            tetajl(j,l)=teta0-delt_y*ddsin*ddsin+eps*ddsin &
                 -delt_z*(1.-ddsin*ddsin)*log(zsig)
+           if (planet_type=="giant") then
+             tetajl(j,l)=teta0+(delt_y*                   &
+                ((sin(rlatu(j)*3.14159*eps+0.0001))**2)   &
+                / ((rlatu(j)*3.14159*eps+0.0001)**2))     &
+                -delt_z*log(zsig)
+           endif
            ! Profil stratospherique isotherme (+vortex)
            w_pv=(1.-tanh((rlatu(j)-phi_pv)/dphi_pv))/2.

LMDZ5/trunk/libf/dyn3dpar/leapfrog_p.F

@@ -977 +977 @@
 !$OMP END DO
 
+!$OMP MASTER
+       if (planet_type.eq."giant") then
+         ! add an intrinsic heat flux at the base of the atmosphere
+         teta(ijb:ije,1) = teta(ijb:ije,1)
+     &        + dtvr * aire(ijb:ije) * ihf / cpp / masse(ijb:ije,1)
+       endif
+!$OMP END MASTER
+!$OMP BARRIER
+
        call Register_Hallo(ucov,ip1jmp1,llm,0,1,1,0,Request_Physic)
        call Register_Hallo(vcov,ip1jm,llm,1,1,1,1,Request_Physic)