Changeset 310 for codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/physiq.F90

Timestamp:

12/15/14 09:22:51 (10 years ago)

Author:

millour

Message:

Add the possibility to have a sponge acting on the temperature
in the physics, using flag "callradsponge=.true.". The sponge
vertical extention is the same as the sponge layer in the dynamics.
EM

File:

• codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/physiq.F90 (modified) (2 diffs)

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/physiq.F90

@@ -280 +280 @@
       real zdvmolvis(ngrid,nlayer)
       real zdqmoldiff(ngrid,nlayer,nq)
+      real dtradsponge(ngrid,nlayer) ! (K/s)
 
 !     Local variables for local calculations:
@@ -1088 +1089 @@
 !-------------------------
 
+        if (callradsponge) then
+          call radsponge(ngrid,nlayer,pplay,pplev,pt,pdt,dtphys,dtradsponge)
+          pdt(1:ngrid,1:nlayer)=pdt(1:ngrid,1:nlayer)+dtradsponge(1:ngrid,1:nlayer)
+        endif
+
       endif ! of if (callrad)