Changeset 1505 for LMDZ5/trunk/libf/dyn3dpar/iniacademic.F90

Timestamp:

Apr 7, 2011, 4:15:05 PM (14 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

File:

• LMDZ5/trunk/libf/dyn3dpar/iniacademic.F90 (modified) (1 diff)

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.