Changeset 310 for codes/icosagcm/branches

Timestamp:

12/15/14 09:22:51 (9 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

Location:

codes/icosagcm/branches/SATURN_DYNAMICO

Files:

• LMDZ.COMMON/libf/phystd/callkeys.h (modified) (2 diffs)
• LMDZ.COMMON/libf/phystd/inifis.F (modified) (1 diff)
• LMDZ.COMMON/libf/phystd/physiq.F90 (modified) (2 diffs)
• LMDZ.COMMON/libf/phystd/radsponge.F90 (added)
• TEST/callphys.def (modified) (1 diff)

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/callkeys.h

@@ -6 +6 @@
 
       COMMON/callkeys_l/callrad,corrk,calldifv,UseTurbDiff,calladj      &
-     &   , co2cond,callsoil                                             &
+     &   , co2cond,callsoil,callradsponge                               &
      &   , season,diurnal,tlocked,rings_shadow,lwrite                   &
      &   , callstats,calleofdump                                        &
@@ -50 +50 @@
      &   , strictboundcorrk                                             
 
+      logical callradsponge
       logical enertest
       logical nonideal

codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/phystd/inifis.F

@@ -204 +204 @@
          if (is_master) write(*,*) "Rmean = ", Rmean
          
+         if (is_master) write(*,*) 
+     &      "Add sponge after radiative transfer?"
+         callradsponge=.false.
+         call getin_p("callradsponge", callradsponge)
+         if (is_master) write(*,*) "callradsponge = ", callradsponge
+
 ! Test of incompatibility:
 ! if tlocked, then diurnal should be false

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)
 

codes/icosagcm/branches/SATURN_DYNAMICO/TEST/callphys.def

@@ -46 +46 @@
 #iradia     = 160
 iradia     = 40
+# Add temperature sponge effect after radiative transfer?
+callradsponge=.true.
 # call multilayer correlated-k radiative transfer ?
 corrk      = .true.