Changeset 559
- Timestamp:
- Nov 9, 2004, 5:55:41 PM (20 years ago)
- Location:
- LMDZ4/trunk/libf
- Files:
-
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/dyn3d/bilan_dyn.F
r541 r559 298 298 299 299 c énergie cinétique 300 ucont(:,:,:)=0 300 301 CALL covcont(llm,ucov,vcov,ucont,vcont) 301 302 CALL enercin(vcov,ucov,vcont,ucont,ecin) -
LMDZ4/trunk/libf/dyn3d/leapfrog.F
r550 r559 4 4 c 5 5 c 6 SUBROUTINE leapfrog(ucov,vcov,teta,ps,masse,phis,nq,q,clesphy0) 6 SUBROUTINE leapfrog(ucov,vcov,teta,ps,masse,phis,nq,q,clesphy0, 7 & time_0) 7 8 8 9 #ifdef INCA -
LMDZ4/trunk/libf/phylmd/clouds_gno.F
r524 r559 58 58 LOGICAL lconv(klon) 59 59 60 60 cym 61 cldf(:,:)=0.0 62 61 63 pi = ACOS(-1.) 62 64 sqrtpi=sqrt(pi) -
LMDZ4/trunk/libf/phylmd/concvl.F
r524 r559 88 88 c 89 89 c 90 cym 91 snow(:)=0 92 90 93 IF (ifrst .EQ. 0) THEN 91 94 ifrst = 1 -
LMDZ4/trunk/libf/phylmd/cv3_routines.F
r524 r559 194 194 real pnk(len), qsnk(len), rh(len), chi(len) 195 195 real A, B ! convect3 196 196 cym 197 plcl=0.0 197 198 c@ !------------------------------------------------------------------- 198 199 c@ ! --- Find level of minimum moist static energy … … 1488 1489 vent(i,k,j)=v(i,j) 1489 1490 elij(i,k,j)=0.0 1490 ment(i,k,j)=0.01491 sij(i,k,j)=0.01491 cym ment(i,k,j)=0.0 1492 cym sij(i,k,j)=0.0 1492 1493 385 continue 1493 1494 390 continue 1494 1495 400 continue 1495 1496 1497 cym 1498 ment(1:ncum,1:nd,1:nd)=0.0 1499 sij(1:ncum,1:nd,1:nd)=0.0 1500 1496 1501 c do k=1,ntra 1497 1502 c do j=1,nd ! instead nlp … … 1620 1625 c===================================================================== 1621 1626 1622 call zilch(asum,ncum*nd) 1623 call zilch(bsum,ncum*nd) 1624 call zilch(csum,ncum*nd) 1627 cym call zilch(asum,ncum*nd) 1628 cym call zilch(bsum,ncum*nd) 1629 cym call zilch(csum,ncum*nd) 1630 call zilch(asum,nloc*nd) 1631 call zilch(csum,nloc*nd) 1632 call zilch(csum,nloc*nd) 1625 1633 1626 1634 do il=1,ncum -
LMDZ4/trunk/libf/phylmd/cv_driver.F
r524 r559 321 321 tp1(i,k)=0.0 322 322 clw1(i,k)=0.0 323 cym 324 clw(i,k)=0.0 323 325 gz1(i,k) = 0. 324 326 -
LMDZ4/trunk/libf/phylmd/fisrtilp.F
r524 r559 123 123 c 124 124 DATA appel1er /.TRUE./ 125 125 cym 126 zdelq=0.0 127 126 128 IF (appel1er) THEN 127 129 c -
LMDZ4/trunk/libf/phylmd/ini_histmth.h
r548 r559 739 739 c endif 740 740 c 741 c CALL histdef(nid_mth, "ozone", "Ozone concentration", "-",742 c . iim,jjmp1,nhori, klev,1,klev,nvert, 32,743 c . "ave(X)", zsto,zout)744 c745 741 if (nqmax.GE.3) THEN 746 742 DO iq=1,nqmax-2 -
LMDZ4/trunk/libf/phylmd/interface_surf.F90
r557 r559 300 300 call abort_gcm(modname,abort_message,1) 301 301 endif 302 cym 303 run_off=0.0 304 cym 305 302 306 !!$PB 303 307 ALLOCATE (tmp_rriv(iim,jjm+1), stat=error) -
LMDZ4/trunk/libf/phylmd/physiq.F
r557 r559 625 625 REAL sens(klon), dsens(klon) ! chaleur sensible et sa derivee 626 626 REAL dlw(klon) ! derivee infra rouge 627 cym 628 SAVE dlw 629 cym 627 630 REAL bils(klon) ! bilan de chaleur au sol 628 631 REAL wfbils(klon,nbsrf) ! bilan de chaleur au sol, pour chaque … … 824 827 c 825 828 INTEGER ibas_con(klon), itop_con(klon) 829 cym 830 SAVE ibas_con,itop_con 831 cym 826 832 REAL rain_con(klon), rain_lsc(klon) 827 833 REAL snow_con(klon), snow_lsc(klon) … … 974 980 LOGICAL ok_ade, ok_aie ! Apply aerosol (in)direct effects or not 975 981 REAL bl95_b0, bl95_b1 ! Parameter in Boucher and Lohmann (1995) 982 cym 983 SAVE ok_ade, ok_aie, bl95_b0, bl95_b1 984 cym 976 985 cjq-end 977 986 c … … 997 1006 bl95_b1=0. 998 1007 !rv 1008 cym 1009 wfbils(:,:)=0 1010 cym 999 1011 modname = 'physiq' 1000 1012 IF (if_ebil.ge.1) THEN … … 2361 2373 CALL aeropt(pplay, paprs, t_seri, sulfate, rhcl, 2362 2374 . tau_ae, piz_ae, cg_ae, aerindex) 2375 cym 2376 ELSE 2377 tau_ae(:,:,:)=0.0 2378 piz_ae(:,:,:)=0.0 2379 cg_ae(:,:,:)=0.0 2380 cym 2363 2381 ENDIF 2364 2382
Note: See TracChangeset
for help on using the changeset viewer.