Changeset 1301 for trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.NewtonCool

Timestamp:

Jun 26, 2014, 12:25:31 PM (11 years ago)

Author:

slebonnois

Message:

SL: many bug corrections in phyvenus, some cleaning, and a new ksi matrix format for Venus IR

File:

• trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.NewtonCool (modified) (4 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.NewtonCool

@@ -2 +2 @@
 ! $Header: /home/cvsroot/LMDZ4/libf/phylmd/radlwsw.F,v 1.2 2004/10/27 10:14:46 lmdzadmin Exp $
 !
-      SUBROUTINE radlwsw(dist, rmu0, fract, 
-     .                  paprs, pplay,tsol, t,
-     .                  heat,cool,radsol,
-     .                  topsw,toplw,solsw,sollw,
-     .                  sollwdown,
-     .                  lwnet, swnet)
-c      
+      SUBROUTINE radlwsw(dist, rmu0, fract, zzlev,
+     .                  paprs, pplay,tsol, pt, nq, nmicro, pq,qaer)
+      
 c======================================================================
 c Auteur(s): Z.X. Li (LMD/CNRS) date: 19960719
@@ -16 +12 @@
 c rmu0-----input-R- cosinus de l'angle zenithal
 c fract----input-R- duree d'ensoleillement normalisee
-c solaire--input-R- constante solaire (W/m**2) (dans clesphys.h)
 c paprs----input-R- pression a inter-couche (Pa)
 c pplay----input-R- pression au milieu de couche (Pa)
 c tsol-----input-R- temperature du sol (en K)
-c t--------input-R- temperature (K)
-c heat-----output-R- echauffement atmospherique (visible) (K/jour)
-c cool-----output-R- refroidissement dans l'IR (K/jour)
-c radsol---output-R- bilan radiatif net au sol (W/m**2) (+ vers le bas)
-c topsw----output-R- flux solaire net au sommet de l'atm. (+ vers le bas)
-c toplw----output-R- ray. IR net au sommet de l'atmosphere (+ vers le haut)
-c solsw----output-R- flux solaire net a la surface (+ vers le bas)
-c sollw----output-R- ray. IR net a la surface (+ vers le bas)
-c sollwdown-output-R- ray. IR descendant a la surface (+ vers le bas)
-c lwnet____output-R- flux IR net (+ vers le haut)
-c swnet____output-R- flux solaire net (+ vers le bas)
+c pt-------input-R- temperature (K)
 c
       
 c   S. Lebonnois    12/04/2007
 c  VERSION NEWTONIAN COOLING pour Venus (no diurnal cycle)
+c  update 01/2014
 
 c======================================================================
       use dimphy
       USE comgeomphy
+      USE phys_state_var_mod, only: falbe,heat,cool,radsol,
+     .      topsw,toplw,solsw,sollw,sollwdown,lwnet,swnet
+      USE write_field_phy
       IMPLICIT none
 #include "dimensions.h"
 #include "YOMCST.h"
 #include "clesphys.h" 
-c
+
+c ARGUMENTS
+      INTEGER nq,nmicro
       real rmu0(klon), fract(klon), dist
-c
-      real paprs(klon,klev+1), pplay(klon,klev)
+ 
+      real zzlev(klon,klev+1),paprs(klon,klev+1), pplay(klon,klev)
       real tsol(klon)
-      real t(klon,klev)
-      real heat(klon,klev), cool(klon,klev)
-      real radsol(klon), topsw(klon), toplw(klon)
-      real solsw(klon), sollw(klon)
-      real sollwdown(klon)
-      REAL swnet(klon,klev+1),lwnet(klon,klev+1)
-c
+      real pt(klon,klev)
+      real pq(klon,klev,nq)
+      REAL qaer(klon,klev,nq)
+ 
+c LOCAL VARIABLES
       INTEGER i,j,k
       integer    nlevCLee,level
@@ -64 +53 @@
       real   ztemp,zdt,fact
       real   dTsdt(klev),zt_eq(klon,klev)
-      real   dureejour
-      parameter (dureejour=10.087e6)
       save   zt_eq
       
@@ -144 +131 @@
 
         DO k = 1, klev
-         heat (j,k) = dTsdt(k)*dureejour ! K/Venusday
+         heat (j,k) = dTsdt(k)    ! K/s
          cool (j,k) = 0.
         ENDDO