Changeset 495 for trunk/LMDZ.TITAN/libf
- Timestamp:
- Jan 10, 2012, 5:59:10 PM (13 years ago)
- Location:
- trunk/LMDZ.TITAN/libf/phytitan
- Files:
-
- 18 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.TITAN/libf/phytitan/cld.F
r3 r495 11 11 COMMON /CLOUD/ RADCLD(NLAYER), XNCLD(NLAYER) 12 12 & , RCLDI(NSPECI), XICLDI(NSPECI), RCLDV(NSPECV), XICLDV(NSPECV) 13 COMMON /PLANT/ CSUBP, RSFI,RSFV,F0PI13 COMMON /PLANT/ CSUBP,F0PI 14 14 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON 15 15 COMMON /CONST/RGAS,RHOP,PI,SIGMA -
trunk/LMDZ.TITAN/libf/phytitan/clesphys.h
r175 r495 20 20 21 21 ! Parametres surface: 22 REAL inertie 22 REAL inertie,emis 23 23 24 24 ! Parametres Chimie: -
trunk/LMDZ.TITAN/libf/phytitan/clmain.F
r102 r495 811 811 zx_alf2(i) = 1.0 - zx_alf1(i) 812 812 zx_coef(i,1) = coef(i,1) 813 . * (1.0+SQRT(u1lay(i)**2+v1lay(i)**2))813 . * SQRT(u1lay(i)**2+v1lay(i)**2) 814 814 . * pplay(i,1)/(RD*t(i,1)) 815 815 zx_coef(i,1) = zx_coef(i,1) * dtime*RG -
trunk/LMDZ.TITAN/libf/phytitan/conf_phys.F90
r474 r495 272 272 !Config Key = inertie 273 273 !Config Desc = 274 !Config Def = 340.275 !Config Help = 276 ! 277 inertie = 340.274 !Config Def = 2000. 275 !Config Help = 276 ! 277 inertie = 2000. 278 278 call getin('inertie',inertie) 279 ! 280 !Config Key = emis 281 !Config Desc = 282 !Config Def = 0.95 283 !Config Help = 284 ! 285 emis = 0.95 286 call getin('emis',emis) 279 287 ! 280 288 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! … … 467 475 write(numout,*)' ok_kzmin = ',ok_kzmin 468 476 write(numout,*)' inertie = ', inertie 477 write(numout,*)' emis = ', emis 469 478 write(numout,*)' iflag_ajs = ', iflag_ajs 470 479 write(numout,*)' chimi = ', chimi -
trunk/LMDZ.TITAN/libf/phytitan/cooling.F
r306 r495 46 46 c spectral ranges. This common is initialized by radtitan. 47 47 c 48 c COMMON /PLANT/ CSUBP, RSFI,RSFV,F0PI48 c COMMON /PLANT/ CSUBP,F0PI 49 49 c This common is initialized by tgmdat. 50 50 c … … 84 84 COMMON /IRTAUS/ dtaui,dtauip 85 85 86 COMMON /PLANT/ CSUBP, RSFI,RSFV,F0PI87 REAL CSUBP, RSFI,RSFV,F0PI86 COMMON /PLANT/ CSUBP,F0PI 87 REAL CSUBP,F0PI 88 88 89 89 c Local: … … 237 237 ENDDO 238 238 c BSURF=PLNCK( WAVEN, TEMP(ig,NL), DW) 239 BSURF=.2*PLNCK* (1-RSFI)240 FUPI(ig,NL) =BSURF * 2.*RPI*UBARI + RSFI*FDI(ig,NL)241 FUPIS(ig,NL,K)=BSURF*2.*RPI*UBARI+ RSFI*FDIS(ig,NL,K)239 BSURF=.2*PLNCK*emis 240 FUPI(ig,NL) =BSURF*2.*RPI*UBARI+(1-emis)*FDI(ig,NL) 241 FUPIS(ig,NL,K)=BSURF*2.*RPI*UBARI+(1-emis)*FDIS(ig,NL,K) 242 242 2310 CONTINUE 243 243 c write(*,*) -
trunk/LMDZ.TITAN/libf/phytitan/heating.F
r175 r495 1 SUBROUTINE heating(dist,rmu0,fract, sol_htg,swnet,icld)1 SUBROUTINE heating(dist,rmu0,fract,falbe,sol_htg,swnet,icld) 2 2 3 3 … … 16 16 c rmu0-----input-R- cosinus de l'angle zenithal 17 17 c fract----input-R- duree d'ensoleillement normalisee 18 c falbe----input-R- surface albedo 18 19 c icld-----input-I- calcul avec nuages. 19 20 c p(klon,nl) pressure (level) … … 46 47 47 48 48 real dist, rmu0(klon), fract(klon) 49 real dist, rmu0(klon), fract(klon), falbe(klon) 49 50 integer icld 50 51 … … 65 66 REAL FNETV(ngrid,NLEVEL),FUPV(ngrid,NLEVEL,NSPECV) 66 67 REAL FDV(ngrid,NLEVEL,NSPECV),FMNETV(ngrid,NLEVEL) 67 REAL CSUBP, RSFI,RSFV,F0PI68 REAL CSUBP,F0PI 68 69 69 70 COMMON /UBARED/ UBARI,UBARV,UBAR0 … … 79 80 & FMNETV 80 81 81 COMMON /PLANT/ CSUBP, RSFI,RSFV,F0PI82 COMMON /PLANT/ CSUBP,F0PI 82 83 83 84 … … 98 99 ubar0=rmu0(ig) 99 100 100 CALL sfluxv(iprint,ig,dist, icld)! #3101 CALL sfluxv(iprint,ig,dist,falbe,icld) ! #3 101 102 102 103 fnetv(ig,:) = fnetv(ig,:) *fract(ig) ! >0 vers le haut -
trunk/LMDZ.TITAN/libf/phytitan/interface_surf.F90
r97 r495 256 256 257 257 zx_coef(i) = coef1lay(i) & 258 & * (1.0+SQRT(u1lay(i)**2+v1lay(i)**2)) &258 & * SQRT(u1lay(i)**2+v1lay(i)**2) & 259 259 & * p1lay(i)/(RD*t1lay(i)) 260 260 -
trunk/LMDZ.TITAN/libf/phytitan/optci.F
r306 r495 66 66 COMMON /IRTAUS/ DTAUP,DTAUPP 67 67 68 COMMON /PLANT/ CSUBP, RSFI,RSFV,F0PI68 COMMON /PLANT/ CSUBP,F0PI 69 69 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON 70 70 COMMON /CONST/RGAS,RHOP,PI,SIGMA -
trunk/LMDZ.TITAN/libf/phytitan/optci_1pt.F
r175 r495 38 38 & DWNI(NSPECI), WLNI(NSPECI) 39 39 40 COMMON /PLANT/ CSUBP, RSFI,RSFV,F0PI40 COMMON /PLANT/ CSUBP,F0PI 41 41 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON 42 42 COMMON /CONST/RGAS,RHOP,PI,SIGMA -
trunk/LMDZ.TITAN/libf/phytitan/optci_1pt_2.F
r306 r495 39 39 & DWNI(NSPECI), WLNI(NSPECI) 40 40 41 COMMON /PLANT/ CSUBP, RSFI,RSFV,F0PI41 COMMON /PLANT/ CSUBP,F0PI 42 42 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON 43 43 COMMON /CONST/RGAS,RHOP,PI,SIGMA -
trunk/LMDZ.TITAN/libf/phytitan/optcv.F
r175 r495 61 61 & ,DWNV(NSPECV),WLNV(NSPECV) 62 62 63 COMMON /PLANT/ CSUBP, RSFI,RSFV,F0PI63 COMMON /PLANT/ CSUBP,F0PI 64 64 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON 65 65 COMMON /CONST/ RGAS,RHOP,PI,SIGMA -
trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt.F
r175 r495 39 39 & ,DWNV(NSPECV),WLNV(NSPECV) 40 40 41 COMMON /PLANT/ CSUBP, RSFI,RSFV,F0PI41 COMMON /PLANT/ CSUBP,F0PI 42 42 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON 43 43 COMMON /CONST/ RGAS,RHOP,PI,SIGMA -
trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt_2.F
r306 r495 38 38 & ,DWNV(NSPECV),WLNV(NSPECV) 39 39 40 COMMON /PLANT/ CSUBP, RSFI,RSFV,F0PI40 COMMON /PLANT/ CSUBP,F0PI 41 41 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON 42 42 COMMON /CONST/ RGAS,RHOP,PI,SIGMA -
trunk/LMDZ.TITAN/libf/phytitan/physiq.F
r306 r495 413 413 414 414 COMMON /TAUD/ TAUHID,TAUCID,TAUGID,TAUHVD,TAUCVD,TAUGVD 415 COMMON /PLANT/ CSUBP, RSFI,RSFV,F0PI416 REAL CSUBP, RSFI,RSFV,F0PI415 COMMON /PLANT/ CSUBP,F0PI 416 REAL CSUBP,F0PI 417 417 418 418 * common relatifs au nuages … … 995 995 c 996 996 DO i = 1, klon 997 dlw(i) = - 4.0* (1-RSFI)*RSIGMA*ftsol(i)**3997 dlw(i) = - 4.0*emis*RSIGMA*ftsol(i)**3 998 998 ENDDO 999 999 … … 1212 1212 c call begintime(tt0) 1213 1213 CALL radlwsw 1214 e (dist, rmu0, fract, dtimerad, zzlev,1214 e (dist, rmu0, fract, falbe, dtimerad, zzlev, 1215 1215 e paprs, pplay,ftsol, t_seri, nqmax, nmicro, 1216 1216 c tr_seri, qaer, -
trunk/LMDZ.TITAN/libf/phytitan/radlwsw.F
r175 r495 1 SUBROUTINE radlwsw(dist, rmu0, fract, dtimerad, zzlev,1 SUBROUTINE radlwsw(dist, rmu0, fract, falbe, dtimerad, zzlev, 2 2 . paprs, pplay,tsol, pt, nq, nmicro, pq, 3 3 . qaer, … … 14 14 c rmu0-----input-R- cosinus de l'angle zenithal 15 15 c fract----input-R- duree d'ensoleillement normalisee 16 c falbe----input-R- surface albedo 16 17 c dtimerad-input-R- intervalle de temps du radiatif 17 18 c zzlev----input-R- altitudes des inter-couches (m) … … 48 49 c ARGUMENTS 49 50 INTEGER nq,nmicro 50 real rmu0(klon), fract(klon), dist, dtimerad51 real rmu0(klon), fract(klon), falbe(klon), dist, dtimerad 51 52 c 52 53 real zzlev(klon,klev+1),paprs(klon,klev+1), pplay(klon,klev) … … 129 130 IF (clouds.eq.1) THEN 130 131 ICLD = 1 ! colonne avec nuages 131 CALL heating(dist,rmu0,fract, zheatc,zswnetc,icld)132 CALL heating(dist,rmu0,fract,falbe,zheatc,zswnetc,icld) 132 133 ELSE 133 134 zheatc = 0. … … 135 136 ENDIF 136 137 ICLD = 0 ! colonne sans nuages 137 CALL heating(dist,rmu0,fract, zheatp,zswnetp,icld)138 CALL heating(dist,rmu0,fract,falbe,zheatp,zswnetp,icld) 138 139 139 140 c inversion de l'axe vertical -
trunk/LMDZ.TITAN/libf/phytitan/sfluxv.F
r175 r495 1 SUBROUTINE SFLUXV(IPRINT,IG,dist_sol, icld)1 SUBROUTINE SFLUXV(IPRINT,IG,dist_sol,falbe,icld) 2 2 3 3 use dimphy … … 10 10 PARAMETER (ngrid=(jjm-1)*iim+2) ! = klon 11 11 c 12 INTEGER NLAYER,NLEVEL,NSPECV,NSPC1V,icld 12 INTEGER IG,IPRINT,icld 13 real dist_sol,falbe(ngrid) 14 15 INTEGER NLAYER,NLEVEL,NSPECV,NSPC1V 13 16 PARAMETER (NLAYER=llm,NLEVEL=NLAYER+1) 14 17 PARAMETER (NSPECV=24,NSPC1V=25) 15 INTEGER IG,IPRINT,NT,NTERM(NSPECV),J,K18 INTEGER NT,NTERM(NSPECV),J,K 16 19 17 20 REAL FUW(NLEVEL),FDW(NLEVEL) … … 35 38 & FMNETV(ngrid,NLEVEL), 36 39 & FMUPV(NLEVEL),FMDV(NLEVEL) 37 REAL CSUBP, RSFI,RSFV,F0PI40 REAL CSUBP,F0PI 38 41 REAL UBARI,UBARV,UBAR0 39 40 real dist_sol41 42 42 43 COMMON /VISGAS/SOLARF,NTERM,PEXPON, … … 60 61 & FMNETV 61 62 62 COMMON /PLANT/ CSUBP, RSFI,RSFV,F0PI63 COMMON /PLANT/ CSUBP,F0PI 63 64 COMMON /UBARED/ UBARI,UBARV,UBAR0 64 65 … … 95 96 DO 912 NT=1,NTERM(K) 96 97 IF (ICLD.eq.1) THEN 97 BSURF=0.+ RSFV*UBAR0*F0PI*EXP(-TAUV(ig,NLEVEL,K,NT)/UBAR0)98 BSURF=0.+ falbe(ig)*UBAR0*F0PI*EXP(-TAUV(ig,NLEVEL,K,NT)/UBAR0) 98 99 ELSE 99 BSURF=0.+ RSFV*UBAR0*F0PI*EXP(-TAUVP(ig,NLEVEL,K,NT)/UBAR0)100 BSURF=0.+ falbe(ig)*UBAR0*F0PI*EXP(-TAUVP(ig,NLEVEL,K,NT)/UBAR0) 100 101 ENDIF 101 102 C … … 139 140 c write(*,*) (CO0(J),J=1,NLAYER) 140 141 c print*,'UBAR0 ',UBAR0 141 c print*,NLEVEL,WNOV(K),F0PI, RSFV,BTOP,BSURF142 c print*,NLEVEL,WNOV(K),F0PI,falbe(ig),BTOP,BSURF 142 143 FUW = 0.0 143 144 FDW = 0.0 … … 146 147 147 148 CALL GFLUXV(NLEVEL,WNOV(K),DT0,T0, 148 & WB0,CO0,F0PI, RSFV,BTOP,BSURF,FUW,FDW,FMUPV,149 & WB0,CO0,F0PI,falbe(ig),BTOP,BSURF,FUW,FDW,FMUPV, 149 150 & FMDV,IPRINT) 150 151 c PRINT*,'sortie gfluxv #: ',ig,K -
trunk/LMDZ.TITAN/libf/phytitan/tgmdat.F
r306 r495 4 4 COMMON /UBARED/ UBARI,UBARV,UBAR0 5 5 COMMON /LAPSE/ DTDP(NLAYER),CONVEQ 6 COMMON /PLANT/ CSUBP, RSFI,RSFV,F0PI6 COMMON /PLANT/ CSUBP,F0PI 7 7 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON 8 8 COMMON /CONST/RGAS,RHOP,PI,SIGMA … … 29 29 C CONSTANT MUST BE CHANGED 30 30 DATA RHOP/1.E4/ 31 C RSF IS THE SURFACE REFLECTANCE FOR VIS AND IR32 DATA RSFV,RSFI/0.15,0.05/33 31 C FOPI IS THE ACTUAL SOLAR FLUX IN ERGS/CM2 34 32 DATA F0PI/1.5E4/ -
trunk/LMDZ.TITAN/libf/phytitan/ustarhb.F
r102 r495 41 41 zxu = u(i,1)*zx_alf1+u(i,2)*zx_alf2 42 42 zxv = v(i,1)*zx_alf1+v(i,2)*zx_alf2 43 zxmod = 1.0+SQRT(zxu**2+zxv**2)43 zxmod = SQRT(zxu**2+zxv**2) 44 44 taux = zxu *zxmod*cd_m(i) 45 45 tauy = zxv *zxmod*cd_m(i)
Note: See TracChangeset
for help on using the changeset viewer.