Changeset 1461 for trunk/LMDZ.TITAN/libf/phytitan
- Timestamp:
- Jun 23, 2015, 10:34:03 AM (9 years ago)
- Location:
- trunk/LMDZ.TITAN/libf/phytitan
- Files:
-
- 15 edited
- 1 moved
-
cld.F (modified) (2 diffs)
-
cooling.F (modified) (2 diffs)
-
gasses.F (modified) (2 diffs)
-
gfluxv.F (modified) (2 diffs)
-
heating.F (modified) (3 diffs)
-
optci.F (modified) (2 diffs)
-
optci_1pt.F (modified) (2 diffs)
-
optci_1pt_2.F (modified) (2 diffs)
-
optci_1pt_3.F (modified) (3 diffs)
-
optcv.F (modified) (2 diffs)
-
optcv_1pt.F (modified) (1 diff)
-
optcv_1pt_2.F (modified) (2 diffs)
-
optcv_1pt_3.F (modified) (3 diffs)
-
radtitan.F (modified) (4 diffs)
-
sfluxv.F (modified) (3 diffs)
-
tgmdat_mod.F90 (moved) (moved from trunk/LMDZ.TITAN/libf/phytitan/tgmdat.F) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.TITAN/libf/phytitan/cld.F
r495 r1461 3 3 C THIS ROUTINE SETS UP THE CLOUD DISTRIBUTION 4 4 C 5 USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, 6 & RCLOUD,FARGON 7 USE TGMDAT_MOD, ONLY: PI 5 8 #include "dimensions.h" 6 9 PARAMETER(NLAYER=llm,NLEVEL=NLAYER+1) … … 11 14 COMMON /CLOUD/ RADCLD(NLAYER), XNCLD(NLAYER) 12 15 & , RCLDI(NSPECI), XICLDI(NSPECI), RCLDV(NSPECV), XICLDV(NSPECV) 13 COMMON /PLANT/ CSUBP,F0PI14 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON15 COMMON /CONST/RGAS,RHOP,PI,SIGMA16 16 TOTALC=0.0 17 17 CCC -
trunk/LMDZ.TITAN/libf/phytitan/cooling.F
r1356 r1461 56 56 57 57 use dimphy 58 use tgmdat_mod, only: CSUBP,F0PI 58 59 IMPLICIT NONE 59 60 #include "dimensions.h" … … 85 86 REAL dtauip(ngrid,NLAYER,NSPECI) 86 87 COMMON /IRTAUS/ dtaui,dtauip 87 88 COMMON /PLANT/ CSUBP,F0PI89 REAL CSUBP,F0PI90 88 91 89 c Local: -
trunk/LMDZ.TITAN/libf/phytitan/gasses.F
r1072 r1461 2 2 C THIS SUBROUTINE SETS UP THE MASS MIXING RATIOS OF THE 3 3 C OPTICALLY ACTIVE GASES: CH4, C2H2, AND C2H6 4 #include "dimensions.h" 4 USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, 5 & RCLOUD,FARGON 6 USE TGMDAT_MOD, ONLY: RGAS,RHOP,PI,SIGMA 7 include "dimensions.h" 5 8 PARAMETER(NLAYER=llm,NLEVEL=NLAYER+1) 6 9 COMMON /ATM/ Z(NLEVEL),PRESS(NLEVEL),DEN(NLEVEL),TEMP(NLEVEL) … … 8 11 & ,XMU(NLEVEL),GAS1(NLAYER),COLDEN(NLAYER) 9 12 COMMON /STRATO/ C2H2(NLAYER),C2H6(NLAYER) 10 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON11 COMMON /CONST/RGAS,RHOP,PI,SIGMA12 13 C* 13 14 C NOW CALCULATE THE LAYER AVERAGE GAS MIXING RATIOS. -
trunk/LMDZ.TITAN/libf/phytitan/gfluxv.F
r3 r1461 2 2 & , F0PI,RSF,BTOP,BSURF,FP,FM,FMIDP,FMIDM,IPRINT) 3 3 c PARAMETER (NL=101) 4 USE TGMDAT_MOD, ONLY: UBARI,UBARV,UBAR0 4 5 IMPLICIT NONE 5 6 INTEGER IPRINT,I,J,K,IDELTA,NAYER,NL,NTODO … … 17 18 C J.A.S., 37, 630-642, 1980. 18 19 REAL B81,B82,R81 19 REAL AP,AM,F0PI,DENOM, UBARI,UBARV,UBAR0,CMMID,TAUMID,RSF,BTOP20 REAL AP,AM,F0PI,DENOM,CMMID,TAUMID,RSF,BTOP 20 21 REAL EM,EP,PI,WAVEN,CPMID,BSURF,G4 21 COMMON /UBARED/ UBARI,UBARV,UBAR022 22 C THIS NEXT ROW OF VARIABLES ARE THOSE ACTUALLY USED IN THE 23 23 C ROUTINE -
trunk/LMDZ.TITAN/libf/phytitan/heating.F
r1356 r1461 33 33 34 34 use dimphy 35 use TGMDAT_MOD, ONLY: UBARI,UBARV,UBAR0 36 use TGMDAT_MOD, ONLY: CSUBP,F0PI 35 37 IMPLICIT NONE 36 #include "dimensions.h"38 include "dimensions.h" 37 39 38 40 INTEGER NLEVEL,NLAYER,NSPECV … … 64 66 c --------------------------------------------- 65 67 66 REAL UBARI,UBARV,UBAR067 68 REAL CH4(NLEVEL),XN2(NLEVEL),H2(NLEVEL),AR(NLEVEL) 68 69 REAL XMU(NLEVEL),GAS1(NLAYER),COLDEN(NLAYER) 69 70 REAL FNETV(ngrid,NLEVEL),FUPV(ngrid,NLEVEL,NSPECV) 70 71 REAL FDV(ngrid,NLEVEL,NSPECV),FMNETV(ngrid,NLEVEL) 71 REAL CSUBP,F0PI72 73 COMMON /UBARED/ UBARI,UBARV,UBAR074 72 75 73 COMMON /GASS/ CH4,XN2 … … 82 80 & FDV, 83 81 & FMNETV 84 85 COMMON /PLANT/ CSUBP,F0PI86 82 87 83 -
trunk/LMDZ.TITAN/libf/phytitan/optci.F
r1126 r1461 3 3 use infotrac 4 4 use common_mod, only:rmcbar,xfbar,ncount,TauHID,TauCID,TauGID 5 USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, 6 & RCLOUD,FARGON 5 7 #include "dimensions.h" 6 8 #include "microtab.h" … … 59 61 COMMON /IRTAUS/ DTAUP,DTAUPP 60 62 61 COMMON /PLANT/ CSUBP,F0PI62 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON63 COMMON /CONST/RGAS,RHOP,PI,SIGMA64 63 COMMON /part/v,rayon,vrat,dr,dv 65 64 -
trunk/LMDZ.TITAN/libf/phytitan/optci_1pt.F
r808 r1461 1 1 SUBROUTINE optci_1pt(zqaer_1pt,rcdb,xfrb,iopti,IPRINT) 2 2 use dimphy 3 USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, 4 & RCLOUD,FARGON 5 USE TGMDAT_MOD, ONLY: RGAS 3 6 #include "dimensions.h" 4 7 #include "microtab.h" … … 38 41 & DWNI(NSPECI), WLNI(NSPECI) 39 42 40 COMMON /PLANT/ CSUBP,F0PI41 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON42 COMMON /CONST/RGAS,RHOP,PI,SIGMA43 43 COMMON /part/v,rayon,vrat,dr,dv 44 44 -
trunk/LMDZ.TITAN/libf/phytitan/optci_1pt_2.F
r808 r1461 2 2 3 3 use dimphy 4 USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, 5 & RCLOUD,FARGON 6 USE TGMDAT_MOD, ONLY: RGAS 4 7 #include "dimensions.h" 5 8 #include "microtab.h" … … 39 42 & DWNI(NSPECI), WLNI(NSPECI) 40 43 41 COMMON /PLANT/ CSUBP,F0PI42 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON43 COMMON /CONST/RGAS,RHOP,PI,SIGMA44 44 COMMON /part/v,rayon,vrat,dr,dv 45 45 -
trunk/LMDZ.TITAN/libf/phytitan/optci_1pt_3.F
r1126 r1461 1 1 SUBROUTINE optci_1pt3(zqaer_1pt,rcdb,xfrb,iopti,IPRINT) 2 2 use dimphy 3 USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, 4 & RCLOUD,FARGON 5 USE TGMDAT_MOD, ONLY: RGAS 3 6 IMPLICIT NONE 4 7 #include "dimensions.h" … … 9 12 integer nlayer, nlevel, nspeci, nspc1i, nspecv, nspc1v 10 13 real z,press, den, temp, ch4, xn2, h2, ar, xmu, gas1, 11 & colden, c2h2, c2h6, hcn, radius, xnumb, reali14 & colden, c2h2, c2h6, hcn, radius, xnumb, 12 15 & ximgi, realv, ximgv, rcldi, xicldi, rcldv, xicldv, rcldi2, 13 & xicldi2, rcldv2, xicldv2,real bwni, wnoi, dwni, wlni, csubp, 14 & f0pi, rhch4, fh2, fhaze, fhvis 15 & reali, ximgi, bwni, fhir, taufac, rcloud, fargon, rgas, rhop, 16 & pi, sigma, prod,reali,fhvis 16 & xicldi2, rcldv2, xicldv2 17 real bwni, wnoi, dwni, wlni, 18 & prod,reali 17 19 18 20 integer k, j,inq,kgas … … 55 57 & DWNI(NSPECI), WLNI(NSPECI) 56 58 57 COMMON /PLANT/ CSUBP,F0PI58 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON59 COMMON /CONST/RGAS,RHOP,PI,SIGMA60 59 COMMON /part/v,rayon,vrat,dr,dv 61 60 -
trunk/LMDZ.TITAN/libf/phytitan/optcv.F
r1126 r1461 4 4 use infotrac 5 5 use common_mod, only:rmcbar,xfbar,ncount,TauHVD,TauCVD,TauGVD 6 USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, 7 & RCLOUD,FARGON 6 8 #include "dimensions.h" 7 9 #include "microtab.h" … … 55 57 & ,DWNV(NSPECV),WLNV(NSPECV) 56 58 57 COMMON /PLANT/ CSUBP,F0PI58 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON59 COMMON /CONST/ RGAS,RHOP,PI,SIGMA60 59 COMMON /part/ v(nrad),rayon(nrad),vrat,dr(nrad),dv(nrad) 61 60 -
trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt.F
r495 r1461 39 39 & ,DWNV(NSPECV),WLNV(NSPECV) 40 40 41 COMMON /PLANT/ CSUBP,F0PI42 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON43 COMMON /CONST/ RGAS,RHOP,PI,SIGMA44 41 * nrad dans microtab.h 45 42 COMMON /part/ v(nrad),rayon(nrad),vrat,dr(nrad),dv(nrad) -
trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt_2.F
r495 r1461 2 2 3 3 use dimphy 4 USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, 5 & RCLOUD,FARGON 4 6 #include "dimensions.h" 5 7 #include "microtab.h" … … 38 40 & ,DWNV(NSPECV),WLNV(NSPECV) 39 41 40 COMMON /PLANT/ CSUBP,F0PI41 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON42 COMMON /CONST/ RGAS,RHOP,PI,SIGMA43 42 * nrad dans microtab.h 44 43 COMMON /part/ v(nrad),rayon(nrad),vrat,dr(nrad),dv(nrad) -
trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt_3.F
r1126 r1461 2 2 3 3 use dimphy 4 USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, 5 & RCLOUD,FARGON 4 6 IMPLICIT NONE 5 7 #include "dimensions.h" … … 9 11 integer nlayer, nlevel, nspeci, nspc1i, nspecv, nspc1v,nterm 10 12 real z,press, den, temp, ch4, xn2, h2, ar, xmu, gas1, 11 & colden, c2h2, c2h6, hcn, radius, xnumb, reali 13 & colden, c2h2, c2h6, hcn, radius, xnumb, reali, 12 14 & ximgi, realv, ximgv, rcldi, xicldi, rcldv, xicldv, rcldi2, 13 & xicldi2, rcldv2, xicldv2,real bwni, wnoi, dwni, wlni, csubp, 14 & f0pi, rhch4, fh2, fhaze, fhvis 15 & reali, ximgi, bwni, fhir, taufac, rcloud, fargon, rgas, rhop, 16 & pi, sigma, prod,reali,fhvis 15 & xicldi2, rcldv2, xicldv2 16 real bwni, wnoi, dwni, wlni,prod 17 17 18 18 integer k, j,inq,nt … … 58 58 & ,DWNV(NSPECV),WLNV(NSPECV) 59 59 60 COMMON /PLANT/ CSUBP,F0PI61 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON62 COMMON /CONST/ RGAS,RHOP,PI,SIGMA63 60 * nrad dans microtab.h 64 61 COMMON /part/ v(nrad),rayon(nrad),vrat,dr(nrad),dv(nrad) -
trunk/LMDZ.TITAN/libf/phytitan/radtitan.F
r1058 r1461 36 36 USE optcld, only : iniqcld 37 37 use moyzon_mod, only:plevmoy 38 USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, 39 & RCLOUD,FARGON 40 USE TGMDAT_MOD, ONLY: RHOP 38 41 IMPLICIT NONE 39 42 #include "dimensions.h" … … 88 91 REAL XMU(NLEVEL),GAS1(NLAYER),COLDEN(NLAYER) 89 92 REAL C2H2(NLAYER),C2H6(NLAYER),HCN(NLAYER) 90 REAL RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON91 REAL RGAS,RHOP,PI,SIGMA92 93 93 94 COMMON /VERTICAL/ DZED … … 104 105 COMMON /STRATO/ C2H2,C2H6 105 106 COMMON /STRAT2/ HCN 106 107 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON108 COMMON /CONST/RGAS,RHOP,PI,SIGMA109 107 110 108 c----------------------------------------------------------------------- … … 179 177 C NOW CALCULATE THE LAYER AVERAGE GAS MIXING RATIOS. 180 178 CALL GASSES(IPRINT) 181 179 182 180 else 183 181 c------------------------------------------------------ -
trunk/LMDZ.TITAN/libf/phytitan/sfluxv.F
r495 r1461 2 2 3 3 use dimphy 4 USE TGMDAT_MOD, ONLY: UBARI,UBARV,UBAR0 5 USE TGMDAT_MOD, ONLY: F0PI 4 6 IMPLICIT NONE 5 7 #include "dimensions.h" … … 38 40 & FMNETV(ngrid,NLEVEL), 39 41 & FMUPV(NLEVEL),FMDV(NLEVEL) 40 REAL CSUBP,F0PI41 REAL UBARI,UBARV,UBAR042 42 43 43 COMMON /VISGAS/SOLARF,NTERM,PEXPON, … … 60 60 & FDV, 61 61 & FMNETV 62 63 COMMON /PLANT/ CSUBP,F0PI64 COMMON /UBARED/ UBARI,UBARV,UBAR065 62 66 63 -
trunk/LMDZ.TITAN/libf/phytitan/tgmdat_mod.F90
r1451 r1461 1 BLOCK DATA TGMDAT 2 #include "dimensions.h" 3 PARAMETER(NLAYER=llm,NLEVEL=NLAYER+1) 4 COMMON /UBARED/ UBARI,UBARV,UBAR0 5 COMMON /LAPSE/ DTDP(NLAYER),CONVEQ 6 COMMON /PLANT/ CSUBP,F0PI 7 COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON 8 COMMON /CONST/RGAS,RHOP,PI,SIGMA 9 C 10 DATA PI/3.14159265358979323846/ 11 C RGAS IS THE UNIVERSAL GAS CONSTANT IN UNITS OF: M SEC-2 AMU K-1 KM 12 DATA RGAS/8.31432/ 13 C SIGMA IS THE STEFAN-BOLTZMAN CONSTATNT IN CGS UNITS 14 DATA SIGMA/5.6677E-5/ 15 C 16 C RT CONSTANTS 17 DATA UBARI,UBARV,UBAR0/0.5,0.5,0.5/ 18 C PLANET SPECIFIC CONSTANTS 19 C 20 C CONVEQ IS THE DRY N2 ADIABATE DLNT/DLNP 21 DATA CONVEQ/0.2994/ 22 C CSUBP IS THE SPECIFIC HEAT AT CONSTANT P OF THE ATMOSPHERE 23 C IN UNITS OF ERGS K-1 G-1 24 DATA CSUBP/1.039E7/ 25 C RHOP IS THE UNITS CONVERSION FROM TO GET MASS UNITS (G CM-2) 26 C FROM PRESSURE (BARS) DEVIDED BY GRAVITY (M SEC-2) 27 C IS EQUAL TO ONE GM CM-2 BARS-1 M SEC-2 28 C IF ONE WHATS TO CHANGE UNITS ON PRESSURE THIS 29 C CONSTANT MUST BE CHANGED 30 DATA RHOP/1.E4/ 31 C FOPI IS THE ACTUAL SOLAR FLUX IN ERGS/CM2 32 DATA F0PI/1.5E4/ 33 C RHCH4 IS THE METHANE RH AT THE SURFACE 34 DATA RHCH4 /0.60/ !! .65 35 C FH2 IS THE CONSTANT MIXING RATIO OF H2 36 DATA FH2/0.04/ !! 0.003 37 C FHAZE IS THE HAZE PRODUCTION SCALING FACTOR 38 DATA FHAZE/0.35/ 39 C FHIR IS THE HAZE INFRARED ABSORPTION SCALE FACTOR 40 DATA FHIR/0.7/ ! anciennement 0.5 41 C FHVIS IS THE HAZE INFRARED ABSORPTION SCALE FACTOR 42 DATA FHVIS/1.1/ !!!! 1.333333333/ 43 C TAUFAC IS THE 200 CM-1 SCALING FACTOR 44 DATA TAUFAC/2.00/ 45 C RCLOUD IS THE PARTICLE SIZE IN THE CLOUD IN MICRONS 46 DATA RCLOUD/60./ !! 100 1 MODULE TGMDAT_MOD 2 IMPLICIT NONE 3 include "dimensions.h" 47 4 48 data FARGON/0./ 49 END 5 INTEGER,PARAMETER :: NLAYER=llm,NLEVEL=NLAYER+1 6 REAL,SAVE :: PI=3.14159265358979323846 7 ! RGAS IS THE UNIVERSAL GAS CONSTANT IN UNITS OF: M SEC-2 AMU K-1 KM 8 REAL,SAVE :: RGAS=8.31432 9 ! SIGMA IS THE STEFAN-BOLTZMAN CONSTATNT IN CGS UNITS 10 REAL,SAVE :: SIGMA=5.6677E-5 11 ! 12 ! RT CONSTANTS 13 REAL,SAVE :: UBARI=0.5,UBARV=0.5,UBAR0=0.5 14 ! PLANET SPECIFIC CONSTANTS 15 ! 16 ! CONVEQ IS THE DRY N2 ADIABATE DLNT/DLNP 17 REAL,SAVE :: CONVEQ=0.2994 18 ! CSUBP IS THE SPECIFIC HEAT AT CONSTANT P OF THE ATMOSPHERE 19 ! IN UNITS OF ERGS K-1 G-1 20 REAL,SAVE :: CSUBP=1.039E7 21 ! RHOP IS THE UNITS CONVERSION FROM TO GET MASS UNITS (G CM-2) 22 ! FROM PRESSURE (BARS) DEVIDED BY GRAVITY (M SEC-2) 23 ! IS EQUAL TO ONE GM CM-2 BARS-1 M SEC-2 24 ! IF ONE WHATS TO CHANGE UNITS ON PRESSURE THIS 25 ! CONSTANT MUST BE CHANGED 26 REAL,SAVE :: RHOP=1.E4 27 ! FOPI IS THE ACTUAL SOLAR FLUX IN ERGS/CM2 28 REAL,SAVE :: F0PI=1.5E4 29 ! RHCH4 IS THE METHANE RH AT THE SURFACE 30 REAL,SAVE :: RHCH4=0.60 !! .65 31 ! FH2 IS THE CONSTANT MIXING RATIO OF H2 32 REAL,SAVE :: FH2=0.04 !! 0.003 33 ! FHAZE IS THE HAZE PRODUCTION SCALING FACTOR 34 REAL,SAVE :: FHAZE=0.35 35 ! FHIR IS THE HAZE INFRARED ABSORPTION SCALE FACTOR 36 REAL,SAVE :: FHIR=0.7 ! anciennement 0.5 37 ! FHVIS IS THE HAZE INFRARED ABSORPTION SCALE FACTOR 38 REAL,SAVE :: FHVIS=1.1 !!!! 1.333333333/ 39 ! TAUFAC IS THE 200 CM-1 SCALING FACTOR 40 REAL,SAVE :: TAUFAC=2.00 41 ! RCLOUD IS THE PARTICLE SIZE IN THE CLOUD IN MICRONS 42 REAL,SAVE ::RCLOUD=60. !! 100 43 44 REAL,SAVE :: FARGON=0. 45 END MODULE TGMDAT_MOD
Note: See TracChangeset
for help on using the changeset viewer.
