- Timestamp:
- Sep 22, 2004, 11:48:34 AM (20 years ago)
- Location:
- LMDZ.3.3/branches/rel-LF/libf
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/dyn3d/bilan_dyn.F
r545 r546 160 160 itau=itau+1 161 161 162 ucont (:,:,:) = 0.0 162 163 if (first) then 163 164 -
LMDZ.3.3/branches/rel-LF/libf/dyn3d/etat0_netcdf.F
r513 r546 41 41 REAL :: solsw(klon), sollw(klon), fder(klon) 42 42 REAL :: deltat(klon), frugs(klon,nbsrf), agesno(klon,nbsrf) 43 REAL :: rugmer(klon) 43 REAL :: rugmer(klon), run_off_lic_0(klon) 44 44 REAL :: zmea(iip1*jjp1), zstd(iip1*jjp1) 45 45 REAL :: zsig(iip1*jjp1), zgam(iip1*jjp1), zthe(iip1*jjp1) … … 671 671 rnebcon = 0.0 672 672 ratqs = 0.0 673 run_off_lic_0 = 0.0 673 674 674 675 cIM call phyredem("startphy.nc",phystep,radpas, co2_ppm, solaire, … … 679 680 $ radsol, frugs, agesno, 680 681 $ zmea, zstd, zsig, zgam, zthe, zpic, zval, rugsrel, 681 $ t_ancien, q_ancien, rnebcon, ratqs, clwcon) 682 $ t_ancien, q_ancien, rnebcon, ratqs, clwcon, 683 $ run_off_lic_0) 682 684 683 685 C Sortie Visu pour les champs dynamiques -
LMDZ.3.3/branches/rel-LF/libf/phylmd/clmain.F
r522 r546 961 961 endif 962 962 963 psref = 0. 963 964 DO i = 1, knon 964 965 psref(i) = paprs(i,1) !pression de reference est celle au sol … … 1051 1052 1052 1053 c initialisation 1053 petAcoef =0.1054 petAcoef =0. 1054 1055 peqAcoef = 0. 1055 1056 petBcoef =0. 1056 1057 peqBcoef = 0. 1057 1058 p1lay =0. 1059 zlev1(:) = 0. 1060 temp_air(:) = 0. 1061 spechum(:) = 0. 1062 epot_air(:) = 0. 1063 tq_cdrag(:) = 0. 1058 1064 1059 1065 c do i = 1, knon -
LMDZ.3.3/branches/rel-LF/libf/phylmd/clouds_gno.F
r495 r546 62 62 ptconv=.false. 63 63 ratqsc=0. 64 64 CLDF = 0. 65 65 66 66 DO 500 K = 1, ND -
LMDZ.3.3/branches/rel-LF/libf/phylmd/cv3_routines.F
r495 r546 193 193 real A, B ! convect3 194 194 195 cLF initialisation 196 plcl(:)=0. 197 tnk(:)=0. 198 qnk(:)=0. 199 gznk(:)=0. 200 pnk(:)=0. 201 qsnk(:)=0. 202 rh(:)=0. 203 chi(:)=0. 195 204 c@ !------------------------------------------------------------------- 196 205 c@ ! --- Find level of minimum moist static energy … … 1245 1254 c ------------------------------------------------------- 1246 1255 1256 m(:,:) = 0. 1257 1247 1258 do k=1,nl 1248 1259 do i=1,ncum … … 1467 1478 c --- INITIALIZE VARIOUS ARRAYS USED IN THE COMPUTATIONS 1468 1479 c===================================================================== 1480 1481 sij(:,:,:) = 0. 1482 ment(:,:,:) = 0. 1483 csum(:,:) = 0. 1469 1484 1470 1485 c ori do 360 i=1,ncum*nlp … … 2158 2173 2159 2174 c inputs: 2160 integer ncum,nd,na,ntra,nloc 2175 integer ncum,nd,na,ntra,nloc, nent(nloc,na) 2161 2176 integer icb(nloc), inb(nloc) 2162 2177 real delt … … 2170 2185 real water(nloc,na), evap(nloc,na), b(nloc,na) 2171 2186 real ment(nloc,na,na), qent(nloc,na,na), uent(nloc,na,na) 2172 real vent(nloc,na,na), nent(nloc,na),elij(nloc,na,na)2187 real vent(nloc,na,na), elij(nloc,na,na) 2173 2188 real traent(nloc,na,na,ntra) 2174 2189 real tv(nloc,nd), tvp(nloc,nd) … … 2212 2227 enddo 2213 2228 2229 qcond(:,:) = 0. 2214 2230 do i=1,nd 2215 2231 do il=1,ncum -
LMDZ.3.3/branches/rel-LF/libf/phylmd/cv_driver.F
r418 r546 309 309 !--------------------------------------------------------------------- 310 310 311 !LF 312 qs(:,:) = 0. 313 qp(:,:) = 0. 314 ep = 0. 315 clw = 0. 316 qcondc = 0. 317 nent = 0 318 311 319 do 20 k=1,nd 312 320 do 10 i=1,len -
LMDZ.3.3/branches/rel-LF/libf/phylmd/fisrtilp.F
r410 r546 374 374 c print*,'ZDQS ',zdqs(i) 375 375 c--Olivier 376 zdelq = min(ratqs(i,k),0.99) * zq(i) 376 377 rhcl(i,k)=(zqs(i)+zq(i)-zdelq)/2./zqs(i) 377 378 IF (rneb(i,k) .LE. 0.0) rhcl(i,k)=zq(i)/zqs(i) -
LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90
r545 r546 1340 1340 INTEGER :: ndexct(iim*(jjm+1)),ndexcs(iim*(jjm+1)) 1341 1341 REAL :: zx_lon(iim,jjm+1), zx_lat(iim,jjm+1), zjulian 1342 integer:: idayref1342 INTEGER,save :: idayref 1343 1343 !med integer :: itau_w 1344 1344 integer,save :: itau_w -
LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F
r544 r546 353 353 REAL nhistoW(kmaxm1,lmaxm1,iwmax,nbregdyn) 354 354 REAL nhistoWt(kmaxm1,lmaxm1,iwmax,nbregdyn) 355 SAVE nhistoWt355 SAVE histoW, nhistoW, nhistoWt 356 356 357 357 INTEGER linv … … 382 382 save ok_regdyn, ecrit_regdyn, nid_regdyn 383 383 384 #define histREGDYN 384 385 #undef histREGDYN 385 #define histREGDYN386 386 #ifdef histREGDYN 387 387 c data ok_regdyn,ecrit_regdyn/.true.,0.125/ … … 611 611 REAL sens(klon), dsens(klon) ! chaleur sensible et sa derivee 612 612 REAL dlw(klon) ! derivee infra rouge 613 SAVE dlw 613 614 REAL bils(klon) ! bilan de chaleur au sol 614 615 REAL wfbils(klon,nbsrf) ! bilan de chaleur au sol, pour chaque … … 811 812 c 812 813 INTEGER ibas_con(klon), itop_con(klon) 814 SAVE itop_con 813 815 REAL rain_con(klon), rain_lsc(klon) 814 816 REAL snow_con(klon), snow_lsc(klon) … … 1217 1219 c 1218 1220 c 1221 cLF Initialisations de variables divers (A.Caubel, M.E. Demory) 1222 clwcon0(:,:)= 0.0 1223 zx_tmp_2d(:,:)= 0.0 1224 paire_ter(:) = 0. 1225 snow_con(:) = 0. 1226 snow_lsc(:) = 0. 1219 1227 ENDIF 1220 1228 c
Note: See TracChangeset
for help on using the changeset viewer.