Index: /trunk/LMDZ.TITAN/libf/phytitan/cld.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/cld.F (revision 1460) +++ /trunk/LMDZ.TITAN/libf/phytitan/cld.F (revision 1461) @@ -3,4 +3,7 @@ C THIS ROUTINE SETS UP THE CLOUD DISTRIBUTION C + USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, + & RCLOUD,FARGON + USE TGMDAT_MOD, ONLY: PI #include "dimensions.h" PARAMETER(NLAYER=llm,NLEVEL=NLAYER+1) @@ -11,7 +14,4 @@ COMMON /CLOUD/ RADCLD(NLAYER), XNCLD(NLAYER) & , RCLDI(NSPECI), XICLDI(NSPECI), RCLDV(NSPECV), XICLDV(NSPECV) - COMMON /PLANT/ CSUBP,F0PI - COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON - COMMON /CONST/RGAS,RHOP,PI,SIGMA TOTALC=0.0 CCC Index: /trunk/LMDZ.TITAN/libf/phytitan/cooling.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/cooling.F (revision 1460) +++ /trunk/LMDZ.TITAN/libf/phytitan/cooling.F (revision 1461) @@ -56,4 +56,5 @@ use dimphy + use tgmdat_mod, only: CSUBP,F0PI IMPLICIT NONE #include "dimensions.h" @@ -85,7 +86,4 @@ REAL dtauip(ngrid,NLAYER,NSPECI) COMMON /IRTAUS/ dtaui,dtauip - - COMMON /PLANT/ CSUBP,F0PI - REAL CSUBP,F0PI c Local: Index: /trunk/LMDZ.TITAN/libf/phytitan/gasses.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/gasses.F (revision 1460) +++ /trunk/LMDZ.TITAN/libf/phytitan/gasses.F (revision 1461) @@ -2,5 +2,8 @@ C THIS SUBROUTINE SETS UP THE MASS MIXING RATIOS OF THE C OPTICALLY ACTIVE GASES: CH4, C2H2, AND C2H6 -#include "dimensions.h" + USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, + & RCLOUD,FARGON + USE TGMDAT_MOD, ONLY: RGAS,RHOP,PI,SIGMA + include "dimensions.h" PARAMETER(NLAYER=llm,NLEVEL=NLAYER+1) COMMON /ATM/ Z(NLEVEL),PRESS(NLEVEL),DEN(NLEVEL),TEMP(NLEVEL) @@ -8,6 +11,4 @@ & ,XMU(NLEVEL),GAS1(NLAYER),COLDEN(NLAYER) COMMON /STRATO/ C2H2(NLAYER),C2H6(NLAYER) - COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON - COMMON /CONST/RGAS,RHOP,PI,SIGMA C* C NOW CALCULATE THE LAYER AVERAGE GAS MIXING RATIOS. Index: /trunk/LMDZ.TITAN/libf/phytitan/gfluxv.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/gfluxv.F (revision 1460) +++ /trunk/LMDZ.TITAN/libf/phytitan/gfluxv.F (revision 1461) @@ -2,4 +2,5 @@ & , F0PI,RSF,BTOP,BSURF,FP,FM,FMIDP,FMIDM,IPRINT) c PARAMETER (NL=101) + USE TGMDAT_MOD, ONLY: UBARI,UBARV,UBAR0 IMPLICIT NONE INTEGER IPRINT,I,J,K,IDELTA,NAYER,NL,NTODO @@ -17,7 +18,6 @@ C J.A.S., 37, 630-642, 1980. REAL B81,B82,R81 - REAL AP,AM,F0PI,DENOM,UBARI,UBARV,UBAR0,CMMID,TAUMID,RSF,BTOP + REAL AP,AM,F0PI,DENOM,CMMID,TAUMID,RSF,BTOP REAL EM,EP,PI,WAVEN,CPMID,BSURF,G4 - COMMON /UBARED/ UBARI,UBARV,UBAR0 C THIS NEXT ROW OF VARIABLES ARE THOSE ACTUALLY USED IN THE C ROUTINE Index: /trunk/LMDZ.TITAN/libf/phytitan/heating.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/heating.F (revision 1460) +++ /trunk/LMDZ.TITAN/libf/phytitan/heating.F (revision 1461) @@ -33,6 +33,8 @@ use dimphy + use TGMDAT_MOD, ONLY: UBARI,UBARV,UBAR0 + use TGMDAT_MOD, ONLY: CSUBP,F0PI IMPLICIT NONE -#include "dimensions.h" + include "dimensions.h" INTEGER NLEVEL,NLAYER,NSPECV @@ -64,12 +66,8 @@ c --------------------------------------------- - REAL UBARI,UBARV,UBAR0 REAL CH4(NLEVEL),XN2(NLEVEL),H2(NLEVEL),AR(NLEVEL) REAL XMU(NLEVEL),GAS1(NLAYER),COLDEN(NLAYER) REAL FNETV(ngrid,NLEVEL),FUPV(ngrid,NLEVEL,NSPECV) REAL FDV(ngrid,NLEVEL,NSPECV),FMNETV(ngrid,NLEVEL) - REAL CSUBP,F0PI - - COMMON /UBARED/ UBARI,UBARV,UBAR0 COMMON /GASS/ CH4,XN2 @@ -82,6 +80,4 @@ & FDV, & FMNETV - - COMMON /PLANT/ CSUBP,F0PI Index: /trunk/LMDZ.TITAN/libf/phytitan/optci.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/optci.F (revision 1460) +++ /trunk/LMDZ.TITAN/libf/phytitan/optci.F (revision 1461) @@ -3,4 +3,6 @@ use infotrac use common_mod, only:rmcbar,xfbar,ncount,TauHID,TauCID,TauGID + USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, + & RCLOUD,FARGON #include "dimensions.h" #include "microtab.h" @@ -59,7 +61,4 @@ COMMON /IRTAUS/ DTAUP,DTAUPP - COMMON /PLANT/ CSUBP,F0PI - COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON - COMMON /CONST/RGAS,RHOP,PI,SIGMA COMMON /part/v,rayon,vrat,dr,dv Index: /trunk/LMDZ.TITAN/libf/phytitan/optci_1pt.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/optci_1pt.F (revision 1460) +++ /trunk/LMDZ.TITAN/libf/phytitan/optci_1pt.F (revision 1461) @@ -1,4 +1,7 @@ SUBROUTINE optci_1pt(zqaer_1pt,rcdb,xfrb,iopti,IPRINT) use dimphy + USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, + & RCLOUD,FARGON + USE TGMDAT_MOD, ONLY: RGAS #include "dimensions.h" #include "microtab.h" @@ -38,7 +41,4 @@ & DWNI(NSPECI), WLNI(NSPECI) - COMMON /PLANT/ CSUBP,F0PI - COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON - COMMON /CONST/RGAS,RHOP,PI,SIGMA COMMON /part/v,rayon,vrat,dr,dv Index: /trunk/LMDZ.TITAN/libf/phytitan/optci_1pt_2.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/optci_1pt_2.F (revision 1460) +++ /trunk/LMDZ.TITAN/libf/phytitan/optci_1pt_2.F (revision 1461) @@ -2,4 +2,7 @@ use dimphy + USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, + & RCLOUD,FARGON + USE TGMDAT_MOD, ONLY: RGAS #include "dimensions.h" #include "microtab.h" @@ -39,7 +42,4 @@ & DWNI(NSPECI), WLNI(NSPECI) - COMMON /PLANT/ CSUBP,F0PI - COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON - COMMON /CONST/RGAS,RHOP,PI,SIGMA COMMON /part/v,rayon,vrat,dr,dv Index: /trunk/LMDZ.TITAN/libf/phytitan/optci_1pt_3.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/optci_1pt_3.F (revision 1460) +++ /trunk/LMDZ.TITAN/libf/phytitan/optci_1pt_3.F (revision 1461) @@ -1,4 +1,7 @@ SUBROUTINE optci_1pt3(zqaer_1pt,rcdb,xfrb,iopti,IPRINT) use dimphy + USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, + & RCLOUD,FARGON + USE TGMDAT_MOD, ONLY: RGAS IMPLICIT NONE #include "dimensions.h" @@ -9,10 +12,9 @@ integer nlayer, nlevel, nspeci, nspc1i, nspecv, nspc1v real z,press, den, temp, ch4, xn2, h2, ar, xmu, gas1, - & colden, c2h2, c2h6, hcn, radius, xnumb, reali + & colden, c2h2, c2h6, hcn, radius, xnumb, & ximgi, realv, ximgv, rcldi, xicldi, rcldv, xicldv, rcldi2, - & xicldi2, rcldv2, xicldv2,real bwni, wnoi, dwni, wlni, csubp, - & f0pi, rhch4, fh2, fhaze, fhvis - & reali, ximgi, bwni, fhir, taufac, rcloud, fargon, rgas, rhop, - & pi, sigma, prod,reali,fhvis + & xicldi2, rcldv2, xicldv2 + real bwni, wnoi, dwni, wlni, + & prod,reali integer k, j,inq,kgas @@ -55,7 +57,4 @@ & DWNI(NSPECI), WLNI(NSPECI) - COMMON /PLANT/ CSUBP,F0PI - COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON - COMMON /CONST/RGAS,RHOP,PI,SIGMA COMMON /part/v,rayon,vrat,dr,dv Index: /trunk/LMDZ.TITAN/libf/phytitan/optcv.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/optcv.F (revision 1460) +++ /trunk/LMDZ.TITAN/libf/phytitan/optcv.F (revision 1461) @@ -4,4 +4,6 @@ use infotrac use common_mod, only:rmcbar,xfbar,ncount,TauHVD,TauCVD,TauGVD + USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, + & RCLOUD,FARGON #include "dimensions.h" #include "microtab.h" @@ -55,7 +57,4 @@ & ,DWNV(NSPECV),WLNV(NSPECV) - COMMON /PLANT/ CSUBP,F0PI - COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON - COMMON /CONST/ RGAS,RHOP,PI,SIGMA COMMON /part/ v(nrad),rayon(nrad),vrat,dr(nrad),dv(nrad) Index: /trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt.F (revision 1460) +++ /trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt.F (revision 1461) @@ -39,7 +39,4 @@ & ,DWNV(NSPECV),WLNV(NSPECV) - COMMON /PLANT/ CSUBP,F0PI - COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON - COMMON /CONST/ RGAS,RHOP,PI,SIGMA * nrad dans microtab.h COMMON /part/ v(nrad),rayon(nrad),vrat,dr(nrad),dv(nrad) Index: /trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt_2.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt_2.F (revision 1460) +++ /trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt_2.F (revision 1461) @@ -2,4 +2,6 @@ use dimphy + USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, + & RCLOUD,FARGON #include "dimensions.h" #include "microtab.h" @@ -38,7 +40,4 @@ & ,DWNV(NSPECV),WLNV(NSPECV) - COMMON /PLANT/ CSUBP,F0PI - COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON - COMMON /CONST/ RGAS,RHOP,PI,SIGMA * nrad dans microtab.h COMMON /part/ v(nrad),rayon(nrad),vrat,dr(nrad),dv(nrad) Index: /trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt_3.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt_3.F (revision 1460) +++ /trunk/LMDZ.TITAN/libf/phytitan/optcv_1pt_3.F (revision 1461) @@ -2,4 +2,6 @@ use dimphy + USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, + & RCLOUD,FARGON IMPLICIT NONE #include "dimensions.h" @@ -9,10 +11,8 @@ integer nlayer, nlevel, nspeci, nspc1i, nspecv, nspc1v,nterm real z,press, den, temp, ch4, xn2, h2, ar, xmu, gas1, - & colden, c2h2, c2h6, hcn, radius, xnumb, reali + & colden, c2h2, c2h6, hcn, radius, xnumb, reali, & ximgi, realv, ximgv, rcldi, xicldi, rcldv, xicldv, rcldi2, - & xicldi2, rcldv2, xicldv2,real bwni, wnoi, dwni, wlni, csubp, - & f0pi, rhch4, fh2, fhaze, fhvis - & reali, ximgi, bwni, fhir, taufac, rcloud, fargon, rgas, rhop, - & pi, sigma, prod,reali,fhvis + & xicldi2, rcldv2, xicldv2 + real bwni, wnoi, dwni, wlni,prod integer k, j,inq,nt @@ -58,7 +58,4 @@ & ,DWNV(NSPECV),WLNV(NSPECV) - COMMON /PLANT/ CSUBP,F0PI - COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON - COMMON /CONST/ RGAS,RHOP,PI,SIGMA * nrad dans microtab.h COMMON /part/ v(nrad),rayon(nrad),vrat,dr(nrad),dv(nrad) Index: /trunk/LMDZ.TITAN/libf/phytitan/radtitan.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/radtitan.F (revision 1460) +++ /trunk/LMDZ.TITAN/libf/phytitan/radtitan.F (revision 1461) @@ -36,4 +36,7 @@ USE optcld, only : iniqcld use moyzon_mod, only:plevmoy + USE TGMDAT_MOD, ONLY: RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC, + & RCLOUD,FARGON + USE TGMDAT_MOD, ONLY: RHOP IMPLICIT NONE #include "dimensions.h" @@ -88,6 +91,4 @@ REAL XMU(NLEVEL),GAS1(NLAYER),COLDEN(NLAYER) REAL C2H2(NLAYER),C2H6(NLAYER),HCN(NLAYER) - REAL RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON - REAL RGAS,RHOP,PI,SIGMA COMMON /VERTICAL/ DZED @@ -104,7 +105,4 @@ COMMON /STRATO/ C2H2,C2H6 COMMON /STRAT2/ HCN - - COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON - COMMON /CONST/RGAS,RHOP,PI,SIGMA c----------------------------------------------------------------------- @@ -179,5 +177,5 @@ C NOW CALCULATE THE LAYER AVERAGE GAS MIXING RATIOS. CALL GASSES(IPRINT) - + else c------------------------------------------------------ Index: /trunk/LMDZ.TITAN/libf/phytitan/sfluxv.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/sfluxv.F (revision 1460) +++ /trunk/LMDZ.TITAN/libf/phytitan/sfluxv.F (revision 1461) @@ -2,4 +2,6 @@ use dimphy + USE TGMDAT_MOD, ONLY: UBARI,UBARV,UBAR0 + USE TGMDAT_MOD, ONLY: F0PI IMPLICIT NONE #include "dimensions.h" @@ -38,6 +40,4 @@ & FMNETV(ngrid,NLEVEL), & FMUPV(NLEVEL),FMDV(NLEVEL) - REAL CSUBP,F0PI - REAL UBARI,UBARV,UBAR0 COMMON /VISGAS/SOLARF,NTERM,PEXPON, @@ -60,7 +60,4 @@ & FDV, & FMNETV - - COMMON /PLANT/ CSUBP,F0PI - COMMON /UBARED/ UBARI,UBARV,UBAR0 Index: unk/LMDZ.TITAN/libf/phytitan/tgmdat.F =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/tgmdat.F (revision 1460) +++ (revision ) @@ -1,49 +1,0 @@ - BLOCK DATA TGMDAT -#include "dimensions.h" - PARAMETER(NLAYER=llm,NLEVEL=NLAYER+1) - COMMON /UBARED/ UBARI,UBARV,UBAR0 - COMMON /LAPSE/ DTDP(NLAYER),CONVEQ - COMMON /PLANT/ CSUBP,F0PI - COMMON /ADJUST/ RHCH4,FH2,FHAZE,FHVIS,FHIR,TAUFAC,RCLOUD,FARGON - COMMON /CONST/RGAS,RHOP,PI,SIGMA -C - DATA PI/3.14159265358979323846/ -C RGAS IS THE UNIVERSAL GAS CONSTANT IN UNITS OF: M SEC-2 AMU K-1 KM - DATA RGAS/8.31432/ -C SIGMA IS THE STEFAN-BOLTZMAN CONSTATNT IN CGS UNITS - DATA SIGMA/5.6677E-5/ -C -C RT CONSTANTS - DATA UBARI,UBARV,UBAR0/0.5,0.5,0.5/ -C PLANET SPECIFIC CONSTANTS -C -C CONVEQ IS THE DRY N2 ADIABATE DLNT/DLNP - DATA CONVEQ/0.2994/ -C CSUBP IS THE SPECIFIC HEAT AT CONSTANT P OF THE ATMOSPHERE -C IN UNITS OF ERGS K-1 G-1 - DATA CSUBP/1.039E7/ -C RHOP IS THE UNITS CONVERSION FROM TO GET MASS UNITS (G CM-2) -C FROM PRESSURE (BARS) DEVIDED BY GRAVITY (M SEC-2) -C IS EQUAL TO ONE GM CM-2 BARS-1 M SEC-2 -C IF ONE WHATS TO CHANGE UNITS ON PRESSURE THIS -C CONSTANT MUST BE CHANGED - DATA RHOP/1.E4/ -C FOPI IS THE ACTUAL SOLAR FLUX IN ERGS/CM2 - DATA F0PI/1.5E4/ -C RHCH4 IS THE METHANE RH AT THE SURFACE - DATA RHCH4 /0.60/ !! .65 -C FH2 IS THE CONSTANT MIXING RATIO OF H2 - DATA FH2/0.04/ !! 0.003 -C FHAZE IS THE HAZE PRODUCTION SCALING FACTOR - DATA FHAZE/0.35/ -C FHIR IS THE HAZE INFRARED ABSORPTION SCALE FACTOR - DATA FHIR/0.7/ ! anciennement 0.5 -C FHVIS IS THE HAZE INFRARED ABSORPTION SCALE FACTOR - DATA FHVIS/1.1/ !!!! 1.333333333/ -C TAUFAC IS THE 200 CM-1 SCALING FACTOR - DATA TAUFAC/2.00/ -C RCLOUD IS THE PARTICLE SIZE IN THE CLOUD IN MICRONS - DATA RCLOUD/60./ !! 100 - - data FARGON/0./ - END Index: /trunk/LMDZ.TITAN/libf/phytitan/tgmdat_mod.F90 =================================================================== --- /trunk/LMDZ.TITAN/libf/phytitan/tgmdat_mod.F90 (revision 1461) +++ /trunk/LMDZ.TITAN/libf/phytitan/tgmdat_mod.F90 (revision 1461) @@ -0,0 +1,45 @@ +MODULE TGMDAT_MOD + IMPLICIT NONE + include "dimensions.h" + + INTEGER,PARAMETER :: NLAYER=llm,NLEVEL=NLAYER+1 + REAL,SAVE :: PI=3.14159265358979323846 +! RGAS IS THE UNIVERSAL GAS CONSTANT IN UNITS OF: M SEC-2 AMU K-1 KM + REAL,SAVE :: RGAS=8.31432 +! SIGMA IS THE STEFAN-BOLTZMAN CONSTATNT IN CGS UNITS + REAL,SAVE :: SIGMA=5.6677E-5 +! +! RT CONSTANTS + REAL,SAVE :: UBARI=0.5,UBARV=0.5,UBAR0=0.5 +! PLANET SPECIFIC CONSTANTS +! +! CONVEQ IS THE DRY N2 ADIABATE DLNT/DLNP + REAL,SAVE :: CONVEQ=0.2994 +! CSUBP IS THE SPECIFIC HEAT AT CONSTANT P OF THE ATMOSPHERE +! IN UNITS OF ERGS K-1 G-1 + REAL,SAVE :: CSUBP=1.039E7 +! RHOP IS THE UNITS CONVERSION FROM TO GET MASS UNITS (G CM-2) +! FROM PRESSURE (BARS) DEVIDED BY GRAVITY (M SEC-2) +! IS EQUAL TO ONE GM CM-2 BARS-1 M SEC-2 +! IF ONE WHATS TO CHANGE UNITS ON PRESSURE THIS +! CONSTANT MUST BE CHANGED + REAL,SAVE :: RHOP=1.E4 +! FOPI IS THE ACTUAL SOLAR FLUX IN ERGS/CM2 + REAL,SAVE :: F0PI=1.5E4 +! RHCH4 IS THE METHANE RH AT THE SURFACE + REAL,SAVE :: RHCH4=0.60 !! .65 +! FH2 IS THE CONSTANT MIXING RATIO OF H2 + REAL,SAVE :: FH2=0.04 !! 0.003 +! FHAZE IS THE HAZE PRODUCTION SCALING FACTOR + REAL,SAVE :: FHAZE=0.35 +! FHIR IS THE HAZE INFRARED ABSORPTION SCALE FACTOR + REAL,SAVE :: FHIR=0.7 ! anciennement 0.5 +! FHVIS IS THE HAZE INFRARED ABSORPTION SCALE FACTOR + REAL,SAVE :: FHVIS=1.1 !!!! 1.333333333/ +! TAUFAC IS THE 200 CM-1 SCALING FACTOR + REAL,SAVE :: TAUFAC=2.00 +! RCLOUD IS THE PARTICLE SIZE IN THE CLOUD IN MICRONS + REAL,SAVE ::RCLOUD=60. !! 100 + + REAL,SAVE :: FARGON=0. +END MODULE TGMDAT_MOD