Changeset 2486 for trunk/LMDZ.VENUS/libf
- Timestamp:
- Mar 26, 2021, 9:31:26 AM (4 years ago)
- Location:
- trunk/LMDZ.VENUS/libf/phyvenus
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.VENUS/libf/phyvenus/physiq_mod.F
r2464 r2486 71 71 USE param_v4_h 72 72 USE compo_hedin83_mod2 73 use radlwsw_newtoncool_mod, only: radlwsw_newtoncool 73 74 ! use ieee_arithmetic 74 75 use time_phylmdz_mod, only: annee_ref, day_ref, itau_phy … … 1501 1502 c ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 1502 1503 if (physideal) then 1503 CALL radlwsw_newtoncool 1504 e (dist, rmu0, fract, zzlev, 1505 e paprs, pplay,ftsol, t_seri) 1504 CALL radlwsw_newtoncool(presnivs,t_seri) 1506 1505 else 1507 1506 CALL radlwsw -
trunk/LMDZ.VENUS/libf/phyvenus/radlwsw_NewtonCool.F
r2135 r2486 2 2 ! $Header: /home/cvsroot/LMDZ4/libf/phylmd/radlwsw.F,v 1.2 2004/10/27 10:14:46 lmdzadmin Exp $ 3 3 ! 4 SUBROUTINE radlwsw_newtoncool(dist, rmu0, fract, zzlev, 5 . paprs, pplay,tsol, pt) 4 MODULE radlwsw_newtoncool_mod 5 6 implicit none 7 8 contains 9 10 SUBROUTINE radlwsw_newtoncool(presnivs,pt) 6 11 7 12 c====================================================================== 8 c Auteur(s): Z.X. Li (LMD/CNRS) date: 199607199 c Objet: interface entre le modele et les rayonnements10 c Arguments:11 c dist-----input-R- distance astronomique terre-soleil12 c rmu0-----input-R- cosinus de l'angle zenithal13 c fract----input-R- duree d'ensoleillement normalisee14 c paprs----input-R- pression a inter-couche (Pa)15 c pplay----input-R- pression au milieu de couche (Pa)16 c tsol-----input-R- temperature du sol (en K)17 c pt-------input-R- temperature (K)18 c19 20 13 c S. Lebonnois 12/04/2007 21 14 c VERSION NEWTONIAN COOLING pour Venus (no diurnal cycle) 22 15 c update 01/2014 23 24 16 c====================================================================== 25 use dimphy 17 use dimphy, only: klon,klev 26 18 USE geometry_mod, ONLY: latitude ! in radians 27 19 USE phys_state_var_mod, only: heat,cool,radsol, 28 20 . topsw,toplw,solsw,sollw,sollwdown,lwnet,swnet,zt_eq 29 USE write_field_phy 21 30 22 IMPLICIT none 31 #include "YOMCST.h" 32 #include "clesphys.h" 23 include "YOMCST.h" 33 24 34 25 c ARGUMENTS 35 real rmu0(klon), fract(klon), dist36 26 37 real zzlev(klon,klev+1),paprs(klon,klev+1), pplay(klon,klev) 38 real tsol(klon) 39 real pt(klon,klev) 27 real,intent(in) :: presnivs(klev) ! approx. pressure of GCM levels (Pa) 28 real,intent(in) :: pt(klon,klev) ! atmospheric temperature (K) 40 29 41 30 c LOCAL VARIABLES … … 71 60 c 72 61 73 logical firstcall 74 data firstcall/.true./ 75 save firstcall 62 logical,save :: firstcall=.true. 76 63 77 64 c Initialisations … … 79 66 80 67 if (firstcall) then 81 68 ! build zt_eq(), reference temperature field towards which to relax. 82 69 PRINT*,"******* ATTENTION, NEWTONIAN COOLING ********" 83 70 … … 90 77 level = 1 91 78 do j=1,nlevCLee 92 if (pressCLee(j).gt.p play(i,k)) level = j79 if (pressCLee(j).gt.presnivs(k)) level = j 93 80 enddo 94 81 95 fact = (log10(p play(i,k))-log10(pressCLee(level)))82 fact = (log10(presnivs(k))-log10(pressCLee(level))) 96 83 . /(log10(pressCLee(level+1))-log10(pressCLee(level))) 97 84 ztemp = tempCLee(level)*(1-fact)+tempCLee(level+1)*fact … … 104 91 ENDDO !i 105 92 93 firstcall = .false. 106 94 endif ! firstcall 107 95 … … 133 121 c+++++++ FIN BOUCLE SUR LA GRILLE +++++++++++++++++++++++++ 134 122 135 firstcall = .false. 136 RETURN 137 END 123 END SUBROUTINE radlwsw_newtoncool 138 124 125 END MODULE radlwsw_newtoncool_mod
Note: See TracChangeset
for help on using the changeset viewer.