Index: /trunk/DOC/000-USERS =================================================================== --- /trunk/DOC/000-USERS (revision 494) +++ /trunk/DOC/000-USERS (revision 495) @@ -28,2 +28,6 @@ Jeremie Burgalat GSMA [URCA] 01 - 06 - 2011 Titan GCM Pascal Rannou GSMA [URCA] 01 - 06 - 2011 Titan GCM +Helen Parish UCLA 01 - 10 - 2011 Venus GCM +David Luz CAAUL [Obs. Lisbonne] 01 - 01 - 2012 Venus GCM +Sarah Camberlain CAAUL [Obs. Lisbonne] 01 - 01 - 2012 Venus GCM + Index: /trunk/LMDZ.COMMON/libf/dyn3d/conf_gcm.F =================================================================== --- /trunk/LMDZ.COMMON/libf/dyn3d/conf_gcm.F (revision 494) +++ /trunk/LMDZ.COMMON/libf/dyn3d/conf_gcm.F (revision 495) @@ -338,4 +338,8 @@ CALL getin('mode_top_bound',mode_top_bound) CALL getin('tau_top_bound',tau_top_bound) + +! FOR TITAN: tidal forces + tidal=.TRUE. + CALL getin('tidal',tidal) !Config Key = coefdis Index: /trunk/LMDZ.COMMON/libf/dyn3d/leapfrog.F =================================================================== --- /trunk/LMDZ.COMMON/libf/dyn3d/leapfrog.F (revision 494) +++ /trunk/LMDZ.COMMON/libf/dyn3d/leapfrog.F (revision 495) @@ -106,4 +106,7 @@ REAL dtetatop(ip1jmp1,llm) REAL dqtop(ip1jmp1,llm,nqtot),dptop(ip1jmp1) + +c TITAN : tendances due au forces de marees */s + REAL dvtidal(ip1jm,llm),dutidal(ip1jmp1,llm) c tendances physiques */s @@ -396,4 +399,13 @@ $ finvmaold ) + IF ((planet_type.eq.titan).and.(tidal)) then +c----------------------------------------------------------------------- +c Marées gravitationnelles causées par Saturne +c B. Charnay (28/10/2010) +c ---------------------------------------------------------- + CALL tidal_forces(rdaym_ini, dutidal, dvtidal) + ucov=ucov+dutidal*dt + vcov=vcov+dvtidal*dt + ENDIF c .P.Le Van (26/04/94 ajout de finvpold dans l'appel d'integrd) Index: /trunk/LMDZ.COMMON/libf/dyn3d/logic.h =================================================================== --- /trunk/LMDZ.COMMON/libf/dyn3d/logic.h (revision 494) +++ /trunk/LMDZ.COMMON/libf/dyn3d/logic.h (revision 495) @@ -10,5 +10,5 @@ COMMON/logicl/ purmats,forward,leapf,apphys, & & statcl,conser,apdiss,apdelq,saison,ecripar,fxyhypb,ysinus & - & ,read_start,ok_guide,ok_strato,ok_gradsfile & + & ,read_start,ok_guide,ok_strato,tidal,ok_gradsfile & & ,ok_limit,ok_etat0,grilles_gcm_netcdf,hybrid @@ -17,5 +17,5 @@ LOGICAL purmats,forward,leapf,apphys,statcl,conser, & & apdiss,apdelq,saison,ecripar,fxyhypb,ysinus & - & ,read_start,ok_guide,ok_strato,ok_gradsfile & + & ,read_start,ok_guide,ok_strato,tidal,ok_gradsfile & & ,ok_limit,ok_etat0,grilles_gcm_netcdf logical hybrid ! vertical coordinate is hybrid if true (sigma otherwise) Index: /trunk/LMDZ.COMMON/libf/dyn3d/tidal_forces.F =================================================================== --- /trunk/LMDZ.COMMON/libf/dyn3d/tidal_forces.F (revision 495) +++ /trunk/LMDZ.COMMON/libf/dyn3d/tidal_forces.F (revision 495) @@ -0,0 +1,109 @@ + SUBROUTINE tidal_forces (t, du, dv) + + IMPLICIT NONE +c +c======================================================================= +c +c Auteur: B. Charnay (10/2010) +c ------- +c +c Objet: +c ------ +c +c ***************************************************************** +c ..... calcul du gradient horizontal du potentiel gravitationnel du aux forces de marees causees par Saturne +c ..... Formule tiree de Tokano 2002 +c ***************************************************************** +c Ces termes sont ajoutes a d(ucov)/dt et a d(vcov)/dt .. +c +c +c du et dv sont des arguments de sortie pour le s-pg .... +c +c======================================================================= +c +#include "dimensions.h" +#include "paramet.h" +#include "logic.h" +#include "comvert.h" +#include "comconst.h" +#include "comgeom.h" +!#include "comorbit.h" + REAL t ! jour de l'annee + REAL du( ip1jmp1,llm ), dv( ip1jm,llm ) + +c variables locales + REAL Vo + PARAMETER (Vo=-4.691e-6) + INTEGER l,ij,i,k + REAL n ! 2pi/periode de rotation siderale (en jours) + REAL a0 ! angle à l'instant initial entre Titan et le perihelie + PARAMETER (a0=0.) + +c cos et sin de la latitude et longitude, calcules au premiers appel + REAL coslonv(ip1jm),sinlonv(ip1jm) + REAL sinlatv(ip1jm),coslatv(ip1jm) + REAL coslonu(ip1jmp1),sinlonu(ip1jmp1) + REAL sinlatu(ip1jmp1),coslatu(ip1jmp1) + + LOGICAl first + + SAVE coslonv,coslonu,sinlonu,sinlonv + SAVE coslatv,coslatu,sinlatu,sinlatv + SAVE first, n + + DATA first /.true./ + +! Calcul des sin et cos aux points consideres + + IF(first) THEN + first=.false. + n=2*3.145!*(1+1/673.) + do i=1,iip1 + do k=1,jjm + coslonv(i+(k-1)*iip1)=cos(rlonv(i)) + sinlonv(i+(k-1)*iip1)=sin(rlonv(i)) + coslatv(i+(k-1)*iip1)=cos(rlatv(k)) + sinlatv(i+(k-1)*iip1)=sin(rlatv(k)) + ENDDO + ENDDO + + + + do i=1,iip1 + do k=1,jjp1 + coslonu(i+(k-1)*iip1)=cos(rlonu(i)) + sinlonu(i+(k-1)*iip1)=sin(rlonu(i)) + coslatu(i+(k-1)*iip1)=cos(rlatu(k)) + sinlatu(i+(k-1)*iip1)=sin(rlatu(k)) + ENDDO + ENDDO + + + + ENDIF + + +! Tendance du aux forces de maree + + DO l = 1,llm + + DO ij = 1, ip1jmp1 + + du(ij,l) = cu(ij)*Vo + $ *(3*sinlonu(ij)*coslonu(ij)*coslatu(ij)*cos(n*t+a0) + $ -2*coslatu(ij)*(2*coslonu(ij)**2-1)*sin(n*t+a0)) + ENDDO + + DO ij = 1, ip1jm + dv(ij,l) = cv(ij)*Vo + $ *(3*sinlatv(ij)*coslatv(ij)*coslonv(ij)**2*cos(n*t+a0) + $ + 4*coslatv(ij)*sinlatv(ij)*sinlonv(ij)*coslonv(ij) + $ *sin(n*t+a0)) + ENDDO + + ENDDO + + +c + RETURN + END Index: /trunk/LMDZ.COMMON/libf/dyn3dpar/conf_gcm.F =================================================================== --- /trunk/LMDZ.COMMON/libf/dyn3dpar/conf_gcm.F (revision 494) +++ /trunk/LMDZ.COMMON/libf/dyn3dpar/conf_gcm.F (revision 495) @@ -352,5 +352,8 @@ CALL getin('tau_top_bound',tau_top_bound) -! +! FOR TITAN: tidal forces + tidal=.TRUE. + CALL getin('tidal',tidal) + !Config Key = coefdis !Config Desc = coefficient pour gamdissip Index: /trunk/LMDZ.COMMON/libf/dyn3dpar/leapfrog_p.F =================================================================== --- /trunk/LMDZ.COMMON/libf/dyn3dpar/leapfrog_p.F (revision 494) +++ /trunk/LMDZ.COMMON/libf/dyn3dpar/leapfrog_p.F (revision 495) @@ -117,4 +117,7 @@ REAL,SAVE :: dptop(ip1jmp1) REAL,DIMENSION(:,:,:),ALLOCATABLE,SAVE :: dqtop + +c TITAN : tendances due au forces de marees */s + REAL,SAVE :: dvtidal(ip1jm,llm),dutidal(ip1jmp1,llm) c variables pour le fichier histoire @@ -696,4 +699,14 @@ $ finvmaold ) + IF ((planet_type.eq.titan).and.(tidal)) then +c----------------------------------------------------------------------- +c Marées gravitationnelles causées par Saturne +c B. Charnay (28/10/2010) +c ---------------------------------------------------------- + CALL tidal_forces(rdaym_ini, dutidal, dvtidal) + ucov=ucov+dutidal*dt + vcov=vcov+dvtidal*dt + ENDIF + ! CALL FTRACE_REGION_END("integrd") c$OMP BARRIER Index: /trunk/LMDZ.COMMON/libf/dyn3dpar/logic.h =================================================================== --- /trunk/LMDZ.COMMON/libf/dyn3dpar/logic.h (revision 494) +++ /trunk/LMDZ.COMMON/libf/dyn3dpar/logic.h (revision 495) @@ -10,5 +10,5 @@ COMMON/logicl/ purmats,forward,leapf,apphys, & & statcl,conser,apdiss,apdelq,saison,ecripar,fxyhypb,ysinus & - & ,read_start,ok_guide,ok_strato,ok_gradsfile & + & ,read_start,ok_guide,ok_strato,tidal,ok_gradsfile & & ,ok_limit,ok_etat0,grilles_gcm_netcdf,hybrid @@ -17,5 +17,5 @@ LOGICAL purmats,forward,leapf,apphys,statcl,conser, & & apdiss,apdelq,saison,ecripar,fxyhypb,ysinus & - & ,read_start,ok_guide,ok_strato,ok_gradsfile & + & ,read_start,ok_guide,ok_strato,tidal,ok_gradsfile & & ,ok_limit,ok_etat0,grilles_gcm_netcdf logical hybrid ! vertical coordinate is hybrid if true (sigma otherwise) Index: /trunk/LMDZ.COMMON/libf/dyn3dpar/tidal_forces.F =================================================================== --- /trunk/LMDZ.COMMON/libf/dyn3dpar/tidal_forces.F (revision 495) +++ /trunk/LMDZ.COMMON/libf/dyn3dpar/tidal_forces.F (revision 495) @@ -0,0 +1,109 @@ + SUBROUTINE tidal_forces (t, du, dv) + + IMPLICIT NONE +c +c======================================================================= +c +c Auteur: B. Charnay (10/2010) +c ------- +c +c Objet: +c ------ +c +c ***************************************************************** +c ..... calcul du gradient horizontal du potentiel gravitationnel du aux forces de marees causees par Saturne +c ..... Formule tiree de Tokano 2002 +c ***************************************************************** +c Ces termes sont ajoutes a d(ucov)/dt et a d(vcov)/dt .. +c +c +c du et dv sont des arguments de sortie pour le s-pg .... +c +c======================================================================= +c +#include "dimensions.h" +#include "paramet.h" +#include "logic.h" +#include "comvert.h" +#include "comconst.h" +#include "comgeom.h" +!#include "comorbit.h" + REAL t ! jour de l'annee + REAL du( ip1jmp1,llm ), dv( ip1jm,llm ) + +c variables locales + REAL Vo + PARAMETER (Vo=-4.691e-6) + INTEGER l,ij,i,k + REAL n ! 2pi/periode de rotation siderale (en jours) + REAL a0 ! angle à l'instant initial entre Titan et le perihelie + PARAMETER (a0=0.) + +c cos et sin de la latitude et longitude, calcules au premiers appel + REAL coslonv(ip1jm),sinlonv(ip1jm) + REAL sinlatv(ip1jm),coslatv(ip1jm) + REAL coslonu(ip1jmp1),sinlonu(ip1jmp1) + REAL sinlatu(ip1jmp1),coslatu(ip1jmp1) + + LOGICAl first + + SAVE coslonv,coslonu,sinlonu,sinlonv + SAVE coslatv,coslatu,sinlatu,sinlatv + SAVE first, n + + DATA first /.true./ + +! Calcul des sin et cos aux points consideres + + IF(first) THEN + first=.false. + n=2*3.145!*(1+1/673.) + do i=1,iip1 + do k=1,jjm + coslonv(i+(k-1)*iip1)=cos(rlonv(i)) + sinlonv(i+(k-1)*iip1)=sin(rlonv(i)) + coslatv(i+(k-1)*iip1)=cos(rlatv(k)) + sinlatv(i+(k-1)*iip1)=sin(rlatv(k)) + ENDDO + ENDDO + + + + do i=1,iip1 + do k=1,jjp1 + coslonu(i+(k-1)*iip1)=cos(rlonu(i)) + sinlonu(i+(k-1)*iip1)=sin(rlonu(i)) + coslatu(i+(k-1)*iip1)=cos(rlatu(k)) + sinlatu(i+(k-1)*iip1)=sin(rlatu(k)) + ENDDO + ENDDO + + + + ENDIF + + +! Tendance du aux forces de maree + + DO l = 1,llm + + DO ij = 1, ip1jmp1 + + du(ij,l) = cu(ij)*Vo + $ *(3*sinlonu(ij)*coslonu(ij)*coslatu(ij)*cos(n*t+a0) + $ -2*coslatu(ij)*(2*coslonu(ij)**2-1)*sin(n*t+a0)) + ENDDO + + DO ij = 1, ip1jm + dv(ij,l) = cv(ij)*Vo + $ *(3*sinlatv(ij)*coslatv(ij)*coslonv(ij)**2*cos(n*t+a0) + $ + 4*coslatv(ij)*sinlatv(ij)*sinlonv(ij)*coslonv(ij) + $ *sin(n*t+a0)) + ENDDO + + ENDDO + + +c + RETURN + END Index: /trunk/LMDZ.TITAN/deftank/gcm.def =================================================================== --- /trunk/LMDZ.TITAN/deftank/gcm.def (revision 494) +++ /trunk/LMDZ.TITAN/deftank/gcm.def (revision 495) @@ -63,5 +63,7 @@ # Coefficient pour la couche eponge (valeur derniere couche) tau_top_bound=4.e-5 -# +## Maree gravitationnelle ou non +tidal=y + ## longitude en degres du centre du zoom clon=0. Index: /trunk/LMDZ.TITAN/deftank/physiq.def =================================================================== --- /trunk/LMDZ.TITAN/deftank/physiq.def (revision 494) +++ /trunk/LMDZ.TITAN/deftank/physiq.def (revision 495) @@ -68,5 +68,6 @@ ok_kzmin=n # -inertie=340. +inertie=2000. +emis=0.95 # # parametres convection seche Index: /trunk/LMDZ.TITAN/libf/phytitan/cld.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/cld.F (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/cld.F (revision 495) @@ -11,5 +11,5 @@ COMMON /CLOUD/ RADCLD(NLAYER), XNCLD(NLAYER) & , RCLDI(NSPECI), XICLDI(NSPECI), RCLDV(NSPECV), XICLDV(NSPECV) - COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI + COMMON /PLANT/ CSUBP,F0PI COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON COMMON /CONST/RGAS,RHOP,PI,SIGMA Index: /trunk/LMDZ.TITAN/libf/phytitan/clesphys.h =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/clesphys.h (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/clesphys.h (revision 495) @@ -20,5 +20,5 @@ ! Parametres surface: - REAL inertie + REAL inertie,emis ! Parametres Chimie: Index: /trunk/LMDZ.TITAN/libf/phytitan/clmain.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/clmain.F (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/clmain.F (revision 495) @@ -811,5 +811,5 @@ zx_alf2(i) = 1.0 - zx_alf1(i) zx_coef(i,1) = coef(i,1) - . * (1.0+SQRT(u1lay(i)**2+v1lay(i)**2)) + . * SQRT(u1lay(i)**2+v1lay(i)**2) . * pplay(i,1)/(RD*t(i,1)) zx_coef(i,1) = zx_coef(i,1) * dtime*RG Index: /trunk/LMDZ.TITAN/libf/phytitan/conf_phys.F90 =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/conf_phys.F90 (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/conf_phys.F90 (revision 495) @@ -272,9 +272,17 @@ !Config Key = inertie !Config Desc = -!Config Def = 340. -!Config Help = -! - inertie = 340. +!Config Def = 2000. +!Config Help = +! + inertie = 2000. call getin('inertie',inertie) +! +!Config Key = emis +!Config Desc = +!Config Def = 0.95 +!Config Help = +! + emis = 0.95 + call getin('emis',emis) ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! @@ -467,4 +475,5 @@ write(numout,*)' ok_kzmin = ',ok_kzmin write(numout,*)' inertie = ', inertie + write(numout,*)' emis = ', emis write(numout,*)' iflag_ajs = ', iflag_ajs write(numout,*)' chimi = ', chimi Index: /trunk/LMDZ.TITAN/libf/phytitan/cooling.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/cooling.F (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/cooling.F (revision 495) @@ -46,5 +46,5 @@ c spectral ranges. This common is initialized by radtitan. c -c COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI +c COMMON /PLANT/ CSUBP,F0PI c This common is initialized by tgmdat. c @@ -84,6 +84,6 @@ COMMON /IRTAUS/ dtaui,dtauip - COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI - REAL CSUBP,RSFI,RSFV,F0PI + COMMON /PLANT/ CSUBP,F0PI + REAL CSUBP,F0PI c Local: @@ -237,7 +237,7 @@ ENDDO c BSURF=PLNCK( WAVEN, TEMP(ig,NL), DW) - BSURF=.2*PLNCK*(1-RSFI) - FUPI(ig,NL)=BSURF * 2.*RPI*UBARI + RSFI*FDI(ig,NL) - FUPIS(ig,NL,K)=BSURF*2.*RPI*UBARI+RSFI*FDIS(ig,NL,K) + BSURF=.2*PLNCK*emis + FUPI(ig,NL) =BSURF*2.*RPI*UBARI+(1-emis)*FDI(ig,NL) + FUPIS(ig,NL,K)=BSURF*2.*RPI*UBARI+(1-emis)*FDIS(ig,NL,K) 2310 CONTINUE c write(*,*) Index: /trunk/LMDZ.TITAN/libf/phytitan/heating.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/heating.F (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/heating.F (revision 495) @@ -1,3 +1,3 @@ - SUBROUTINE heating(dist,rmu0,fract,sol_htg,swnet,icld) + SUBROUTINE heating(dist,rmu0,fract,falbe,sol_htg,swnet,icld) @@ -16,4 +16,5 @@ c rmu0-----input-R- cosinus de l'angle zenithal c fract----input-R- duree d'ensoleillement normalisee +c falbe----input-R- surface albedo c icld-----input-I- calcul avec nuages. c p(klon,nl) pressure (level) @@ -46,5 +47,5 @@ - real dist, rmu0(klon), fract(klon) + real dist, rmu0(klon), fract(klon), falbe(klon) integer icld @@ -65,5 +66,5 @@ REAL FNETV(ngrid,NLEVEL),FUPV(ngrid,NLEVEL,NSPECV) REAL FDV(ngrid,NLEVEL,NSPECV),FMNETV(ngrid,NLEVEL) - REAL CSUBP,RSFI,RSFV,F0PI + REAL CSUBP,F0PI COMMON /UBARED/ UBARI,UBARV,UBAR0 @@ -79,5 +80,5 @@ & FMNETV - COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI + COMMON /PLANT/ CSUBP,F0PI @@ -98,5 +99,5 @@ ubar0=rmu0(ig) - CALL sfluxv(iprint,ig,dist,icld) ! #3 + CALL sfluxv(iprint,ig,dist,falbe,icld) ! #3 fnetv(ig,:) = fnetv(ig,:) *fract(ig) ! >0 vers le haut Index: /trunk/LMDZ.TITAN/libf/phytitan/interface_surf.F90 =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/interface_surf.F90 (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/interface_surf.F90 (revision 495) @@ -256,5 +256,5 @@ zx_coef(i) = coef1lay(i) & - & * (1.0+SQRT(u1lay(i)**2+v1lay(i)**2)) & + & * SQRT(u1lay(i)**2+v1lay(i)**2) & & * p1lay(i)/(RD*t1lay(i)) Index: /trunk/LMDZ.TITAN/libf/phytitan/optci.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/optci.F (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/optci.F (revision 495) @@ -66,5 +66,5 @@ COMMON /IRTAUS/ DTAUP,DTAUPP - COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI + COMMON /PLANT/ CSUBP,F0PI COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON COMMON /CONST/RGAS,RHOP,PI,SIGMA Index: /trunk/LMDZ.TITAN/libf/phytitan/optci_1pt.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/optci_1pt.F (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/optci_1pt.F (revision 495) @@ -38,5 +38,5 @@ & DWNI(NSPECI), WLNI(NSPECI) - COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI + COMMON /PLANT/ CSUBP,F0PI COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON COMMON /CONST/RGAS,RHOP,PI,SIGMA Index: /trunk/LMDZ.TITAN/libf/phytitan/optci_1pt_2.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/optci_1pt_2.F (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/optci_1pt_2.F (revision 495) @@ -39,5 +39,5 @@ & DWNI(NSPECI), WLNI(NSPECI) - COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI + COMMON /PLANT/ CSUBP,F0PI COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON COMMON /CONST/RGAS,RHOP,PI,SIGMA Index: /trunk/LMDZ.TITAN/libf/phytitan/optcv.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/optcv.F (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/optcv.F (revision 495) @@ -61,5 +61,5 @@ & ,DWNV(NSPECV),WLNV(NSPECV) - COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI + COMMON /PLANT/ CSUBP,F0PI COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON COMMON /CONST/ RGAS,RHOP,PI,SIGMA Index: /trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt.F (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt.F (revision 495) @@ -39,5 +39,5 @@ & ,DWNV(NSPECV),WLNV(NSPECV) - COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI + COMMON /PLANT/ CSUBP,F0PI COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON COMMON /CONST/ RGAS,RHOP,PI,SIGMA Index: /trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt_2.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt_2.F (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt_2.F (revision 495) @@ -38,5 +38,5 @@ & ,DWNV(NSPECV),WLNV(NSPECV) - COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI + COMMON /PLANT/ CSUBP,F0PI COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON COMMON /CONST/ RGAS,RHOP,PI,SIGMA Index: /trunk/LMDZ.TITAN/libf/phytitan/physiq.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/physiq.F (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/physiq.F (revision 495) @@ -413,6 +413,6 @@ COMMON /TAUD/ TAUHID,TAUCID,TAUGID,TAUHVD,TAUCVD,TAUGVD - COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI - REAL CSUBP,RSFI,RSFV,F0PI + COMMON /PLANT/ CSUBP,F0PI + REAL CSUBP,F0PI * common relatifs au nuages @@ -995,5 +995,5 @@ c DO i = 1, klon - dlw(i) = - 4.0*(1-RSFI)*RSIGMA*ftsol(i)**3 + dlw(i) = - 4.0*emis*RSIGMA*ftsol(i)**3 ENDDO @@ -1212,5 +1212,5 @@ c call begintime(tt0) CALL radlwsw - e (dist, rmu0, fract, dtimerad, zzlev, + e (dist, rmu0, fract, falbe, dtimerad, zzlev, e paprs, pplay,ftsol, t_seri, nqmax, nmicro, c tr_seri, qaer, Index: /trunk/LMDZ.TITAN/libf/phytitan/radlwsw.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/radlwsw.F (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/radlwsw.F (revision 495) @@ -1,3 +1,3 @@ - SUBROUTINE radlwsw(dist, rmu0, fract, dtimerad, zzlev, + SUBROUTINE radlwsw(dist, rmu0, fract, falbe, dtimerad, zzlev, . paprs, pplay,tsol, pt, nq, nmicro, pq, . qaer, @@ -14,4 +14,5 @@ c rmu0-----input-R- cosinus de l'angle zenithal c fract----input-R- duree d'ensoleillement normalisee +c falbe----input-R- surface albedo c dtimerad-input-R- intervalle de temps du radiatif c zzlev----input-R- altitudes des inter-couches (m) @@ -48,5 +49,5 @@ c ARGUMENTS INTEGER nq,nmicro - real rmu0(klon), fract(klon), dist, dtimerad + real rmu0(klon), fract(klon), falbe(klon), dist, dtimerad c real zzlev(klon,klev+1),paprs(klon,klev+1), pplay(klon,klev) @@ -129,5 +130,5 @@ IF (clouds.eq.1) THEN ICLD = 1 ! colonne avec nuages - CALL heating(dist,rmu0,fract,zheatc,zswnetc,icld) + CALL heating(dist,rmu0,fract,falbe,zheatc,zswnetc,icld) ELSE zheatc = 0. @@ -135,5 +136,5 @@ ENDIF ICLD = 0 ! colonne sans nuages - CALL heating(dist,rmu0,fract,zheatp,zswnetp,icld) + CALL heating(dist,rmu0,fract,falbe,zheatp,zswnetp,icld) c inversion de l'axe vertical Index: /trunk/LMDZ.TITAN/libf/phytitan/sfluxv.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/sfluxv.F (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/sfluxv.F (revision 495) @@ -1,3 +1,3 @@ - SUBROUTINE SFLUXV(IPRINT,IG,dist_sol,icld) + SUBROUTINE SFLUXV(IPRINT,IG,dist_sol,falbe,icld) use dimphy @@ -10,8 +10,11 @@ PARAMETER (ngrid=(jjm-1)*iim+2) ! = klon c - INTEGER NLAYER,NLEVEL,NSPECV,NSPC1V,icld + INTEGER IG,IPRINT,icld + real dist_sol,falbe(ngrid) + + INTEGER NLAYER,NLEVEL,NSPECV,NSPC1V PARAMETER (NLAYER=llm,NLEVEL=NLAYER+1) PARAMETER (NSPECV=24,NSPC1V=25) - INTEGER IG,IPRINT,NT,NTERM(NSPECV),J,K + INTEGER NT,NTERM(NSPECV),J,K REAL FUW(NLEVEL),FDW(NLEVEL) @@ -35,8 +38,6 @@ & FMNETV(ngrid,NLEVEL), & FMUPV(NLEVEL),FMDV(NLEVEL) - REAL CSUBP,RSFI,RSFV,F0PI + REAL CSUBP,F0PI REAL UBARI,UBARV,UBAR0 - - real dist_sol COMMON /VISGAS/SOLARF,NTERM,PEXPON, @@ -60,5 +61,5 @@ & FMNETV - COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI + COMMON /PLANT/ CSUBP,F0PI COMMON /UBARED/ UBARI,UBARV,UBAR0 @@ -95,7 +96,7 @@ DO 912 NT=1,NTERM(K) IF (ICLD.eq.1) THEN - BSURF=0.+ RSFV*UBAR0*F0PI*EXP(-TAUV(ig,NLEVEL,K,NT)/UBAR0) + BSURF=0.+ falbe(ig)*UBAR0*F0PI*EXP(-TAUV(ig,NLEVEL,K,NT)/UBAR0) ELSE - BSURF=0.+ RSFV*UBAR0*F0PI*EXP(-TAUVP(ig,NLEVEL,K,NT)/UBAR0) + BSURF=0.+ falbe(ig)*UBAR0*F0PI*EXP(-TAUVP(ig,NLEVEL,K,NT)/UBAR0) ENDIF C @@ -139,5 +140,5 @@ c write(*,*) (CO0(J),J=1,NLAYER) c print*,'UBAR0 ',UBAR0 -c print*,NLEVEL,WNOV(K),F0PI,RSFV,BTOP,BSURF +c print*,NLEVEL,WNOV(K),F0PI,falbe(ig),BTOP,BSURF FUW = 0.0 FDW = 0.0 @@ -146,5 +147,5 @@ CALL GFLUXV(NLEVEL,WNOV(K),DT0,T0, - & WB0,CO0,F0PI,RSFV,BTOP,BSURF,FUW,FDW,FMUPV, + & WB0,CO0,F0PI,falbe(ig),BTOP,BSURF,FUW,FDW,FMUPV, & FMDV,IPRINT) c PRINT*,'sortie gfluxv #: ',ig,K Index: /trunk/LMDZ.TITAN/libf/phytitan/tgmdat.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/tgmdat.F (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/tgmdat.F (revision 495) @@ -4,5 +4,5 @@ COMMON /UBARED/ UBARI,UBARV,UBAR0 COMMON /LAPSE/ DTDP(NLAYER),CONVEQ - COMMON /PLANT/ CSUBP,RSFI,RSFV,F0PI + COMMON /PLANT/ CSUBP,F0PI COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON COMMON /CONST/RGAS,RHOP,PI,SIGMA @@ -29,6 +29,4 @@ C CONSTANT MUST BE CHANGED DATA RHOP/1.E4/ -C RSF IS THE SURFACE REFLECTANCE FOR VIS AND IR - DATA RSFV,RSFI/0.15,0.05/ C FOPI IS THE ACTUAL SOLAR FLUX IN ERGS/CM2 DATA F0PI/1.5E4/ Index: /trunk/LMDZ.TITAN/libf/phytitan/ustarhb.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/ustarhb.F (revision 494) +++ /trunk/LMDZ.TITAN/libf/phytitan/ustarhb.F (revision 495) @@ -41,5 +41,5 @@ zxu = u(i,1)*zx_alf1+u(i,2)*zx_alf2 zxv = v(i,1)*zx_alf1+v(i,2)*zx_alf2 - zxmod = 1.0+SQRT(zxu**2+zxv**2) + zxmod = SQRT(zxu**2+zxv**2) taux = zxu *zxmod*cd_m(i) tauy = zxv *zxmod*cd_m(i) Index: /trunk/LMDZ.VENUS/deftank/gcm.def =================================================================== --- /trunk/LMDZ.VENUS/deftank/gcm.def (revision 494) +++ /trunk/LMDZ.VENUS/deftank/gcm.def (revision 495) @@ -45,5 +45,5 @@ ## Dissipation horizontale dissip_fac_mid=2. -dissip_fac_up=2. +dissip_fac_up=10. # deltaz et hdelta en km dissip_deltaz=10.