source: trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.NewtonCool @ 1723

Last change on this file since 1723 was 1545, checked in by emillour, 9 years ago

Venus and Titan GCMs:

Adaptation wrt previous changes for Titan and Venus where
longitude and latitude arrays (in phycommon/geometry_mod) were overwritten
with values from startphy.nc files, where values are given in degrees.
For the sake of homegeneity with other physics package, revert to "default"
behaviour: longitude/latitude are in radians and longitude_deg/latitude_deg
are in degrees.
Also added checking in phyetat0 that the longitude/latitude read in the
restartphy.nc files match the ones provided by the dynamics.

EM
File size: 4.5 KB
Line 
1!
2! $Header: /home/cvsroot/LMDZ4/libf/phylmd/radlwsw.F,v 1.2 2004/10/27 10:14:46 lmdzadmin Exp $
3!
4      SUBROUTINE radlwsw(dist, rmu0, fract, zzlev,
5     .                  paprs, pplay,tsol, pt)
6     
7c======================================================================
8c Auteur(s): Z.X. Li (LMD/CNRS) date: 19960719
9c Objet: interface entre le modele et les rayonnements
10c Arguments:
11c dist-----input-R- distance astronomique terre-soleil
12c rmu0-----input-R- cosinus de l'angle zenithal
13c fract----input-R- duree d'ensoleillement normalisee
14c paprs----input-R- pression a inter-couche (Pa)
15c pplay----input-R- pression au milieu de couche (Pa)
16c tsol-----input-R- temperature du sol (en K)
17c pt-------input-R- temperature (K)
18c
19     
20c   S. Lebonnois    12/04/2007
21c  VERSION NEWTONIAN COOLING pour Venus (no diurnal cycle)
22c  update 01/2014
23
24c======================================================================
25      use dimphy
26      USE geometry_mod, ONLY: latitude ! in radians
27      USE phys_state_var_mod, only: heat,cool,radsol,
28     .      topsw,toplw,solsw,sollw,sollwdown,lwnet,swnet
29      USE write_field_phy
30      IMPLICIT none
31#include "YOMCST.h"
32#include "clesphys.h"
33
34c ARGUMENTS
35      real rmu0(klon), fract(klon), dist
36 
37      real zzlev(klon,klev+1),paprs(klon,klev+1), pplay(klon,klev)
38      real tsol(klon)
39      real pt(klon,klev)
40 
41c LOCAL VARIABLES
42      INTEGER i,j,k
43      integer    nlevCLee,level
44      parameter (nlevCLee=30)
45      REAL   pressCLee(nlevCLee+1),tempCLee(nlevCLee+1)
46      real   dt_epCLee(nlevCLee+1),etaCLee(nlevCLee+1)
47      real   tauCLee
48      parameter (tauCLee=25*86400) ! en s
49      real   ztemp,zdt,fact
50      real   dTsdt(klev),zt_eq(klon,klev)
51      save   zt_eq
52     
53      data     etaCLee/9.602e-1,8.679e-1,7.577e-1,6.420e-1,5.299e-1,
54     .                 4.273e-1,3.373e-1,2.610e-1,1.979e-1,1.472e-1,
55     .                 1.074e-1,7.672e-2,5.361e-2,3.657e-2,2.430e-2,
56     .                 1.569e-2,9.814e-3,5.929e-3,3.454e-3,1.934e-3,
57     .                 1.043e-3,5.400e-4,2.710e-4,1.324e-4,6.355e-5,
58     .                 3.070e-5,1.525e-5,7.950e-6,4.500e-6,2.925e-6,
59     .                 2.265e-6/
60      data   tempCLee/728.187,715.129,697.876,677.284,654.078,628.885,
61     .                602.225,574.542,546.104,517.339,488.560,459.932,
62     .                431.741,404.202,377.555,352.042,327.887,305.313,
63     .                284.556,265.697,248.844,233.771,220.368,208.247,
64     .                197.127,187.104,178.489,171.800,167.598,165.899,
65     .                165.676/
66      data   dt_epCLee/6.101 , 6.136 , 6.176 , 6.410 , 6.634 , 6.678 ,
67     .                 6.719 , 6.762 , 7.167 , 7.524 , 9.840 ,14.948 ,
68     .                21.370 ,28.746 ,36.373 ,43.315 ,48.534 ,51.175 ,
69     .                50.757 ,47.342 ,41.536 ,34.295 ,26.758 ,19.807 ,
70     .                14.001 , 9.599 , 6.504 , 4.439 , 3.126 , 2.370 ,
71     .                2.000/
72c
73
74      logical firstcall
75      data    firstcall/.true./
76      save    firstcall
77     
78c  Initialisations
79c-----------------
80
81      if (firstcall) then
82
83        PRINT*,"******* ATTENTION, NEWTONIAN COOLING ********"
84
85        pressCLee = etaCLee * 9.2e6
86
87        DO i = 1, klon
88       
89          do k = 1,klev
90         
91            level = 1
92            do j=1,nlevCLee
93              if (pressCLee(j).gt.pplay(i,k)) level = j
94            enddo
95           
96            fact  = (log10(pplay(i,k))-log10(pressCLee(level)))
97     .        /(log10(pressCLee(level+1))-log10(pressCLee(level)))
98            ztemp = tempCLee(level)*(1-fact)+tempCLee(level+1)*fact
99            zdt   = dt_epCLee(level)*(1-fact)+dt_epCLee(level+1)*fact
100c           zt_eq(i,k) = ztemp + zdt*(cos(latitude(i))-2./RPI)
101            zt_eq(i,k) = ztemp + zdt*(cos(latitude(i))-RPI/4.)
102           
103          enddo
104         
105        ENDDO !i
106
107      endif ! firstcall
108     
109c+++++++ BOUCLE SUR LA GRILLE +++++++++++++++++++++++++
110      DO j = 1,klon
111 
112          do k = 1,klev
113             dTsdt(k) = -(t(j,k)-zt_eq(j,k))/tauCLee   ! en K/s
114          enddo
115       
116         radsol(j) = 0.           ! + vers bas
117         topsw(j) = 0.            ! + vers bas
118         toplw(j) = 0.            ! + vers haut
119         solsw(j) = 0.            ! + vers bas
120         sollw(j) = 0.            ! + vers bas
121         sollwdown(j) = 0.        ! + vers bas
122
123        DO k = 1, klev+1
124         lwnet  (j,k)   = 0.
125         swnet  (j,k)   = 0.
126        ENDDO
127
128        DO k = 1, klev
129         heat (j,k) = dTsdt(k)    ! K/s
130         cool (j,k) = 0.
131        ENDDO
132c
133      ENDDO !j
134c+++++++ FIN BOUCLE SUR LA GRILLE +++++++++++++++++++++++++
135
136      firstcall = .false.
137      RETURN
138      END
139
Note: See TracBrowser for help on using the repository browser.