Changeset 1106 for LMDZ4/branches/LMDZ4-dev
- Timestamp:
- Feb 24, 2009, 3:46:23 PM (16 years ago)
- Location:
- LMDZ4/branches/LMDZ4-dev/libf
- Files:
-
- 9 edited
-
dyn3d/abort_gcm.F (modified) (1 diff)
-
dyn3d/comconst.h (modified) (2 diffs)
-
dyn3d/comgeom2.h (modified) (1 diff)
-
dyn3dpar/comconst.h (modified) (2 diffs)
-
dyn3dpar/comgeom2.h (modified) (1 diff)
-
phylmd/calcul_fluxs_mod.F90 (modified) (1 diff)
-
phylmd/oasis.F90 (modified) (1 diff)
-
phylmd/radiation_AR4.F (modified) (3 diffs)
-
phylmd/screenp.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4-dev/libf/dyn3d/abort_gcm.F
r1004 r1106 19 19 C ierr = severity of situation ( = 0 normal ) 20 20 21 character *(*) modname21 character(len=*) modname 22 22 integer ierr 23 character *(*) message23 character(len=*) message 24 24 25 25 ! write(lunout,*) 'in abort_gcm' -
LMDZ4/branches/LMDZ4-dev/libf/dyn3d/comconst.h
r524 r1106 2 2 ! $Header$ 3 3 ! 4 c-----------------------------------------------------------------------5 cINCLUDE comconst.h4 !----------------------------------------------------------------------- 5 ! INCLUDE comconst.h 6 6 7 COMMON/comconst/im,jm,lllm,imp1,jmp1,lllmm1,lllmp1,lcl, 8 * dtvr,daysec,9 *pi,dtphys,dtdiss,rad,r,cpp,kappa,cotot,unsim,g,omeg7 COMMON/comconst/im,jm,lllm,imp1,jmp1,lllmm1,lllmp1,lcl, & 8 & dtvr,daysec, & 9 & pi,dtphys,dtdiss,rad,r,cpp,kappa,cotot,unsim,g,omeg 10 10 11 11 INTEGER im,jm,lllm,imp1,jmp1,lllmm1,lllmp1,lcl … … 14 14 REAL cotot,unsim,g,omeg 15 15 16 c-----------------------------------------------------------------------16 !----------------------------------------------------------------------- -
LMDZ4/branches/LMDZ4-dev/libf/dyn3d/comgeom2.h
r524 r1106 2 2 ! $Header$ 3 3 ! 4 *CDK comgeom25 COMMON/comgeom/ 6 1 cu(iip1,jjp1),cv(iip1,jjm),unscu2(iip1,jjp1),unscv2(iip1,jjm) ,7 2 aire(iip1,jjp1),airesurg(iip1,jjp1),aireu(iip1,jjp1) ,8 3 airev(iip1,jjm),unsaire(iip1,jjp1),apoln,apols ,9 4 unsairez(iip1,jjm),airuscv2(iip1,jjm),airvscu2(iip1,jjm) ,10 5 aireij1(iip1,jjp1),aireij2(iip1,jjp1),aireij3(iip1,jjp1) ,11 6 aireij4(iip1,jjp1),alpha1(iip1,jjp1),alpha2(iip1,jjp1) ,12 7 alpha3(iip1,jjp1),alpha4(iip1,jjp1),alpha1p2(iip1,jjp1) ,13 8 alpha1p4(iip1,jjp1),alpha2p3(iip1,jjp1),alpha3p4(iip1,jjp1) ,14 9 fext(iip1,jjm),constang(iip1,jjp1), rlatu(jjp1),rlatv(jjm),15 1 rlonu(iip1),rlonv(iip1),cuvsurcv(iip1,jjm),cvsurcuv(iip1,jjm) ,16 1 cvusurcu(iip1,jjp1),cusurcvu(iip1,jjp1) ,17 2 cuvscvgam1(iip1,jjm),cuvscvgam2(iip1,jjm),cvuscugam1(iip1,jjp1),18 3 cvuscugam2(iip1,jjp1),cvscuvgam(iip1,jjm),cuscvugam(iip1,jjp1) ,19 4 unsapolnga1,unsapolnga2,unsapolsga1,unsapolsga2 ,20 5 unsair_gam1(iip1,jjp1),unsair_gam2(iip1,jjp1) ,21 6 unsairz_gam(iip1,jjm),aivscu2gam(iip1,jjm),aiuscv2gam(iip1,jjm)22 7, xprimu(iip1),xprimv(iip1)4 !CDK comgeom2 5 COMMON/comgeom/ & 6 & cu(iip1,jjp1),cv(iip1,jjm),unscu2(iip1,jjp1),unscv2(iip1,jjm) , & 7 & aire(iip1,jjp1),airesurg(iip1,jjp1),aireu(iip1,jjp1) , & 8 & airev(iip1,jjm),unsaire(iip1,jjp1),apoln,apols , & 9 & unsairez(iip1,jjm),airuscv2(iip1,jjm),airvscu2(iip1,jjm) , & 10 & aireij1(iip1,jjp1),aireij2(iip1,jjp1),aireij3(iip1,jjp1) , & 11 & aireij4(iip1,jjp1),alpha1(iip1,jjp1),alpha2(iip1,jjp1) , & 12 & alpha3(iip1,jjp1),alpha4(iip1,jjp1),alpha1p2(iip1,jjp1) , & 13 & alpha1p4(iip1,jjp1),alpha2p3(iip1,jjp1),alpha3p4(iip1,jjp1) , & 14 & fext(iip1,jjm),constang(iip1,jjp1), rlatu(jjp1),rlatv(jjm), & 15 & rlonu(iip1),rlonv(iip1),cuvsurcv(iip1,jjm),cvsurcuv(iip1,jjm) , & 16 & cvusurcu(iip1,jjp1),cusurcvu(iip1,jjp1) , & 17 & cuvscvgam1(iip1,jjm),cuvscvgam2(iip1,jjm),cvuscugam1(iip1,jjp1), & 18 & cvuscugam2(iip1,jjp1),cvscuvgam(iip1,jjm),cuscvugam(iip1,jjp1) , & 19 & unsapolnga1,unsapolnga2,unsapolsga1,unsapolsga2 , & 20 & unsair_gam1(iip1,jjp1),unsair_gam2(iip1,jjp1) , & 21 & unsairz_gam(iip1,jjm),aivscu2gam(iip1,jjm),aiuscv2gam(iip1,jjm) & 22 & , xprimu(iip1),xprimv(iip1) 23 23 24 c 25 REAL 26 1 cu,cv,unscu2,unscv2,aire,airesurg,aireu,airev,apoln,apols,unsaire 27 2 ,unsairez,airuscv2,airvscu2,aireij1,aireij2,aireij3,aireij4 , 28 3 alpha1,alpha2,alpha3,alpha4,alpha1p2,alpha1p4,alpha2p3,alpha3p4 , 29 4 fext,constang,rlatu,rlatv,rlonu,rlonv,cuvscvgam1,cuvscvgam2 , 30 5 cvuscugam1,cvuscugam2,cvscuvgam,cuscvugam,unsapolnga1 , 31 6 unsapolnga2,unsapolsga1,unsapolsga2,unsair_gam1,unsair_gam2 , 32 7 unsairz_gam,aivscu2gam,aiuscv2gam,cuvsurcv,cvsurcuv,cvusurcu , 33 8 cusurcvu,xprimu,xprimv 34 c 24 25 REAL & 26 & cu,cv,unscu2,unscv2,aire,airesurg,aireu,airev,apoln,apols,unsaire & 27 & ,unsairez,airuscv2,airvscu2,aireij1,aireij2,aireij3,aireij4 , & 28 & alpha1,alpha2,alpha3,alpha4,alpha1p2,alpha1p4,alpha2p3,alpha3p4 , & 29 & fext,constang,rlatu,rlatv,rlonu,rlonv,cuvscvgam1,cuvscvgam2 , & 30 & cvuscugam1,cvuscugam2,cvscuvgam,cuscvugam,unsapolnga1 , & 31 & unsapolnga2,unsapolsga1,unsapolsga2,unsair_gam1,unsair_gam2 , & 32 & unsairz_gam,aivscu2gam,aiuscv2gam,cuvsurcv,cvsurcuv,cvusurcu , & 33 & cusurcvu,xprimu,xprimv -
LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/comconst.h
r774 r1106 2 2 ! $Header$ 3 3 ! 4 c-----------------------------------------------------------------------5 cINCLUDE comconst.h4 !----------------------------------------------------------------------- 5 ! INCLUDE comconst.h 6 6 7 COMMON/comconst/im,jm,lllm,imp1,jmp1,lllmm1,lllmp1,lcl, 8 * dtvr,daysec,9 *pi,dtphys,dtdiss,rad,r,cpp,kappa,cotot,unsim,g,omeg7 COMMON/comconst/im,jm,lllm,imp1,jmp1,lllmm1,lllmp1,lcl, & 8 & dtvr,daysec, & 9 & pi,dtphys,dtdiss,rad,r,cpp,kappa,cotot,unsim,g,omeg 10 10 11 11 INTEGER im,jm,lllm,imp1,jmp1,lllmm1,lllmp1,lcl … … 14 14 REAL cotot,unsim,g,omeg 15 15 16 c-----------------------------------------------------------------------16 !----------------------------------------------------------------------- -
LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/comgeom2.h
r774 r1106 2 2 ! $Header$ 3 3 ! 4 *CDK comgeom25 COMMON/comgeom/ 6 1 cu(iip1,jjp1),cv(iip1,jjm),unscu2(iip1,jjp1),unscv2(iip1,jjm) ,7 2 aire(iip1,jjp1),airesurg(iip1,jjp1),aireu(iip1,jjp1) ,8 3 airev(iip1,jjm),unsaire(iip1,jjp1),apoln,apols ,9 4 unsairez(iip1,jjm),airuscv2(iip1,jjm),airvscu2(iip1,jjm) ,10 5 aireij1(iip1,jjp1),aireij2(iip1,jjp1),aireij3(iip1,jjp1) ,11 6 aireij4(iip1,jjp1),alpha1(iip1,jjp1),alpha2(iip1,jjp1) ,12 7 alpha3(iip1,jjp1),alpha4(iip1,jjp1),alpha1p2(iip1,jjp1) ,13 8 alpha1p4(iip1,jjp1),alpha2p3(iip1,jjp1),alpha3p4(iip1,jjp1) ,14 9 fext(iip1,jjm),constang(iip1,jjp1), rlatu(jjp1),rlatv(jjm),15 1 rlonu(iip1),rlonv(iip1),cuvsurcv(iip1,jjm),cvsurcuv(iip1,jjm) ,16 1 cvusurcu(iip1,jjp1),cusurcvu(iip1,jjp1) ,17 2 cuvscvgam1(iip1,jjm),cuvscvgam2(iip1,jjm),cvuscugam1(iip1,jjp1),18 3 cvuscugam2(iip1,jjp1),cvscuvgam(iip1,jjm),cuscvugam(iip1,jjp1) ,19 4 unsapolnga1,unsapolnga2,unsapolsga1,unsapolsga2 ,20 5 unsair_gam1(iip1,jjp1),unsair_gam2(iip1,jjp1) ,21 6 unsairz_gam(iip1,jjm),aivscu2gam(iip1,jjm),aiuscv2gam(iip1,jjm)22 7, xprimu(iip1),xprimv(iip1)4 !CDK comgeom2 5 COMMON/comgeom/ & 6 & cu(iip1,jjp1),cv(iip1,jjm),unscu2(iip1,jjp1),unscv2(iip1,jjm) , & 7 & aire(iip1,jjp1),airesurg(iip1,jjp1),aireu(iip1,jjp1) , & 8 & airev(iip1,jjm),unsaire(iip1,jjp1),apoln,apols , & 9 & unsairez(iip1,jjm),airuscv2(iip1,jjm),airvscu2(iip1,jjm) , & 10 & aireij1(iip1,jjp1),aireij2(iip1,jjp1),aireij3(iip1,jjp1) , & 11 & aireij4(iip1,jjp1),alpha1(iip1,jjp1),alpha2(iip1,jjp1) , & 12 & alpha3(iip1,jjp1),alpha4(iip1,jjp1),alpha1p2(iip1,jjp1) , & 13 & alpha1p4(iip1,jjp1),alpha2p3(iip1,jjp1),alpha3p4(iip1,jjp1) , & 14 & fext(iip1,jjm),constang(iip1,jjp1), rlatu(jjp1),rlatv(jjm), & 15 & rlonu(iip1),rlonv(iip1),cuvsurcv(iip1,jjm),cvsurcuv(iip1,jjm) , & 16 & cvusurcu(iip1,jjp1),cusurcvu(iip1,jjp1) , & 17 & cuvscvgam1(iip1,jjm),cuvscvgam2(iip1,jjm),cvuscugam1(iip1,jjp1), & 18 & cvuscugam2(iip1,jjp1),cvscuvgam(iip1,jjm),cuscvugam(iip1,jjp1) , & 19 & unsapolnga1,unsapolnga2,unsapolsga1,unsapolsga2 , & 20 & unsair_gam1(iip1,jjp1),unsair_gam2(iip1,jjp1) , & 21 & unsairz_gam(iip1,jjm),aivscu2gam(iip1,jjm),aiuscv2gam(iip1,jjm) & 22 & , xprimu(iip1),xprimv(iip1) 23 23 24 c 25 REAL 26 1 cu,cv,unscu2,unscv2,aire,airesurg,aireu,airev,apoln,apols,unsaire 27 2 ,unsairez,airuscv2,airvscu2,aireij1,aireij2,aireij3,aireij4 , 28 3 alpha1,alpha2,alpha3,alpha4,alpha1p2,alpha1p4,alpha2p3,alpha3p4 , 29 4 fext,constang,rlatu,rlatv,rlonu,rlonv,cuvscvgam1,cuvscvgam2 , 30 5 cvuscugam1,cvuscugam2,cvscuvgam,cuscvugam,unsapolnga1 , 31 6 unsapolnga2,unsapolsga1,unsapolsga2,unsair_gam1,unsair_gam2 , 32 7 unsairz_gam,aivscu2gam,aiuscv2gam,cuvsurcv,cvsurcuv,cvusurcu , 33 8 cusurcvu,xprimu,xprimv 34 c 24 25 REAL & 26 & cu,cv,unscu2,unscv2,aire,airesurg,aireu,airev,apoln,apols,unsaire & 27 & ,unsairez,airuscv2,airvscu2,aireij1,aireij2,aireij3,aireij4 , & 28 & alpha1,alpha2,alpha3,alpha4,alpha1p2,alpha1p4,alpha2p3,alpha3p4 , & 29 & fext,constang,rlatu,rlatv,rlonu,rlonv,cuvscvgam1,cuvscvgam2 , & 30 & cvuscugam1,cvuscugam2,cvscuvgam,cuscvugam,unsapolnga1 , & 31 & unsapolnga2,unsapolsga1,unsapolsga2,unsair_gam1,unsair_gam2 , & 32 & unsairz_gam,aivscu2gam,aiuscv2gam,cuvsurcv,cvsurcuv,cvusurcu , & 33 & cusurcvu,xprimu,xprimv -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/calcul_fluxs_mod.F90
r1067 r1106 99 99 WRITE(*,*)' radsol (min, max)', & 100 100 MINVAL(radsol(1:knon)), MAXVAL(radsol(1:knon)) 101 CALL flush(6)102 101 ENDIF 103 102 -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/oasis.F90
r1079 r1106 280 280 WRITE (nuout,*) 'Fromcpl: Reading fields from CPL, ktime=',ktime 281 281 WRITE (nuout,*) ' ' 282 CALL flush (nuout)283 282 284 283 istart=ii_begin -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/radiation_AR4.F
r998 r1106 2123 2123 ZALPHA = 3.*ZWCP*ZRM2*(1.+ZGP*ZWM)/ZX2 2124 2124 ZBETA = 3.*ZWCP* PRMUZ(JL) *(1.+3.*ZGP*ZRM2*ZWM)/ZX2 2125 CMAF ZARG=MIN(ZTOP/PRMUZ(JL),200.) 2126 ZARG=MIN(ZTOP/PRMUZ(JL),2.0d+2) 2125 ZARG=MIN(ZTOP/PRMUZ(JL),200._8) 2127 2126 ZEXMU0=EXP(-ZARG) 2128 CMAF ZARG2=MIN(ZRK*ZTOP,200.) 2129 ZARG2=MIN(ZRK*ZTOP,2.0d+2) 2127 ZARG2=MIN(ZRK*ZTOP,200._8) 2130 2128 ZEXKP=EXP(ZARG2) 2131 2129 ZEXKM = 1./ZEXKP … … 2898 2896 ZTX2=ZTX*ZTX 2899 2897 ZZABLY = ZABLY(JL,6,JAE1)+ZABLY(JL,6,JAE2)+ZABLY(JL,6,JAE3) 2900 CMAF ZUP=MIN( MAX( 0.5*R10E*LOG( ZZABLY ) + 5., 0.), 6.0) 2901 ZUP=MIN( MAX( 0.5*R10E*LOG( ZZABLY ) + 5., 0.d+0), 6.d+0) 2898 ZUP=MIN( MAX( 0.5*R10E*LOG( ZZABLY ) + 5., 0._8), 6._8) 2902 2899 ZCAH1=AT(1,1)+ZUP*(AT(1,2)+ZUP*(AT(1,3))) 2903 2900 ZCBH1=BT(1,1)+ZUP*(BT(1,2)+ZUP*(BT(1,3))) … … 2930 2927 ZZABLY = ZABLY(JL,9,JAE1)+ZABLY(JL,9,JAE2)+ZABLY(JL,9,JAE3) 2931 2928 ZALUP = R10E * LOG ( ZZABLY ) 2932 CMAF ZUP = MAX( 0.0 , 5.0 + 0.5 * ZALUP ) 2933 ZUP = MAX( 0.d+0 , 5.0 + 0.5 * ZALUP ) 2929 ZUP = MAX( 0._8, 5.0 + 0.5 * ZALUP ) 2934 2930 ZPSC2(JL) = (ZTAVI(JL)/TREF) ** ZUP 2935 2931 ZCAC8=AT(8,1)+ZUP*(AT(8,2)+ZUP*(AT(8,3))) -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/screenp.F90
r902 r1106 64 64 delu(i) = (ustar(i)/RKAR)* & 65 65 (log(zref/(rugos(i))+1.) + & 66 min(5. 0, 5.0 *(zref - rugos(i))/lmon(i)))66 min(5.d0, 5.0 *(zref - rugos(i))/lmon(i))) 67 67 delte(i) = (testar(i)/RKAR)* & 68 68 (log(zref/(rugos(i))+1.) + & 69 min(5. 0, 5.0 * (zref - rugos(i))/lmon(i)))69 min(5.d0, 5.0 * (zref - rugos(i))/lmon(i))) 70 70 delq(i) = (qstar(i)/RKAR)* & 71 71 (log(zref/(rugos(i))+1.) + & 72 min(5. 0, 5.0 * (zref - rugos(i))/lmon(i)))72 min(5.d0, 5.0 * (zref - rugos(i))/lmon(i))) 73 73 ELSE 74 74 delu(i) = 0.1 * speed(i)
Note: See TracChangeset
for help on using the changeset viewer.
